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.

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.

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

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.

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 have observed the central region of the nearby starburst galaxy NGC 253 with the Kyoto Tridimensional Spectrograph II (Kyoto3DII) Fabry-Perot mode in order to investigate the properties of its galactic wind. Since this galaxy has a large inclination, it is easy to observe its galactic wind. We produced the Ha, [N II]6583, and [S II]6716,6731 images, as well as those line ratio maps. The [N II]/Ha ratio in the galactic wind region is larger than those in H II regions in the galactic disk. The [N II]/Ha ratio in the southeastern filament, a part of the galactic wind, is the largest and reaches about 1.5. These large [N II]/Ha ratios are explained by shock ionization/excitation. Using the [S II]/Ha ratio map, we spatially separate the galactic wind region from the starburst region. The kinetic energy of the galactic wind can be sufficiently supplied by supernovae in a starburst region in the galactic center. The shape of the galactic wind and the line ratio maps are non-axisymmetric about the galactic minor axis, which is also seen in M82. In the [N II]6583/[S II]6716,6731 map, the positions with large ratios coincide with the positions of star clusters found in the Hubble Space Telescope (HST) observation. This means that intense star formation causes strong nitrogen enrichment in these regions. Our unique data of the line ratio maps including [S II] lines have demonstrated their effectiveness for clearly distinguishing between shocked gas regions and starburst regions, determining the extent of galactic wind and its mass and kinetic energy, and discovering regions with enhanced nitrogen abundance.

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 performed H-alpha imaging observations of 22 luminous infrared galaxies to investigate how the distribution of star-forming regions in these galaxies is related to galaxy interactions. Based on correlation diagrams between H-alpha flux and continuum emission for individual galaxies, a sequence for the distribution of star-forming regions was found: very compact (~100 pc) nuclear starbursts with almost no star-forming activity in the outer regions (type 1), dominant nuclear starbursts < 1 kpc in size and a negligible contribution from the outer regions (type 2), nuclear starbursts > 1 kpc in size and a significant contribution from the outer regions (type 3), and extended starbursts with relatively faint nuclei (type 4). These classes of star-forming region were found to be strongly related to global star-forming properties such as star-formation efficiency, far-infrared color, and dust extinction. There was a clear tendency for the objects with more compact distributions of star-forming regions to show a higher star-formation efficiency and hotter far-infrared color. An appreciable fraction of the sample objects were dominated by extended starbursts (type 4), which is unexpected in the standard scenario of interaction-induced starburst galaxies. We also found that the distribution of star-forming regions was weakly but clearly related to galaxy morphology: severely disturbed objects had a more concentrated distribution of star-forming regions. This suggests that the properties of galaxy interactions, such as dynamical phase and orbital parameters, play a more important role than the internal properties of progenitor galaxies, such as dynamical structure or gas mass fraction. We also discuss the evolution of the distribution of star-forming regions in interacting galaxies.