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.

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.

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.

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.

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.

Aims. We present a spectroscopic study of the properties of 64 Balmer break galaxies that show signs of star formation. The studied sample of star-forming galaxies spans a redshift range from 0.094 to 1.475 with stellar masses in the range 10$^{8}-$10$^{12}$ $M_{\odot}$. The sample also includes eight broad emission line galaxies with redshifts between 1.5 $<z<$ 3.0. Methods. We derived star formation rates (SFRs) from emission line luminosities and investigated the dependence of the SFR and specific SFR (SSFR) on the stellar mass and color. Furthermore, we investigated the evolution of these relations with the redshift. Results. We found that the SFR correlates with the stellar mass, our data is consistent with previous results from other authors in that there is a break in the correlation, which reveals the presence of massive galaxies with lower SFR values (i.e., decreasing star formation). We also note an anticorrelation for the SSFR with the stellar mass. Again in this case, our data is also consistent with a break in the correlation, revealing the presence of massive star-forming galaxies with lower SSFR values, thereby increasing the anticorrelation. These results might suggest a characteristic mass ($M_{0}$) at which the red sequence could mostly be assembled. In addition, at a given stellar mass, high-redshift galaxies have on average higher SFR and SSFR values than local galaxies. Finally, we explored whether a similar trend could be observed with redshift in the SSFR$-(u-B)$ color diagram, and we hypothesize that a possible $(u-B)_{0}$ break color may define a characteristic color for the formation of the red sequence.

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

We report the multi-wavelength identification of the X-ray sources found in the Subaru-XMM-Newton Deep Survey (SXDS) using deep imaging data covering the wavelength range between the far-UV to the mid-IR. We select a primary counterpart of each X-ray source by applying the likelihood ratio method to R-band, 3.6micron, near-UV, and 24micron source catalogs as well as matching catalogs of AGN candidates selected in 1.4GHz radio and i'-band variability surveys. Once candidates of Galactic stars, ultra-luminous X-ray sources in a nearby galaxy, and clusters of galaxies are removed there are 896 AGN candidates in the sample. We conduct spectroscopic observations of the primary counterparts with multi-object spectrographs in the optical and NIR; 65\% of the X-ray AGN candidates are spectroscopically-identified. For the remaining X-ray AGN candidates, we evaluate their photometric redshift with photometric data in 15 bands. Utilising the multi-wavelength photometric data of the large sample of X-ray selected AGNs, we evaluate the stellar masses, M*, of the host galaxies of the narrow-line AGNs. The distribution of the stellar mass is remarkably constant from z=0.1 to 4.0. The relation between M* and 2--10 keV luminosity can be explained with strong cosmological evolution of the relationship between the black hole mass and M*. We also evaluate the scatter of the UV-MIR spectral energy distribution (SED) of the X-ray AGNs as a function of X-ray luminosity and absorption to the nucleus. The scatter is compared with galaxies which have redshift and stellar mass distribution matched with the X-ray AGN. The UV-NIR SEDs of obscured X-ray AGNs are similar to those of the galaxies in the matched sample. In the NIR-MIR range, the median SEDs of X-ray AGNs are redder, but the scatter of the SEDs of the X-ray AGN broadly overlaps that of the galaxies in the matched sample.

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.

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 present a spectroscopic study with the derivation of the physical properties of 37 Balmer break galaxies, which have the necessary lines to locate them in star-forming-AGN diagnostic diagrams. These galaxies span a redshift range from 0.045 to 0.93 and are somewhat less massive than similar samples of previous works. The studied sample has multiwavelength photometric data coverage from the ultraviolet to MIR Spitzer bands. We investigate the connection between star formation and AGN activity via optical, mass-excitation (MEx) and MIR diagnostic diagrams. Through optical diagrams, 31 (84%) star-forming galaxies, 2 (5%) composite galaxies and 3 (8%) AGNs were classified, whereas from the MEx diagram only one galaxy was classified as AGN. A total of 19 galaxies have photometry available in all the IRAC/Spitzer bands. Of these, 3 AGN candidates were not classified as AGN in the optical diagrams, suggesting they are dusty/obscured AGNs, or that nuclear star formation has diluted their contributions. Furthermore, the relationship between SFR surface density (\Sigma_{SFR}) and stellar mass surface density per time unit (\Sigma_{M_{\ast}/\tau}) as a function of redshift was investigated using the [OII] \lambda3727, 3729, H\alpha \lambda6563 luminosities, which revealed that both quantities are larger for higher redshift galaxies. We also studied the SFR and SSFR versus stellar mass and color relations, with the more massive galaxies having higher SFR values but lower SSFR values than less massive galaxies. These results are consistent with previous ones showing that, at a given mass, high-redshift galaxies have on average larger SFR and SSFR values than low-redshift galaxies. Finally, bluer galaxies have larger SSFR values than redder galaxies and for a given color the SSFR is larger for higher redshift galaxies.

In order to investigate the growth of super-massive black holes (SMBHs), we construct the black hole mass function (BHMF) and Eddington ratio distribution function (ERDF) of X-ray-selected broad-line AGNs at z~1.4 in the Subaru XMM-Newton Deep Survey field. In this redshift range, a significant part of the accretion growth of SMBHs is thought to be taking place. Black hole masses of X-ray-selected broad-line AGNs are estimated using the width of the broad MgII line and the 3000A monochromatic luminosity. We supplement the MgII FWHM values with the Ha FWHM obtained from our NIR spectroscopic survey. Using the black hole masses of broad-line AGNs at redshifts between 1.18 and 1.68, the binned broad-line AGN BHMF and ERDF are calculated using the Vmax method. To properly account for selection effects that impact the binned estimates, we derive the corrected broad-line AGN BHMF and ERDF by applying the Maximum Likelihood method, assuming that the ERDF is constant regardless of the black hole mass. We do not correct for the non-negligible uncertainties in virial BH mass estimates. If we compare the corrected broad-line AGN BHMF with that in the local Universe, the corrected BHMF at z~1.4 has a higher number density above 10^8 Msolar but a lower number density below that mass range. The evolution may be indicative of a down-sizing trend of accretion activity among the SMBH population. The evolution of broad-line AGN ERDF from z=1.4 to 0 indicates that the fraction of broad-line AGNs with accretion rate close to the Eddington-limit is higher at higher redshifts.

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.

Context: It is crucial to develop a method for classifying objects detected in deep surveys at infrared wavelengths. We specifically need a method to separate galaxies from stars using only the infrared information to study the properties of galaxies, e.g., to estimate the angular correlation function, without introducing any additional bias. Aims. We aim to separate stars and galaxies in the data from the AKARI North Ecliptic Pole (NEP) Deep survey collected in nine AKARI / IRC bands from 2 to 24 {\mu}m that cover the near- and mid-infrared wavelengths (hereafter NIR and MIR). We plan to estimate the correlation function for NIR and MIR galaxies from a sample selected according to our criteria in future research. Methods: We used support vector machines (SVM) to study the distribution of stars and galaxies in the AKARIs multicolor space. We defined the training samples of these objects by calculating their infrared stellarity parameter (sgc). We created the most efficient classifier and then tested it on the whole sample. We confirmed the developed separation with auxiliary optical data obtained by the Subaru telescope and by creating Euclidean normalized number count plots. Results: We obtain a 90% accuracy in pinpointing galaxies and 98% accuracy for stars in infrared multicolor space with the infrared SVM classifier. The source counts and comparison with the optical data (with a consistency of 65% for selecting stars and 96% for galaxies) confirm that our star/galaxy separation methods are reliable. Conclusions: The infrared classifier derived with the SVM method based on infrared sgc- selected training samples proves to be very efficient and accurate in selecting stars and galaxies in deep surveys at infrared wavelengths carried out without any previous target object selection.

We present a new catalogue of mid-IR sources using the AKARI NEP-Deep survey. The InfraRed Camera (IRC) onboard AKARI has a comprehensive mid-IR wavelength coverage with 9 photometric bands at 2 - 24 micron. We utilized all of these bands to cover a nearly circular area adjacent to the North Ecliptic Pole (NEP). We designed the catalogue to include most of sources detected in 7, 9, 11, 15 and 18 micron bands, and found 7284 sources in a 0.67 deg^2 area. From our simulations, we estimate that the catalogue is ~80 per cent complete to 200 micro Jy at 15 - 18 micron, and ~10 per cent of sources are missed, owing to source blending. Star-galaxy separation is conducted using only AKARI photometry, as a result of which 10 per cent of catalogued sources are found to be stars. The number counts at 11, 15, 18, and 24 micron are presented for both stars and galaxies. A drastic increase in the source density is found in between 11 and 15 micron at the flux level of ~300 micro Jy. This is likely due to the redshifted PAH emission at 8 micron, given our rough estimate of redshifts from an AKARI colour-colour plot. Along with the mid-IR source catalogue, we present optical-NIR photometry for sources falling inside a Subaru/Sprime-cam image covering part of the AKARI NEP-Deep field, which is deep enough to detect most of AKARI mid-IR sources, and useful to study optical characteristics of a complete mid-IR source sample.

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.

The Westerbork Radio Synthesis Telescope, WSRT, has been used to make a deep radio survey of an ~ 1.7 sq degree field coinciding with the AKARI North Ecliptic Pole Deep Field. The observations, data reduction and source count analysis are presented, along with a description of the overall scientific objectives. The survey consisted of 10 pointings, mosaiced with enough overlap to maintain a similar sensitivity across the central region that reached as low as 21 microJy per beam at 1.4 GHz. A catalogue containing 462 sources detected with a resolution of 17"x15" is presented. The differential source counts calculated from the WSRT data have been compared with those from the shallow VLA-NEP survey of Kollgaard et al 1994, and show a pronounced excess for sources fainter than ~ 1 mJy, consistent with the presence of a population of star forming galaxies at sub-mJy flux levels. The AKARI North Ecliptic Pole Deep field is the focus of a major observing campaign conducted across the entire spectral region. The combination of these data sets, along with the deep nature of the radio observations will allow unique studies of a large range of topics including the redshift evolution of the luminosity function of radio sources, the clustering environment of radio galaxies, the nature of obscured radio-loud active galactic nuclei, and the radio/far-infrared correlation for distant galaxies. This catalogue provides the basic data set for a future series of paper dealing with source identifications, morphologies, and the associated properties of the identified radio sources.

We report the detection and measurement of the absolute brightness and spatial fluctuations of the cosmic infrared background (CIB) with the AKARI satellite. We have carried out observations at 65, 90, 140 and 160 um as a cosmological survey in AKARI Deep Field South (ADF-S), which is one of the lowest cirrus regions with contiguous area on the sky. After removing bright galaxies and subtracting zodiacal and Galactic foregrounds from the measured sky brightness, we have successfully measured the CIB brightness and its fluctuations across a wide range of angular scales from arcminutes to degrees. The measured CIB brightness is consistent with previous results reported from COBE data but significantly higher than the lower limits at 70 and 160 um obtained with the Spitzer satellite from the stacking analysis of 24-um selected sources. The discrepancy with the Spitzer result is possibly due to a new galaxy population at high redshift obscured by hot dust. From power spectrum analysis at 90 um, three components are identified: shot noise due to individual galaxies; Galactic cirrus emission dominating at the largest angular scales of a few degrees; and an additional component at an intermediate angular scale of 10-30 arcminutes, possibly due to galaxy clustering. The spectral shape of the clustering component at 90 um is very similar to that at longer wavelengths as observed by Spitzer and BLAST. Moreover, the color of the fluctuations indicates that the clustering component is as red as Ultra-luminous infrared galaxies (ULIRGs) at high redshift, These galaxies are not likely to be the majority of the CIB emission at 90 um, but responsible for the clustering component. Our results provide new constraints on the evolution and clustering properties of distant infrared galaxies.

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.

We present an 880 micron Submillimeter Array (SMA) detection of the submillimeter galaxy SXDF850.6. SXDF850.6 is a bright source (S(850 micron) = 8 mJy) detected in the SCUBA Half Degree Extragalactic Survey (SHADES), and has multiple possible radio counterparts in its deep radio image obtained at the VLA. Our new SMA detection finds that the submm emission coincides with the brightest radio emission that is found ~8" north of the coordinates determined from SCUBA. Despite the lack of detectable counterparts in deep UV/optical images, we find a source at the SMA position in near-infrared and longer wavelength images. We perform SED model fits to UV-optical-IR photometry (u, B, V, R, i', z', J, H, K, 3.6 micron, 4.5 micron, 5.8 micron, and 8.0 micron) and to submm-radio photometry (850 micron, 880 micron, 1100 micron, and 21 cm) independently, and we find both are well described by starburst templates at a redshift of z ~= 2.2 (+/- 0.3). The best-fit parameters from the UV-optical-IR SED fit are a redshift of z = 1.87 (+0.15/-0.07), a stellar mass of M_star = 2.5 +2.2/-0.3 x 10^11 M_sun, an extinction of A_V = 3.0 (+0.3/-1.0) mag, and an age of 720 (+1880/-210) Myr. The submm-radio SED fit provides a consistent redshift of z ~ 1.8-2.5, an IR luminosity of L_IR = (7-26) x 10^12 L_sun, and a star formation rate of 1300-4500 M_sun/yr. These results suggest that SXDF850.6 is a mature system already having a massive amount of old stellar population constructed before its submm bright phase and is experiencing a dusty starburst, possibly induced by major mergers.

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.

There is a huge gap between properties of red-sequence selected massive galaxy clusters at z<1 and Lyman-break selected proto-clusters at z>3. It is important to understand when and how the z>3 proto-clusters evolve into passive clusters at z<1. We aim to fill this cluster desert by using the space-based N4(4um) imaging with the AKARI. The z'-N4 color is a powerful separator of cluster galaxies at z>1, taking advantage of the 4000A break and the 1.6um bump. We carefully selected 16 promising cluster candidates at 0.9<z<1.7, which all show obvious over-density of galaxies and a prominent red-sequence. At this redshift range, the mid-infrared S15um/S9um flux ratio is an extinction-free indicator of galaxy star formation activity due to the redshifted PAH emission lines (6.2,7.7 and 8.6um). We show statistically that the cluster galaxies have a lower S15um/S9um flux ratio than field galaxies, i.e., cluster galaxies already have lower star-formation activity at 0.9<z<1.7, pushing the formation epoch of these galaxy clusters to a higher redshift.

We investigate the relationship between two massive star-forming galaxy populations at redshift z~ 2; i.e. submillimetre galaxies (SMGs) and BzK-selected galaxies (BzKs). Out of 60 SMGs found in the Subaru/XMM-Newton deep field, we collect optical--NIR photometry of 28 radio counterparts for 24 SMGs, based on refined sky positions with a radio map for 35 SMGs (Ivison et al. 2007). We find a correlation between their K-band magnitudes and BzK [= (z-K)-(B-z)] colours: almost all of the K-faint (K_AB > 21.3) radio-detected SMGs have BzK>-0.2, and therefore BzKs. This result gives strong support to perform direct optical identification of SMGs by searching for BzKs around SMGs. We calculate the formal significance (P' value) for each of the BzK associations around radio-undetected SMGs, and find 6 new robust identifications, including one double identification. From this analysis, we obtain the current best estimate on the surface density of BzK-selected SMGs, which indicates that only ~1 per cent of BzKs are SMGs. If BzKs are normal disk-like galaxies at z~ 2 as indicated by the correlation between their star formation rate (SFR) and stellar mass and also by dynamical properties, SMGs are likely to be merging BzKs. In this case, a typical enhancement of SFR due to merging is only a factor of ~ 3, which is an order of magnitude lower than that of local ULIRGs. This may indicate that most of the merging BzKs could be observed as SMGs. Considering a possible high fraction of mergers at z~ 2 (at least it would be higher than the fraction at z> 1 of ~ 10 per cent), it is rather puzzling to find such a low fraction of SMGs in the progenitor population, i.e. BzKs.

We present the number counts of Ks-band selected high redshift galaxy populations such as extremely red objects (EROs), B, z & K -band selected galaxies (BzKs) and distant red galaxies (DRGs) in the AKARI NEP field. The final catalogue contains 308 EROs (Ks<19.0 ; 54 percent are dusty star-forming EROs and the rest are passive old EROs), 137 star-forming BzKs and 38 passive old BzKs (Ks<19.0) and 64 DRGs (Ks<18.6). We also produce individual component source counts for both the dusty star-forming and passive populations. We compare the observed number counts of the high redshift passively evolving galaxy population with a backward pure luminosity evolution (PLE) model allowing different degrees of number density evolution. We find that the PLE model without density evolution fails to explain the observed counts at faint magnitudes, while the model incorporating negative density evolution is consistent with the observed counts of the passively evolving population. We also compare our observed counts of dusty star-forming EROs with a phenomenological evolutionary model postulating that the near-infrared EROs can be explained by the source densities of the far-infrared - submillimetre populations. Our model predicts that the dusty ERO source counts can be explained assuming a 25 percent contribution of submillimetre star-forming galaxies with the majority of brighter Ks -band detected dusty EROs having luminous (rather than HR10 type ultra-luminous) submillimetre counterparts. We propose that the fainter Ks>19.5 population is dominated by the sub-millijansky submillimetre population. We also predict a turnover in in dusty ERO counts around 19<Ks<20.

We present the results of optical identifications for 257 mid-infrared sources detected with a deep 15um survey over approximately 80 arcmin^2 area in the AKARI performance verification field near the North Ecliptic Pole. The 15um fluxes of the sources range from 1 mJy down to 40 uJy, approximately a half of which are below 100 uJy. Optical counterparts were searched for within a 2-3 arcsec radius in both the BVRi'z' catalog generated by using the deep Subaru/Suprime-cam field which covers one-third of the performance verification field, and the g'r'i'z' catalog based on observations made with MegaCam at CFHT. We found B-R and R-z' colours of sources with successful optical identifications are systematically redder than that of the entire optical sample in the same field. Moreover, approximately 40% of the 15um sources show colours R-L15>5, which cannot be explained by the spectral energy distribution (SED) of normal quiescent spiral galaxies, but are consistent with SEDs of redshifted (z>1) starburst or ultraluminous infrared galaxies. This result indicates that the fraction of the ultraluminous infrared galaxies in our faint 15um sample is much larger than that in our brighter 15um sources, which is consistent with the evolving mid-infrared luminosity function derived by recent studies based on the Spitzer 24um deep surveys. Based on an SED fitting technique, the nature of the faint 15um sources is further discussed for a selected number of sources with available K_s-band data.

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 have studied the submillimetre (submm) properties of the following classes of near-infrared (NIR)-selected massive galaxies at high redshifts: BzK-selected star-forming galaxies (BzKs); distant red galaxies (DRGs); and extremely red objects (EROs). We used the SCUBA HAlf Degree Extragalactic Survey (SHADES), the largest uniform submm survey to date. Partial overlap of SIRIUS/NIR images and SHADES in SXDF has allowed us to identify 4 submm-bright NIR-selected galaxies, which are detected in the mid-infrared, 24 micron, and the radio, 1.4 GHz. We find that all of our submm-bright NIR-selected galaxies satisfy the BzK selection criteria, except for one galaxy whose B-z and z-K colours are however close to the BzK colour boundary. Two of the submm-bright NIR-selected galaxies satisfy all of the selection criteria we considered, i.e. they belong to the BzK-DRG-ERO overlapping population, or `extremely red' BzKs. Although these extremely red BzKs are rare (0.25 arcmin^{-2}), up to 20 % of this population could be submm galaxies. This fraction is significantly higher than that found for other galaxy populations studied here. Via a stacking analysis, we have detected the 850 micron flux of submm-faint BzKs and EROs in our SCUBA maps. While the contribution of z~2 BzKs to the submm background is about 10--15 % and similar to that from EROs typically at z~1, BzKs have a higher fraction (~30 %) of submm flux in resolved sources compared with EROs and submm sources as a whole. From the SED fitting analysis for both submm-bright and submm-faint BzKs, we found no clear signature that submm-bright BzKs are experiencing a specifically luminous evolutionary phase, compared with submm-faint BzKs. An alternative explanation might be that submm-bright BzKs are more massive than submm-faint ones.

We present an initial analysis of AKARI 18um-selected galaxies using all 9 photometric bands at 2-24um available in the InfraRed Camera (IRC), in order to demonstrate new capabilities of AKARI cosmological surveys. We detected 72 sources at 18um in an area of 50.2 arcmin^2 in the AKARI/IRC monitor field towards the North Ecliptic Pole (NEP). From this sample, 25 galaxies with probable redshifts z>~ 0.5 are selected with a single colour cut (N2-N3>0.1) for a detailed SED analysis with ground-based BVRi'z'JK data. Using an SED radiative transfer model of starbursts covering the wavelength range UV -- submm, we derive photometric redshifts from the optical-MIR SEDs of 18um-selected galaxies. From the best-fit SED models, we show that the IRC all-band photometry is capable of tracing the steep rise in flux at the blue side of the PAH 6.2um emission feature. This indicates that the IRC all-band photometry is useful to constrain the redshift of infrared galaxies, specifically for dusty galaxies with a less prominent 4000A break. Also, we find that the flux dip between the PAH 7.7 and 11.2um emission feature is recognizable in the observed SEDs of galaxies at z~1. By using such a colour anomaly due to the PAH and silicate absorption features, unique samples of ULIRGs at z~1, `silicate-break' galaxies, can be constructed from large cosmological surveys of AKARI towards the NEP, i.e. the NEP-Deep and NEP-Wide survey. This pilot study suggests the possibility of detecting many interesting galaxy properties in the NEP-Deep and Wide surveys, such as a systematic difference in SEDs between high- and low-z ULIRGs, and a large variation of the PAH inter-band strength ratio in galaxies at high redshifts. [abridged]

We present a five (u*,g',r',i',z') band optical photometry catalog of the sources in the North Ecliptic Pole (NEP) region based on deep observations made with MegaCam at CFHT. The source catalog covers about 2 square degree area centered at the NEP and reaches depths of about 26 mag for u*, g', r' bands, about 25 mag for i' band, and about 24 mag for z' band (4 sigma detection over an 1 arcsec aperture). The total number of cataloged sources brighter than r'= 23 mag is about 56,000 including both point sources and extended sources. From the investigation of photometric properties using the color-magnitude diagrams and color-color diagrams, we have found that the colors of extended sources are mostly (u*-r') < 3.0 and (g'-z') > 0.5. This can be used to separate the extended sources from the point sources reliably, even for the faint source domain where typical morphological classification schemes hardly work efficiently. We have derived an empirical color-redshift relation of the red sequence galaxies using the Sloan Digital Sky Survey data. By applying this relation to our photometry catalog and searching for any spatial overdensities, we have found two galaxy clusters and one nearby galaxy group.

We present the J- and Ks-band galaxy counts and galaxy colors covering 750 square arcminutes in the deep AKARI North Ecliptic Pole (NEP) field, using the FLoridA Multi-object Imaging Near-ir Grism Observational Spectrometer (FLAMINGOS) on the Kitt Peak National Observatory (KPNO) 2.1m telescope. The limiting magnitudes with a signal-to-noise ratio of three in the deepest regions are 21.85 and 20.15 in the J- and Ks-bands respectively in the Vega magnitude system. The J- and Ks-band galaxy counts in the AKARI NEP field are broadly in good agreement with those of other results in the literature, however we find some indication of a change in the galaxy number count slope at J~19.5 and over the magnitude range 18.0 < Ks < 19.5. We interpret this feature as a change in the dominant population at these magnitudes because we also find an associated change in the B - Ks color distribution at these magnitudes where the number of blue samples in the magnitude range 18.5 < Ks < 19.5 is significantly larger than that of Ks < 17.5.

AKARI (formerly ASTRO-F) is an infrared space telescope designed for an all-sky survey at 10-180 (mu)m, and deep pointed surveys of selected areas at 2-180 (mu)m. The deep pointed surveys with AKARI will significantly advance our understanding of galaxy evolution, the structure formation of the Universe, the nature of the buried AGNs, and the cosmic infrared background. Here we describe the important characteristics of the AKARI mission: the orbit, and the attitude control system, and investigate the optimum survey area based on the updated pre-flight sensitivities of AKARI, taking into account the cirrus confusion noise as well as the surface density of bright stars. The North Ecliptic Pole (NEP) is concluded to be the best area for 2-26 (mu)m deep surveys, while the low-cirrus noise regions around the South Ecliptic Pole (SEP) are worth considering for 50-180 (mu)m pointed surveys to high sensitivities limited by the galaxy confusion noise. Current observational plans of these pointed surveys are described in detail. Comparing these surveys with the deep surveys with the Spitzer Space Telescope, the AKARI deep surveys are particularly unique in respect of their continuous wavelength coverage over the 2-26 (mu)m range in broad-band deep imaging, and their slitless spectroscopy mode over the same wavelength range.

We analyse the spectral energy distributions (SEDs) of 23 submm galaxies and 3 ISO-detected EROs, all of which have the spectroscopic redshifts, by using an evolutionary SED model of starbursts. This SED model allows us to investigate intrinsic properties of starbursts, such as the starburst age and the mean stellar metallicity, as it takes into account the chemical evolution. Also, the intrinsic size of the starburst region is estimated from observed SEDs. Using this SED model, we predict colour, magnitude, and size of present-day descendants of submm galaxies, and derive scaling relations, such as the present-day colour-magnitude and size-magnitude relations. We argue that submm galaxies are the progenitors of present-day elliptical galaxies, provided that the initial mass function (IMF) of submm galaxies is slightly flatter than the Salpeter IMF. In this case, we find that 1) the mean present-day magnitude of submm galaxies is similar to that of L* elliptical galaxies, 2) the present-day colour-magnitude relation is consistent with that of elliptical galaxies, 3) the present-day size-magnitude relation of elliptical galaxies can be reproduced if massive submm galaxies consist of multiple starburst regions. We estimate the effect of feedback in submm galaxies as a function of starburst age. It is found that starburst regions in submm galaxies are likely to be self-regulated; i.e. the effect of feedback is nearly balanced with the self-gravity of starburst regions.

We construct an evolutionary spectral energy distribution (SED) model of a starburst region, from the ultraviolet to submillimetre wavelengths. This model allows us to derive the star formation rate, optical depth by dust and apparent effective radius of starburst regions at various wavelengths; as a result, the intrinsic surface brightness of starburst regions can be derived. Using this SED model, we analyse 16 UV-selected starburst galaxies and 10 ultraluminous infrared galaxies. The derived star formation rates and optical depths are compared with emission line measurements and found to be consistent. The derived apparent effective radii are also consistent with observations. From the SED analysis, we find a bimodal property of the star formation rate with the optical depth and the compactness of stellar distributions. While mild starbursts have a limiting intrinsic surface brightness L(bol) x r_e^{-2} ~ 10^{12} Lsun kpc^{-2}, intense starbursts tend to be more heavily obscured and concentrated within a characteristic scale of r_e ~ 0.3 kpc. We suggest that the mild starbursts can be triggered by a self-gravitating disc instability in which feedback is effective in the shallow gravitational potential. On the other hand, the intense starbursts can be induced via an external dynamical perturbation like galaxy merging, in which feedback is less effective due to the deep gravitational potential attained by the large gas concentration within the central starburst region.

In order to study the statistics of the objects with hierarchical merging, we propose the skeleton tree formalism, which can analytically distinguish the episodic merging and the continuous accretion in the mass growth processes. The distinction was not clear in extended Press-Schechter (PS) formalism. The skeleton tree formalism is a natural extension of the peak theory which is an alternative formalism for the statistics of the bound objects. The fluctuation field smoothing with Gaussian filter produces the landscape with adding the extra-dimension of the filter resolution scale to the spatial coordinate of the original fluctuation. In the landscape, some smoothing peaks are nesting into the neighboring peaks at a type of critical points called sloping saddles appears, which can be interpreted as merging events of the objects in the context of the hierarchical structure formation. The topological properties of the landscape can be abstracted in skeleton trees, which consist of line process of the smoothing peaks and the point process of the sloping saddles. According to this abstract topological picture, in this paper, we present the concept and the basic results of the skeleton tree formalism to describe (1) the distinction between the accretion and the merger in the hierarchical structure formation from various initial random Gaussian fields; (2) the instantaneous number density of the sloping saddles which gives the instantaneous scale function of the objects with the destruction and reformation in the mergers; (3) the rates of the destruction, the reformation, and the relative accretion growth; (4) the self-consistency of the formalism for the statistics of the mass growth processes of the objects; (5) the mean growth history of the objects at the fixed mass.