{'fraction_non_alphanumeric': 0.04169955632978579, 'fraction_numerical': 0.02534091837012439, 'mean_word_length': 4.318471337579618, 'pattern_counts': {'":': 0, '<': 18, '<?xml version=': 0, '>': 13, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Polycyclic aromatic hydrocarbon (PAH) luminous galaxies at z~1

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.


Introduction
In recent studies of galaxy formation and evolution, the most massive and luminous galaxies have played a leading role. Most luminous galaxies, such as submillimetre galaxies and ultraluminous infrared galaxies (ULIRGs 1 ) at high redshifts are now believed to be the progenitors of massive spheroids, and highlight the early formation of massive galaxies (e.g. Smail et al. 2004;Chapman et al. 2005;Tacconi et al. 2008). This supports the scenario of anti-hierarchical or 'down-sizing' galaxy formation (e.g. Heavens et al. 2004;Cimatti et al. 2006).
Mid-infrared (IR) surveys are a powerful tool for studying the most luminous galaxies at high redshifts, by providing the infrared luminosity function and the star formation rate in a large cosmic volume as a function of redshift (e.g Pérez-González et al. 2005;Le Floc'h et al. 2005;Babbedge et al. 2006;Caputi et al. 2007;Goto et al. 2010). However, this work can suffer uncertainty from the k-correction, 1 Galaxies with a total infrared (8 -1000 µm) luminosity of 10 12−13 L ⊙ . which requires knowing the appropriate mid-IR spectral energy distribution (SED). In fact, sensitive IRS observations with the Spitzer Space Telescope reveal that the mid-IR spectra of submillimetre galaxies resembles that of over 2 orders of magnitude less luminous starburst galaxies (Valiante et al. 2007;Pope et al. 2008;Farrah et al. 2008), which have prominent emission features of polycyclic aromatic hydrocarbons (PAHs). This raises questions about the nature of most luminous galaxies, and how they differ systematically from low-z to high-z.
Multi-wavelength mid-IR surveys with AKARI provide a unique measure of the mid-IR properties of infrared galaxies at high redshifts with a statistically valid sample. Takagi et al. (2007a) demonstrate that AKARI/IRC all-band photometry is capable of identifying the approximate spectral shape of the PAH emission, specifically the steep rise of flux at the blue side of the PAH 6.2 µm feature. Such a steep flux rise can only be accounted for by PAH emission. Using this fact, we propose a new photometric selection method for galaxies whose mid-IR emission is dominated by PAH emission. This is possible if photometric observations are sufficiently comprehensive enough to trace the PAH emission feature in the SED, as in the AKARI/IRC all-band photometric survey. We call galaxies selected in this method PAH-selected galaxies.
In this paper, we utilize the multi-wavelength mid-IR coverage of the NEP-Deep survey with AKARI, in order to identify galaxies with strong PAH emission. Sections 2 and 3 describe the data and sample selection, respectively. Our SEDfitting method using the optical to mid-IR wavelength range, and its results are presented in section 4. We discuss the properties of PAH-selected galaxies in section 5 and summarise the results in section 6. Throughout this paper, a flat cosmology with H 0 = 70 km s −1 Mpc −1 and Ω Λ = 0.7 is used. All magnitudes are given in the AB system, unless otherwise stated.
We performed an optical spectroscopic survey of 242 mid-IR sources with R < ∼ 24 mag in the S-cam field using Keck/DEIMOS (Takagi et al. in prep). Here we use the spectroscopic redshifts for the calibration of photometric redshifts. The optical emission line diagnostics of mid-IR sources will be given elsewhere.

All-band-detected sources
In order to evaluate and to maximize the new diagnostic capabilities becoming available from the NEP-Deep survey, we have constructed a source catalogue of objects, detected in all nine IRC bands (hereafter all-band-detected sources). We have combined the S 7, S 9W, S 11-band images for source detection and used SExtractor (Bertin & Arnouts 1996) to make an initial catalogue. We then follow Takagi et al. (2007a) for photometry with IRC images, in which the aperture photometry with 2 pixel and 3 pixel radii are used for near-IR (N2, N3, N4) and mid-IR images (S 7, S 9W, S 11, L15, L18W, L24), respectively. Appropriate aperture corrections are then applied. We also follow Takagi et al. (2007a) for band merging. Starting from the centroid position of S 11 sources, we searched for the centroid in the other images. If the angular distance between centroids 1 Corrected from the values in Wada et al. (2008), in which the adopted conversion factors for N2, N3, and N4 bands are not correct. The correct sensitivity is better by a factor of 1.47. A colour-colour plot of all-band-detected sources with 15-to-9 µm and 11-to-7 µm flux ratios. Dashed lines indicate the colour cuts for selecting PAH-selected galaxies. Points with error bars indicate the PAH-selected galaxies.  2. 11-to-7 µm flux ratio as a function of redshift. Solid circles with error bars indicates the spectroscopic sample with S 7 , S 9 , S 11 > 30 µJy and S 15 > 100 µJy. a) Flux ratio estimated from the SED template of Polletta et al. (2007). The names of galaxies in the template are given in the legend. b) Same as a), but with the SBURT model with different starburst age and compactness of the starburst region (Θ=1.4, 1.6, 1.8, 2.0) , which controls the optical depth. The MW dust model is adopted.  Figure 2b), but for 15-to-9 µm flux ratio. is less than 3 ′′ , the objects are considered to be the same (see Wada et al. 2008). Since we have not adopted a specific signalto-noise ratio for the detection at all of the IRC bands other than S 11, the resulting catalogue includes less significant detections (< 3 σ) in some cases. Nevertheless, we include these sources in our all-band-detected sample, since the less significant detections are rare, and the resulting SED is still well constrained.
There are 1100 all-band-detected sources in the entire NEP-Deep field, and 630 in the S-cam field. We searched for optical identification using sky positions in the N2 band, where the astrometry is based on the Two Micron All Sky Survey (2MASS) and therefore most reliable. The resulting optical identifications were all visually inspected. We found that about 10 % of sources have ambiguous optical identifications. Excluding these sources, we obtained 568 all-band-detected sources in the S-cam field with unambiguous optical identification, of which 113 have optical spectra with DEIMOS.

PAH-selected galaxies
At the mid-IR wavelength region, emission from galaxies is usually dominated by small and hot dust grains experiencing the temperature fluctuation. PAHs are believed to be the smallest of those, consisting of ∼ 100 carbon atoms grouped into ∼10 Å in dimensions, and having emission bands. PAH emission features have sometimes been used as star formation indicators, which are prominent in starburst galaxies (Brandl et al. 2006), but weak or absent in AGNs (Weedman et al. 2005).
Using the very steep flux rise on the blue side of the 6.2 µm PAH feature, we attempt to identify galaxies whose mid-IR emission is dominated by PAHs, based on AKARI photometric observations. The steep flux rise due to PAH emission can be characterized by the flux ratio at the rest-frame wavelengths of approximately 7 to 4 µm. For galaxies at z ∼ 0.5 and 1, we could use 11-to-7 µm flux ratio and 15-to-9 µm flux ratio, respectively, in order to identify galaxies with prominent PAH emission. For convenience, we hereafter use the term, the PAH-to-continuum flux ratio, for these flux ratios. Figure 1 shows a colour-colour plot of all-band-detected sources using these flux ratios. There are galaxies with the PAH-to-continuum flux ratio as high as 10 for both ratios. We hereafter call galaxies with the flux ratio of > 8 as PAH-selected galaxies.
For star-forming galaxies, it is likely that the flux at the restframe 4 µm is dominated by the stellar continuum, while PAH 7.7 µm feature dominate the flux at the rest-frame 7 µm. The PAH luminosity is believed to be a good tracer of total infrared luminosity of star-forming galaxies, and hence the SFR (e.g. Genzel et al. 1998;Rigopoulou et al. 1999;Farrah et al. 2007). Thus, these flux ratios are roughly proportional to the massnormalized SFRs, i.e. specific SFR. We hence expect that PAHselected galaxies are the best candidates of starburst galaxies.
In Figure 2a, we show the 11-to-7 µm flux ratio as a function of redshift, using the SED template of various galaxy typesthe SWIRE template library (Polletta et al. 2007). This template confirms that the 11-to-7 µm flux ratio has a peak at z ∼ 0.5, although no SEDs in the template can reproduce the flux ratios greater than 8. We note that the SWIRE template library has only a few SEDs of starburst-dominated galaxies. In order to explain the high flux ratio of PAH-selected galaxies, we may need more comprehensive SED templates of starburst galaxies.
For such SED template of starburst galaxies, we adopt an evolutionary starburst SED model with radiative transfer (Takagi et al. 2003a, SBURT). This is an SED model for starburst regions, assumed to be spherical and have a concentrated stellar distribution and homogeneous dust distribution. SBURT adopts a stellar population synthesis model of Kodama & Arimoto (1997) as input. The amount of dust is calculated, based on the chemical evolution assuming gas infall. The time-scale of gas infall and of star formation is assumed to be 0.1 Gyr 2 . Thus, the SBURT model can be regarded as a population synthesis model with physically consistent treatment of dust extinction and re-emission. The optical depth of the system is controlled by the compactness (Θ) of the starburst region. In the SBURT model, there are three choices of the extinction curve, corresponding to the dust model for the Milky Way (MW), Large or Small Magellanic Clouds (LMC or SMC). The fraction of PAHs in dust is 5, 1 and 0.5 % in the MW, LMC, and SMC dust model, respectively. Thus, models with the MW dust has the largest PAH-to-continuum flux ratio. See Takagi et al. (2003a,b) for more details of SBURT.
In Figure 2b and 3, we show the 11-to-7 µm and 15-to-9 µm flux ratios estimated from the SBURT model, respectively. We found that the MW dust model is required, in order to reproduce high flux ratios of PAH-selected galaxies. According to the model, younger galaxies may have larger flux ratios. This is because the contribution from the stellar component to the 4 µm continuum is less for younger starbursts. The flux ratios of greater than 8 are satisfied with the model younger than 0.4 Gyr. However, this is only one possible interpretation, since the dust model adopted has only three choices (MW, LMC or SMC), and limited capability of reproducing large variations of mid-IR properties. Specifically, this model does not take the ionization status of PAHs into account, which significantly affects the optical constants of PAHs (e.g. Li & Draine 2001).
Among all-band-detected sources, we identify 56 (39) PAHselected galaxies from 11-to-7 µm (15-to-9 µm) flux ratio in the entire NEP-Deep field, of which 38 (18) lie in the S-cam field with unambiguous optical identification. We include sources with less significant detection at the rest-frame 4 µm, if the 2.5 σ upper limit is consistent with the selection criterion. In Table  A.1, we tabulate the IRC photometry of PAH-selected galaxies.

Method
We fit the optical-to-MIR SEDs of all-band-detected sources with the SBURT model. We use the model with the starburst age of 0.01 -0.6 Gyr, which is enough to represent the starburst phase. The SBURT model originally covers the wavelength range from UV to submillimetre. Takagi et al. (2007b) empirically added the radio component to reproduce the observed radio-IR correlation. This model does not include any emission from AGN. Thus, it is expected that the SBURT model underpredicts the mid-IR fluxes from AGN hosts, which usually have hot dust excess compared to pure starburst galaxies (Hunt et al. 1997). We adopt the same SED fitting method as in Takagi et al. (2007a), using a standard χ 2 minimization technique. Takagi et al. (2007a) adopted relatively large minimum photometric errors, 20 -30 % depending on the band. Since the photometric calibration of IRC is recently improved (Tanabé et al. 2008), we here set the minimum photometric errors to 15 % for all IRC bands. As in Takagi et al. (2007a), we quadratically added an additional 20 % error for data at the rest-frame UV wavelengths, < 4000 Å, in order to account for the uncertainty of the extinction curve.
We apply the SBURT model for all-band-detected sources. By nature, all-band-detected sources are heterogeneous sample, even including stars. Since we use a starburst SED template, SBURT, for this sample, it is important to discriminate starburst candidates from other type of sources. To do this, we use the goodness of SBURT fitting and reject sources with no acceptable SBURT model in further analyses. We reject the best-fitting SBURT model if the value of resulting χ 2 is so large that the corresponding probability is less than 1 %.

All-band-detected sources
We performed SED-fitting for 568 all-band-detected sources using the SBURT model. We obtained good fits for ∼ 40 % (211/568) of all-band-detected sources. On the other hand, using the SWIRE template library, instead of SBURT, we obtained good fits for only ∼10 % (69/568) of sources, although it covers a wide range of SED type. Note that the SWIRE library has only a few SEDs of starbursts. We think that a large variety of starburst SEDs, as in the SBURT model, is necessary to reproduce a good fraction of observed UV-mid-IR SEDs of IR-selected galaxies.
In order to show what kind of sources are rejected in the SBURT fitting, we plot both accepted and rejected samples in a colour-colour plot of N3 − S 7 vs N2 − N3 in Figure 4. Takagi et al. (2007a) use this colour-colour plot to discriminate AGNs from normal star-forming galaxies. AGNs identified with optical spectroscopy have red colours in both N2 − N3 and N3 − S 7. The accepted fits avoid this AGN colour region, because the SBURT model does not have any AGN component. The accepted sample also avoids the clump of objects with blue colours at N2 − N3 ≃ 0.6 and N3 − S 7 ≃ −1.5, which are stars. The rest of the rejected sample have similar N2−N3 and N3−S 7 colours to those of the accepted sample. These include quiescent spiral galaxies, which have systematically redder optical-near-IR colours than any of the SBURT models. As we discuss in Section 5.1, limitations on the dust model may also cause rejection of the best-fitting model.
In order to derive any firm conclusions, photometric redshifts of the good-fit sample should be accurate enough. Figure  5 shows the comparison of photometric redshifts to spectroscopic redshifts. Solid circles and crosses represent galaxies with acceptable and rejected fits, respectively. Defining a catastrophic error by ∆z (1+z) > 0.2, we find no catastrophic errors and σ ∆z (1+z) = 0.034 for the good-fit sample, and 21/68 catastrophic errors for the bad-fit sample. Thus, the goodness of the SBURT fit is actually a good measure of the accuracy of photometric redshifts.
In Figure 5, we also plot photometric redshifts derived from ground-based optical-NIR (BVRi ′ z ′ JK) bands only, using a widely-used photometric redshift code, Hyperz (Bolzonella et al. 2000). In Hyperz, we adopt a population synthesis model of Bruzual & Charlot (2003), which is distributed along with the code, and the extinction curve of starbursts (Calzetti et al. 2000) for dust reddening. In order to compare the photometric redshifts with AKARI bands to those without AKARI, we show the good-fit sample (with SBURT fitting) only in Figure 5. We find one catastrophic error in ground-based photometric redshifts, while this galaxy has a reasonable photometric redshift with AKARI bands. Furthermore, ground-based photometric redshifts seem to have systematic errors; i.e. at z ∼ 0.5, most of photometric redshifts are lower than the spectroscopic redshifts. Such errors can be removed by using AKARI bands.
In order to check the statistical reliability of our photometric redshifts, we use a correlation between the NIR colour N2 − N3 and redshift (Takagi et al. 2007a). With this test, we could detect systematic effects, if any, due to the usage of mid-IR fluxes for deriving photometric redshifts. In Figure 6, we show the correlations for both spectroscopic and photometric samples. For the photometric sample, we obtain the correlation of (N2 − N3) = (2.37 ± 0.12) log(1 + z) − 0.56, which is consistent with that of the spectroscopic sample, i.e. (N2 − N3) = (2.43 ± 0.05) log(1 + z) − 0.59. Therefore, we find no systematic errors in our photometric redshifts.

PAH-selected galaxies
We here show more details of SED fitting results specifically for PAH-selected galaxies. In Figure A.1 -A.4, we show the best-fitting (both accepted and rejected) SED models for PAHselected galaxies at z ∼ 0.5 and 1. We obtained acceptable SED fits for 22/38 and 10/18 of the PAH-selected galaxies at z ∼ 0.5 and 1, respectively. Despite the fact that PAH-selected galaxies are starbursts, the SBURT model does not provide acceptable fits for 40 % of PAH-selected galaxies, although the success rate of the fitting is better than that of whole sample. This topic is discussed in detail in the next section.
We show the redshift distributions of PAH-selected galaxies in Figure 7. They are consistent with the expectation from the model templates, shown in Figure 2 and 3. There are no systematic differences in photometric redshifts between the good-and bad-fit samples. This may indicate that the resulting large χ 2 is not due to an incorrect redshift. The total-IR luminosity from the SED model is indicated in each panel of Figure A.1 -A.4. We note that PAH-selected galaxies at z ∼ 1 include ULIRGs. This type of ULIRGs can be found only at high redshifts (see Section 5.2).

Fig. 7.
Redshift distribution of PAH-selected galaxies. Vertically and diagonally shaded histograms represent 11-to-7 µm and 15-to-9 µm selected galaxies with good SED fits, respectively. Histograms with thin dashed and solid lines indicate the redshift distribution including the bad-fit sample. Spectroscopic redshifis are adopted if available. All-band-detected sources and PAH-selected galaxies are represented by small dots and open circles, respectively. We show galaxies with good SED fits only.

Discussion
5.1. Nature of PAH-selected galaxies In Figure 8, we show 8 µm luminosities and stellar masses of PAH-selected galaxies, along with those of all-band-detected sources. The 8 µm luminosity is derived from the observed flux in the IRC bands closest to the reset-frame 8 µm and the luminosity distance from photometric redshifts. Some of PAHselected galaxies at z ∼ 1 have L 8 µm > 10 11 L ⊙ , which corre- Fig. 9. Ratio of observed to model fluxes as a function of restframe wavelengths for PAH-selected galaxies at 0.4 < z < 0.5. Open and solid circles indicate individual galaxies and the weighted average for each photometric band, assuming the mean redshift of 0.46 and 0.45 for good and bad fit sample, respectively.
spond to the luminosity of ULIRGs for a typical starburst SED (Weedman & Houck 2008). Stellar masses are derived from the best-fitting SED model when it is accepted. The adopted initial mass function (IMF) is a top-heavy one with the power law index of x = 1.10 (the Salpeter IMF has x = 1.35), and the lower and upper mass limits are 0.1 and 60 M ⊙ (see also Takagi et al. 2004). For the Salpeter IMF, stellar mass increases by a factor of 2. Stellar masses thus derived should be regarded as a lower limit, since the contribution from evolved stellar populations is not included. Most of PAHselected galaxies at z ∼ 1 have stellar masses of > 3 × 10 10 M ⊙ . Even some fraction of PAH-selected galaxies at z ∼ 0.5 have similar stellar masses. What are descendants of such massive PAH-selected, i.e. starburst galaxies? At redshift z ∼ 2, massive starburst galaxies have been found mainly as submillimetre galaxies (SMGs; e.g. Borys et al. 2005) and they are good candidates of massive spheroids Smail et al. 2004;Bouché et al. 2007;Tacconi et al. 2008). However, PAHselected galaxies are not as massive as SMGs. Their luminosity is typically an order of magnitude less than that of SMGs, corresponding to luminous infrared galaxies (LIRGs 3 ). Recent high resolution images of LIRGs at z ∼ 1 show that the majority of LIRGs have disk morphology (Elbaz et al. 2007;Melbourne et al. 2008), although the sample size is still small. Clearly, we need more comprehensive study to identify what 3 Galaxies with a total infrared (8 -1000 µm) luminosity of 10 11−12 L ⊙ . triggers starburst activity of PAH-selected galaxies, in order to constrain their evolutionary path.
What is the main cause of a large PAH-to-continuum flux ratio? As shown in Figure 2 and 3, the SBURT model predicts that galaxies with large PAH-to-continuum flux ratios are young starbursts with the starburst age of < ∼ 0.4 Gyr. For the whole sample of PAH-selected galaxies with good SED fits, we derive the mean starburst age of 0.4 Gyr. We note that, in the SED fitting, the starburst age depends not only on the PAH-to-continuum flux ratio, but also on the optical colours (Takagi et al. 1999). The lack of very young starbursts in the sample may be accounted for by a selection effect, i.e. such young systems are hard to be observed due to a short life time.
Indeed, star-forming regions, i.e. young stellar systems in a galaxy have a large PAH-to-continuum flux ratio, as shown in Spitzer/IRAC images of nearby galaxies in the literature. Wang et al. (2004) study the Antennae galaxies (NGC4038/4039) and show the map of the 8.0 4 -to-4.5 µm flux ratio, which is close to the PAH-to-continuum flux ratio in our definition. This flux ratio is quite high only at the overlap region of the two disks where the most of IR emission come from. They estimate that such active star-forming region comprise only ∼10 % of the total stellar mass of the entire system; i.e. the activity is localized rather than global. On the other hand, PAHselected galaxies should have global star-forming activity, since the flux ratio for the entire system is high.
For fair discussion, we should also pay attention to the case with no acceptable fits and clarify the cause of poor fits. In Figure 9, we show the observed-to-model flux ratio as a function of rest-frame wavelength for PAH-selected galaxies at 0.4 < z < 0.5. All of the rejected models, but one in Figure 9 have a starburst age of 0.6 Gyr, i.e. the oldest and reddest model in the template. This means that the model-fitting fails for galaxies with relatively red optical colours, but with strong PAH emission (see also Ohyama et al. in preparation). Since the SBURT model is a rather simple model of starbursts, which have a single starburst region and starburst stellar population only, it is not surprising to find its limitation. Figure 9 shows that a large discrepancy between observations and the model can be found at the 5 -7 µm range with sharp rising trend towards longer wavelength. We can see the same trend for the good-fit sample as well, although it is less significant. This indicates that the dust model adopted in the SBURT may need to be modified. One possibility is the ionization of PAHs. The optical constant of PAHs in the SBURT model is for neutral PAHs, which may not be true for every PAHs. Li & Draine (2001) suggest that ionized PAHs have systematically higher absorption cross section at 4 -10 µm, compared to neutral PAHs. If some fraction of PAHs are ionized, the discrepancy found in the SED fitting may be recovered.
Another possibility is that stars dominating optical light is not the same as those responsible for strong PAH emission. As starbursts age, dust heating would be mostly dominated by newly formed stars in molecular clouds, while most of optical light come from stars having already exited from molecular clouds; i.e. older stars are less attenuated and dominate optical light, while younger stars embedded in dense molecular clouds dominate infrared light. Such age-dependent dust attenuation would be more effective in the later stage of starbursts. This may explain why we find larger discrepancy in optically red, probably old, starbursts.
Above two possibilities could be discriminated with the far-IR photometry. The former, the case of ionized PAHs, would change only mid-IR luminosity and does not affect the far-IR luminosity. If the latter (age-dependent attenuation) is the case, the current model would underestimate not only mid-IR luminosity but also far-IR luminosity. Thus, expected PAH-to-total-IR luminosity depends on the scenario. Sensitive far-IR photometry of red PAH-selected galaxies, e.g. with the Herschel Space Observatory, could have diagnostic power to identify the main cause of the discrepancy.

PAH-to-total IR luminosity relation
The PAH luminosity of local ULIRGs is less than that expected from the correlation of PAH and total IR luminosity for less luminous starburst galaxies. In Figure 10, we show this luminosity relation for galaxies at z < 0.2 with crosses. The PAH luminosity is represented with the PAH 7.7 µm peak luminosity taken from Weedman & Houck (2008). The total infrared luminosity is taken from the references in Weedman & Houck (2008), and converted to the cosmology adopted here. The 7.7 µm peak luminosity of ULIRGs is systematically below the correlation for local starbursts. A most plausible explanation of this systematic effect is an AGN contribution to the total IR luminosity. Imanishi et al. (2007) claim that 30-50 % of optical non-Seyfert ULIRGs at z < 0.15 have dust-enshrouded AGN. Combined with optical Seyfert ULIRGs, their results indicate that ∼ 50 % of ULIRGs at z < 0.15 harbor AGN (see also Genzel et al. 1998;Lutz et al. 1998). However, Figure 10 shows that none of z < 0.2 ULIRGs are consistent with the pure starburst case. This may indicate another explanation for the lower 7.7 µm luminosity of ULIRGs. If the starburst region is heavily obscured, the absorption of PAH emissions by silicate dust may not be negligible.
High redshift ULIRGs may behave differently, owing to possible evolutionary effects. We investigate the PAH-to-total IR luminosity relation of the all-band-detected sources, including PAH-selected galaxies. The 7.7 µm peak luminosities are estimated from the IRC photometry by the following method. For a local starburst spectral template of Brandl et al. (2006), filterconvolved flux at the rest-frame 8 µm is approximately half of the 7.7 µm peak flux; 0.50 in the L15 filter at z = 1 and 0.42 for the S 11 filter at z = 0.5 for example. Here we adopt 0.50 for the conversion factor, since z ∼ 1 sources are our prime targets. This estimate would be accurate within a 30 % uncertainty, as far as the 8 µm emission is dominated by starbursts. We note that intrinsic spectral variation at 5 -8 µm is found to be small in starbursts (Nardini et al. 2008). The rest-frame 8 µm luminosities are estimated from the observed fluxes at the band closest to the rest-frame 8 µm and the luminosity distance derived from the redshift estimates. We have not done any interpolation or k-correction in this step, considering the uncertainty of the photometric redshift. The total IR luminosities are derived from the best-fitting SED model. Therefore, we consider the good-fit sample only. A well-known tight correlation between the global far-IR and radio emission allows us to check the reliability of estimated total-IR luminosity from the SED fitting. Among allband-detected sources, we detected five galaxies in our 1.4 GHz observation with WSRT (White et al. 2010 in prep), which covers a part of the NEP-Deep field. Comparing the observed radio fluxes to those predicted from the best-fitting SBURT models based on the far-IR/radio correlation, we estimate that the derived total-IR luminosities have uncertainties of 50 %. Since the derivations of both the 7.7 µm and total IR luminosity assume that starbursts dominate the infrared emission, the results for Fig. 11. Same as Figure 9 but for PAH-selected galaxies at 1.1 < z < 1.4. Open and solid circles indicate individual galaxies and the average for each photometric band, assuming a mean redshift of 1.22. Dashed lines bracket the 3.3 µm region.
PAH-selected galaxies, the most plausible starburst candidates in our sample, should be most reliable.
As shown in Figure 10, some PAH-selected galaxies at z ∼ 1 are ULIRGs. Their 7.7 µm peak luminosities are larger than local ULIRGs with a similar luminosity, and consistent with the PAHto-total IR luminosity ratio of less luminous starbursts. They have no local counterparts and we specifically call these PAHselected ULIRGs "PAH-luminous galaxies". Total IR luminosity of individual galaxies can be found in the legend of Figure  A.1 -A.4. In our sample, we find no PAH-luminous galaxies at z ∼ 0.5.
Normal starbursts, which have unabsorbed PAH emission features and no AGN contribution, can reach higher luminosity at higher redshifts. Weedman & Houck (2008) reach the same conclusion by compiling various Spitzer IRS spectroscopic samples. Rigby et al. (2008) investigate the rest-frame 8 µm luminosity of galaxies at z ∼ 2 as a function of total IR luminosity, and find a similar trend. Using Spitzer/IRS, Huang et al. (2009) also report ULIRGs with high PAH-to-total IR luminosity ratio at z ∼ 2. Our results show that there are extreme starbursts with ULIRG luminosity even at z ∼ 1.

Photometric detection of PAH 3.3 µm feature?
For galaxies at z ∼ 1, the PAH 3.3 µm feature falls into S 7 band. PAH-selected galaxies at z ∼ 1 are good candidates for detecting this feature, because of both redshift and PAH prominence. Magnelli et al. (2008) made a similar attempt using Spitzer IRAC photometry for galaxies at 0.6 < z < 0.9, and detected an excess over the stellar continuum, which may originate from the PAH 3.3 µm feature.
In Figure A.3, we can see that some objects, i.e. ID2711, 2836, 2932, 3033, 4956, have blue S 7 − S 9W colours, which might be caused by PAH 3.3 µm feature. Surprisingly, this excess is notable even if compared with the model flux which already includes the contribution from 3.3 µm feature. We show the comparison of observed fluxes to model fluxes in Figure 11. Here we use the sample of PAH-selected galaxies at 1.1 < z < 1.4 with good SED fits. The average observed-to-model flux ratios range from 0.75 to 1.23, depending on the wavelength. At PAHdominated wavelengths, i.e. > 7 µm, the scatter is large and the ratios are marginally consistent with the model. At the rest-frame 3 µm band, including PAH 3.3 µm feature, the flux ratio is systematically high by a factor of 1.23 on average. If this excess is Spitzer IRS local sample All-band detected starbursts PAH-selected galaxies Houck et al. (2007) galaxies at z>1 Fig. 10. 7.7 µm peak luminosity versus infrared luminosity. Solid and dashed lines represent the correlation between these luminosities found for less luminous starburst galaxies and its 1 σ uncertainty (Houck et al. 2007). Legend describes symbols for local galaxies from Spitzer IRS sample from Weedman & Houck (2008), all-band-detected sources, and PAH-selected galaxies at both z ∼ 0.5 and 1. Solid circles indicate galaxies with z > 1. Typical uncertainty of both luminosities is indicated at the bottom right corner.
accounted for by a systematically strong 3.3 µm feature only, this feature must be unusually strong by a factor of 3 -7, compared with the model. In the literature, no galaxies with such strong 3.3 µm feature can be found to our knowledge. This anomaly could originate from the continuum emission, rather than PAH 3.3 µm feature. However, we stress that, for galaxies at z ∼ 0.5, the SBURT model reproduces the rest-frame 1 -3 µm continuum emission well, as shown in Figure 9.
It may be possible that the model lacks additional very hot dust component which contribute to the rest frame 3-4 µm emission. Such a near-IR excess with a colour temperature of ∼10 3 K is reported in study of nearby galaxies (e.g. Lu et al. 2003;Flagey et al. 2006). However, we note that galaxies with large 3 µm flux tend to have small 4 µm flux as well, compared with the model. If the additional component has a gray body emission, it is difficult to explain such a trend. Furthermore, PAH-selected galaxies at z ∼ 0.5 show no systematic near-IR excess, compared to the model, as shown in Figure 9.
We need more careful investigation of the origin of the flux anomaly at this wavelength range. In fact, the model does not take into account the absorption features due to Hydrogenated Amorphous Carbon (HAC) and/or H 2 O. Unfortunately, we have to wait for the launch of JWST (Clampin 2008) and SPICA (Nakagawa 2004) to obtain mid-IR spectra of these galaxies and reveal the origin of this flux anomaly.

Summary
Using a multi-wavelength mid-IR survey, the AKARI NEP-Deep, we have constructed a catalogue of all-band-detected sources. From this catalogue, we photometrically identified galaxies whose mid-IR emission is dominated by PAH emission. These PAH-selected galaxies have large PAH-to-continuum flux ratio, i.e. 11-to-7 µm and 15-to-9 µm flux ratios at z ∼ 0.5 and 1, respectively. PAH-selected galaxies are the best candidates of starburst galaxies at these redshifts. Some of PAH-selected galaxies have stellar masses of > 3 × 10 10 M ⊙ , and could be progenitors of present-day massive galaxies, which are mostly spheroids.
An evolutionary SED model of starbursts, SBURT, reproduces observed optical-to-mid-IR SEDs of more than half of PAH-selected galaxies. According to the SED model, the PAHto-continuum flux ratio of > 8 can be explained with the starburst age of 0.4 Gyr or younger. The average starburst age from the SED fitting is 0.4 Gyr, which is marginally consistent with the PAH-to-continuum flux ratio. The lack of young starbursts may be the selection effect due to a short life time of young starbursts. On the other hand, SBURT has a difficulty to reproduce large PAH fluxes of optically red PAH-selected galaxies, which are probably evolved starbursts. This may require the SED model to include the age-dependent extinction, and/or to adopt improved optical properties of PAHs.
At z ∼ 1, the infrared luminosity of some PAH-selected galaxies corresponds to that of ULIRGs, and we call them PAHluminous galaxies. They have large PAH luminosity, compared with local ULIRGs, but the PAH-to-total IR luminosity ratio is comparable to that of less luminous starbursts. PAH-luminous galaxies seem to be a unique galaxy population at high redshifts. The number density of PAH-luminous galaxies would be given elsewhere, using whole area of the NEP-Deep survey.
There is a hint that the PAH 3.3 µm feature is quite strong in PAH-luminous galaxies by a factor of 3 -7, but this needs to be confirmed with next generation infrared telescopes, such as JWST and SPICA, with infrared spectroscopy.    (3) and (4)