We present the physical properties of AKARI sources without optical counterparts in optical images from the Hyper Suprime-Cam (HSC) on the Subaru telescope. Using the AKARI infrared (IR) source catalog and HSC optical catalog, we select 583 objects that do not have HSC counterparts in the AKARI North Ecliptic Pole (NEP) wide survey field ($\sim 5$ deg$^{2}$). Because the HSC limiting magnitude is deep ($g_{\rm AB}$ $\sim 28.6$), these are good candidates for extremely red star-forming galaxies (SFGs) and/or active galactic nuclei (AGNs), possibly at high redshifts. We compile multi-wavelength data out to 500 $\mu$m and use it for Spectral Energy Distribution (SED) fitting with CIGALE to investigate the physical properties of AKARI galaxies without optical counterparts. We also compare their physical quantities with AKARI mid-IR selected galaxies with HSC counterparts. The estimated redshifts of AKARI objects without HSC counterparts range up to $z\sim 4$, significantly higher than that of AKARI objects with HSC counterparts. We find that: (i) 3.6 $-$ 4.5 $\mu$m color, (ii) AGN luminosity, (iii) stellar mass, (iv) star formation rate, and (v) $V$-band dust attenuation in the interstellar medium of AKARI objects without HSC counterparts are systematically larger than those of AKARI objects with counterparts. These results suggest that our sample includes luminous, heavily dust-obscured SFGs/AGNs at $z\sim 1-4$ that are missed by previous optical surveys, providing very interesting targets for the coming James Webb Space Telescope era.

Observations of the molecular gas in galaxies are vital to understanding the evolution and star-forming histories of galaxies. However, galaxies with molecular gas maps of their whole discs having sufficient resolution to distinguish galactic structures are severely lacking. Millimeter wavelength studies at a high angular resolution across multiple lines and transitions are particularly needed, severely limiting our ability to infer the universal properties of molecular gas in galaxies. Hence, we conducted a legacy project with the 45 m telescope of the Nobeyama Radio Observatory, called the CO Multi-line Imaging of Nearby Galaxies (COMING), which simultaneously observed 147 galaxies with high far-infrared flux in $^{12}$CO, $^{13}$CO, and C$^{18}$O $J=1-0$ lines. The total molecular gas mass was derived using the standard CO-to-H$_2$ conversion factor and found to be positively correlated with the total stellar mass derived from the WISE $3.4 \mu$m band data. The fraction of the total molecular gas mass to the total stellar mass in galaxies does not depend on their Hubble types nor the existence of a galactic bar, although when galaxies in individual morphological types are investigated separately, the fraction seems to decrease with the total stellar mass in early-type galaxies and vice versa in late-type galaxies. No differences in the distribution of the total molecular gas mass, stellar mass, and the total molecular gas to stellar mass ratio was observed between barred and non-barred galaxies, which is likely the result of our sample selection criteria, in that we prioritized observing FIR bright (and thus molecular gas-rich) galaxies.

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.

We present the results of $^{12}$CO($J$=1-0) and $^{13}$CO($J$=1-0) simultaneous mappings toward the nearby barred spiral galaxy NGC 4303 as a part of the CO Multi-line Imaging of Nearby Galaxies (COMING) project. Barred spiral galaxies often show lower star-formation efficiency (SFE) in their bar region compared to the spiral arms. In this paper, we examine the relation between the SFEs and the volume densities of molecular gas $n(\rm{H}_2)$ in the eight different regions within the galactic disk with CO data combined with archival far-ultraviolet and 24 $\mu$m data. We confirmed that SFE in the bar region is lower by 39% than that in the spiral arms. Moreover, velocity-alignment stacking analysis was performed for the spectra in the individual regions. The integrated intensity ratios of $^{12}$CO to $^{13}$CO ($R_{12/13}$) range from 10 to 17 as the results of stacking. Fixing a kinetic temperature of molecular gas, $n(\rm{H}_2)$ was derived from $R_{12/13}$ via non-local thermodynamic equilibrium (non-LTE) analysis. The density $n(\rm{H}_2)$ in the bar is lower by 31-37% than that in the arms and there is a rather tight positive correlation between SFEs and $n(\rm{H}_2)$, with a correlation coefficient of $\sim 0.8$. Furthermore, we found a dependence of $n(\rm{H}_2)$ on the velocity dispersion of inter-molecular clouds ($\Delta V/ \sin i$). Specifically, $n(\rm{H}_2)$ increases as $\Delta V/ \sin i$ increases when $\Delta V/ \sin i < 100$ km s$^{-1}$. On the other hand, $n(\rm{H}_2)$ decreases as $\Delta V/ \sin i$ increases when $\Delta V/ \sin i > 100$ km s$^{-1}$. These relations indicate that the variations of SFE could be caused by the volume densities of molecular gas, and the volume densities could be governed by the dynamical influence such as cloud-cloud collisions, shear and enhanced inner-cloud turbulence.

Much of the star formation is obscured by dust. For the complete understanding of the cosmic star formation history (CSFH), infrared (IR) census is indispensable. AKARI carried out deep mid-infrared observations using its continuous 9-band filters in the North Ecliptic Pole (NEP) field (5.4 deg$^2$). This took significant amount of satellite's lifetime, $\sim$10\% of the entire pointed observations. By combining archival Spitzer (5 bands) and WISE (4 bands) mid-IR photometry, we have, in total, 18 band mid-IR photometry, which is the most comprehensive photometric coverage in mid-IR for thousands of galaxies. However previously, we only had shallow optical imaging ($\sim$25.9ABmag) in a small area of 1.0 deg$^2$. As a result, there remained thousands of AKARI's infrared sources undetected in optical. Using the new Hyper Suprime-Cam on Subaru telescope, we obtained deep enough optical images of the entire AKARI NEP field in 5 broad bands ($g\sim$27.5mag). These provided photometric redshift, and thereby IR luminosity for the previously undetected faint AKARI IR sources. Combined with the accurate mid-IR luminosity measurement, we constructed mid-IR LFs, and thereby performed a census of dust-obscured CSFH in the entire AKARI NEP field. We have measured restframe 8$\mu$m, 12$\mu$m luminosity functions (LFs), and estimated total infrared LFs at 0.35$<$z$<$2.2. Our results are consistent with our previous work, but with much reduced statistical errors thanks to the large area coverage of the new data. We have possibly witnessed the turnover of CSFH at $z\sim$2.

Revealing what fraction of galaxies harbor AGN is central in understanding black hole accretion history of the Universe. However, optical and soft X-ray surveys miss the most highly obscured AGNs. Infrared (IR), instead, is more robust against absorption. Previous IR photometric surveys, however, only had 4 or 5 filters in mid-IR. Our AKARI North Ecliptic Pole (NEP) wide field sample has 18 filters in mid-IR (9 from AKARI, 4 from WISE, and 5 from Spitzer), for the first time, allowing a sophisticated mid-IR SED fitting diagnosis for a statistical number of sources (89178 over 5.4 deg$^2$). By using a SED fitting technique, we investigate the evolution of AGN fraction as a function of redshift and IR (8-1000 $\mu$m) luminosity in an extinction-free way. We found that the AGN fraction (F$_{\rm AGN}$) shows no sign of strong redshift evolution. Instead, F$_{\rm AGN}$ increases with increasing IR luminosity in all redshifts bins ($0<z<2$).

We examined radial variations in molecular-gas based star formation efficiency (SFE), which is defined as star formation rate per unit molecular gas mass, for 80 galaxies selected from the CO Multi-line Imaging of Nearby Galaxies project (Sorai et al. 2019). The radial variations in SFE for individual galaxies are typically a factor of 2 -- 3, which suggests that SFE is nearly constant along galactocentric radius. We found the averaged SFE in 80 galaxies of $(1.69 \pm 1.1) \times 10^{-9}$ yr$^{-1}$, which is consistent with Leroy et al. 2008 if we consider the contribution of helium to the molecular gas mass evaluation and the difference in the assumed initial mass function between two studies. We compared SFE among different morphological (i.e., SA, SAB, and SB) types, and found that SFE within the inner radii ($r/r_{25} < 0.3$, where $r_{25}$ is $B$-band isophotal radius at 25 mag arcsec$^{-2}$) of SB galaxies is slightly higher than that of SA and SAB galaxies. This trend can be partly explained by the dependence of SFE on global stellar mass, which probably relates to the CO-to-H$_2$ conversion factor through the metallicity. For two representative SB galaxies in our sample, NGC 3367 and NGC 7479, the ellipse of $r/r_{25}$ = 0.3 seems to cover not only the central region but also the inner part of the disk, mainly the bar. These two galaxies show higher SFE in the bar than in spiral arms. However, we found an opposite trend in NGC 4303; SFE is lower in the bar than in spiral arms, which is consistent with earlier studies (e.g., Momose et al. 2010). These results suggest diversity of star formation activities in the bar.

The $^{12}$CO $(J=1\rightarrow0)$ velocity fields of a sample of 20 nearby spiral galaxies, selected from the CO Multi-line Imaging of Nearby Galaxies (COMING) legacy project of Nobeyama Radio Observatory, have been analyzed by Fourier decomposition to determine their basic kinematic properties, such as circular and noncircular velocities. On average, the investigated barred (SAB and SB) galaxies exhibit a ratio of noncircular to circular velocities of molecular gas larger by a factor of 1.5-2 than non-barred (SA) spiral galaxies at radii within the bar semimajor axis $a_\mathrm{b}$ at 1 kpc resolution, with a maximum at a radius of $R/a_\mathrm{b}\sim0.3$. Residual velocity field images, created by subtracting model velocity fields from the data, reveal that this trend is caused by kpc-scale streaming motions of molecular gas in the bar region. Applying a new method based on radial velocity reversal, we estimated the corotation radius $R_\mathrm{CR}$ and bar pattern speed $\Omega_\mathrm{b}$ in seven SAB and SB systems. The ratio of the corotation to bar radius is found to be in a range of $\mathcal{R}\equiv R_\mathrm{CR}/a_\mathrm{b}\sim0.8\mathrm{-}1.6$, suggesting that intermediate (SBb-SBc), luminous barred spiral galaxies host fast and slow rotator bars. Tentative negative correlations are found for $\Omega_\mathrm{b}$ vs. $a_\mathrm{b}$ and $\Omega_\mathrm{b}$ vs. total stellar mass $M_\ast$, indicating that bars in massive disks are larger and rotate slower, possibly a consequence of angular momentum transfer. The kinematic properties of SAB and SB galaxies, derived from Fourier decomposition, are compared with recent numerical simulations that incorporate various rotation curve models and galaxy interactions.

Understanding infrared (IR) luminosity is fundamental to understanding the cosmic star formation history and AGN evolution. Japanese infrared satellite, AKARI, provided unique data sets to probe this both at low and high redshift; the AKARI all sky survey in 6 bands (9-160 $\mu$m), and the AKARI NEP survey in 9 bands (2-24$\mu$m). The AKARI performed all sky survey in 6 IR bands (9, 18, 65, 90, 140, and 160 $\mu$m) with 3-10 times better sensitivity than IRAS, covering the crucial far-IR wavelengths across the peak of the dust emission. Combined with a better spatial resolution, we measure the total infrared luminosity ($L_{TIR}$) of individual galaxies, and thus, the total infrared luminosity density of the local Universe much more precisely than previous work. In the AKARI NEP wide field, AKARI has obtained deep images in the mid-infrared (IR), covering 5.4 deg$^2$. However, our previous work was limited to the central area of 0.25 deg$^2$ due to the lack of deep optical coverage. To rectify the situation, we used the newly advent Subaru telescope's Hyper Suprime-Cam to obtain deep optical images over the entire 5.4 deg$^2$ of the AKARI NEP wide field. With this deep and wide optical data, we, for the first time, can use the entire AKARI NEP wide data to construct restframe 8$\mu$m, 12$\mu$m, and total infrared (TIR) luminosity functions (LFs) at 0.15$<z<$2.2. A continuous 9-band filter coverage in the mid-IR wavelength (2.4, 3.2, 4.1, 7, 9, 11, 15, 18, and 24$\mu$m) by the AKARI satellite allowed us to estimate restframe 8$\mu$m and 12$\mu$m luminosities without using a large extrapolation based on a SED fit, which was the largest uncertainty in previous work. By combining these two results, we reveal dust-hidden cosmic star formation history and AGN evolution from z=0 to z=2.2, all probed by the AKARI satellite.

The AKARI North Ecliptic Pole (NEP) survey consists of two survey projects: NEP-Deep (0.5 sq.deg) and NEP-Wide (5.4 sq.deg), providing with tens of thousands of galaxies. A continuous filter coverage in the mid-infrared wavelengths (7, 9, 11, 15, 18 and 24 $\mu$m) is unique to diagnose the contributions from dusty star-formation activity and AGNs. Here we present current status focused on the newly obtained optical images and near-future prospects with a new X-ray telescope. Hyper Suprime-Cam on Subaru telescope is a gigantic optical camera with huge Field of View (FoV). Thanks to the wide FoV, we successfully obtained deep optical images at g, r, i, z and Y-bands covering most of the NEP-Wide field. Using the deep optical images, we identified over 5000 optical counterparts of the mid-IR sources, presumably deeply obscured galaxies in NEP-Wide field. We also investigated properties of these infrared sources with SED-fitting. eROSITA, to be launched early 2018, is a new all-sky X-ray survey telescope, and expected to conduct ultra deep 2-10 keV imaging toward NEP. We expect unprecedentedly numerous Compton-thick AGN candidates when combined with the multi-wavelength data in NEP region.

In this research, we provide a new, efficient method to select infrared (IR) active galatic nucleus (AGN). In the past, AGN selection in IR had been established by many studies using color-color diagrams. However, those methods have a problem in common that the number of bands is limited. The AKARI North Ecliptic Pole (NEP) survey was carried out by the AKARI Infrared Camera (IRC), which has 9 filters in mid-IR with a continuous wavelength coverage from 2 to 24$\mu$m$^{-1}$. Based on the intrinsic different mid-IR features of AGN and star-forming galaxies (SFGs), we performed SED fitting to separate these two populations by the best-fitting model. In the X-ray AGN sample, our method by SED fitting selects 50$\%$ AGNs, while the previous method by colour criteria recovers only 30$\%$ of them, which is a significant improvement. Furthermore, in the whole NEP deep sample, SED fitting selects two times more AGNs than the color selection. This may imply that the black hole accretion history could be more stronger than people expected before.

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.

We built a 8um selected sample of galaxies in the NEP-AKARI field by defining 4 redshift bins with the four AKARI bands at 11, 15, 18 and 24 microns (0.15<z<0.49, 0.75<z<1.34, 1.34<z<1.7 and 1.7<z<2.05) . Our sample contains 4079 sources, 599 are securely detected with Herschel/PACS. Also adding ultraviolet (UV) data from GALEX, we fit the spectral energy distributions using the physically motivated code CIGALE to extract the star formation rate, stellar mass, dust attenuation and the AGN contribution to the total infrared luminosity (L_{IR}). We discuss the impact of the adopted attenuation curve and that of the wavelength coverage to estimate these physical parameters. We focus on galaxies with a luminosity close the characteristic L_{IR}^* in the different redshift bins to study the evolution with redshift of the dust attenuation in these galaxies.

We have developed an efficient Active Galactic Nucleus (AGN) selection method using 18-band Spectral Energy Distribution (SED) fitting in mid-infrared (mid-IR). AGNs are often obscured by gas and dust, and those obscured AGNs tend to be missed in optical, UV and soft X-ray observations. Mid-IR light can help us to recover them in an obscuration free way using their thermal emission. On the other hand, Star-Forming Galaxies (SFG) also have strong PAH emission features in mid-IR. Hence, establishing an accurate method to separate populations of AGN and SFG is important. However, in previous mid-IR surveys, only 3 or 4 filters were available, and thus the selection was limited. We combined AKARI's continuous 9 mid-IR bands with WISE and Spitzer data to create 18 mid-IR bands for AGN selection. Among 4682 galaxies in the AKARI NEP deep field, 1388 are selected to be AGN hosts, which implies an AGN fraction of 29.6$\pm$0.8$\%$ (among them 47$\%$ are Seyfert 1.8 and 2). Comparing the result from SED fitting into WISE and Spitzer colour-colour diagram reveals that Seyferts are often missed by previous studies. Our result has been tested by stacking median magnitude for each sample. Using X-ray data from Chandra, we compared the result of our SED fitting with WISE's colour box selection. We recovered more X-ray detected AGN than previous methods by 20$\%$.

Mass, metallicity, and star formation rate (SFR) of a galaxy are crucial parameters in understanding galaxy formation and evolution. However, the relation among these is still a matter of debate for luminous infrared galaxies, which carry a bulk of SFR budget of the universe at $z\sim1$. We have investigated the relation among stellar mass, gas-phase oxygen abundance, and SFR of AKARI-detected mid-IR galaxies at $z\sim0.88$ in the AKARI NEP deep field. We observed about 350 AKARI sources with Subaru/FMOS NIR spectrograph, and detected secure and expected H$\alpha$ emission lines from 25 and 44 galaxies, respectively. The SFR of our sample is almost constant ($\sim 25M_{\odot}/yr$) over the stellar mass range of our sample. Compared with main-sequence (MS) galaxies at a similar redshift range, the average SFR of our detected sample is comparable for massive galaxies ($\sim10^{10.58}~M_{\odot}$), while higher by $\sim$0.6dex for less massive galaxies ($\sim 10^{10.05}~M_{\odot}$). We measure metallicities from the [NII]/H$\alpha$ emission line ratio. We find that the mass-metallicity relation of our individually measured sources agrees with that for optical-selected star-forming galaxies at $z\sim0.1$, while metallicities of stacked spectra agree with that of MS galaxies at $z\sim0.78$. Considering high SFR of individually measured sources, FMR of the IR galaxies is different from that at $z\sim0.1$. However, on the mass-metallicity plane, they are consistent with the MS galaxies, highlighting higher SFR of the IR galaxies. This suggests the evolutionary path of our IR galaxies is different from that of MS galaxies. A possible physical interpretation includes that the star-formation activities of IR galaxies at $z\sim0.88$ in our sample are enhanced by interaction and/or merger of galaxies, but the inflow of metal-poor gas is not yet induced, keeping the metallicity intact.

We performed near-diffraction-limited (~0.4 FWHM) N-band imaging of one of the nearest Active Galactic Nucleus (AGN) in M51 with 8.2m Subaru telescope to study the nuclear structure and spectral energy distribution (SED) at 8-13 um. We found that the nucleus is composed of an unresolved core (at ~13 pc resolution, or intrinsic size corrected for the instrumental effect of <6 pc) and an extended halo (at a few tens pc scale), and each of their SEDs is almost flat. We examined the SED by comparing with the archival Spitzer IRS spectrum processed to mimic our chopping observation of the nucleus, and the published radiative-transfer model SEDs of the AGN clumpy dusty torus. The halo SED is likely due to circumnuclear star formation showing little Polycyclic Aromatic Hydrocarbon (PAH) emission due to the AGN. The core SED is likely dominated by the AGN because of the following two reasons. Firstly, the clumpy torus model SEDs can reproduce the red mid-infrared continuum with apparently moderate silicate 9.7 um absorption. Secondly, the core 12 um luminosity and the absorption-corrected X-ray luminosity at 2-10 keV in the literature follow the mid-infrared-X-ray luminosity correlation known for the nearby AGNs including the Compton-thick ones.

We present the results of new infrared spectroscopic observations of 37 quasars at z~3, selected based on the optical r'-band magnitude and the availability of nearby bright stars for future imaging follow-up with Adaptive Optics system. The supermassive black hole (SMBH) masses (M_BH) were successfully estimated in 28 out of 37 observed objects from the combination of the H_beta emission linewidth and continuum luminosity at rest-frame 5100A. Comparing these results with those from previous studies of quasars with similar redshift, our sample exhibited slightly lower (~ -0.11 dex in median) Eddington ratios; and, the SMBH masses are slightly (~ 0.38 dex in median) higher. The SMBH growth time, t_grow, was calculated by dividing the estimated SMBH mass by the mass accretion rate measured using optical luminosity. We found, given reasonable assumptions, that t_grow was smaller than the age of the universe at the redshift of individual quasars for a large fraction of observed sources, suggesting that the SMBHs in many of our observed quasars are in growing phase with high accretion rates. A comparison of the SMBH masses derived from our H_beta data and archived CIV data indicated considerable scattering, as indicated in previous studies. All quasars with measured SMBH masses have at least one nearby bright star, such that they are suitable targets for adaptive optics observations to study the mass relationship between SMBHs and host galaxies' stellar component at high redshift.

Using star-forming galaxies sample in the nearby Universe (0.02<z<0.10) selected from the SDSS (DR7) and GALEX all-sky survey (GR5), we present a new empirical calibration for predicting dust extinction of galaxies from H-alpha-to-FUV flux ratio. We find that the H-alpha dust extinction (A(Ha)) derived with H-alpha/H-beta ratio (Balmer decrement) increases with increasing H-alpha/UV ratio as expected, but there remains a considerable scatter around the relation, which is largely dependent on stellar mass and/or H-alpha equivalent width (EW(Ha)). At fixed H-alpha/UV ratio, galaxies with higher stellar mass (or galaxies with lower EW(Ha)) tend to be more highly obscured by dust. We quantify this trend and establish an empirical calibration for predicting A(Ha) with a combination of H-alpha/UV ratio, stellar mass and EW(Ha), with which we can successfully reduce the systematic uncertainties accompanying the simple H-alpha/UV approach by ~15-30%. The new recipes proposed in this study will provide a convenient tool for predicting dust extinction level of galaxies particularly when Balmer decrement is not available. By comparing A(Ha) (derived with Balmer decrement) and A(UV) (derived with IR/UV luminosity ratio) for a subsample of galaxies for which AKARI FIR photometry is available, we demonstrate that more massive galaxies tend to have higher extra extinction towards the nebular regions compared to the stellar continuum light. Considering recent studies reporting smaller extra extinction towards nebular regions for high-redshift galaxies, we argue that the dust geometry within high-redshift galaxies resemble more like low-mass galaxies in the nearby Universe.

We present infrared galaxy luminosity functions (LFs) in the AKARI North Ecliptic Pole (NEP) deep field using recently-obtained, wider CFHT optical/near-IR images. AKARI has obtained deep images in the mid-infrared (IR), covering 0.6 deg$^2$ of the NEP deep field. However, our previous work was limited to the central area of 0.25 deg$^2$ due to the lack of optical coverage of the full AKARI NEP survey. To rectify the situation, we recently obtained CFHT optical and near-IR images over the entire AKARI NEP deep field. These new CFHT images are used to derive accurate photometric redshifts, allowing us to fully exploit the whole AKARI NEP deep field. AKARI's deep, continuous filter coverage in the mid-IR wavelengths (2.4, 3.2, 4.1, 7, 9, 11, 15, 18, and 24$\mu$m) exists nowhere else, due to filter gaps of other space telescopes. It allows us to estimate restframe 8$\mu$m and 12$\mu$m luminosities without using a large extrapolation based on spectral energy distribution (SED) fitting, which was the largest uncertainty in previous studies. Total infrared luminosity (TIR) is also obtained more reliably due to the superior filter coverage. The resulting restframe 8$\mu$m, 12$\mu$m, and TIR LFs at $0.15<z<2.2$ are consistent with previous works, but with reduced uncertainties, especially at the high luminosity-end, thanks to the wide field coverage. In terms of cosmic infrared luminosity density ($\Omega_{\mathrm{IR}}$), we found that the $\Omega_{\mathrm{IR}}$ evolves as $\propto (1+z)^{4.2\pm 0.4}$.

We present the photometric properties of a sample of infrared (IR) bright dust obscured galaxies (DOGs). Combining wide and deep optical images obtained with the Hyper Suprime-Cam (HSC) on the Subaru Telescope and all-sky mid-IR (MIR) images taken with Wide-Field Infrared Survey Explorer (WISE), we discovered 48 DOGs with $i - K_\mathrm{s} > 1.2$ and $i - [22] > 7.0$, where $i$, $K_\mathrm{s}$, and [22] represent AB magnitude in the $i$-band, $K_\mathrm{s}$-band, and 22 $\mu$m, respectively, in the GAMA 14hr field ($\sim$ 9 deg$^2$). Among these objects, 31 ($\sim$ 65 %) show power-law spectral energy distributions (SEDs) in the near-IR (NIR) and MIR regime, while the remainder show a NIR bump in their SEDs. Assuming that the redshift distribution for our DOGs sample is Gaussian, with mean and sigma $z$ = 1.99 $\pm$ 0.45, we calculated their total IR luminosity using an empirical relation between 22 $\mu$m luminosity and total IR luminosity. The average value of the total IR luminosity is (3.5 $\pm$ 1.1) $\times$ $10^{13}$ L$_{\odot}$, which classifies them as hyper-luminous infrared galaxies (HyLIRGs). We also derived the total IR luminosity function (LF) and IR luminosity density (LD) for a flux-limited subsample of 18 DOGs with 22 $\mu$m flux greater than 3.0 mJy and with $i$-band magnitude brighter than 24 AB magnitude. The derived space density for this subsample is log $\phi$ = -6.59 $\pm$ 0.11 [Mpc$^{-3}$]. The IR LF for DOGs including data obtained from the literature is well fitted by a double-power law. The derived lower limit for the IR LD for our sample is $\rho_{\mathrm{IR}}$ $\sim$ 3.8 $\times$ 10$^7$ [L$_{\odot}$ Mpc$^{-3}$] and its contributions to the total IR LD, IR LD of all ultra-luminous infrared galaxies (ULIRGs), and that of all DOGs are $>$ 3 %, $>$ 9 %, and $>$ 15 %, respectively.

(Abridged) We aim to study the evolution of dust attenuation in galaxies selected in the IR in the redshift range in which they are known to dominate the star formation activity in the universe. The comparison with other measurements of dust attenuation in samples selected using different criteria will give us a global picture of the attenuation at work in star-forming galaxies and its evolution with redshift. Using multiple filters of IRC instrument, we selected more than 4000 galaxies from their rest-frame emission at 8 microns, from z~0.2 to 2$. We built SEDs from the rest-frame UV to the far-IR by adding data in the optical-NIR and from GALEX and Herschel surveys. We fit SEDs with the physically-motivated code CIGALE. We test different templates for AGNs and recipes for dust attenuation and estimate stellar masses, SFRs, amount of dust attenuation, and AGN contribution to the total IR luminosity. The AGN contribution to the total IR luminosity is found to be on average approximately 10% with a slight increase with redshift. Dust attenuation in galaxies dominating the IR luminosity function is found to increase from z=0 to z=1 and to remain almost constant from z=1 to z=1.5. Conversely, when galaxies are selected at a fixed IR luminosity, their dust attenuation slightly decreases as redshift increases but with a large dispersion. The attenuation in our mid-IR selected sample is found ~ 2 mag higher than that found globally in the universe or in UV and Halpha line selections in the same redshift range. This difference is well explained by an increase of dust attenuation with the stellar mass, in global agreement with other recent studies. Starbursting galaxies do not systematically exhibit a high attenuation

In this work, we investigate the dependence of the covering factor (CF) of active galactic nuclei (AGNs) (i) on the mid-infrared (MIR) luminosity and (ii) on the redshift. We constructed 12- and 22-micron luminosity functions (LFs) at 0.006 < z < 0.3 using the Wide-field Infrared Survey Explorer} (WISE) data. Combining the WISE catalog with the Sloan Digital Sky Survey (SDSS) spectroscopic data, we selected 223,982 galaxies at 12 micron and 25,721 galaxies at 22 micron for spectroscopic classification. We then identified 16,355 AGNs at 12 micron and 4,683 AGNs at 22 micron by their optical emission lines and cataloged classifications in the SDSS. Following that, we estimated the CF as the fraction of type 2 AGN in all AGNs whose MIR emissions are dominated by the active nucleus (not their host galaxies) based on their MIR colors. We found that (i) the CF decreased with increasing MIR luminosity, regardless of the choice of type 2 AGN classification criteria, and (ii) the CF did not change significantly with the redshift for z < 0.2. Furthermore, we carried out various tests to determine the influence of selection bias and confirmed similar dependences exist even when taking these uncertainties into account. The luminosity dependence of the CF can be explained by the receding torus model, but the "modified" receding torus model gives a slightly better fit, as suggested by Simpson.

Aims. We present an 8-band (u*, g', r', i', z', Y, J, Ks) optical to NIR deep photometric catalog based on the observations made with MegaCam and WIRCam at CFHT, and compute photometric redshifts, zp in the North Ecliptic Pole (NEP) region. Our catalog provides us to identify the counterparts, and zp for AKARI NIR/MIR sources. Results. The estimated 4sigma detection limits within an 1" aperture radius are 26.7, 25.9, 25.1, and 24.1 mag [AB] for g', r', i', and z'-bands and 23.4, 23.0, and 22.7 mag for Y, J, and Ks-bands, respectively. There are a total of 85797 sources in the band-merged catalog. An astrometric accuracy of this catalog determined by examining coordinate offsets with regard to 2MASS is 0.013" with a root mean square offset of 0.32". We distinguish 5441 secure stars from extended sources using u*-J vs. g'-Ks colors, combined with the SExtractor stellarity index of the images. Comparing with galaxy spectroscopic redshifts, we find a photometric redshift dispersion, sigma_(dz/(1+z)), of 0.032 and catastrophic failure rate, dz/(1+z)>0.15, of 5.8% at z<1, while a dispersion of 0.117 and a catastrophic failures rate of 16.6% at z>1. We extend estimate of the zp uncertainty over the full magnitude/redshift space with a redshift probability distribution function and find that our redshift are highly accurate with z'<22 at zp<2.5 and for fainter sources with z'<24 at z<1. From the investigation of photometric properties of AKARI infrared sources (23354 sources) using the g'z'Ks diagram, <5% of AKARI sources with optical counterparts are classified as high-z (1.4<z<2.5) star-forming galaxies. Among the high-z star-forming galaxies, AKARI MIR detected sources seem to be affected by stronger dust extinction compared with sources with non-detections in the AKARI MIR bands. The full, electronic version of our catalog with zp will be available at the CDS.

Several lines of argument support the existence of a link between activity at the nuclei of galaxies, in the form of an accreting supermassive black hole, and star-formation activity in these galaxies. The exact nature of this link is still under debate. Radio jets have long been argued to be an ideal mechanism that allows AGN to interact with their host galaxy and regulate star-formation. In this context, we are using a sample of radio sources in the North Ecliptic Pole (NEP) field to study the nature of the putative link between AGN activity and star-formation. This is done by means of spectral energy distribution (SED) fitting. We use the excellent spectral coverage of the AKARI infrared space telescope together with the rich ancillary data available in the NEP to build SEDs extending from UV to far-IR wavelengths. Through SED fitting we constrain both the AGN and host galaxy components. We find a significant AGN component in our sample of relatively faint radio-sources ($<$mJy), that increases in power with increasing radio-luminosity. At the highest radio-luminosities, the presence of powerful jets dominates the radio emission of these sources. A positive correlation is found between the luminosity of the AGN component and that of star-formation in the host galaxy, independent of the radio luminosity. By contrast, for a given redshift and AGN luminosity, we find that increasing radio-luminosity leads to a decrease in the specific star-formation rate. The most radio-loud AGN are found to lie on the main sequence of star-formation for their respective redshifts. For the first time, such a two-sided feedback process is seen in the same sample. We conclude that radio jets do suppress star-formation in their host galaxies but appear not to totally quench it. Our results therefore support the maintenance nature of "radio-mode" feedback from radio-AGN jets.

We present high-resolution mid-infrared (MIR) imaging, nuclear spectral energy distributions (SEDs) and archival Spitzer spectra for 22 low-luminosity active galactic nuclei (LLAGN; Lbol \lesssim 10^42 erg/sec). Infrared (IR) observations may advance our understanding of the accretion flows in LLAGN, the fate of the obscuring torus at low accretion rates, and, perhaps, the star formation histories of these objects. However, while comprehensively studied in higher-luminosity Seyferts and quasars, the nuclear IR properties of LLAGN have not yet been well-determined. We separate the present LLAGN sample into three categories depending on their Eddington ratio and radio emission, finding different IR characteristics for each class. (I) At the low-luminosity, low-Eddington ratio (log Lbol/LEdd < -4.6) end of the sample, we identify "host-dominated" galaxies with strong polycyclic aromatic hydrocarbon bands that may indicate active (circum-)nuclear star formation. (II) Some very radio-loud objects are also present at these low Eddington ratios. The IR emission in these nuclei is dominated by synchrotron radiation, and some are likely to be unobscured type 2 AGN that genuinely lack a broad line region. (III) At higher Eddington ratios, strong, compact nuclear sources are visible in the MIR images. The nuclear SEDs of these galaxies are diverse; some resemble typical Seyfert nuclei, while others lack a well-defined MIR "dust bump". Strong silicate emission is present in many of these objects. We speculate that this, together with high ratios of silicate strength to hydrogen column density, could suggest optically thin dust and low dust-to-gas ratios, in accordance with model predictions that LLAGN do not host a Seyfert-like obscuring torus.

We have obtained high-resolution mid-infrared (MIR) imaging, nuclear spectral energy distributions (SEDs) and archival Spitzer spectra for 22 low-luminosity active galactic nuclei (LLAGN; L_bol < 5 x 10^42 erg/s). Infrared (IR) observations may advance our understanding of the accretion flows in LLAGN, the fate of the obscuring torus at low accretion rates, and, perhaps, the star formation histories of these objects. However, while comprehensively studied in higher-luminosity Seyferts and quasars, the nuclear IR properties of LLAGN have not yet been well-determined. In these proceedings we summarise the results for the LLAGN at the relatively high-luminosity, high-Eddington ratio end of the sample. Strong, compact nuclear sources are visible in the MIR images of these objects, with luminosities consistent with or slightly in execss of that predicted by the standard MIR/X-ray relation. Their broadband nuclear SEDs are diverse; some resemble typical Seyfert nuclei, while others possess less of a well-defined MIR ``dust bump''. Strong silicate emission is present in many of these objects. We speculate that this, together with high ratios of silicate strength to hydrogen column density, could suggest optically thin dust and low dust-to-gas ratios, in accordance with model predictions that LLAGN do not host a Seyfert-like obscuring torus.

We present the discovery and mass measurement of the cold, low-mass planet MOA-2009-BLG-266Lb, made with the gravitational microlensing method. This planet has a mass of m_p = 10.4 +- 1.7 Earth masses and orbits a star of mass M_* = 0.56 +- 0.09 Solar masses at a semi-major axis of a = 3.2 (+1.9 -0.5) AU and an orbital period of P = 7.6 (+7.7 -1.5} yrs. The planet and host star mass measurements are enabled by the measurement of the microlensing parallax effect, which is seen primarily in the light curve distortion due to the orbital motion of the Earth. But, the analysis also demonstrates the capability to measure microlensing parallax with the Deep Impact (or EPOXI) spacecraft in a Heliocentric orbit. The planet mass and orbital distance are similar to predictions for the critical core mass needed to accrete a substantial gaseous envelope, and thus may indicate that this planet is a "failed" gas giant. This and future microlensing detections will test planet formation theory predictions regarding the prevalence and masses of such planets.

We report high resolution 12CO(J = 2-1), 13CO(J = 2-1), and 12CO(J = 3-2) imaging of the Seyfert 1/starburst ring galaxy NGC 1097 with the Submillimeter Array to study the physical and kinematic properties of the 1-kpc circumnuclear starburst ring. Individual star clusters as detected in the HST map of Pa_alpha line emission have been used to determine the star formation rate, and are compared with the properties of the molecular gas. The molecular ring has been resolved into individual clumps at GMA-scale of 200--300 pc in all three CO lines. The intersection between the dust lanes and the starburst ring, which is associated with orbit-crowding region, is resolved into two physically/kinematically distinct features in the 1.5x1.0 (105x70 pc) 12CO(J = 2-1) map. The clumps associated with the dust lanes have broader line width, higher surface gas density, and lower star formation rate, while the narrow line clumps associated with the starburst ring have opposite characteristics. Toomre-Q value under unity at the radius of the ring suggests that the molecular ring is gravitationally unstable to fragment at the scale of the GMA. The line widths and surface density of gas mass of the clumps show an azimuthal variation related to the large scale dynamics. The star formation rate, on the other hand, is not significantly affected by the dynamics, but has a correlation with the intensity ratio of 12CO (J = 3-2) and 12CO(J = 2-1), which traces the denser gas associated with star formation. Our resolved CO map, especially in the orbit-crowding region, for the first time demonstrates observationally that the physical/kinematic properties of the GMAs are affected by the large scale bar-potential dynamics in NGC 1097.

X-ray and near-infrared ($J$-$H$-$K_{\rm s}$) observations of the Galactic black hole binary GX 339--4 in the low/hard state were performed with Suzaku and IRSF in 2009 March. The spectrum in the 0.5--300 keV band is dominated by thermal Comptonization of multicolor disk photons, with a small contribution from a direct disk component, indicating that the inner disk is almost fully covered by hot corona with an electron temperature of $\approx$175 keV. The Comptonizing corona has at least two optical depths, $\tau \approx 1,0.4$. Analysis of the iron-K line profile yields an inner disk radius of $(13.3^{+6.4}_{-6.0}) R_{\rm g}$ ($R_{\rm g} $ represents the gravitational radius $GM/c^2$), with the best-fit inclination angle of $\approx50^\circ$. This radius is consistent with that estimated from the continuum fit by assuming the conservation of photon numbers in Comptonization. Our results suggest that the standard disk of GX 339--4 is likely truncated before reaching the innermost stable circular orbit (for a non rotating black hole) in the low/hard state at $\sim$1% of the Eddington luminosity. The one-day averaged near-infrared light curves are found to be correlated with hard X-ray flux with $F_{\rm Ks} \propto F_{\rm X}^{0.45}$. The flatter near infrared $\nu F_{\nu}$ spectrum than the radio one suggests that the optically thin synchrotron radiation from the compact jets dominates the near-infrared flux. Based on a simple analysis, we estimate the magnetic field and size of the jet base to be $5\times10^4$ G and $6\times 10^8$ cm, respectively. The synchrotron self Compton component is estimated to be approximately 0.4% of the total X-ray flux.

We present the results of a ground-based, high spatial resolution infrared 18 micron imaging study of nearby luminous infrared galaxies (LIRGs), using the Subaru 8.2-m and Gemini South 8.1-m telescopes. The diffraction-limited images routinely achieved with these telescopes in the Q-band (17-23 micron) allow us to investigate the detailed spatial distribution of infrared emission in these LIRGs. We then investigate whether the emission surface brightnesses are modest, as observed in starbursts, or are so high that luminous active galactic nuclei (AGNs; high emission surface brightness energy sources) are indicated. The sample consists of 18 luminous buried AGN candidates and starburst-classified LIRGs identified in earlier infrared spectroscopy. We find that the infrared 18 micron emission from the buried AGN candidates is generally compact, and the estimated emission surface brightnesses are high, sometimes exceeding the maximum value observed in and theoretically predicted for a starburst phenomenon. The starburst-classified LIRGs usually display spatially extended 18 micron emission and the estimated emission surface brightnesses are modest, within the range sustained by a starburst phenomenon. The general agreement between infrared spectroscopic and imaging energy diagnostic methods suggests that both are useful tools for understanding the hidden energy sources of the dusty LIRG population.

We present the results of infrared L-band (3-4 micron) slit spectroscopy of 30 PG QSOs at z < 0.17, the representative sample of local high-luminosity, optically selected AGNs. The 3.3 micron polycyclic aromatic hydrocarbon (PAH) emission feature is used to probe nuclear (< a few kpc) starburst activity and to investigate the connections between AGNs and nuclear starbursts in PG QSOs. The 3.3 micron PAH emission is detected in the individual spectra of 5/30 of the observed PG QSOs. We construct a composite spectrum of PAH-undetected PG QSOs and discern the presence of the 3.3 micron PAH emission therein. We estimate the nuclear-starburst and AGN luminosities from the observed 3.3 micron PAH emission and 3.35 micron continuum luminosities, respectively, and find that the nuclear-starburst-to-AGN luminosity ratios in PG QSOs are similar to those of previously studied AGN populations with lower luminosities, suggesting that AGN-nuclear starburst connections are valid over the wide luminosity range of AGNs in the local universe. The observed nuclear-starburst-to-AGN luminosity ratios in PG QSOs with available supermassive black hole masses are comparable to a theoretical prediction based on the assumption that the growth of a supermassive black hole is controlled by starburst-induced turbulence.

We present near-infrared K- (1.9-2.5 micrometer) and L- (2.8-4.2 micrometer) band spectroscopy of 22 Seyfert nuclei. We use two methods to investigate the presence of nuclear starbursts: (1) the Polycyclic Aromatic Hydrocarbon (PAH) emission feature at lambda_rest = 3.3 micrometer in the rest frame of L-band spectrum (a starburst indicator) and (2) the CO absorption feature at lambda_rest = 2.3-2.4 micrometer in the rest frame of the K-band spectrum, originating in the CO molecule. We clearly detected the 3.3 micrometer PAH emission features in five objects and the CO absorption features in 17 objects. Seyfert 2 galaxies tend to show bluer K-L colors compared with Seyfert 1 galaxies. We interpret the discrepancy as resulting from relative strength of stellar emission because AGN emission is affected by dust extinction. The 3.3 micrometer PAH emission luminosity (L_3.3PAH) distributions for the Seyfert 1s and Seyfert 2s are very similar when normalized to the AGN power. Star-formation rates estimated from L_3.3PAH could be large enough to inflate the dusty torus by supernova explosion. We find that L_3.3PAH positively correlates with N-band luminosity with small aperture over a wide AGN luminosity range, and is independent of physical area we probed. The results suggest that nuclear region has a concentration of star formation and the star formation would control AGN activity.

We report on the results of systematic infrared 2.5-5 micron spectroscopy of 45 nearby ultraluminous infrared galaxies (ULIRGs) at z < 0.3 using IRC onboard the AKARI satellite. This paper investigates whether the luminosities of these ULIRGs are dominated by starburst activity, or optically elusive buried AGNs are energetically important. Our criteria include the strengths of the 3.3 micron polycyclic aromatic hydrocarbon (PAH) emission features and the optical depths of absorption features at 3.1 micron due to ice-covered dust grains and at 3.4 micron from bare carbonaceous dust grains. Because of the AKARI IRC's spectroscopic capability in the full 2.5-5 micron wavelength range, unaffected by Earth's atmosphere, we can apply this energy diagnostic method to ULIRGs at z > 0.15. We estimate the intrinsic luminosities of extended (several kpc), modestly obscured (Av < 15 mag) starburst activity based on the 3.3 micron PAH emission luminosities measured in AKARI IRC slitless spectra, and confirm that such starbursts are energetically unimportant in nearby ULIRGs. In roughly half of the observed ULIRGs classified optically as non-Seyferts, we find signatures of luminous energy sources that produce no PAH emission and/or are more centrally concentrated than the surrounding dust. We interpret these energy sources as buried AGNs. The fraction of ULIRGs with detectable buried AGN signatures increases with increasing infrared luminosity. Our overall results support the scenario that luminous buried AGNs are important in many ULIRGs at z < 0.3 classified optically as non-Seyferts, and that the optical undetectability of such buried AGNs occurs merely because of a large amount of nuclear dust, which can make the sightline of even the lowest dust column density opaque to the ionizing radiation of the AGNs.

We present the results on millimeter interferometric observations of four luminous infrared galaxies (LIRGs), Arp 220, Mrk 231, IRAS 08572+3915, and VV 114, and one Wolf-Rayet galaxy, He 2-10, using the Nobeyama Millimeter Array (NMA). Both the HCN(1-0) and HCO+(1-0) molecular lines were observed simultaneously and their brightness-temperature ratios were derived. High-quality infrared L-band (2.8-4.1 micron) spectra were also obtained for the four LIRGs to better constrain their energy sources deeply buried in dust and molecular gas. When combined with other LIRGs we have previously observed with NMA, the final sample comprised nine LIRGs (12 LIRGs' nuclei) with available interferometric HCN(1-0) and HCO+(1-0) data-sufficient to investigate the overall trend in comparison with known AGNs and starburst galaxies. We found that LIRGs with luminous buried AGN signatures at other wavelengths tend to show high HCN(1-0)/HCO+(1-0) brightness-temperature ratios as seen in AGN-dominated galaxies, while the Wolf-Rayet galaxy He 2-10 displays a small ratio. An enhanced HCN abundance in the interstellar gas surrounding a strongly X-ray-emitting AGN, as predicted by some chemical calculations, is a natural explanation of our results.