The need for a method to construct multidimensional distribution function is increasing recently, in the era of huge multiwavelength surveys. We have proposed a systematic method to build a bivariate luminosity or mass function of galaxies by using a copula. It allows us to construct a distribution function when only its marginal distributions are known, and we have to estimate the dependence structure from data. A typical example is the situation that we have univariate luminosity functions at some wavelengths for a survey, but the joint distribution is unknown. Main limitation of the copula method is that it is not easy to extend a joint function to higher dimensions ($d > 2$), except some special cases like multidimensional Gaussian. Even if we find such a multivariate analytic function in some fortunate case, it would often be inflexible and impractical. In this work, we show a systematic method to extend the copula method to unlimitedly higher dimensions by a vine copula. This is based on the pair-copula decomposition of a general multivariate distribution. We show how the vine copula construction is flexible and extendable. We also present an example of the construction of an stellar mass--atomic gas--molecular gas 3-dimensional mass function. We demonstrate the maximum likelihood estimation of the best functional form for this function, as well as a proper model selection via vine copula.

(abridged) The scale of the acoustic oscillation of baryons at the baryon-photon decoupling is imprinted on the spatial distribution of galaxies in the Universe, known as the baryon acoustic oscillation (BAO). The correlation functions and power spectrum are used as a central tool for the studies on the BAO analysis. In this work, we analyzed the spatial distribution of galaxies with a method from the topological data analysis (TDA), in order to detect and examine the BAO signal in the galaxy distribution. The TDA provides a method to treat various types of "holes" in point set data, by constructing the persistent homology (PH) group from the geometric structure of data points and handling the topological information of the dataset. We can obtain the information on the size, position, and statistical significance of the holes in the data. A particularly strong point of the persistent homology is that it can classify the holes by their spatial dimension, i.e., a 0-dim separation, 1-dim loop, 2-dim shell, etc. We first analyzed the simulation datasets with and without the baryon physics to examine the performance of the PH method. We found that the PH is indeed able to detect the BAO signal: simulation data with baryon physics present a prominent signal from the BAO, while data without baryon physics does not show this signal. Then, we applied the PH to a quasar sample at $z <1.0$ from extended Baryon Oscillation Spectroscopic Survey in Sloan Digital Sky Survey Data Release 14. We discovered a characteristic hole (a hollow shell) at a scaler $\sim150\ [{\rm Mpc}]$. This exactly corresponds to the BAO signature imprinted in the galaxy/quasar distribution. We performed this analysis on a small subsample of 2000 quasars. This clearly demonstrates that the PH analysis is very efficient in finding this type of topological structures even if the sampling is very sparse.

We discuss the existence of an acceleration scale in galaxies and galaxy clusters. The presence of the same acceleration scale found at very different scales and in very different astrophysical objects strongly supports the existence of a fundamental acceleration scale governing the observed gravitational physics. We also comment on the implication of such a fundamental acceleration scale on the problem of dark matter. We discuss the relevance of the fundamental acceleration for the nature of dark matter as well as for structure formation to be explored in future numerical simulations.

Photometric surveys have provided incredible amounts of astronomical information in the form of images. However, astronomical images often contain artifacts that can critically hinder scientific analysis by misrepresenting intensities or contaminating catalogs as artificial objects. These affected pixels need to be masked and dealt with in any data reduction pipeline. In this paper, we present a flexible, iterative algorithm to recover (unmask) astronomical images where some pixels are lacking. We demonstrate the application of the method on some intensity calibration source images in CO Multi-line Imaging of Nearby Galaxies (COMING) Project conducted using the 45m telescope at Nobeyama Radio Observatory (NRO). The proposed algorithm restored artifacts due to a detector error in the intensity calibration source images. The restored images were used to calibrate 11 out of 147 observed galaxy maps in the survey. The tests show that the algorithm can restore measured intensities at sub 1% error even for noisy images (SNR = 2.4), despite lacking a significant part of the image. We present the formulation of the reconstruction algorithm, discuss its possibilities and limitations for extensions to other astronomical signals and the results of the COMING application.

The first generation of stars were born a few hundred million years after the big bang. These stars synthesized elements heavier than H and He, that are later expelled into the interstellar medium, initiating the rise of metals. Within this enriched medium, the first dust grains formed. This event is cosmological crucial for molecule formation as dust plays a major role by cooling low-metallicity star-forming clouds which can fragment to create lower mass stars. Collecting information on these first dust grains is difficult because of the negative alliance of large distances and low dust masses. We combine the observational information from galaxies at redshifts 5 < z < 10 to constrain their dust emission and theoretically understand the first evolutionary phases of the dust cycle. Spectral energy distributions (SEDs) are fitted with CIGALE and the physical parameters and their evolution are modelled. From this SED fitting, we build a dust emission template for this population of galaxies in the epoch of reionization. Our new models explain why some early galaxies are observed and others are not. We follow in time the formation of the first grains by supernovae later destroyed by other supernova blasts and expelled in the circumgalactic and intergalactic media. We have found evidence for the first dust grains formed in the universe. But, above all, this letter underlines the need to collect more data and to develop new facilities to further constrain the dust cycle in galaxies in the epoch of reionization.

The Borexino detector measures solar neutrino fluxes via neutrino-electron elastic scattering. Observed spectra are determined by the solar-$\nu_{e}$ survival probability $P_{ee}(E)$, and the chiral couplings of the neutrino and electron. Some theories of physics beyond the Standard Model postulate the existence of Non-Standard Interactions (NSI's) which modify the chiral couplings and $P_{ee}(E)$. In this paper, we search for such NSI's, in particular, flavor-diagonal neutral current interactions that modify the $\nu_e e$ and $\nu_\tau e$ couplings using Borexino Phase II data. Standard Solar Model predictions of the solar neutrino fluxes for both high- and low-metallicity assumptions are considered. No indication of new physics is found at the level of sensitivity of the detector and constraints on the parameters of the NSI's are placed. In addition, with the same dataset the value of $\sin^2\theta_W$ is obtained with a precision comparable to that achieved in reactor antineutrino experiments.

We study the spectrum of semiclassical rotating strings on de Sitter space and its consistency. Even though a naive extrapolation of the linear Regge trajectory on flat space implies a violation of the Higuchi bound (a unitarity bound on the mass of higher-spin particles in de Sitter space), the curved space effects turn out to modify the trajectory to respect the bound. Interestingly, there exists a maximum spin for each Regge trajectory as a consequence of accelerated expansion, which is helpful to make the spectrum consistent with the Higuchi bound, but at the same time it could be an obstruction to stringy UV completion based on an infinite higher-spin tower. By pushing further this observation, we demonstrate that the vacuum energy $V$ inflating the universe has to be bounded by the string scale $M_s$ as $V\lesssim M_s^4$, if UV completion is achieved by the leading Regge trajectory with higher spin states up to the 4D Planck scale. Its application to inflation at the early universe implies an upper bound on the tensor-to-scalar ratio, $r\lesssim 0.01\times(M_s/10^{16} \text{GeV})^{4}$, which is within the scope of the near future CMB experiments. We also discuss another possibility that UV completion is achieved by multiple Regge trajectories.

An important aspect of quenching star formation is the removal of the cold interstellar medium (ISM; non-ionised gas and dust) from a galaxy. In addition, dust grains can be destroyed in a hot or turbulent medium. The adopted timescale of dust removal usually relies on uncertain theoretical estimates. It is tricky to track the dust removal, because usually dust is constantly replenished by consecutive generations of stars. Our objective is to measure observationally the timescale of dust removal. We here explore an approach to select galaxies which do have detectable amounts of dust and cold ISM but exhibit a low current dust production rate. Any decrease of the dust and gas content as a function of the age of such galaxies therefore must be attributed to processes governing the ISM removal. We used a sample of galaxies detected by Herschel in the far-infrared with visually assigned early-type morphology or spirals with red colours. We also obtained JCMT/SCUBA-2 observations for five of them. We discovered an exponential decline of the dust-to-stellar mass ratio with age, which we interpret as an evolutionary trend of dust removal from these galaxies. For the first time we directly measure the dust removal timescale in such galaxies to be tau=(2.5+-0.4) Gyr (the corresponding half-life time is (1.75+-0.25) Gyr). This quantity may be used in models in which it must be assumed a priori and cannot be derived. Any process which removes dust in these galaxies, such as dust grain destruction, cannot happen on shorter timescales. The timescale is comparable to the quenching timescales found in simulations for galaxies with similar stellar masses. The dust is likely of internal, not external origin. It was either formed in the past directly by supernovae, or from seeds produced by SNe and with grain growth in the ISM contributing substantially to the dust mass accumulation.

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 Atacama Large Millimetre/submillimeter Array (ALMA) Cycle-5 follow-up observations of two candidate [CII] emitters at z ~ 6 in the ALMA deep field in SSA22 (ADF22). The candidates were detected blindly in a Cycle-2 ALMA survey covering ~ 5 square arcmins, with a single tuning, along with two CO lines associated with galaxies at lower redshifts. Various tests suggested at least one of the two > 6-sigma [CII] candidates should be robust (Hayatsu et al. 2017). Nevertheless, our new, deeper observations recover neither candidate, demonstrating a higher contamination rate than expected. The cause of the spurious detections is under investigation but at present it remains unclarified.

We present the results of $^{12}$CO($J$=1-0) mosaicing observations of the cD galaxy NGC 1316 at kpc-resolution performed with the Morita Array of the Atacama Large Millimeter/submillimeter Array (ALMA). We reveal the detailed distribution of the molecular gas in the central region for the first time: a shell structure in the northwest, a barely resolved blob in the southeast of the center and some clumps between them. The total molecular gas mass obtained with a standard Milky-Way CO-to-H$_2$ conversion factor is $(5.62 \pm 0.53)\times10^8$ M$_\odot$, which is consistent with previous studies. The disturbed velocity field of the molecular gas suggests that the molecular gas is injected very recently ($<1$ Gyr) if it has an external origin and is in the process of settling into a rotating disk. Assuming that a low-mass gas-rich galaxy has accreted, the gas-to-dust ratio and H$_2$-to-HI ratio are unusually low ($\sim 28$) and high ($\sim 5.6$), respectively. To explain these ratios, additional processes should be taken into accounts such as an effective dust formation and conversion from atomic to molecular gas during the interaction. We also discuss the interaction between the nuclear jet and the molecular gas.

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.

We present the Atacama Large Millimeter/submillimeter Array (ALMA) detection of the [O III] 88 $\mu$m line and rest-frame 90 $\mu$m dust continuum emission in a Y-dropout Lyman break galaxy (LBG), MACS0416_Y1, lying behind the Frontier Field cluster MACS J0416.1-2403. This [O III] detection confirms the LBG with a spectroscopic redshift of $z = 8.3118 \pm 0.0003$, making this object one of the furthest galaxies ever identified spectroscopically. The observed 850 $\mu$m flux density of $137 \pm 26$ $\mu$Jy corresponds to a de-lensed total infrared (IR) luminosity of $L_{\rm IR} = (1.7 \pm 0.3) \times 10^{11} L_{\odot}$ if assuming a dust temperature of $T_{\rm dust} = 50$ K and an emissivity index of $\beta = 1.5$, yielding a large dust mass of $4 \times 10^6 M_{\odot}$. The ultraviolet-to-far IR spectral energy distribution modeling where the [O III] emissivity model is incorporated suggests the presence of a young ($\tau_{\rm age} \approx 4$ Myr), star-forming (SFR $\approx 60 M_{\odot}$ yr$^{-1}$), moderately metal-polluted ($Z \approx 0.2 Z_{\odot}$) stellar component with a mass of $M_{\rm star} = 3 \times 10^8 M_{\odot}$. An analytic dust mass evolution model with a single episode of star-formation does not reproduce the metallicity and dust mass in $\tau_{\rm age} \approx 4$ Myr, suggesting a pre-existing evolved stellar component with $M_{\rm star} \sim 3 \times 10^9 M_{\odot}$ and $\tau_{\rm age} \sim 0.3$ Gyr as the origin of the dust mass.

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.

We use the new release of the AKARI Far-Infrared all sky Survey matched with the NVSS radio database to investigate the local ($z<0.25$) far infrared-radio correlation (FIRC) of different types of extragalactic sources. To obtain the redshift information for the AKARI FIS sources we crossmatch the catalogue with the SDSS DR8. This also allows us to use emission line properties to divide sources into four categories: i) star-forming galaxies (SFGs), ii) composite galaxies (displaying both star-formation and active nucleus components), iii) Seyfert galaxies, and iv) low-ionization nuclear emission-line region (LINER) galaxies. We find that the Seyfert galaxies have the lowest FIR/radio flux ratios and display excess radio emission when compared to the SFGs. We conclude that FIRC can be used to separate SFGs and AGNs only for the most radio-loud objects.

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.

Various galaxy classification schemes have been developed so far to constrain the main physical processes regulating evolution of different galaxy types. In the era of a deluge of astrophysical information and recent progress in machine learning, a new approach to galaxy classification becomes imperative. We employ a Fisher Expectation-Maximization unsupervised algorithm working in a parameter space of 12 rest-frame magnitudes and spectroscopic redshift. The model (DBk) and the number of classes (12) were established based on the joint analysis of standard statistical criteria and confirmed by the analysis of the galaxy distribution with respect to a number of classes and their properties. This new approach allows us to classify galaxies based just on their redshifts and UV-NIR spectral energy distributions. The FEM unsupervised algorithm has automatically distinguished 12 classes: 11 classes of VIPERS galaxies and an additional class of broad-line AGNs. After a first broad division into blue, green and red categories we obtained a further sub-division into three red, three green, and five blue galaxy classes. The FEM classes follow the galaxy sequence from the earliest to the latest types that is reflected in their colours (which are constructed from rest-frame magnitudes used in classification procedure) but also their morphological, physical, and spectroscopic properties (not included in the classification scheme). We demonstrate that the members of each class share similar physical and spectral properties. In particular, we are able to find three different classes of red passive galaxy populations. Thus, we demonstrate the potential of an unsupervised approach to galaxy classification and we retrieve the complexity of galaxy populations at z~0.7, a task that usual simpler colour-based approaches cannot fulfil.

Techniques to classify galaxies solely based on photometry will be necessary for future large cosmology missions, such as Euclid or LSST. However, the precision of classification is always lower in photometric surveys and can be systematically biased with respect to classifications based upon spectroscopic data. We verified how precisely the detailed classification scheme introduced by Siudek et al. (2018, hereafter: S1) for galaxies at z~0.7 could be reproduced if only photometric data are available. We applied the Fisher Expectation-Maximization (FEM) unsupervised clustering algorithm to 54,293 VIPERS galaxies working in a parameter space of reliable photometric redshifts and 12 corresponding rest-frame magnitudes. The FEM algorithm distinguishes four main groups: (1) red, (2) green, (3) blue, and (4) outliers. Each group is further divided into 3, 3, 4, and 2 subclasses, respectively. The accuracy of reproducing galaxy classes using spectroscopic data is high: 92%, 84%, 96% for red, green, and blue classes, respectively, except for dusty star-forming galaxies. The presented verification of the photometric classification demonstrates that large photometric samples can be used to distinguish different galaxy classes at z > 0.5 with an accuracy provided so far only by spectroscopic data except for particular galaxy classes.

We want to investigate how planet formation is imprinted on stellar surface composition using up-to-date stellar evolution models. We simulate the evolution of pre-main-sequence stars as a function of the efficiency of heat injection during accretion, the deuterium mass fraction, and the stellar mass. For simplicity, we assume that planet formation leads to the late accretion of zero-metallicity gas, diluting the surface stellar composition as a function of the mass of the stellar outer convective zone. We adopt $150\,{\mathrm{M}_\oplus}(M_\star/\mathrm{M}_\odot)(Z/\mathrm{Z}_\odot)$ as an uncertain but plausible estimate of the mass of heavy elements that is not accreted by stars with giant planets, including our Sun. By combining our stellar evolution models to these estimates, we evaluate the consequences of planet formation on stellar surface composition. We show that after the first $\sim0.1$ Myr, the evolution of the convective zone follows classical evolutionary tracks within a factor of two in age. We find that planet formation should lead to a scatter in stellar surface composition that is larger for high-mass stars than for low-mass stars. We predict a spread in [Fe/H] of approximately $0.02$ dex for stars with $T_\mathrm{eff}\sim 5500\,$K, marginally compatible with differences in metallicities observed in some binary stars with planets. Stars with $T_\mathrm{eff}\geq 7000\,$K may show much larger [Fe/H] deficits, by 0.6 dex or more, compatible with the existence of refractory-poor $\lambda$ Boo stars. We also find that planet formation may explain the lack of refractory elements seen in the Sun as compared to solar twins, but only if the ice-to-rock ratio in the solar-system planets is less than $\approx0.4$ and planet formation began less than $\approx1.3$ Myr after the beginning of the formation of the Sun. (abbreviated)

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.

We present the dust mass function (DMF) of 15,750 galaxies with redshift $z< 0.1$, drawn from the overlapping area of the GAMA and {\it H-}ATLAS surveys. The DMF is derived using the density corrected $V_{\rm max}$ method, where we estimate $V_{\rm max}$ using: (i) the normal photometric selection limit ($pV_{\rm max}$) and (ii) a bivariate brightness distribution (BBD) technique, which accounts for two selection effects. We fit the data with a Schechter function, and find $M^{*}=(4.65\pm0.18)\times 10^{7}\,h^2_{70}\, M_{\odot}$, $\alpha=(1.22\pm 0.01)$, $\phi^{*}=(6.26\pm 0.28)\times 10^{-3}\,h^3_{70}\,\rm Mpc^{-3}\,dex^{-1}$. The resulting dust mass density parameter integrated down to $10^4\,M_{\odot}$ is $\Omega_{\rm d}=(1.11 \pm0.02)\times 10^{-6}$ which implies the mass fraction of baryons in dust is $f_{m_b}=(2.40\pm0.04)\times 10^{-5}$; cosmic variance adds an extra 7-17\,per\,cent uncertainty to the quoted statistical errors. Our measurements have fewer galaxies with high dust mass than predicted by semi-analytic models. This is because the models include too much dust in high stellar mass galaxies. Conversely, our measurements find more galaxies with high dust mass than predicted by hydrodynamical cosmological simulations. This is likely to be from the long timescales for grain growth assumed in the models. We calculate DMFs split by galaxy type and find dust mass densities of $\Omega_{\rm d}=(0.88\pm0.03)\times 10^{-6}$ and $\Omega_{\rm d}=(0.060\pm0.005)\times 10^{-6}$ for late-types and early-types respectively. Comparing to the equivalent galaxy stellar mass functions (GSMF) we find that the DMF for late-types is well matched by the GMSF scaled by $(8.07\pm0.35) \times 10^{-4}$.

Modified Dark Matter (MDM) is a phenomenological model of dark matter, inspired by gravitational thermodynamics, that naturally accounts for the universal acceleration constant observed in galactic rotation curve data; a critical acceleration related to the cosmological constant, $\Lambda$, appears as a phenomenological manifestation of MDM. We show that the resulting mass profiles, which are sensitve to $\Lambda$, are consistent with observations at the galactic and galaxy cluster scales. Our results suggest that dark matter mass profiles contain information about the cosmological constant in a non-trivial way.

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.

Modified dark matter (MDM) is a phenomenological model of dark matter, inspired by gravitational thermodynamics. For an accelerating Universe with positive cosmological constant ($\Lambda$), such phenomenological considerations lead to the emergence of a critical acceleration parameter related to $\Lambda$. Such a critical acceleration is an effective phenomenological manifestation of MDM, and it is found in correlations between dark matter and baryonic matter in galaxy rotation curves. The resulting MDM mass profiles, which are sensitive to $\Lambda$, are consistent with observational data at both the galactic and cluster scales. In particular, the same critical acceleration appears both in the galactic and cluster data fits based on MDM. Furthermore, using some robust qualitative arguments, MDM appears to work well on cosmological scales, even though quantitative studies are still lacking. Finally, we comment on certain non-local aspects of the quanta of modified dark matter, which may lead to novel non-particle phenomenology and which may explain why, so far, dark matter detection experiments have failed to detect dark matter particles.

Using the AKARI, Wide-field Infrared Survey Explorer (WISE), Infrared Astronomical Satellite (IRAS), Sloan Digital Sky Survey (SDSS) and Hubble Space Telescope (HST) data, we investigated the relation of polycyclic aromatic hydrocarbon (PAH) mass ($M_{\rm PAH}$), very small grain mass ($M_{\rm VSG}$), big grain mass ($M_{\rm BG}$) and stellar mass ($M_{\rm star}$) with galaxy merger for 55 star-forming galaxies at redshift $z<0.2$. Using the SDSS image at $z<0.1$ and the HST image at $z>0.1$, we divided the galaxies into merger galaxies and non-merger galaxies with the morphological parameter asymmetry $A$, and quantified merging stages of galaxies based on the morphological indicators, the second-order momentum of the brightest 20$\%$ region $M_{20}$ and the Gini coefficient. We find that $M_{\rm PAH}/M_{\rm BG}$ of merger galaxies tend to be lower than that of non-merger galaxies and there are no systematic differences of $M_{\rm VSG}/M_{\rm BG}$ and $M_{\rm BG}/M_{\rm star}$ between merger galaxies and non-merger galaxies. We find that galaxies with very low $M_{\rm PAH}/M_{\rm BG}$ seem to be merger galaxies at late stages. These results suggest that PAHs are partly destroyed at late stages of merging processes. Furthermore, we investigated $M_{\rm PAH}/M_{\rm BG}$ variations in radiation field intensity strength $G_0$ and the emission line ratio of $[{\rm O\,{\scriptsize I}}]\lambda 6300/{\rm H}\alpha$ which is a shock tracer for merger galaxies and find that $M_{\rm PAH}/M_{\rm BG}$ decreases with increasing both $G_0$ and $[{\rm O\,{\scriptsize I}}]/{\rm H}\alpha$. PAH destruction is likely to be caused by two processes; strong radiation fields and large-scale shocks during merging processes of galaxies.

Recent theoretical work has shown that the pre-main-sequence (PMS) evolution of stars is much more complex than previously envisioned. Instead of the traditional steady, one-dimensional solution, accretion may be episodic and not necessarily symmetrical, thereby affecting the energy deposited inside the star and its interior structure. Given this new framework, we want to understand what controls the evolution of accreting stars. We use the MESA stellar evolution code with various sets of conditions. In particular, we account for the (unknown) efficiency of accretion in burying gravitational energy into the protostar through a parameter, $\xi$, and we vary the amount of deuterium present. We confirm the findings of previous works that the evolution changes significantly with the amount of energy that is lost during accretion. We find that deuterium burning also regulates the PMS evolution. In the low-entropy accretion scenario, the evolutionary tracks in the H-R diagram are significantly different from the classical tracks and are sensitive to the deuterium content. A comparison of theoretical evolutionary tracks and observations allows us to exclude some cold accretion models ($\xi\sim 0$) with low deuterium abundances. We confirm that the luminosity spread seen in clusters can be explained by models with a somewhat inefficient injection of accretion heat. The resulting evolutionary tracks then become sensitive to the accretion heat efficiency, initial core entropy, and deuterium content. In this context, we predict that clusters with a higher D/H ratio should have less scatter in luminosity than clusters with a smaller D/H. Future work on this issue should include radiation-hydrodynamic simulations to determine the efficiency of accretion heating and further observations to investigate the deuterium content in star-forming regions. (abbrev.)

We report the identification of four millimeter line emitting galaxies with the Atacama Large Milli/submillimeter Array (ALMA) in SSA22 Field (ADF22). We analyze the ALMA 1.1 mm survey data, with an effective survey area of 5 arcmin$^2$, a frequency range of 253.1--256.8 and 269.1--272.8 GHz, angular resolution of 0".7 and RMS noise of 0.8 mJy beam$^{-1}$ at 36 km s$^{-1}$ velocity resolution. We detect four line emitter candidates with significance levels above $6 \sigma$. We identify one of the four sources as a CO(9-8) emitter at $z = 3.1$ in a member of the proto-cluster known in this field. Another line emitter with an optical counterpart is likely a CO(4-3) emitter at $z = 0.7$. The other two sources without any millimeter continuum or optical/near-infrared counterpart are likely to be [CII] emitter candidates at $z = 6.0$ and $6.5$. The equivalent widths of the [CII] candidates are consistent with those of confirmed high-redshift [CII] emitters and candidates, and are a factor of 10 times larger than that of the CO(9-8) emitter detected in this search. The [CII] luminosity of the candidates are $4-7 \times 10^8~\rm L_\odot$. The star formation rates (SFRs) of these sources are estimated to be $10-20~\rm M_{\odot}~yr^{-1}$ if we adopt an empirical [CII] luminosity - SFR relation. One of them has a relatively low-S/N ratio, but shows features characteristic of emission lines. Assuming that at least one of the two candidates is a [CII] emitter, we derive a lower limit of [CII]-based star formation rate density (SFRD) at $z~\sim~6$. The resulting value of $> 10^{-2}$ $\rm M_\odot yr^{-1} Mpc^{-3}$ is consistent with the dust-uncorrected UV-based SFRD. Future millimeter/submillimeter surveys can be used to detect a number of high redshift line emitters, with which to study the star formation history in the early Universe.

In the Solar System, interplanetary dust particles (IDPs) originating mainly from asteroid collisions and cometary activities drift to the Earth orbit due to the Poynting-Robertson drag. We analyzed the thermal emission from IDPs that was observed by the first Japanese infrared astronomical satellite, AKARI. The observed surface brightness in the trailing direction of the Earth orbit is 3.7% greater than that in the leading direction in the $9{\rm \mu m}$ band and 3.0% in the $18{\rm \mu m}$ band. In order to reveal dust properties causing the leading-trailing surface brightness asymmetry, we numerically integrated orbits of the Sun, the Earth, and a dust particle as a restricted three-body problem including radiation from the Sun. The initial orbits of particles are determined according to the orbits of main-belt asteroids or Jupiter-family comets. The orbital trapping in mean motion resonances results in a significant leading-trailing asymmetry so that intermediate sized dust (~10-100${\rm \mu m}$) produces a greater asymmetry than the zodiacal light has. The leading-trailing surface brightness difference integrated over the size distribution of the asteroidal dust is obtained to be the values of 27.7% and 25.3% in the $9{\rm \mu m}$ and $18{\rm \mu m}$ bands, respectively. In contrast, the brightness difference for cometary dust is calculated as the values of 3.6% and 3.1% in the $9{\rm \mu m}$ and $18{\rm \mu m}$ bands, respectively, if the maximum dust radius is set to be $s_{\rm max} = 3000{\rm \mu m}$. Taking into account these values and their errors, we conclude that the contribution of asteroidal dust to the zodiacal infrared emission is less than ~10%, while cometary dust of the order of 1 mm mainly accounts for the zodiacal light in infrared.

The aim of this work is to characterize physical properties of Ultra Luminous Infrared Galaxies (ULIRGs) and Luminous Infrared Galaxies (LIRGs) detected in the far-infrared (FIR) 90um band in the AKARI Deep Field-South (ADF-S) survey. In particular, we want to estimate the AGN contribution to the [U]LIRGs' infrared emission and which types of AGNs are related to their activity. We examined 69 galaxies at z>0.05 detected at 90um by the AKARI satellite in the ADF-S, with optical counterparts and spectral coverage from the ultraviolet to the FIR. We used two independent spectral energy distribution fitting codes: one fitting the SED from FIR to FUV (CIGALE) and gray-body + power spectrum fit for the infrared part of the spectra (CMCIRSED) in order to identify a subsample of [U]LIRGs, and to estimate their properties. Based on the CIGALE SED fitting, we have found that [U]LIRGs selected at the 90um AKARI band compose ~56% of our sample (we found 17 ULIRGs and 22 LIRGs, spanning over the redshift range 0.06<z<1.23). Their physical parameters, such as stellar mass, star formation rate (SFR), and specific SFR are consistent with the ones found for other samples selected at IR wavelengths. We have detected a significant AGN contribution to the MIR luminosity for 63% of LIRGs and ULIRGs. Our LIRGs contain Type 1, Type 2, and intermediate types of AGN, whereas for ULIRGs, a majority (more than 50%) of AGN emission originates from Type 2 AGNs. The temperature--luminosity and temperature--mass relations for the dust component of ADF--S LIRGs and ULIRGs indicate that these relations are shaped by the dust mass and not by the increased dust heating. We conclude that LIRGs contain Type 1, Type 2, and intermediate types of AGNs, with an AGN contribution to the MIR emission at the median level of 13+/-3%, whereas the majority of ULIRGs contain Type 2 AGNs, with a median AGN fraction equal to 19+/-8%.

We trace the evolution and the star formation history of passive galaxies, using a subset of the VIMOS Public Extragalactic Redshift Survey (VIPERS). We extracted from the VIPERS survey a sample of passive galaxies in the redshift range 0.4<z<1.0 and stellar mass range 10<$log(M_{star}/M_{\odot})$<12. The sample was selected using an evolving cut in the rest-frame U-V color distribution and additional quality-ensuring cuts. We use the stacked spectra to measure the 4000$\AA$ break (D4000) and the $H\delta$ Lick index ($H\delta_{A}$) with high precision. We compare the results with a grid of synthetic spectra to constrain the star formation epochs of these galaxies. We characterize the formation redshift-stellar mass relation for intermediate-redshift passive galaxies. We find that at $z\sim1$ stellar populations in low-mass passive galaxies are younger than in high-mass passive galaxies, similarly to what is observed at the present epoch. Over the full analyzed redshift and stellar mass range, the $D4000$ index increases with redshift, while $H\delta_{A}$ gets lower. This implies that the stellar populations are getting older with increasing stellar mass. Comparison to the spectra of passive galaxies in the SDSS survey shows that the shape of the relations of $D4000$, and $H\delta_{A}$ with stellar mass has not changed significantly with redshift. Assuming a single burst formation, this implies that high-mass passive galaxies formed their stars at $z_{form}\sim2$, while low-mass galaxies formed their main stellar population more recently, at $z_{form}\sim1$. The consistency of these results, obtained using two independent estimator of the formation redshift ($D4000$ and $H\delta_{A}$), further strengthens a scenario in which star formation proceeds from higher- to lower-mass systems as time passes, i.e. what has become known as the 'downsizing' picture.

The DUNE/LBNF program aims to address key questions in neutrino physics and astroparticle physics. Realizing DUNE's potential to reconstruct low-energy particles in the 10-100 MeV energy range will bring significant benefits for all DUNE's science goals. In neutrino physics, low-energy sensitivity will improve neutrino energy reconstruction in the GeV range relevant for the kinematics of DUNE's long-baseline oscillation program. In astroparticle physics, low-energy capabilities will make DUNE's far detectors the world's best apparatus for studying the electron-neutrino flux from a supernova. This will open a new window to unrivaled studies of the dynamics and neutronization of a star's central core in real time, the potential discovery of the neutrino mass hierarchy, provide new sensitivity to physics beyond the Standard Model, and evidence of neutrino quantum-coherence effects. The same capabilities will also provide new sensitivity to `boosted dark matter' models that are not observable in traditional direct dark matter detectors.

We aim to use the a new and improved version of AKARI all sky survey catalogue of far-infrared sources to recalibrate the formula to derive the total infrared luminosity. We cross-match the faint source catalogue (FSC) of IRAS with the new AKARI-FIS and obtained a sample of 2430 objects. Then we calculate the total infrared (TIR) luminosity $L_{\textrm{TIR}}$ from the Sanders at al. (1996) formula and compare it with total infrared luminosity from AKARI FIS bands to obtain new coefficients for the general relation to convert FIR luminosity from AKARI bands to the TIR luminosity.

We explore the evolution of the statistical distribution of galaxy morphological properties and colours over the redshift range $0.5<z<1$, combining high-quality imaging data from the CFHT Legacy Survey with the large number of redshifts and extended photometry from the VIPERS survey. Galaxy structural parameters are measured by fitting S\'ersic profiles to $i$-band images and then combined with absolute magnitudes, colours and redshifts, to trace the evolution in a multi-parameter space. We analyse, using a new method, the combination of colours and structural parameters of early- and late-type galaxies in luminosity--redshift space. We found that both the rest-frame colour distributions in the (U-B) vs. (B-V) plane and the S\'ersic index distributions are well fitted by a sum of two Gaussians, with a remarkable consistency of red-spheroidal and blue-disky galaxy populations, over the explored redshift ($0.5<z<1$) and luminosity ($-1.5<B-B_*<1.0$) ranges. The combination of the UBV rest-frame colour and S\'ersic index $n$ as a function of redshift and luminosity allows us to present the structure of early- and late-type galaxies and their evolution. We found that early type galaxies display only a slow change of their concentrations since $z\sim1$; it is already established by $z\sim1$ and depends much more strongly on their luminosities. In contrast, late-type galaxies get clearly more concentrated with cosmic time since $z\sim1$, with only little evolution in colour, which remains dependent mainly on their luminosity. This flipped luminosity (mass) and redshift dependence likely reflects different evolutionary tracks of early- and late-type galaxies before and after $z\sim1$. The combination of rest-frame colours and S\'ersic index $n$ as a function of redshift and luminosity leads to a precise statistical description of the structure of galaxies and their evolution.

Protoplanetary disks with non-axisymmetric structures have been observed. The Rossby wave instability (RWI) is considered as one of the origins of the non-axisymmetric structures. We perform linear stability analyses of the RWI in barotropic flow using four representative types of the background flow on a wide parameter space. We find that the co-rotation radius is located at the background vortensity minimum with large concavity if the system is marginally stable to the RWI, and this allows us to check the stability against the RWI easily. We newly derive the necessary and sufficient condition for the onset of the RWI in semi-analytic form. We discuss the applicability of the new condition in realistic systems and the physical nature of the RWI.

A giant planet creates a gap in a protoplanetary disk, which might explain the observed gaps in protoplanetary disks. The width and depth of the gaps depend on the planet mass and disk properties. We have performed two--dimensional hydrodynamic simulations for various planet masses, disk aspect ratios and viscosities, to obtain an empirical formula for the gap width. The gap width is proportional to the square root of the planet mass, -3/4 power of the disk aspect ratio and -1/4 power of the viscosity. This empirical formula enables us to estimate the mass of a planet embedded in the disk from the width of an observed gap. We have applied the empirical formula for the gap width to the disk around HL~Tau, assuming that each gap observed by ALMA observations is produced by planets, and discussed the planet masses within the gaps. The estimate of planet masses from the gap widths is less affected by the observational resolution and dust filtration than that from the gap depth.

Formation and evolution of galaxies have been a central driving force in the studies of galaxies and cosmology. Recent studies provided a global picture of cosmic star formation history. However, what drives the evolution of star formation activities in galaxies has long been a matter of debate. The key factor of the star formation is the transition of hydrogen from atomic to molecular state, since the star formation is associated with the molecular phase. This transition is also strongly coupled with chemical evolution, because dust grains, i.e., tiny solid particles of heavy elements, play a critical role in molecular formation. Therefore, a comprehensive understanding of neutral-molecular gas transition, star formation and chemical enrichment is necessary to clarify the galaxy formation and evolution. Here we present the activity of SKA-JP galaxy evolution sub-science working group (subSWG) Our activity is focused on three epochs: z \sim 0, 1, and z > 3. At z \sim 0, we try to construct a unified picture of atomic and molecular hydrogen through nearby galaxies in terms of metallicity and other various ISM properties. Up to intermediate redshifts z \sim 1, we explore scaling relations including gas and star formation properties, like the main sequence and the Kennicutt-Schmidt law of star forming galaxies. To connect the global studies with spatially-resolved investigations, such relations will be plausibly a viable way. For high redshift objects, the absorption lines of HI 21-cm line will be a very promising observable to explore the properties of gas in galaxies. By these studies, we will surely witness a real revolution in the studies of galaxies by SKA.

We present the Spitzer-IRAC/MIPS Extragalactic survey (SIMES) in the South Ecliptic Pole (SEP) field. The large area covered (7.7 deg$^2$), together with one of the lowest Galactic cirrus emissions in the entire sky and a very extensive coverage by Spitzer, Herschel, Akari, and GALEX, make the SIMES field ideal for extragalactic studies. The elongated geometry of the SIMES area ($\approx$4:1), allowing for a significant cosmic variance reduction, further improves the quality of statistical studies in this field. Here we present the reduction and photometric measurements of the Spitzer/IRAC data. The survey reaches a depth of 1.93 and 1.75 $\mu$Jy (1$\sigma$) at 3.6 and 4.5 $\mu$m, respectively. We discuss the multiwavelength IRAC--based catalog, completed with optical, mid-- and far--IR observations. We detect 341,000 sources with F$_{3.6\mu m} \geq 3\sigma$. Of these, 10% have an associated 24 $\mu$m counterpart, while 2.7% have an associated SPIRE source. We release the catalog through the NASA/IPAC Infrared Science Archive (IRSA). Two scientific applications of these IRAC data are presented in this paper: first we compute integral number counts at 3.6 $\mu$m. Second, we use the [3.6]--[4.5] color index to identify galaxy clusters at z$>$1.3. We select 27 clusters in the full area, a result consistent with previous studies at similar depth.

Modified dark matter (MDM, formerly known as MoNDian dark matter) is a phenomenological model of dark matter, inspired by quantum gravity. We review the construction of MDM by generalizing entropic gravity to de-Sitter space as is appropriate for an accelerating universe (in accordance with the Lambda-CDM model). Unlike cold dark matter models, the MDM mass profile depends on the baryonic mass. We successfully fit the rotation curves to a sample of 30 local spiral galaxies with a single free parameter (viz., the mass-to-light ratio for each galaxy). We show that dynamical and observed masses agree in a sample of 93 galactic clusters. We also comment on strong gravitational lensing in the context of MDM.

We discuss the possibility that the cold dark matter mass profiles contain information on the cosmological constant, and that such information constrains the nature of cold dark matter (CDM). We call this approach Modified Dark Matter (MDM). In particular, we examine the ability of MDM to explain the observed mass profiles of 13 galaxy clusters. Using general arguments from gravitational thermodynamics, we provide a theoretical justification for our MDM mass profile and successfully compare it to the NFW mass profiles both on cluster and galactic scales. Our results suggest that indeed the CDM mass profiles contain information about the cosmological constant in a non-trivial way.

Layered accretion is one of the inevitable ingredients in protoplanetary disks when disk turbulence is excited by magnetorotational instabilities (MRIs). In the accretion, disk surfaces where MRIs fully operate have a high value of disk accretion rate ($\dot{M}$), while the disk midplane where MRIs are generally quenched ends up with a low value of $\dot{M}$. Significant progress on understanding MRIs has recently been made by a number of dedicated MHD simulations, which requires improvement of the classical treatment of $\alpha$ in 1D disk models. To this end, we obtain a new expression of $\alpha$ by utilizing an empirical formula that is derived from recent MHD simulations of stratified disks with Ohmic diffusion. It is interesting that this new formulation can be regarded as a general extension of the classical $\alpha$. Armed with the new $\alpha$, we perform a linear stability analysis of protoplanetary disks that undergo layered accretion, and find that a viscous instability can occur around the outer edge of dead zones. Disks become stable in using the classical $\alpha$. We identify that the difference arises from $\Sigma-$dependence of $\dot{M}$; whereas $\Sigma$ is uniquely determined for a given value of $\dot{M}$ in the classical approach, the new approach leads to $\dot{M}$ that is a multi-valued function of $\Sigma$. We confirm our finding both by exploring a parameter space as well as by performing the 1D, viscous evolution of disks. We finally discuss other non-ideal MHD effects that are not included in our analysis, but may affect our results.

We investigate the dust and gas distribution in the disk around HD 142527 based on ALMA observations of dust continuum, 13CO(3-2), and C18O(3-2) emission. The disk shows strong azimuthal asymmetry in the dust continuum emission, while gas emission is more symmetric. In this paper, we investigate how gas and dust are distributed in the dust-bright northern part of the disk and in the dust-faint southern part. We construct two axisymmetric disk models. One reproduces the radial profiles of the continuum and the velocity moments 0 and 1 of CO lines in the north and the other reproduces those in the south. We have found that the dust is concentrated in a narrow ring having ~50AU width (in FWHM; w_d=30AU in our parameter definition) located at ~170-200AU from the central star. The dust particles are strongly concentrated in the north. We have found that the dust surface density contrast between the north and south amounts to ~70. Compared to the dust, the gas distribution is more extended in the radial direction. We find that the gas component extends at least from ~100AU to ~250AU from the central star, and there should also be tenuous gas remaining inside and outside of these radii. The azimuthal asymmetry of gas distribution is much smaller than dust. The gas surface density differs only by a factor of ~3-10 between the north and south. Hence, gas-to-dust ratio strongly depends on the location of the disk: ~30 at the location of the peak of dust distribution in the south and ~3 at the location of the peak of dust distribution in the north. Despite large uncertainties, the overall gas-to-dust ratio is inferred to be ~10-30, indicating that the gas depletion may have already been under way.

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

Concerted effort is currently ongoing to open up the Epoch of Reionization (EoR) ($z\sim$15-6) for studies with IR and radio telescopes. Whereas IR detections have been made of sources (Lyman-$\alpha$ emitters, quasars and drop-outs) in this redshift regime in relatively small fields of view, no direct detection of neutral hydrogen, via the redshifted 21-cm line, has yet been established. Such a direct detection is expected in the coming years, with ongoing surveys, and could open up the entire universe from $z\sim$6-200 for astrophysical and cosmological studies, opening not only the EoR, but also its preceding Cosmic Dawn ($z\sim$30-15) and possibly even the later phases of the Dark Ages ($z\sim$200-30). All currently ongoing experiments attempt statistical detections of the 21-cm signal during the EoR, with limited signal-to-noise. Direct imaging, except maybe on the largest (degree) scales at lower redshifts, as well as higher redshifts will remain out of reach. The Square Kilometre Array(SKA) will revolutionize the field, allowing direct imaging of neutral hydrogen from scales of arc-minutes to degrees over most of the redshift range $z\sim$6-28 with SKA1-LOW, and possibly even higher redshifts with the SKA2-LOW. In this SKA will be unique, and in parallel provide enormous potential of synergy with other upcoming facilities (e.g. JWST). In this chapter we summarize the physics of 21-cm emission, the different phases the universe is thought to go through, and the observables that the SKA can probe, referring where needed to detailed chapters in this volume (Abridged).

Coronagraphic imagery of the circumstellar disk around HD 169142 in H-band polarized intensity (PI) with Subaru/HiCIAO is presented. The emission scattered by dust particles at the disk surface in 0.2" <= r <= 1.2", or 29 <= r <= 174 AU, is successfully detected. The azimuthally-averaged radial profile of the PI shows a double power-law distribution, in which the PIs in r=29-52 AU and r=81.2-145 AU respectively show r^{-3}-dependence. These two power-law regions are connected smoothly with a transition zone (TZ), exhibiting an apparent gap in r=40-70 AU. The PI in the inner power-law region shows a deep minimum whose location seems to coincide with the point source at \lambda = 7 mm. This can be regarded as another sign of a protoplanet in TZ. The observed radial profile of the PI is reproduced by a minimally flaring disk with an irregular surface density distribution or with an irregular temperature distribution or with the combination of both. The depletion factor of surface density in the inner power-law region (r< 50 AU) is derived to be <= 0.16 from a simple model calculation. The obtained PI image also shows small scale asymmetries in the outer power-law region. Possible origins for these asymmetries include corrugation of the scattering surface in the outer region, and shadowing effect by a puffed up structure in the inner power-law region.

A giant planet embedded in a protoplanetary disk forms a gap. An analytic relationship among the gap depth, planet mass $M_{p}$, disk aspect ratio $h_p$, and viscosity $\alpha$ has been found recently, and the gap depth can be written in terms of a single parameter $K= (M_{p}/M_{\ast})^2 h_p^{-5} \alpha^{-1}$. We discuss how observed gap features can be used to constrain the disk and/or planet parameters based on the analytic formula for the gap depth. The constraint on the disk aspect ratio is critical in determining the planet mass so the combination of the observations of the temperature and the image can provide a constraint on the planet mass. We apply the formula for the gap depth to observations of HL~Tau and HD~169142. In the case of HL~Tau, we propose that a planet with $\gtrsim 0.3$ is responsible for the observed gap at $30$~AU from the central star based on the estimate that the gap depth is $\lesssim 1/3$. In the case of HD~169142, the planet mass that causes the gap structure recently found by VLA is $\gtrsim 0.4 M_J$. We also argue that the spiral structure, if observed, can be used to estimate the lower limit of the disk aspect ratio and the planet mass.

Current constraints on dark matter density profiles from weak lensing are typically limited to radial scales greater than 50-100 kpc. In this paper, we explore the possibility of probing the very inner regions of galaxy/halo density profiles by measuring stacked weak lensing on scales of only a few tens of kpc. Our forecasts focus on scales smaller than the equality radius (Req) where the stellar component and the dark matter component contribute equally to the lensing signal. We compute the evolution of Req as a function of lens stellar mass and redshift and show that Req=7-34 kpc for galaxies with the stellar mass of 10^{9.5}-10^{11.5} solar masses. Unbiased shear measurements will be challenging on these scales. We introduce a simple metric to quantify how many source galaxies overlap with their neighbours and for which shear measurements will be challenging. Rejecting source galaxies with close-by companions results in about a 20 per cent decrease in the overall source density. Despite this decrease, we show that Euclid and WFIRST will be able to constrain galaxy/halo density profiles at Req with signal-to-noise ratio >20 for the stellar mass of >10^{10} solar masses. Weak lensing measurements at Req, in combination with stellar kinematics on smaller scales, will be a powerful means by which to constrain both the inner slope of the dark matter density profile as well as the mass and redshift dependence of the stellar initial mass function.

The gap formation induced by a giant planet is important in the evolution of the planet and the protoplanetary disc. We examine the gap formation by a planet with a new formulation of one-dimensional viscous discs which takes into account the deviation from Keplerian disc rotation due to the steep gradient of the surface density. This formulation enables us to naturally include the Rayleigh stable condition for the disc rotation. It is found that the derivation from Keplerian disc rotation promotes the radial angular momentum transfer and makes the gap shallower than in the Keplerian case. For deep gaps, this shallowing effect becomes significant due to the Rayleigh condition. In our model, we also take into account the propagation of the density waves excited by the planet, which widens the range of the angular momentum deposition to the disc. The effect of the wave propagation makes the gap wider and shallower than the case with instantaneous wave damping. With these shallowing effects, our one-dimensional gap model is consistent with the recent hydrodynamic simulations.

When and how galaxy morphology such as disk and bulge seen in the present-day universe emerged is still not clear. In the universe at $z\gtrsim 2$, galaxies with various morphology are seen, and star-forming galaxies at $z\sim2$ show an intrinsic shape of bar-like structure. Then, when did round disk structure form? Here we take a simple and straightforward approach to see the epoch when a round disk galaxy population emerged by constraining the intrinsic shape statistically based on apparent axial ratio distribution of galaxies. We derived the distributions of the apparent axial ratios in the rest-frame optical light ($\sim 5000$ \AA) of star-forming main sequence galaxies at $2.5>z>1.4$, $1.4>z>0.85$, and $0.85>z>0.5$, and found that the apparent axial ratios of them show peaky distributions at $z\gtrsim0.85$, while a rather flat distribution at the lower redshift. By using a tri-axial model ($A>B>C$) for the intrinsic shape, we found the best-fit models give the peaks of the $B/A$ distribution of $0.81\pm0.04$, $0.84\pm0.04$, and $0.92\pm0.05$ at $2.5>z>1.4$, $1.4>z>0.85$, and $0.85>z>0.5$, respectively. The last value is close to the local value of 0.95. Thickness ($C/A$) is $\sim0.25$ at all the redshifts and is close to the local value (0.21). The results indicate the shape of the star-forming galaxies in the main sequence changes gradually, and the round disk is established at around $z\sim0.9$. Establishment of the round disk may be due to a cease of violent interaction of galaxies or a growth of a bulge and/or a super-massive black hole resides at the center of a galaxy which dissolves the bar structure.

The discoveries of huge amounts of dust and unusual extinction curves in high-redshift quasars (z > 4) cast challenging issues on the origin and properties of dust in the early universe. In this Letter, we investigate the evolutions of dust content and extinction curve in a high-z quasar, based on the dust evolution model taking account of grain size distribution. First, we show that the Milky-Way extinction curve is reproduced by introducing a moderate fraction (~0.2) of dense molecular-cloud phases in the interstellar medium for a graphite-silicate dust model. Then we show that the peculier extinction curves in high-z quasars can be explained by taking a much higher molecular-cloud fraction (>0.5), which leads to more efficient grain growth and coagulation, and by assuming amorphous carbon instead of graphite. The large dust content in high-z quasar hosts is also found to be a natural consequence of the enhanced dust growth. These results indicate that grain growth and coagulation in molecular clouds are key processes that can increase the dust mass and change the size distribution of dust in galaxies, and that, along with a different dust composition, can contribute to shape the extinction curve.

We study the time evolution of a large-scale magnetic flux threading an accretion disk. Induction equation of the mean poloidal field is solved under the standard viscous disk model. Magnetic flux evolution is controlled by the two timescales: One is the timescale of the inward advection of the magnetic flux, tau_{adv}. This is induced by the dragging of the flux by the accreting gas. The other is the outward diffusion timescale of the magnetic flux tau_{dif}. We consider diffusion due to the Ohmic resistivity. These timescales can be significantly different from the disk viscous timescale tau_{disk}. The behaviors of the magnetic flux evolution is quite different depending on the magnitude relationship of the timescales tau_{adv}, \tau_{dif}, and tau_{disk}. The most interesting phenomena occurs when tau_{adv} << tau_{dif}, tau_{disk}. In such a case, the magnetic flux distribution approaches a quasi-steady profile much faster than the viscous evolution of the gas disk, and also the magnetic flux has been tightly bundled to the inner part of the disk. In the inner part, although the poloidal magnetic field becomes much stronger than the interstellar magnetic field, the field strength is limited to the maximum value that is analytically given by our previous work (Okuzumi et al. 2014, ApJ, 785, 127). We also find a condition for that the initial large magnetic flux, which is a fossil of the magnetic field dragging during the early phase of star formation, survives for a duration in which significant gas disk evolution proceeds.

Galaxies' rest-frame ultraviolet (UV) properties are often used to directly infer the degree to which dust obscuration affects the measurement of star formation rates. While much recent work has focused on calibrating dust attenuation in galaxies selected at rest-frame ultraviolet wavelengths, locally and at high-$z$, here we investigate attenuation in dusty, star-forming galaxies (DSFGs) selected at far-infrared wavelengths. By combining multiwavelength coverage across 0.15--500\,$\mu$m in the COSMOS field, in particular making use of {\it Herschel} imaging, and a rich dataset on local galaxies, we find a empirical variation in the relationship between rest-frame UV slope ($\beta$) and ratio of infrared-to-ultraviolet emission ($L_{\rm IR}/L_{\rm UV}\equiv\,IRX$) as a function of infrared luminosity, or total star formation rate, SFR. Both locally and at high-$z$, galaxies above SFR$\gt$50\,M$_\odot$\,yr$^{-1}$ deviate from the nominal $IRX-\beta$ relation towards bluer colors by a factor proportional to their increasing IR luminosity. We also estimate contamination rates of DSFGs on high-$z$ dropout searches of $\ll1$\%\ at $z\lt4-10$, providing independent verification that contamination from very dusty foreground galaxies is low in LBG searches. Overall, our results are consistent with the physical interpretation that DSFGs, e.g. galaxies with $>50$\,M$_\odot$\,yr$^{-1}$, are dominated at all epochs by short-lived, extreme burst events, producing many young O and B stars that are primarily, yet not entirely, enshrouded in thick dust cocoons. The blue rest-frame UV slopes of DSFGs are inconsistent with the suggestion that most DSFGs at $z\sim2$ exhibit steady-state star formation in secular disks.

MONDian dark matter (MDM) is a new form of dark matter quantum that naturally accounts for Milgrom's scaling, usually associated with modified Newtonian dynamics (MOND), and theoretically behaves like cold dark matter (CDM) at cluster and cosmic scales. In this paper, we provide the first observational test of MDM by fitting rotation curves to a sample of 30 local spiral galaxies (z approximately 0.003). For comparison, we also fit the galactic rotation curves using MOND, and CDM. We find that all three models fit the data well. The rotation curves predicted by MDM and MOND are virtually indistinguishable over the range of observed radii (~1 to 30 kpc). The best-fit MDM and CDM density profiles are compared. We also compare with MDM the dark matter density profiles arising from MOND if Milgrom's formula is interpreted as Newtonian gravity with an extra source term instead of as a modification of inertia. We find that discrepancies between MDM and MOND will occur near the center of a typical spiral galaxy. In these regions, instead of continuing to rise sharply, the MDM mass density turns over and drops as we approach the center of the galaxy. Our results show that MDM, which restricts the nature of the dark matter quantum by accounting for Milgrom's scaling, accurately reproduces observed rotation curves.

We analytically calculate the marginally stable surface density profile for rotational instability of protoplanetary disks. The derived profile can be utilized for considering the region in a rotating disk where radial pressure gradient force is comparable to the gravitational force, such as an inner edge, steep gaps or bumps and an outer region of the disk. In this paper we especially focus on the rotational instability in the outer region of disks. We find an protoplanetary disk with a surface density profile of similarity solution becomes rotationally unstable at a certain radius, depending on its temperature profile and a mass of the central star. If the temperature is relatively low and the mass of the central star is high, disks have rotationally stable similarity profiles. Otherwise, deviation from the similarity profiles of surface density could be observable, using facilities with high sensitivity, such as ALMA.

We carried out 12CO(J=2-1) observations toward three star-forming galaxies on the main sequence at z~1.4 with the Nobeyama 45m radio telescope. These galaxies are detected with Spitzer/MIPS in 24 um, Herschel/SPIRE in 250 um and 350 um, and their gas metallicity, derived from optical emission line ratios based on near infrared spectroscopic observations, is close to the solar metallicity. Although weak signal-like features of CO were seen, we could not detect significant CO emission. The dust mass and the upper limits on the molecular gas mass are (3.4-6.7) x 10^{8} Msun and (9.7-14) x 10^{10} Msun, respectively. The upper limits on the gas-to-dust ratios at z~1.4 are 150-410 which are comparable to the gas-to-dust ratios in local galaxies with similar gas metallicity. A line stacking analysis enables us to detect a significant CO emission and to derive an average molecular gas mass of 1.3 x 10^{11} Msun and gas-to-dust ratio of 250. This gas-to-dust ratio is also near that in local galaxies with solar metallicity. These results suggest that the gas-to-dust ratio in star-forming galaxies with solar metallicity does not evolve significantly up to z~1.4. By comparing to a theoretical calculation, a rapid increase of the dust mass in an earlier epoch of galaxy evolution is suggested.

Large-scale magnetic fields are key ingredients of magnetically driven disk accretion. We study how large-scale poloidal fields evolve in accretion disks, with the primary aim of quantifying the viability of magnetic accretion mechanisms in protoplanetary disks. We employ a kinematic mean-field model for poloidal field transport and focus on steady states where inward advection of a field balances with outward diffusion due to effective resistivities. We analytically derive the steady-state radial distribution of poloidal fields in highly conducting accretion disks. The analytic solution reveals an upper limit on the strength of large-scale vertical fields attainable in steady states. Any excess poloidal field will be diffused away within a finite time, and we demonstrate this with time-dependent numerical calculations of the mean-field equations. We apply this upper limit to large-scale vertical fields threading protoplanetary disks. We find that the maximum attainable strength is about 0.1 G at 1 AU, and about 1 mG at 10 AU from the central star. When combined with recent magnetic accretion models, the maximum field strength translates into the maximum steady-state accretion rate of $\sim 10^{-7} M_\odot {\rm yr}^{-1}$, in agreement with observations. We also find that the maximum field strength is ~ 1 kG at the surface of the central star provided that the disk extends down to the stellar surface. This implies that any excess stellar poloidal field of strength >~ kG can be transported to the surrounding disk. This might in part resolve the magnetic flux problem in star formation.

We present a new version of a semi-analytic model of cosmological galaxy formation, incorporating a star formation law with a feedback depending on the galaxy-scale mean dust opacity and metallicity, motivated by recent observations of star formation in nearby galaxies and theoretical considerations. This new model is used to investigate the effect of such a feedback on shaping the galaxy luminosity function and its evolution. Star formation activity is significantly suppressed in dwarf galaxies by the new feedback effect, and the faint-end slope of local luminosity functions can be reproduced with a reasonable strength of supernova feedback, which is in contrast to the previous models that require a rather extreme strength of supernova feedback. Our model can also reproduce the early appearance of massive galaxies manifested in the bright-end of high redshift K-band luminosity functions. Though some of the previous models also succeeded in reproducing this, they assumed a star formation law depending on the galaxy-scale dynamical time, which is not supported by observations. We argue that the feedback depending on dust opacity (or metal column density) is essential, rather than that simply depending on gas column density, to get these results.

We investigate the evolution of extinction curves in galaxies based on our evolution model of grain size distribution. In this model, we considered various processes: dust formation by SNe II and AGB stars, dust destruction by SN shocks in the ISM, metal accretion onto the surface of grains (referred to as grain growth), shattering and coagulation. We find that the extinction curve is flat in the earliest stage of galaxy evolution. As the galaxy is enriched with dust, shattering becomes effective to produce a large abundance of small grains ($a \la 0.01\;\mu$m). Then, grain growth becomes effective at small grain radii, forming a bump at $a \sim 10^{-3}${--}$10^{-2}\;\mu$m on the grain size distribution. Consequently, the extinction curve at ultraviolet (UV) wavelengths becomes steep, and a bump at $1/\lambda \sim 4.5\;\mu{\rm m}^{-1}\;(\lambda: \mbox{wavelength})$ on the extinction curve becomes prominent. Once coagulation becomes effective, the extinction curves become flatter, but the UV extinction remains overproduced when compared with the Milky Way extinction curve. This discrepancy can be resolved by introducing a stronger contribution of coagulation. Consequently, an interplay between shattering and coagulation could be important to reproduce the Milky Way extinction curve.

This paper analyses the behaviour of the gas-to-dust mass ratio (G/D) of local Universe galaxies over a large metallicity range. We combine three samples: the Dwarf Galaxy Survey, the KINGFISH survey and a subsample from Galametz et al. (2011) totalling 126 galaxies, covering a 2 dex metallicity range, with 30% of the sample with 12+log(O/H) < 8.0. The dust masses are homogeneously determined with a semi-empirical dust model, including submm constraints. The atomic and molecular gas masses are compiled from the literature. Two XCO are used to estimate molecular gas masses: the Galactic XCO, and a XCO depending on the metallicity (as Z^{-2}). Correlations with morphological types, stellar masses, star formation rates and specific star formation rates are discussed. The trend between G/D and metallicity is empirically modelled using power-laws (slope of -1 and free) and a broken power-law. We compare the evolution of the G/D with predictions from chemical evolution models. We find that out of the five tested galactic parameters, metallicity is the galactic property driving the observed G/D. The G/D versus metallicity relation cannot be represented by a power-law with a slope of -1 over the whole metallicity range. The observed trend is steeper for metallicities lower than ~ 8.0. A large scatter is observed in the G/D for a given metallicity, with a dispersion of 0.37 dex in metallicity bins of ~0.1 dex. The broken power-law reproduces best the observed G/D and provides estimates of the G/D that are accurate to a factor of 1.6. The good agreement of the G/D and its scatter with the three tested chemical evolution models shows that the scatter is intrinsic to galactic properties, reflecting the different star formation histories, dust destruction efficiencies, dust grain size distributions and chemical compositions across the sample. (abriged)

The main aim of this work is the characterization of physical properties of galaxies detected in the far infrared (FIR) in the AKARI Deep Field-South (ADF-S) survey. Starting from a catalog of the 1 000 brightest ADF-S sources in the WIDE-S (90$\mu$m) AKARI band, we constructed a subsample of galaxies with spectral coverage from the ultraviolet to the far infrared. We then analyzed the multiwavelength properties of this 90$\mu$m selected sample of galaxies. For galaxies without known spectroscopic redshifts we computed photometric redshifts using the codes Photometric Analysis for Redshift Estimate (Le PHARE) and Code Investigating GALaxy Emission (CIGALE), tested these photometric redshifts using spectroscopic redshifts, and compared the performances of both codes. To test the reliability of parameters obtained by fitting spectral energy distributions, a mock cataloge was generated. We built a large multiwavelength catalog of more than 500 ADF-S galaxies. We successfully fitted Spectral Energy Distributions of 186 galaxies with $\rm{\chi^2_{min}<4}$, and analyzed the output parameters of the fits. We conclude that our sample consists mostly of nearby actively star-forming galaxies, and all our galaxies have a relatively high metallicity. We estimated photometric redshifts for 113 galaxies from the whole ADF-S sample. Comparing the performance of Le PHARE and CIGALE, we found that CIGALE gives more reliable redshift estimates for our galaxies, which implies that including the IR photometry allows for substantial improvement of photometric redshift estimation.

We present the results of the analysis of multiwavelength Spectral Energy Distributions (SEDs) of far-infrared galaxies detected in the AKARI Deep Field-South (ADF--S) Survey. The analysis uses a carefully selected sample of 186 sources detected at the 90 $\mu$m AKARI band, identified as galaxies with cross-identification in public catalogues. For sources without known spectroscopic redshifts, we estimate photometric redshifts after a test of two independent methods: one based on using mainly the optical -- mid infrared range, and one based on the whole range of ultraviolet -- far infrared data. We observe a vast improvement in the estimation of photometric redshifts when far infrared data are included, compared with an approach based mainly on the optical -- mid infrared range. We discuss the physical properties of our far-infrared-selected sample. We conclude that this sample consists mostly of rich in dust and young stars nearby galaxies, and, furthermore, that almost 25% of these sources are (Ultra)Luminous Infrared Galaxies. Average SEDs normalized at 90 $\mu$m for normal galaxies (138 sources), LIRGs (30 sources), and ULIRGs (18 galaxies) a the significant shift in the peak wavelength of the dust emission, and an increasing ratio between their bolometric and dust luminosities which varies from 0.39 to 0.73.

We report ALMA observations of dust continuum, 13CO J=3--2, and C18O J=3--2 line emission toward a gapped protoplanetary disk around HD 142527. The outer horseshoe-shaped disk shows the strong azimuthal asymmetry in dust continuum with the contrast of about 30 at 336 GHz between the northern peak and the southwestern minimum. In addition, the maximum brightness temperature of 24 K at its northern area is exceptionally high at 160 AU from a star. To evaluate the surface density in this region, the grain temperature needs to be constrained and was estimated from the optically thick 13CO J=3--2 emission. The lower limit of the peak surface density was then calculated to be 28 g cm-2 by assuming a canonical gas-to-dust mass ratio of 100. This finding implies that the region is locally too massive to withstand self-gravity since Toomre's Q <~1--2, and thus, it may collapse into a gaseous protoplanet. Another possibility is that the gas mass is low enough to be gravitationally stable and only dust grains are accumulated. In this case, lower gas-to-dust ratio by at least 1 order of magnitude is required, implying possible formation of a rocky planetary core.

The evolution of galaxies and the history of star formation in the Universe are among the most important topics in today's astrophysics. Especially, the role of small, irregular galaxies in the star-formation history of the Universe is not yet clear. Using the data from the AKARI IRC survey of the Large Magellanic Cloud at 3.2, 7, 11, 15, and 24 {\mu}m wavelengths, i.e., at the mid- and near-infrared, we have constructed a multiwavelength catalog containing data from a cross-correlation with a number of other databases at different wavelengths. We present the separation of different classes of stars in the LMC in color-color, and color-magnitude, diagrams, and analyze their contribution to the total LMC flux, related to point sources at different infrared wavelengths.

We propose a co-ordinated multi-observatory survey at the North Ecliptic Pole. This field is the natural extragalactic deep field location for most space observatories (e.g. containing the deepest Planck, WISE and eROSITA data), is in the continuous viewing zones for e.g. Herschel, HST, JWST, and is a natural high-visibility field for the L2 halo orbit of SPICA with deep and wide-field legacy surveys already planned. The field is also a likely deep survey location for the forthcoming Euclid mission. It is already a multi-wavelength legacy field in its own right (e.g. AKARI, LOFAR, SCUBA-2): the outstanding and unparalleled continuous mid-IR photometric coverage in this field and nowhere else enables a wide range of galaxy evolution diagnostics unachievable in any other survey field, by spanning the wavelengths of redshifted PAH and silicate features and the peak energy output of AGN hot dust. We argue from the science needs of Euclid and JWST, and from the comparative multiwavelength depths, that the logical approach is (1) a deep (H-UDF) UV/optical tile in the NEP over ~10 square arcminutes, and (2) an overlapping wide-field UV/optical HST survey tier covering >100 square arcminutes, with co-ordinated submm SPIRE mapping up to or beyond the submm point source confusion limit over a wider area and PACS data over the shallower HST tier.

We construct a dust evolution model taking into account the grain size distribution, and investigate what kind of dust processes determine the grain size distribution at each stage of galaxy evolution. In addition to the dust production by type II supernovae (SNeII) and asymptotic giant branch (AGB) stars, we consider three processes in the ISM: (i) dust destruction by SN shocks, (ii) metal accretion onto the surface of preexisting grains in the cold neutral medium (CNM) (called grain growth), and (iii) grain-grain collisions (shattering and coagulation) in the warm neutral medium (WNM) and CNM. We found that the grain size distribution in galaxies is controlled by stellar sources in the early stage of galaxy evolution, and that afterwards the main processes that govern the size distribution changes to those in the ISM. Since shattering produces a large abundance of small grains (consequently, the surface-to-volume ratio of grains increases), it enhances the efficiency of grain growth, contributing to the significant increase of the total dust mass. Grain growth creates a large bump in the grain size distribution around a ~ 0.01 \mu m. Coagulation occurs effectively after the number of small grains is enhanced by shattering, and the grain size distribution is deformed to have a bump at a ~ 0.03 - 0.05 \mu m at t ~ 10 Gyr. We conclude that the evolutions of the total dust mass and the grain size distribution in galaxies are closely related to each other, and the grain size distribution changes considerably through the galaxy evolution because the dominant dust processes which regulate the grain size distribution change.

Fast Independent Component Analysis (FastICA) is a component separation algorithm based on the levels of non-Gaussianity. Here we apply the FastICA to the component separation problem of the microwave background including carbon monoxide (CO) line emissions that are found to contaminate the PLANCK High Frequency Instrument (HFI) data. Specifically we prepare 100GHz, 143GHz, and 217GHz mock microwave sky maps including galactic thermal dust, NANTEN CO line, and the Cosmic Microwave Background (CMB) emissions, and then estimate the independent components based on the kurtosis. We find that the FastICA can successfully estimate the CO component as the first independent component in our deflection algorithm as its distribution has the largest degree of non-Gaussianity among the components. By subtracting the CO and the dust components from the original sky maps, we will be able to make an unbiased estimate of the cosmological CMB angular power spectrum.

Merger processes play an important role in galaxy formation and evolution. To study the influence of merger processes on the evolution of dust properties and cosmic star formation rate, we investigate a local sample of major merger galaxies and a control sample of isolated galaxies using GALEX ultraviolet (UV) and Spitzer infrared (IR) images. Through a statistical study, we find that dust attenuation in merger galaxies is enhanced with respect to isolated galaxies. We find this enhancement is contributed mainly by spiral galaxies in spiral-spiral (S-S) pairs, and increases with the increasing stellar mass of a galaxy. Combining the IR and UV parts of star formation rates (SFRs), we then calculated the total SFRs and specific star formation rates (SSFRs). We find the SSFRs to be enhanced in merger galaxies. This enhancement depends on galaxy stellar mass and the companion's morphology, but depends little on whether the galaxy is a primary or secondary component or on the separation between two components. These results are consistent with a previous study based only on IR images. In addition, we investigate the nuclear contributions to SFRs. SFRs in paired galaxies are more concentrated in the central part of the galaxies than in isolate galaxies. Our studies of dust attenuation show that the nuclear parts of pairs most resemble ULIRGs. Including UV data in the present work not only provides reliable information on dust attenuation, but also refines analyses of SFRs.

The Leonids show meteor storms in a period of 33 years, and known as one of the most active meteor showers. It has recently shown a meteor stream consisting of several narrow dust trails made by meteoroids ejected from a parent comet. Hence, an analysis of the temporal behavior of the meteor flux is important to study the structure of the trails. However, statistical inference for the count data is not an easy task, because of its Poisson characteristics. We carried out a wide-field video observation of the Leonid meteor storm in 2001. We formulated a state-of-the-art statistical analysis, which is called a self-organizing state space model, to infer the true behavior of the dust density of the trails properly from the meteor count data. {}From this analysis, we found that the trails have a fairly smooth spatial structure, with small and dense clumps that cause a temporal burst of meteor flux. We also proved that the time behavior (trend) of the fluxes of bright meteors and that of faint meteors are significantly different. In addition we comment on some other application of the self-organizing state-space model in fields related to astronomy and astrophysics.

We present the first measurement of the angular two-point correlation function for AKARI 90-$\mu$m point sources, detected outside of the Milky Way plane and other regions characterized by high Galactic extinction, and categorized as extragalactic sources according to our far-infrared-color based criterion (Pollo et al. 2010). This is the first measurement of the large-scale angular clustering of galaxies selected in the far-infrared after IRAS measurements. Although a full description of clustering properties of these galaxies will be obtained by more detailed studies, using either spatial correlation function, or better information about properties and at least photometric redshifts of these galaxies, the angular correlation function remains the first diagnostics to establish the clustering properties of the catalog and observed galaxy population. We find a non-zero clustering signal in both hemispheres extending up to $\sim 40$ degrees, without any significant fluctuations at larger scales. The observed correlation function is well fitted by a power law function. The notable differences between a northern and southern hemisphere are found, which can be probably attributed to the photometry problems and point out to a necessity of performing a better calibration in the data from southern hemisphere.

An accurate estimation of the star formation-related properties of galaxies is crucial for understanding the evolution of galaxies. In galaxies, ultraviolet (UV) light emitted by recently formed massive stars is attenuated by dust, which is also produced by star formation (SF) activity, and is reemitted at mid- and far- infrared (IR) wavelengths. In this study, we investigate the star formation rate (SFR) and dust extinction using UV and IR data. We selected local galaxies which are detected at AKARI FIS 90 um and matched the IRAS IIFSCz 60 um select catalog. We measured FUV and NUV flux densities from GALEX images. We examined the SF and extinction of Local galaxies using four bands of AKARI. Then, we calculated FUV and total IR luminosities, and obtained the SF luminosity, L_{SF}, the total luminosity related to star formation activity, and the SFR. We find that in most galaxies, L_{SF} is dominated by L_{dust}. We also find that galaxies with higher SF activity have a higher fraction of their SF hidden by dust. In fact, the SF of galaxies with SFRs >20 M_{sun}/yr is almost completely hidden by dust. Our results boast a significantly higher precision with respect to previously published works, due to the use of much larger object samples from the AKARI and GALEX all sky surveys.

We present a photometric catalog of infrared (IR) sources based on the North Ecliptic PoleWide field (NEP-Wide) survey of AKARI, which is an infrared space telescope launched by Japan. The NEP-Wide survey covered 5.4 deg2 area, a nearly circular shape centered on the North Ecliptic Pole, using nine photometric filter-bands from 2 - 25 {\mu}m of the Infrared Camera (IRC). Extensive efforts were made to reduce possible false objects due to cosmic ray hits, multiplexer bleeding phenomena around bright sources, and other artifacts. The number of detected sources varied depending on the filter band: with about 109,000 sources being cataloged in the near-IR bands at 2 - 5 {\mu}m, about 20,000 sources in the shorter parts of the mid-IR bands between 7 - 11 {\mu}m, and about 16,000 sources in the longer parts of the mid-IR bands, with \sim 4,000 sources at 24 {\mu}m. The estimated 5? detection limits are approximately 21 magnitude (mag) in the 2 - 5 {\mu}m bands, 19.5 - 19 mag in the 7 - 11 {\mu}m, and 18.8 - 18.5 mag in the 15 - 24 {\mu}m bands in the AB magnitude scale. The completenesses for those bands were evaluated as a function of magnitude: the 50% completeness limits are about 19.8 mag at 3 {\mu}m, 18.6 mag at 9 {\mu}m, and 18 mag at 18 {\mu}m band, respectively. To construct a reliable source catalog, all of the detected sources were examined by matching them with those in other wavelength data, including optical and ground-based near-IR bands. The final band-merged catalog contains about 114,800 sources detected in the IRC filter bands. The properties of the sources are presented in terms of the distributions in various color-color diagrams.

Far-ultraviolet (FUV) and far-infrared (FIR) luminosity functions (LFs) of galaxies show a strong evolution from $z = 0$ to $z = 1$, but the FIR LF evolves much stronger than the FUV one. The FUV is dominantly radiated from newly formed short-lived OB stars, while the FIR is emitted by dust grains heated by the FUV radiation field. It is known that dust is always associated with star formation activity. Thus, both FUV and FIR are tightly related to the star formation in galaxies, but in a very complicated manner. In order to disentangle the relation between FUV and FIR emissions, we estimate the UV-IR bivariate LF (BLF) of galaxies with {\sl GALEX} and {\sl AKARI} All-Sky Survey datasets. Recently we invented a new mathematical method to construct the BLF with given marginals and prescribed correlation coefficient. This method makes use of a tool from mathematical statistics, so called "copula". The copula enables us to construct a bivariate distribution function from given marginal distributions with prescribed correlation and/or dependence structure. With this new formulation and FUV and FIR univariate LFs, we analyze various FUV and FIR data with {\sl GALEX}, {\sl Spitzer}, and {\sl AKARI} to estimate the UV-IR BLF. The obtained BLFs naturally explain the nonlinear complicated relation between FUV and FIR emission from star-forming galaxies. Though the faint-end of the BLF was not well constrained for high-$z$ samples, the estimated linear correlation coefficient $\rho$ was found to be very high, and is remarkably stable with redshifts (from 0.95 at $z = 0$ to 0.85 at $z = 1.0$). This implies the evolution of the UV-IR BLF is mainly due to the different evolution of the univariate LFs, and may not be controlled by the dependence structure.

The relation between the ratio of infrared (IR) and ultraviolet (UV) flux densities (the infrared excess: IRX) and the slope of the UV spectrum (\beta) of galaxies plays a fundamental role in the evaluation of the dust attenuation of star forming galaxies especially at high redshifts. Many authors, however, pointed out that there is a significant dispersion and/or deviation from the originally proposed IRX-\beta relation depending on sample selection. We reexamined the IRX-\beta relation by measuring the far- and near-UV flux densities of the original sample galaxies with GALEX and AKARI imaging data, and constructed a revised formula. We found that the newly obtained IRX values were lower than the original relation because of the significant underestimation of the UV flux densities of the galaxies, caused by the small aperture of IUE, Further, since the original relation was based on IRAS data which covered a wavelength range of \lambda = 42--122\mum, using the data from AKARI which has wider wavelength coverage toward longer wavelengths, we obtained an appropriate IRX-\beta relation with total dust emission (TIR): \log(L_{\rm TIR}/L_{\rm FUV}) = \log [10^{0.4(3.06+1.58\beta)}-1] +0.22. This new relation is consistent with most of the preceding results for samples selected at optical and UV, though there is a significant scatter around it. We also found that even the quiescent class of IR galaxies follows this new relation, though luminous and ultraluminous IR galaxies distribute completely differently as well known before.

We present H-band polarimetric imagery of UX Tau A taken with HiCIAO/AO188 on the Subaru Telescope. UX Tau A has been classified as a pre-transitional disk object, with a gap structure separating its inner and outer disks. Our imagery taken with the 0.15 (21 AU) radius coronagraphic mask has revealed a strongly polarized circumstellar disk surrounding UX Tau A which extends to 120 AU, at a spatial resolution of 0.1 (14 AU). It is inclined by 46 \pm 2 degree as the west side is nearest. Although SED modeling and sub-millimeter imagery suggested the presence of a gap in the disk, with the inner edge of the outer disk estimated to be located at 25 - 30 AU, we detect no evidence of a gap at the limit of our inner working angle (23 AU) at the near-infrared wavelength. We attribute the observed strong polarization (up to 66 %) to light scattering by dust grains in the disk. However, neither polarization models of the circumstellar disk based on Rayleigh scattering nor Mie scattering approximations were consistent with the observed azimuthal profile of the polarization degrees of the disk. Instead, a geometric optics model of the disk with nonspherical grains with the radii of 30 micron meter is consistent with the observed profile. We suggest that the dust grains have experienced frequent collisional coagulations and have grown in the circumstellar disk of UX Tau A.

This paper investigate what is the main driver of the dust mass growth in the interstellar medium (ISM) by using a chemical evolution model of galaxy with metals (elements heavier than helium) in dust phase in addition to the total amount of metals. We consider asymptotic giant branch (AGB) stars, type II supernovae (SNe II) and the dust mass growth in the ISM as the sources of dust, and SN shocks as the destruction mechanism of dust. Further, to describe the dust evolution precisely, our model takes into account the age and metallicity (the ratio of metal mass to ISM mass) dependence of the sources of dust. We particularly focused on the dust mass growth, and found that the dust mass growth in the ISM is regulated by the metallicity. To quantify this aspect, we introduce a "critical metallicity", which is a metallicity at which the contribution of stars (AGB stars and SNe II) equals that of the dust mass growth in the ISM. If the star formation timescale is shorter, the value of the critical metallicity is higher, but the galactic age at which the metallicity reaches the critical metallicity is shorter. From observations, it was expected that the dust mass growth was the dominant source of dust in the Milky Way and dusty QSOs at high redshifts. By introducing the critical metallicity, it is clearly shown that the dust mass growth is the main source of dust in such galaxies with various star formation timescales and ages. The dust mass growth in the ISM is regulated by metallicity, and we stress that the critical metallicity works as an indicator to judge whether the grain growth in the ISM is dominant source of dust in a galaxy, especially because of the strong and nonlinear dependence on the metallicity.

This white paper addresses the hypothesis of light sterile neutrinos based on recent anomalies observed in neutrino experiments and the latest astrophysical data.

The standard models of inflation predict statistically homogeneous and isotropic primordial fluc- tuations, which should be tested by observations. In this paper we test the statistical isotropy of the mean of the CMB temperature fluctuations measured by the WMAP in the spherical harmonic space. A classical method to test a mean, like the simple student's-t test, is not appropriate for this purpose because the WMAP data contain anisotropic instrumental noise and suffer from the effect of the mask for the foreground emissions which breaks the statistical independence. Here we perform a band-power analysis with Monte-Carlo simulations in which we take into account the anisotropic noise and the mask. We find evidence of a non-zero mean at 99.93 % confidence level in a particular range of multipoles. The evidence against the zero-mean assumption as a whole is still significant at the 99 % confidence level even if the fact is taken into account that we have tested multiple ranges.

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 high-resolution, H-band, imaging observations, collected with Subaru/HiCIAO, of the scattered light from the transitional disk around SAO 206462 (HD 135344B). Although previous sub-mm imagery suggested the existence of the dust-depleted cavity at r~46AU, our observations reveal the presence of scattered light components as close as 0.2" (~28AU) from the star. Moreover, we have discovered two small-scale spiral structures lying within 0.5" (~70AU). We present models for the spiral structures using the spiral density wave theory, and derive a disk aspect ratio of h~0.1, which is consistent with previous sub-mm observations. This model can potentially give estimates of the temperature and rotation profiles of the disk based on dynamical processes, independently from sub-mm observations. It also predicts the evolution of the spiral structures, which can be observable on timescales of 10-20 years, providing conclusive tests of the model. While we cannot uniquely identify the origin of these spirals, planets embedded in the disk may be capable of exciting the observed morphology. Assuming that this is the case, we can make predictions on the locations and, possibly, the masses of the unseen planets. Such planets may be detected by future multi-wavelengths observations.

We estimate minimum dust abundances required for secular gravitational instability (SGI) to operate at the midplane dust layer of protoplanetary disks. For SGI to be a viable process, the growth time of the instability T_grow must be shorter than the radial drift time of the dust T_drift. The growth time depends on the turbulent diffusion parameter alpha, because the modes with short wavelengths are stabilized by turbulent diffusion. Assuming that turbulence is excited via the Kelvin-Helmholtz or streaming instabilities in the dust layer, and that its strength is controlled by the energy supply rate from dust accretion, we estimate the diffusion parameter and the growth time of the instability. The condition T_grow < T_drift requires that the dust abundance must be greater than a critical abundance Z_min, which is a function of the Toomre parameter Q_g and aspect ratio h_g / r of the gas disk. For a wide range of parameter space, the required dust abundance is less than 0.1. A slight increase in dust abundance opens a possible route for the dust to directly collapse to planetesimals.

We have studied the properties of more than 1600 low-redshift galaxies by utilizing high-quality infrared flux measurements of the AKARI All-Sky Survey and physical quantities based on optical and 21-cm observations. Our goal is to understand the physics determining the infrared spectral energy distribution (SED). The ratio of the total infrared luminosity L_TIR, to the star-formation rate (SFR) is tightly correlated by a power-law to specific SFR (SSFR), and L_TIR is a good SFR indicator only for galaxies with the largest SSFR. We discovered a tight linear correlation for normal galaxies between the radiation field strength of dust heating, estimated by infrared SED fits (U_h), and that of galactic-scale infrared emission (U_TIR ~ L_TIR/R^2), where R is the optical size of a galaxy. The dispersion of U_h along this relation is 0.3 dex, corresponding to 13% dispersion in the dust temperature. This scaling and the U_h/U_TIR ratio can be explained physically by a thin layer of heating sources embedded in a thicker, optically-thick dust screen. The data also indicate that the heated fraction of the total dust mass is anti-correlated to the dust column density, supporting this interpretation. In the large U_TIR limit, the data of circumnuclear starbursts indicate the existence of an upper limit on U_h, corresponding to the maximum SFR per gas mass of ~ 10 Gyr^{-1}. We find that the number of galaxies sharply drops when they become optically thin against dust-heating radiation, suggesting that a feedback process to galaxy formation (likely by the photoelectric heating) is working when dust-heating radiation is not self-shielded on a galactic scale. Implications are discussed for the M_HI-size relation, the Kennicutt-Schmidt relation, and galaxy formation in the cosmological context.

AKARI/Infrared Camera (IRC) Point Source Catalog provides a large amount of flux data at {\it S9W} ($9\ {\rm \mu m}$) and {\it L18W} ($18\ {\rm \mu m}$) bands. With the goal of constructing Star-Formation Rate(SFR) calculations using IRC data, we analyzed an IR selected GALEX-SDSS-2MASS-AKARI(IRC/Far-Infrared Surveyor) sample of 153 nearby galaxies. The far-infrared fluxes were obtained from AKARI diffuse maps to correct the underestimation for extended sources raised by the point-spread function photometry. SFRs of these galaxies were derived by the spectral energy distribution fitting program CIGALE. In spite of complicated features contained in these bands, both the {\it S9W} and {\it L18W} emission correlate with the SFR of galaxies. The SFR calibrations using {\it S9W} and {\it L18W} are presented for the first time. These calibrations agree well with previous works based on Spitzer data within the scatters, and should be applicable to dust-rich galaxies.

We study the turbulence induced in the dust layer of a protoplanetary disk based on the energetics of dust accretion due to gas drag. We estimate turbulence strength from the energy supplied by dust accretion, using the radial drift velocity of the dust particles in a laminar disk. Our estimate of the turbulence strength agrees with previous analytical and numerical research on the turbulence induced by Kelvin-Helmholtz and/or streaming instabilities for particles whose stopping time is less than the Keplerian time. For such small particles, the strongest turbulence is expected to occur when the dust-to-gas ratio of the disk is ~C_eff^(1/2) (h_g / r) ~ 10^(-2), where C_eff ~ 0.2 represents the energy supply efficiency to turbulence and h_g / r ~ 5 x 10^(-2) is the aspect ratio of the gas disk. The maximum viscosity parameter is alpha_max ~ C_eff T_s (h_g / r)^2 ~ 10^(-4) T_s, where T_s (<1) is the non-dimensional stopping time of the dust particles. Modification in the dust-to-gas ratio from the standard value, 10^(-2), by any process, results in weaker turbulence and a thinner dust layer, and consequently may accelerate the growth process of the dust particles.

We present a new analytic approach to the disk-planet interaction that is especially useful for planets with eccentricity larger than the disk aspect ratio. We make use of the dynamical friction formula to calculate the force exerted on the planet by the disk, and the force is averaged over the period of the planet. The resulting migration and eccentricity damping timescale agrees very well with the previous works in which the planet eccentricity is moderately larger than the disk aspect ratio. The advantage of this approach is that it is possible to apply this formulation to arbitrary large eccentricity. We have found that the timescale of the orbital evolution depends largely on the adopted disk model in the case of highly eccentric planets. We discuss the possible implication of our results to the theory of planet formation.

The nearby universe remains the best laboratory to understand physical properties of galaxies and is a reference for any comparison with high redshift observations. The all sky (or very large) surveys performed from the ultraviolet (UV) to the far-infrared (far-IR) gives us large datasets with a very large wavelength coverage to perform a reference study. We want to investigate dust attenuation characteristics as well as star formation rate (SFR) calibrations on a sample of nearby galaxies observed over 13 bands from 0.15 to 160 microns. A sample of 363 galaxies is built from the AKARI /FIS all sky survey cross-correlated with SDSS and GALEX surveys. Broad band spectral energy distributions are fitted with the CIGALE code optimized to analyze variations in dust attenuation curves and SFR measurements and based on an energetic budget between the stellar and dust emission. Our galaxy sample is primarily selected in far-IR and mostly constituted of massive, actively star forming galaxies. There is some evidence for a dust attenuation law slightly steeper than that used for starburst galaxies but we are unable to constrain the presence or not of a bump at 220 nm. We confirm that a time dependent dust attenuation is necessary to perform the best fits. Various calibrations of the dust attenuation in the UV as a function of UV-optical colours are discussed. A calibration of the current star formation rate combining UV and total IR emissions is proposed with an accurate estimation of dust heating by old stars: for the whole sample 17 % of the total dust luminosity is not related to the recent star formation

Coagulation of submicron-sized dust grains into porous aggregates is the initial step of dust evolution in protoplanetary disks. Recently, it has been pointed out that negative charging of dust in the weakly ionized disks could significantly slow down the coagulation process. In this paper, we apply the growth criteria obtained in Paper I to finding out a location ("frozen" zone) where the charging stalls dust growth at the fractal growth stage. For low-turbulence disks, we find that the frozen zone can cover the major part of the disks at a few to 100 AU from the central star. The maximum mass of the aggregates is approximately 10^{-7} g at 1 AU and as small as a few monomer masses at 100 AU. Strong turbulence can significantly reduce the size of the frozen zone, but such turbulence will cause the fragmentation of macroscopic aggregates at later stages. We examine a possibility that complete freezeout of dust evolution in low-turbulence disks could be prevented by global transport of dust in the disks. Our simple estimation shows that global dust transport can lead to the supply of macroscopic aggregates and the removal of frozen aggregates on a timescale of 10^6 yr. This overturns the usual understanding that tiny dust particles get depleted on much shorter timescales unless collisional fragmentation is effective. The frozen zone together with global dust transport might explain "slow" (\sim 10^6 yr) dust evolution suggested by infrared observation of T Tauri stars and by radioactive dating of chondrites.

Collisional growth of submicron-sized dust grains into macroscopic aggregates is the first step of planet formation in protoplanetary disks. These grains are expected to carry nonzero negative charges in the weakly ionized disks, but its effect on their collisional growth has not been fully understood so far. In this paper, we investigate how the charging affects the evolution of the dust size distribution properly taking into account the charging mechanism in a weakly ionized gas as well as porosity evolution through low-energy collisions. To clarify the role of the size distribution, we divide our analysis into two steps. First, we analyze the collisional growth of charged aggregates assuming a monodisperse (i.e., narrow) size distribution. We show that the monodisperse growth stalls due to the electrostatic repulsion when a certain condition is met, as is already expected in the previous work. Second, we numerically simulate dust coagulation using Smoluchowski's method to see how the outcome changes when the size distribution is allowed to freely evolve. We find that, under certain conditions, the dust undergoes bimodal growth where only a limited number of aggregates continue to grow carrying the major part of the dust mass in the system. This occurs because remaining small aggregates efficiently sweep up free electrons to prevent the larger aggregates from being strongly charged. We obtain a set of simple criteria that allows us to predict how the size distribution evolves for a given condition. In Paper II (arXiv:1009.3101), we apply these criteria to dust growth in protoplanetary disks.

By cross-correlating AKARI infrared (IR) sources with the SDSS galaxies, we identified 2357 infrared galaxies with a spectroscopic redshift. This is not just one of the largest samples of local IR galaxies, but AKARI provides crucial FIR bands (9, 18, 65, 90, 140, and 160um) in accurately measuring galaxy SED across the peak of the dust emission at ~100um. By fitting modern IR SED models to the AKARI photometry, we measured the total infrared luminosity (L_IR) of individual galaxies more accurately. Using this L_IR, we constructed luminosity functions of infrared galaxies at a median redshift of z=0.031, with 4 times larger sample than previous work. The LF agrees well with that at z=0.0082 (RBGS), showing smooth and continuous evolution toward higher redshift LFs measured in the AKARI NEP deep field. The derived local cosmic IR luminosity density is Omega_IR=3.8x10^8 LsunMpc^-3. We separate galaxies into AGN, star-forming, and composite by using the [NII]/Ha vs [OIII]/Hb line ratios. The fraction of AGN shows a continuous increase with increasing L_IR from 25% to 90% at 9<log L_IR<12.5. The SFR_Ha and L_[OIII] show good correlations with L_IR for SFG (star-forming galaxies) and AGN, respectively. The self-absorption corrected Ha/Hb ratio shows a weak increase with L_IR with a substantial scatter. When we separate IR LFs into contributions from AGN and SFG, the AGN contribution becomes dominant at L_IR>10^11Lsun, coinciding the break of the both SFG and AGN IR LFs. At L_IR<10^11Lsun, SFG dominates IR Lfs. Only 1.1% of Omega_IR is produced by LIRG, and only 0.03% is by ULIRG in the local Universe. This work also provides the most accurate infrared luminosity density of the local Universe to date. Compared with high redshift results from the AKARI NEP deep survey, we observed a strong evolution of Omega_IR^SFG ~(1+z)^4.1+-0.4 and Omega_IR^AGN ~(1+z)^4.1+-0.5 (abridged).

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.

To estimate cosmological parameters from a given dataset, we need to construct a likelihood function, which sometimes has a complicated functional form. We introduce the copula, a mathematical tool to construct an arbitrary multivariate distribution function from one-dimensional marginal distribution functions with any given dependence structure. It is shown that a likelihood function constructed by the so-called Gaussian copula can reproduce very well the n-dimensional probability distribution of the cosmic shear power spectrum obtained from a large number of ray-tracing simulations. This suggests that the Copula likelihood will be a powerful tool for future weak lensing analyses, instead of the conventional multivariate Gaussian likelihood.

The estimation of cosmological parameters from a given data set requires a construction of a likelihood function which, in general, has a complicated functional form. We adopt a Gaussian copula and constructed a copula likelihood function for the convergence power spectrum from a weak lensing survey. We show that the parameter estimation based on the Gaussian likelihood erroneously introduces a systematic shift in the confidence region, in particular for a parameter of the dark energy equation of state w. Thus, the copula likelihood should be used in future cosmological observations.

We present results of a 1.1 mm deep survey of the AKARI Deep Field South (ADF-S) with AzTEC mounted on the Atacama Submillimetre Telescope Experiment (ASTE). We obtained a map of 0.25 sq. deg area with an rms noise level of 0.32-0.71 mJy. This is one of the deepest and widest maps thus far at millimetre and submillimetre wavelengths. We uncovered 198 sources with a significance of 3.5-15.6 sigma, providing the largest catalog of 1.1 mm sources in a contiguous region. Most of the sources are not detected in the far-infrared bands of the AKARI satellite, suggesting that they are mostly at z ~ 1.5 given the detection limits. We constructed differential and cumulative number counts in the ADF-S, the Subaru/XMM Newton Deep Field (SXDF), and the SSA 22 field surveyed by AzTEC/ASTE, which provide currently the tightest constraints on the faint end. The integration of the best-fit number counts in the ADF-S find that the contribution of 1.1 mm sources with fluxes >=1 mJy to the cosmic infrared background (CIB) at 1.1 mm is 12-16%, suggesting that the large fraction of the CIB originates from faint sources of which the number counts are not yet constrained. We estimate the cosmic star-formation rate density contributed by 1.1 mm sources with >=1 mJy using the best-fit number counts in the ADF-S and find that it is lower by about a factor of 5-10 compared to those derived from UV/optically-selected galaxies at z ~ 2-3. The fraction of stellar mass of the present-day universe produced by 1.1 mm sources with >=1 mJy at z >= 1 is ~20%, calculated by the time integration of the star-formation rate density. If we consider the recycled fraction of >0.4, which is the fraction of materials forming stars returned to the interstellar medium, the fraction of stellar mass produced by 1.1 mm sources decrease to <~10%.

Infrared (IR) luminosity is fundamental to understanding the cosmic star formation history and AGN evolution. The AKARI IR space telescope performed all sky survey in 6 IR bands (9, 18, 65, 90, 140, and 160um) 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, AKARI can much more precisely measure the total infrared luminosity (L_TIR) of individual galaxies, and thus, the total infrared luminosity density in the local Universe. By fitting IR SED models, we have re-measured L_TIR of the IRAS Revised Bright Galaxy Sample. We present mid-IR monochromatic luminosity to L_TIR conversions for Spitzer 8,24um, AKARI 9,18um, IRAS 12um, WISE 12,22um, and ISO 15um filters, with scatter ranging 13-44%. The resulting AKARI IR luminosity function (LF) agrees well with that from the IRAS. We integrate the LF weighted by L_TIR to obtain a cosmic IR luminosity density of Omega_TIR= (8.5^{+1.5}_{-2.3})x 10^7 L Mpc^-3, of which 7+-1% is produced by LIRGs, and only 0.4+-0.1% is from ULIRGs in the local Universe. Once IR contributions from AGN and star-forming galaxies (SFG) are separated, SFG IR LF shows a steep decline at the bright-end. Compared with high-redshift results from the AKARI NEP deep survey, these data show a strong evolution of Omega_TIRSF propto (1+z)^4.0+-0.5, and Omega_TIRAGN propto (1+z)^4.4+-0.4. For Omega_TIRAGN, the ULIRG contribution exceeds that from LIRG already by z~1. A rapid evolution in both Omega_TIRAGN and Omega_TIRSFG suggests the correlation between star formation and black hole accretion rate continues up to higher redshifts. We compare the evolution of Omega_TIRAGN to that of X-ray luminosity density. The Omega_TIRAGN/Omega_X-rayAGN ratio shows a possible increase at z>1, suggesting an increase of obscured AGN at z>1.

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

We show an analytic method to construct a bivariate distribution function (DF) with given marginal distributions and correlation coefficient. We introduce a convenient mathematical tool, called a copula, to connect two DFs with any prescribed dependence structure. If the correlation of two variables is weak (Pearson's correlation coefficient $|\rho| <1/3 $), the Farlie-Gumbel-Morgenstern (FGM) copula provides an intuitive and natural way for constructing such a bivariate DF. When the linear correlation is stronger, the FGM copula cannot work anymore. In this case, we propose to use a Gaussian copula, which connects two given marginals and directly related to the linear correlation coefficient between two variables. Using the copulas, we constructed the BLFs and discuss its statistical properties. Especially, we focused on the FUV--FIR BLF, since these two luminosities are related to the star formation (SF) activity. Though both the FUV and FIR are related to the SF activity, the univariate LFs have a very different functional form: former is well described by the Schechter function whilst the latter has a much more extended power-law like luminous end. We constructed the FUV-FIR BLFs by the FGM and Gaussian copulas with different strength of correlation, and examined their statistical properties. Then, we discuss some further possible applications of the BLF: the problem of a multiband flux-limited sample selection, the construction of the SF rate (SFR) function, and the construction of the stellar mass of galaxies ($M_*$)--specific SFR ($SFR/M_*$) relation. The copulas turned out to be a very useful tool to investigate all these issues, especially for including the complicated selection effects.

To separate stars and galaxies in the far infrared AKARI All-Sky Survey data, we have selected a sample with the complete color information available in the low extinction regions of the sky and constructed color-color plots for these data. We looked for the method to separate stars and galaxies using the color information. We performed an extensive search for the counterparts of these selected All-Sky Survey sources in the NED and SIMBAD databases. Among 5176 objects, we found 4272 galaxies, 382 other extragalactic objects, 349 Milky Way stars, 50 other Galactic objects, and 101 sources detected before in various wavelengths but of an unknown origin. 22 sources were left unidentified. Then, we checked colors of stars and galaxies in the far-infrared flux-color and color-color plots. In the resulting diagrams, stars form two clearly separated clouds. One of them is easy to be distinguished from galaxies and allows for a simple method of excluding a large part of stars using the far-infrared data. The other smaller branch, overplotting galaxies, consists of stars known to have an infrared excess, like Vega and some fainter stars discovered by IRAS or 2MASS. The color properties of these objects in any case make them very difficult to distinguish from galaxies. We conclude that the FIR color-color diagrams allow for a high-quality star-galaxy separation. With the proposed simple method we can select more that 95 % of galaxies rejecting at least 80 % of stars.

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

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.

The AKARI All-Sky Survey provided the first bright point source catalog detected at 90um. Starting from this catalog, we selected galaxies by matching AKARI sources with those in the IRAS PSCz. Next, we have measured total GALEX FUV and NUV flux densities. Then, we have matched this sample with SDSS and 2MASS galaxies. By this procedure, we obtained the final sample which consists of 607 galaxies. If we sort the sample with respect to 90um, their average SED shows a coherent trend: the more luminous at 90um, the redder the global SED becomes. The M_r--NUV-r color-magnitude relation of our sample does not show bimodality, and the distribution is centered on the green valley between the blue cloud and red sequence seen in optical surveys. We have established formulae to convert FIR luminosity from AKARI bands to the total infrared (IR) luminosity L_TIR. With these formulae, we calculated the star formation directly visible with FUV and hidden by dust. The luminosity related to star formation activity (L_SF) is dominated by L_TIR even if we take into account the far-infrared (FIR) emission from dust heated by old stars. At high star formation rate (SFR) (> 20 Msun yr^-1), the fraction of directly visible SFR, SFR_FUV, decreases. We also estimated the FUV attenuation A_FUV from FUV-to-total IR (TIR) luminosity ratio. We also examined the L_TIR/L_FUV-UV slope (FUV- NUV) relation. The majority of the sample has L_TIR/L_FUV ratios 5 to 10 times lower than expected from the local starburst relation, while some LIRGs and all the ULIRGs of this sample have higher L_TIR/L_FUV ratios. We found that the attenuation indicator L_TIR/L_FUV is correlated to the stellar mass of galaxies, M*, but there is no correlation with specific SFR (SSFR), SFR/M*, and dust attenuation L_TIR/L_FUV. (abridged)

To investigate the nature and properties of far-infrared (FIR) sources from the AKARI Deep Field South (ADF-S), we performed an extensive search for the counterparts of 1000 ADF-S objects brighter than 0.0301 Jy in the WIDE-S (90 $\mu$m) AKARI band in the public databases (NED and SIMBAD). We analyzed the properties of the resulting sample: statistic of the identified objects, number counts, redshift distribution and morphological types. We also made a crude analysis of the clustering properties of the sources and constructed spectral energy distributions (SEDs) of 47 selected objects with the best photometry. Among 1000 investigated ADF-S sources, 545 were identified at other wavelengths. From them, 518 are known galaxies, and 343 of them were not known previously as infra-red sources. We found redshifts of 48 extragalactic objects and morphological types of 77 galaxies. We conclude that the bright FIR point sources observed in the ADF-S are mostly nearby galaxies.Their properties are very similar to properties of the local population of optically bright galaxies, except an unusually high ratio of peculiar or interacting objects and a smaller percentage of elliptical galaxies. The percentage of lenticular galaxies is the same as in the optically bright population which suggests they may frequently contain a significant amount of cool dust. It is possible that source confusion plays a significant role in more than 34% of measurements. The SEDs display a variety of galaxy types, from very actively star forming to very quiescent. Thanks to the AKARI long wavelength bands it was revealed for the first time that these galaxies form a population of objects with very cool dust.

The Herschel ATLAS is the largest open-time key project that will be carried out on the Herschel Space Observatory. It will survey 510 square degrees of the extragalactic sky, four times larger than all the other Herschel surveys combined, in five far-infrared and submillimetre bands. We describe the survey, the complementary multi-wavelength datasets that will be combined with the Herschel data, and the six major science programmes we are undertaking. Using new models based on a previous submillimetre survey of galaxies, we present predictions of the properties of the ATLAS sources in other wavebands.

We select galaxies in UV rest-frame at z=0, z~0.7 and z~1 together with a sample of LBGs at z~1, the samples are built in order to sample the same range of luminosity at any redshift. The evolution of the IR and UV luminosities with z is analysed for individual galaxies as well as in terms of luminosity functions. The L_IR/L_UV ratio is used to measure dust attenuation. This ratio does not exhibit a strong evolution with z for the bulk of our sample galaxies but some trends are found for galaxies with a strong dust attenuation and for UV luminous sources: galaxies with L_IR/L_UV>10 are more frequent at z>0 than at z=0 and the largest values of L_IR/L_UV are found for UV faint objects; conversely the most luminous galaxies of our samples (L_UV> 2 10^{10} L_sun$), detected at z=1, exhibit a lower dust attenuation than the fainter ones. L_IR/L_UV increases with the K rest-frame luminosity of the galaxies at all the redshifts considered and shows a residual anti-correlation with L_UV. The most massive and UV luminous galaxies exhibit quite large specific star formation rates. LBGs exhibit systematically lower dust attenuation than UV selected galaxies of same luminosity but similar specific star formation rates. The analysis of the UV+IR luminosity functions leads to the conclusion that up to z = 1 most of the star formation activity of UV selected galaxies is emitted in IR. Whereas we are able to retrieve all the star formation from our UV selection at z=0.7, at z = 1 we miss a large fraction of galaxies more luminous than ~ 10^{11} L_sun. The effect is found larger for Lyman Break Galaxies.

We investigated the relationship between supermassive black hole (SMBH) mass and host starburst luminosity in Seyfert galaxies and Palomar-Green QSOs, focusing on the host galaxy morphology. Host starburst luminosity was derived from the 11.3 micron polycyclic aromatic hydrocarbon luminosity. We found that the SMBH masses of elliptical-dominated host galaxies are more massive than those of disk-dominated host galaxies statistically. We also found that the SMBH masses of disk-dominated host galaxies seem to be suppressed even under increasing starburst luminosity. These findings imply that final SMBH mass is strongly regulated by host galaxy morphology. This can be understood by considering the radiation drag model as the SMBH growth mechanism, taking into account the radiation efficiency of the host galaxy.

In order to find counterparts of the detected objects in the AKARI Deep Field South (ADFS) in all available wavelengths, we searched public databases (NED, SIMBAD and others). Checking 500 sources brighter than 0.0482 Jy in the AKARI Wide-S band, we found 114 sources with possible counterparts, among which 78 were known galaxies. We present these sources as well as our first attempt to construct spectral energy distributions (SEDs) for the most secure and most interesting sources among them, taking into account all the known data together with the AKARI measurements in four bands.

We present recent advances in theoretical studies of the formation and evolution of dust in primordial supernovae (SNe) that are considered to be the main sources of dust in the early universe. Being combined with the results of calculations of dust formation in the ejecta of Population III SNe, the investigations of the evolution of newly formed dust within supernova remnants (SNRs) show that smaller grains are predominantly destroyed by sputtering in the shocked gas, while larger grains are injected into the ambient medium. The mass of dust grains surviving the destruction in SNRs reaches up to 0.1--15 $M_\odot$, which is high enough to account for the content of dust observed for the host galaxies of quasars at $z > 5$. In addition, the transport of dust formed in the ejecta causes the formation of low-mass stars in the dense shell of primordial SNRs and affects the elemental composition of those stars. We also show that the flat extinction curve is expected in the high-redshift universe where SNe are the possible sources of dust.

Deep observations of the CDFS have been secured at 15um with AKARI/IRC infrared space telescope (ESA open time). From these observations, we define a sample of MIR-selected galaxies at 15um and we also obtain 15um flux densities for a sample of LBGs at z=1 already observed at 24um with Spitzer/MIPS. Number counts for the MIR-selected sample show a bump around a 15um flux density of 0.2mJy that can be attributed to galaxies at z>0.4 and at z>0.8 for the fainter part of the bump. This bump seems to be shifted as compared to other works and a possible origin can be the Cosmic variance. Thanks to this dataset, we have tested, on the two above samples at z=1, the validity of the conversions from monochromatic luminosities nu.f(nu) at a rest-frame wavelength of 8um by a comparison with total dust luminosities estimated from Spitzer rest-frame 12um data that we use as a reference. We find that the 8um dust luminosities are not all consistent and that some of them are better when compared to L(dust) evaluated from longer wavelength luminosities. We also find that the rest-frame 8um luminosities provide globally good estimates of L(dust). By comparing our data for the two samples to several libraries of SEDs, we find that models can explain the diversity of the observed f(24)/f(15) ratio quite reasonably for the MIR-selected sample and better for the LBG sample which are less dispersed than the MIR selection. However, when we analyse the luminosity dependence of this ratio, we find important discrepancies. Finally, we revisit the evolution of L(dust)/L(UV) ratio with the redshift z by re-calibrating previous L(dust) at z=2 based on our results and added new data points at higher redshifts. The decreasing trend is amplified as compared to the previous estimate.

We investigate whether the mean star formation activity of star-forming galaxies from z=0 to z=0.7 in the GOODS-S field can be reproduced by simple evolution models of these systems. In this case, such models might be used as first order references for studies at higher z to decipher when and to what extent a secular evolution is sufficient to explain the star formation history in galaxies. We selected star-forming galaxies at z=0 and at z=0.7 in IR and in UV to have access to all the recent star formation. We focused on galaxies with a stellar mass ranging between 10^{10} and 10^{11} M_sun for which the results are not biased by the selections. We compared the data to chemical evolution models developed for spiral galaxies and originally built to reproduce the main characteristics of the Milky Way and nearby spirals without fine-tuning them for the present analysis. We find a shallow decrease in the specific star formation rate (SSFR) when the stellar mass increases. The evolution of the SSFR characterizing both UV and IR selected galaxies from z=0 to z=0.7 is consistent with the models built to reproduce the present spiral galaxies. There is no need to strongly modify of the physical conditions in galaxies to explain the average evolution of their star formation from z=0 to z=0.7. We use the models to predict the evolution of the star formation rate and the metallicity on a wider range of redshift and we compare these predictions with the results of semi-analytical models.

We investigate the extinction curves of young galaxies in which dust is supplied from Type II supernovae (SNe II) and/or pair instability supernovae (PISNe). Since at high redshift (z>5), low-mass stars cannot be dominant sources for dust grains, SNe II and PISNe, whose progenitors are massive stars with short lifetimes, should govern the dust production. Here, we theoretically investigate the extinction curves of dust produced by SNe II and PISNe, taking into account reverse shock destruction induced by collision with ambient interstellar medium. We find that the extinction curve is sensitive to the ambient gas density around a SN, since the efficiency of reverse shock destruction strongly depends on it. The destruction is particularly efficient for small-sized grains, leading to a flat extinction curve in the optical and ultraviolet wavelengths. Such a large ambient density as n_H > 1 cm^{-3} produces too flat an extinction curve to be consistent with the observed extinction curve for SDSS J104845.05+463718.3 at z=6.2. Although the extinction curve is highly sensitive to the ambient density, the hypothesis that the dust is predominantly formed by SNe at z~6 is still allowed by the current observational constraints. For further quantification, the ambient density should be obtained by some other methods. Finally we also discuss the importance of our results for observations of high-z galaxies, stressing a possibility of flat extinction curves.

We study dust accumulation by photophoresis in optically thin gas disks. Using formulae of the photophoretic force that are applicable for the free molecular regime and for the slip-flow regime, we calculate dust accumulation distances as a function of the particle size. It is found that photophoresis pushes particles (smaller than 10 cm) outward. For a Sun-like star, these particles are transported to 0.1-100 AU, depending on the particle size, and forms an inner disk. Radiation pressure pushes out small particles (< 1 mm) further and forms an extended outer disk. Consequently, an inner hole opens inside ~0.1 AU. The radius of the inner hole is determined by the condition that the mean free path of the gas molecules equals the maximum size of the particles that photophoresis effectively works on (100 micron - 10 cm, depending on the dust property). The dust disk structure formed by photophoresis can be distinguished from the structure of gas-free dust disk models, because the particle sizes of the outer disks are larger, and the inner hole radius depends on the gas density.

Ultraviolet (UV) galaxies have been selected from GALEX. The presence of a FUV-dropout in their spectral energy distributions proved to be a very complete (83.3%) but not very efficient (21.4%) tool for identifying Lyman Break Galaxies (LBGs) at z~1. We divide the LBG sample into two sub-classes: red LBGs (RLBGs) detected at 24 micron which are mainly Luminous IR Galaxies (LIRGs) and blue LBGs (BLBGs) undetected at 24 microns down to 83 microJy. Two of the RLBGs are also detected at 70 micron. The median SED of the RLBGs is similar (above lambda~1 micron) to the dusty starburst HR10. However, unlike local (U)LIRGs, RLBGs are UV bright objects. We suggest that these objects contain a large amount of dust but that some bare stellar populations are also directly visible. The median SED of the BLBGs is consistent with their containing the same stellar population as the RLBGs but with a lower dust content. The luminosity function of our LBG sample at z~1 is similar to the luminosity function of NUV-selected galaxies at the same redshift. The integrated luminosity densities of z~1 LBGs and NUV-selected galaxies are very consistent. We show that star formation rates (SFRs) estimated from UV measurements and corrected using the IRX-beta method provide average total SFR_TOT in agreement with SFR_UV + SFR_dust. However, IRX-beta-based SFR_TOT shows a large dispersion. Summing up the detected UV (1150A rest-frame) and IR-based star formation rates of the detected objects, we find that only one third of the total (i.e. UV + dust) LBG SFR resides in BLBGs and two thirds in RLBGs, even though most LBGs at z~1 are BLBGs. On the other hand, the total SFR of LBGs accounts for only 11% of the total SFR at z~1. Finally, we observe a regular decrease of L_TIR / L_FUV from z=0 to z~2 for UV-selected samples.

AKARI, the first Japanese satellite dedicated to infrared astronomy, was launched on 2006 February 21, and started observations in May of the same year. AKARI has a 68.5 cm cooled telescope, together with two focal-plane instruments, which survey the sky in six wavelength bands from the mid- to far-infrared. The instruments also have the capability for imaging and spectroscopy in the wavelength range 2 - 180 micron in the pointed observation mode, occasionally inserted into the continuous survey operation. The in-orbit cryogen lifetime is expected to be one and a half years. The All-Sky Survey will cover more than 90 percent of the whole sky with higher spatial resolution and wider wavelength coverage than that of the previous IRAS all-sky survey. Point source catalogues of the All-Sky Survey will be released to the astronomical community. The pointed observations will be used for deep surveys of selected sky areas and systematic observations of important astronomical targets. These will become an additional future heritage of this mission.

We study the ultraviolet to far-infrared (hereafter UV-to-IR) SEDs of a sample of intermediate redshift (0.2 < z < 0.7) UV-selected galaxies from the ELAIS-N1 and ELAIS-N2 fields by fitting a multi-wavelength dataset to a library of GRASIL templates. Star formation related properties of the galaxies are derived from the library of models by using the Bayesian statistics. We find a decreasing presence of galaxies with low attenuation and low total luminosity as redshift decreases, which does not hold for high total luminosity galaxies. In addition the dust attenuation of low mass galaxies increases as redshift decreases, and this trend seems to disappear for galaxies with M* > 10^11 M_sun. This result is consistent with a mass dependent evolution of the dust to gas ratio, which could be driven by a mass dependent efficiency of star formation in star forming galaxies. The specific star formation rates (SSFR) decrease with increasing stellar mass at all redshifts, and for a given stellar mass the SSFR decreases with decreasing redshift. The differences in the slope of the M*--SSFR relation found between this work and others at similar redshift could be explained by the adopted selection criteria of the samples which, for a UV selected sample, favours blue, star forming galaxies.

The TIR and FUV luminosity functions of galaxies and the related luminosity densities rho(TIR) and rho(FUV) are known to evolve from z=0 to z~1 but with a different rate: the galaxy populations appear brighter in the past at both wavelengths but the evolution in TIR is larger than in FUV. It leads to an increase of the ratio of TIR to FUV luminosity densities rho(TIR)/rho(FUV) which can be interpreted as a global increase of the dust attenuation from z=0 to z~1. Our aim is to understand the origin of this increase. We focus on infrared galaxies with L(TIR>10^{11} Lsun at z~0.7 observed by SPITZER/MIPS and we measure their UV emission at 2310 A from GALEX. These Luminous InfraRed Galaxies (LIRGs) represent the bulk of the TIR luminosity density at intermediate redshift. Some evolution of L(TIR)/L(FUV) and therefore of dust attenuation is found for these galaxies: LIRGs at z=0.7 span a larger range of L(TIR)/L(FUV) ratios than at z=0 and their mean dust attenuation at FUV is found ~0.5 mag lower than for their local counterparts. The decrease of dust attenuation is found lower than that reported in other studies for bright galaxies selected in UV rest-frame at z=1 and 2. From a semi-quantitative analysis we find that the slight decrease of dust attenuation for LIRGs at z=0.7 remains consistent with the increase of rho(TIR)/rho(FUV) with redshift.

We study dust attenuation and stellar mass of $\rm z\sim 0.6$ star-forming galaxies using new SWIRE observations in IR and GALEX observations in UV. Two samples are selected from the SWIRE and GALEX source catalogs in the SWIRE/GALEX field ELAIS-N1-00 ($\Omega = 0.8$ deg$^2$). The UV selected sample has 600 galaxies with photometric redshift (hereafter photo-z) $0.5 \leq z \leq 0.7$ and NUV$\leq 23.5$ (corresponding to $\rm L_{FUV} \geq 10^{9.6} L_\sun$). The IR selected sample contains 430 galaxies with $f_{24\mu m} \geq 0.2$ mJy ($\rm L_{dust} \geq 10^{10.8} L_\sun$) in the same photo-z range. It is found that the mean $\rm L_{dust}/L_{FUV}$ ratios of the z=0.6 UV galaxies are consistent with that of their z=0 counterparts of the same $\rm L_{FUV}$. For IR galaxies, the mean $\rm L_{dust}/L_{FUV}$ ratios of the z=0.6 LIRGs ($\rm L_{dust} \sim 10^{11} L_\sun$) are about a factor of 2 lower than local LIRGs, whereas z=0.6 ULIRGs ($\rm L_{dust} \sim 10^{12} L_\sun$) have the same mean $\rm L_{dust}/L_{FUV}$ ratios as their local counterparts. This is consistent with the hypothesis that the dominant component of LIRG population has changed from large, gas rich spirals at z$>0.5$ to major-mergers at z=0. The stellar mass of z=0.6 UV galaxies of $\rm L_{FUV} \leq 10^{10.2} L_\sun$ is about a factor 2 less than their local counterparts of the same luminosity, indicating growth of these galaxies. The mass of z=0.6 UV lunmous galaxies (UVLGs: $\rm L_{FUV} > 10^{10.2} L_\sun$) and IR selected galaxies, which are nearly exclusively LIRGs and ULIRGs, is the same as their local counterparts.

We use deep GALEX images of CDFS in UV to define the first large sample of 420 Lyman Break Galaxies at z~1. We use a PSF fitting to estimate UV magnitudes on these deep crowded images. Deep Spitzer IRAC and MIPS provide the first detection of a large sample of Lyman Break Galaxies in the mid- to far-infrared range. We are therefore able to study and compare the UV and TIR emission of Lyman Break Galaxies. We find that about 15% of the LBG sample are strong emitters at 24 microns (Red LBGs). Most of them are Luminous IR Galaxies (LIRGs) while the rest (Blue LBGs) are undetected at the 83 microJy level of MIPS GTO image. We find that Blue LBGs have a Spectral Energy Distribution similar to high redshift Lyman Break Galaxies. Finally, the dust-to-FUV ratio of this sample is compared with similar ratios at z=1 and z~2. This work suggests an evolution (decrease) of the dust-to-FUV ratio with the redshift.

Based on GALEX and IRAS/Spitzer datasets, we have found that both FUV and FIR luminosity functions (LFs) show a strong evolution from z=0 to z=1, but the FIR LF evolves much stronger than the FUV one. Consequently, the FIR/FUV luminosity density ratio increases from 4 (z=0) to 15 (z=1). It means that more than 80% of the star-forming activity in the Universe is hidden by dust at z=1. To explore this issue further, we have performed a combined analysis of the galaxy sample in FUV and FIR. For the Local Universe we used GALEX-IRAS sample, whereas at z=1 we used the Lyman-break galaxy sample selected by GALEX bands constructed by Burgarella et al. (2005), which is known to be representative of visible (i.e., non-obscured) star-forming galaxies at z=1. From these datasets, we constructed the LFs of the FUV-selected galaxies by the survival analysis to, take into account the upper-limit data properly. We discovered that the FIR LF of the Lyman-break galaxies show a significant evolution comparing with the local FIR LF, but it is a factor of 2-3 lower than the global FIR LF (Le Floc'h et al. 2005). This indicates that the evolution of visible galaxies is not strong enough to explain the drastic evolution of the FIR LF. Namely, a FIR-luminous, rapidly diminishing population of galaxies is required.

We select far-infrared (FIR-60 microns) and far-ultraviolet (FUV-1530 A) samples of nearby galaxies in order to discuss the biases encountered by monochromatic surveys (FIR or FUV). Very different volumes are sampled by each selection and much care is taken to apply volume corrections to all the analyses. The distributions of the bolometric luminosity of young stars are compared for both samples: they are found to be consistent with each other for galaxies of intermediate luminosities but some differences are found for high (>5 10^{10} L_sun) luminosities. The shallowness of the IRAS survey prevents us from securing comparison at low luminosities (<2 10^9 L_sun). The ratio of the total infrared (TIR) luminosity to the FUV luminosity is found to increase with the bolometric luminosity in a similar way for both samples up to 5 10^{10} L_sun. Brighter galaxies are found to have a different behavior according to their selection: the L_TIR/L_FUV ratio of the FUV-selected galaxies brighter than 5 10^{10} L_sun reaches a plateau whereas L_TIR/L_FUV continues to increase with the luminosity of bright galaxies selected in FIR. The volume-averaged specific star formation rate (SFR per unit galaxy stellar mass, SSFR) is found to decrease toward massive galaxies within each selection. The SSFR is found to be larger than that measured for optical and NIR-selected sample over the whole mass range for the FIR selection, and for masses larger than 10^{10} M_sun for the FUV selection. Luminous and massive galaxies selected in FIR appear as active as galaxies with similar characteristics detected at z ~ 0.7.

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

We discuss dust properties in the interstellar medium (ISM) of nearby normal galaxies, by comparing observations in the ultraviolet (UV) with simulations by a radiative transfer model. The observed UV colours of nearby galaxies show a reddening relative to their expected intrinsic colours. Some authors argued that the Milky Way dust cannot reproduce the reddening because of the prominent 2175 \AA absorption bump. Other authors proposed a reduction mechanism of the bump strength in an {\it attenuation law} derived from the ratio of the observed intensity to the intrinsic one through an age-selective attenuation (i.e., young stars are more attenuated selectively). We newly find that the wavelength dependence of the scattering albedo also has a strong effect on the UV colour; an albedo decreasing toward shorter wavelengths (except for the absorption bump range) produces a significant UV reddening. After comparing the observed UV colours of nearby normal galaxies with those expected from radiative transfer simulations assumed several dust models, we find two sorts of dust suitable for these galaxies: (1) dust with a bump and a smaller albedo for a shorter wavelength (except for the bump range), and (2) dust without any bump but with an almost constant albedo. If very small carbonaceous grains responsible for the common unidentified infrared emission band are also the bump carrier, the former dust is favorable. Finally, we derive mean attenuation laws of various dust models as a function of the UV attenuation, and derive some relations between the UV attenuation and observable/theoretical quantities.

We study two samples of local galaxies, one is UV (GALEX) selected and the other FIR (IRAS) selected, to address the question whether UV and FIR surveys see the two sides ('bright' and 'dark') of the star formation of the same population of galaxies or two different populations of star forming galaxies. No significant difference between the L$_{tot}$ ($=L_{60}+L_{FUV}$) luminosity functions of the UV and FIR samples is found. Also, after the correction for the `Malmquist bias' (bias for flux limited samples), the FIR-to-UV ratio v.s. L$_{tot}$ relations of the two samples are consistent with each other. In the range of $9 \la \log(L_{tot}/L_\sun) \la 12$, both can be approximated by a simple linear relation of $\log (L_{60}/L_{FUV})=\log(L_{tot}/L_\sun)-9.66$. These are consistent with the hypothesis that the two samples represent the same population of star forming galaxies, and their well documented differences in L$_{tot}$ and in FIR-to-UV ratio are due only to the selection effect. A comparison between the UV luminosity functions shows marginal evidence for a population of faint UV galaxies missing in the FIR selected sample. The contribution from these 'FIR-quiet' galaxies to the overall UV population is insignificant, given that the K-band luminosity functions (i.e. the stellar mass functions) of the two samples do not show any significant difference.

We investigate the possibility to search for metal-poor high column density (> 10^{23} cm^{-2}) clouds at high redshift (z) by using gamma ray burst (GRB) afterglows. Such clouds could be related to primeval galaxies which may cause a burst of star formation. We show that a large part of hydrogen is in molecular form in such a high column density environment. Therefore, hydrogen molecules (H2) rather than hydrogen atoms should be searched for. Then we show that infrared H2 lines are detectable for metal-poor (< 0.01 solar metallicity) high column density (log NH [cm^{-2}] > 23.5) systems at high-z without suffering dust extinction. Optical properties of dust in infrared could also be constrained by observations of high column density systems. Some possible scenarios of producing high column density systems are finally discussed in the context of galaxy evolution.

Dust plays crucial roles in galaxy formation and evolution. In the early epoch of galaxy evolution dust is only supplied by supernovae (SNe). With the aid of a new physical model of dust production by SNe, we constructed a model of dust emission from forming galaxies. We show the evolution of the spectral energy distribution (SED). Then we adopt this model to a local starbursting dwarf galaxy SBS 0335$-$052. Further we discuss the SEDs of high redshift galaxies, and consider their observational feasibility.

This work presents the main ultraviolet (UV) and far-infrared (FIR) properties of two samples of nearby galaxies selected from the GALEX ($\lambda = 2315$\AA, hereafter NUV) and IRAS ($\lambda = 60\mu$m) surveys respectively. They are built in order to get detection at both wavelengths for most of the galaxies. Star formation rate (SFR) estimators based on the UV and FIR emissions are compared. Systematic differences are found between the SFR estimators for individual galaxies based on the NUV fluxes corrected for dust attenuation and on the total IR luminosity. A combined estimator based on NUV and IR luminosities seems to be the best proxy over the whole range of values of SFR. Although both samples present similar average values of the birthrate parameter b, their star-formation-related properties are substantially different: NUV-selected galaxies tend to show larger values of $b$ for lower masses, SFRs and dust attenuations, supporting previous scenarios for the star formation history (SFH). Conversely, about 20% of the FIR-selected galaxies show high values of $b$, SFR and NUV attenuation. These galaxies, most of them being LIRGs and ULIRGs, break down the downsizing picture for the SFH, however their relative contribution per unit volume is small in the local Universe. Finally, the cosmic SFR density of the local Universe is estimated in a consistent way from the NUV and IR luminosities.

We present the first large, unbiased sample of Lyman Break Galaxies (LBGs) at z ~ 1. Far ultraviolet-dropout (1530 A) galaxies in the Chandra Deep Field South have been selected using GALEX data. This first large sample in the z ~ 1 universe provides us with a high quality reference sample of LBGs. We analyzed the sample from the UV to the IR using GALEX, SPITZER, ESO and HST data. The morphology (obtained from GOODS data) of 75 % of our LBGs is consistent with a disk. The vast majority of LBGs with an IR detection are also Luminous Infrared Galaxies (LIRGs). As a class, the galaxies not detected at 24 microns are an order of magnitude fainter relative to the UV compared with those detected individually, suggesting that there may be two types of behavior within the sample. For the IR-bright galaxies, there is an apparent upper limit for the UV dust attenuation and this upper limit is anti-correlated with the observed UV luminosity. Previous estimates of dust attenuations based on the ultraviolet slope are compared to new ones based on the FIR/UV ratio (for LBGs detected at 24 microns), which is usually a more reliable estimator. Depending on the calibration we use to estimate the total IR luminosity, beta-based attenuations A_{FUV} are larger by 0.2 to 0.6 mag. than the ones estimated from FIR/UV ratio. Finally, for IR-bright LBGs, median estimated beta-based SFRs are 2-3 times larger than the total SFRs estimated as SFR_{TOT} = SFR_{UV} + SFR_{IR} while IR-based SFRs provide values below SFR_{TOT} by 15 - 20 %. We use a stacking method to statistically constrain the 24 microns flux of LBGs non individually detected. The results suggest that these LBGs do not contain large amounts of dust.

We constructed a local luminosity function (LF) of galaxies using a flux-limited sample (S_170 > 0.195Jy) of 55 galaxies at z < 0.3 taken from the ISO FIRBACK survey at 170um. The overall shape of the 170-um LF is found to be different from that of the total 60-um LF (Takeuchi et al. 2003): the bright end of the LF declines more steeply than that of the 60-um LF. This behavior is quantitatively similar to the LF of the cool subsample of the IRAS PSCz galaxies. We also estimated the strength of the evolution of the LF by assuming the pure luminosity evolution (PLE): L(z) \propto (1+z)^Q. We obtained Q=5.0^{+2.5}_{-0.5} which is similar to the value obtained by recent Spitzer observations, in spite of the limited sample size. Then, integrating over the 170-um LF, we obtained the local luminosity density at 170um, \rho_L(170um). A direct integration of the LF gives \rho_L(170um) = 1.1 \times 10^8 h Lsun Mpc^{-3}, whilst if we assume a strong PLE with Q=5, the value is 5.2 \times 10^7 h Lsun Mpc^{-3}. This is a considerable contribution to the local FIR luminosity density. By summing up with other available infrared data, we obtained the total dust luminosity density in the Local Universe, \rho_L(dust)=1.1 \times 10^8 h Lsun Mpc^{-3}. Using this value, we estimated the cosmic star formation rate (SFR) density hidden by dust in the Local Universe. We obtained \rho_SFR(dust) \simeq 1.1-1.2 h \times 10^{-2} Msun yr^{-1} Mpc^{-3}$, which means that 58.5% of the star formation is obscured by dust in the Local Universe.

The interplay between dust and molecules is of fundamental importance in early galaxy evolution. First we present the prediction for the dust emission from forming galaxies. Then we discuss the observational strategy for molecules in early galaxies by infrared absorption lines of a bright continuum source behind the clouds. By combining these two approaches, we will be able to have a coherent picture of the very early stage of galaxy evolution.

The first large sample of Lyman Break Galaxies (LBGs) at z~1 is selected from GALEX Far-UV and Near-UV images of the Chandra Deep Field South. This multi-wavelength study is based on the wealth of data available in this field: SPITZER-MIPS to estimate total IR luminosities and dust attenuations, redshifts from COMBO-17, EIS and GOODS to check the validity of dropout identifications and GOODS to derive the morphology of about 40 of our LBGs. The advantages of this "nearby" LBG sample are that images reach fainter fluxes and surface brightnesses, allowing a deeper detection at all wavelengths and a better morphological analysis. Moreover, the 24-mic fluxes (i.e. 12-mic rest-frame) permit a good estimation of dust attenuations and total IR luminosities. The main results are that the vast majority of our LBGs are also Luminous Infrared Galaxies (LIRGs). The morphology of 75% of our LBGs is consistent with a disk. Previous estimates of dust attenuations based on the ultraviolet slope are likely to be too large by a factor of about 2, which implies that estimated star formation rates were too large by the same amount. This sample in the z~1 universe provides us with a high quality reference sample of LBGs.

The ratio between far-ultraviolet (FUV) and infrared (IR) luminosity densities from z=0 to z=1 is discussed by using the luminosity functions (LFs) of both wavelengths. The FUV LF (z=0-1) based on GALEX has been reported by Arnouts et al. (2005), whilst for the IR LF, we used the IRAS PSCz 60-um LF for the local universe and the Spitzer 15-um LF at higher-z as used by Le Floch et al. (2005). Both luminosity densities show a significant evolutionary trend, but the IR evolves much faster than the FUV. Consequently, the ratio $\rho_{dust}/\rho_{FUV}$ increases toward higher-z, from 4 (local) to 15 (z = 1). It is also shown that more than 70% of the star formation activity in the universe is obscured by dust at 0.5 < z < 1.2.

In the early epoch of galaxy evolution, dust is only supplied by supernovae (SNe). With the aid of a new physical model of dust production by SNe developed by Nozawa et al. (2003) (N03), we constructed a model of dust emission from forming galaxies on the basis of the theoretical framework of Takeuchi et al. (2003) (T03). N03 showed that the produced dust species depends strongly on the mixing within SNe. We treated both unmixed and mixed cases and calculated the infrared (IR) spectral energy distribution (SED) of forming galaxies for both cases. Our model SED is less luminous than the SED of T03 model by a factor of 2-3. The difference is due to our improved treatment of UV photon absorption cross section, as well as different grain size and species newly adopted in this work. The SED for the unmixed case is found to have an enhanced near to mid-IR (N-MIR) continuum radiation in its early phase of the evolution (age < 10^{7.25} yr) compared with that for the mixed case. The strong N--MIR continuum is due to the emission from Si grains, which only exist in the species of the unmixed dust production. We also calculated the IR extinction curves for forming galaxies. Then we calculated the SED of a local starbursting dwarf galaxy SBS 0335-052. Our present model SED naturally reproduced the strong N--MIR continuum and the lack of cold FIR emission of SBS 0335-052. We found that only the SED of unmixed case can reproduce the NIR continuum of this galaxy. We then made a prediction for the SED of another typical star-forming dwarf, I Zw 18. We also presented the evolution of the SED of LBGs. Finally, we discussed the possibility of observing forming galaxies at z > 5.

We construct a protostellar disk model that takes into account the combined effect of viscous evolution, photoevaporation and the differential radial motion of dust grains and gas. For T Tauri disks, the lifetimes of dust disks that are mainly composed of millimeter sized grains are always shorter than the gas disks' lifetimes, and become similar only when the grains are fluffy (density < 0.1 g cm^{-3}). If grain growth during the classical T Tauri phase produces plenty of millimeter sized grains, such grains completely accrete onto the star in 10^7 yr, before photoevaporation begins to drain the inner gas disk and the star evolves to the weak line T Tauri phase. In the weak line phase, only dust-poor gas disks remain at large radii (> 10 AU), without strong signs of gas accretion nor of millimeter thermal emission from the dust. For Herbig AeBe stars, the strong photoevaporation clears the inner disks in 10^6 yr, before the dust grains in the outer disk migrate to the inner region. In this case, the grains left behind in the outer gas disk accumulate at the disk inner edge (at 10-100 AU from the star). The dust grains remain there even after the entire gas disk has been photoevaporated, and form a gas-poor dust ring similar to that observed around HR 4796A. Hence, depending on the strength of the stellar ionizing flux, our model predicts opposite types of products around young stars. For low mass stars with a low photoevaporation rate, dust-poor gas disks with an inner hole would form, whereas for high mass stars with a high photoevaporation rate, gas-poor dust rings would form. This prediction should be examined by observations of gas and dust around weak line T Tauri stars and evolved Herbig AeBe stars.

From millimeter observations of classical T Tauri stars, it is suggested that dust grains in circumstellar disks have grown to millimeter size or larger. However, gas drag on such large grains induces rapid accretion of the dust. We examine the evolution of dust disks composed of millimeter sized grains, and show that rapid accretion of the dust disk causes attenuation of millimeter continuum emission. If a dust disk is composed mainly of grains of 1 cm to 1 m, its millimeter emission goes off within 10^6 yr. Hence, grains in this size range cannot be a main population of the dust. Considering our results together with grain growth suggested by the millimeter continuum observations, we expect that the millimeter continuum emission of disks comes mainly from grains in a narrow size range of [1 mm -1 cm]. This suggests either that growth of millimeter sized grains to centimeter size takes more than 10^6 yr, or that millimeter sized grains are continuously replenished. In the former case, planet formation is probably difficult, especially in the outer disks. In the latter case, reservoirs of millimeter grains are possibly large (> 10 m) bodies, which can reside in the disk more than 10^6 yr. Constraints on the grain growth time-scale are discussed for the above two cases.

The infrared (IR) emission plays a crucial role for understanding the star formation in galaxies hidden by dust. We first examined four estimators of the IR luminosity of galaxies, L_fir (Helou et al. 1988), L_tir (Dale et al. 2001), revised version of L_tir (Dale & Helou 2002) (we denote L_tir2), and L_ir (Sanders & Mirabel 1996) by using the observed SEDs of well-known galaxies. We found that L_ir provides excellent estimates of the total IR luminosity for a variety of galaxy SEDs. The performance of L_tir2 was also found to be very good. Using L_ir, we then statistically analyzed the IRAS PSCz galaxy sample (Saunders et al. 2000) and found useful formulae relating the MIR monochromatic luminosities [L(12um) and L(25um)], and L_ir. For this purpose we constructed a subsample of 1420 galaxies with all IRAS four band (12, 25, 60, and 100um) flux densities. We found linear relations between L_ir and MIR luminosities, L(12um) and L(25um). The prediction error with 95-% confidence level is a factor of 4-5. Hence, these formulae are useful for the estimation of the total IR luminosity only from 12um or 25um observations. We further tried to make an `interpolation' formula for galaxies at 0<z<1. For this purpose we construct the formula of the relation between 15-um luminosity and the total IR luminosity. We conclude that the 15-um formula can be used as an estimator of the total IR luminosity from 24um observation of galaxies at z \simeq 0.6.

We investigate the extinction curves of young galaxies in which dust is supplied from Type II supernovae (SNe II) and/or pair instability supernovae (PISNe). We adopt Nozawa et al. (2003) for compositions and size distribution of grains formed in SNe II and PISNe. We find that the extinction curve is quite sensitive to internal metal mixing in supernovae (SNe). The extinction curves predicted from the mixed SNe are dominated by SiO2 and is characterised by steep rise from infrared to ultraviolet (UV). The dust from unmixed SNe shows shallower extinction curve, because of the contribution from large-sized (~ 0.1 um) Si grains. However, the progenitor mass is important in unmixed SNe II: If the progenitor mass is smaller than ~ 20 Msun, the extinction curve is flat in UV; otherwise, the extinction curve rises toward the short wavelength. The extinction curve observed in a high-redshift quasar (z=6.2) favours the dust production by unmixed SNe II. We also provide some useful observational quantities, so that our model might be compared with future high-z extinction curves.

Astrometric detection of a stellar wobble on the plane of the sky will provide us a next breakthrough in searching extrasolar planets. The Space Interferometry Mission (SIM) is expected to achieve a high-precision astrometry as accurate as 1 micro-as, which is precise enough to discover a new-born Jupiter mass planet around a pre-main-sequence (PMS) star in the Taurus-Auriga star forming region. PMS stars, however, have circum-stellar disks that may be obstacles to the precise measurement of the stellar position. We present results on disk influences to the stellar wobble. The density waves excited by a planet move both of the disk's mass center and the photo-center. The motion of the disk mass center induces an additional wobble of the stellar position, and the motion of the disk photo-center causes a contamination in the measurement of the stellar position. We show that the additional stellar motion dynamically caused by the disk's gravity is always negligible, but that the contamination of the disk light can interfere with the precise measurement of the stellar position, if the planet's mass is smaller than ~10 Jupiter mass. The motion of the disk photo-center is sensitive to a slight change in the wave pattern and the disk properties. Measurements by interferometers are generally insensitive to extended sources such as disks. Because of this property SIM will not suffer significant contaminations of the disk light, even if the planet's mass is as small as 1 Jupiter mass.

Lyman-break galaxies (LBGs) contain a non-negligible amount of dust. Takeuchi et al. (2003)(T03) constructed a model of the infrared spectral energy distribution (SED) for very young galaxies by taking into account the dust size distribution in the early stage of galaxy evolution, which can be different from that of present-day evolved galaxies. We applied T03 model to LBGs and constructed their expected SED. In order to examine the grain size distribution of dust, we calculated the SEDs based on two distinct type of the distribution models: a single-sized distribution and a power-law distribution with a slope of $dN/da \propto a^{-3.5}$. We found that the single-sized and power-law dust size distributions yield a very similar detectability of LBGs at the submillimetre (submm). We also found that galaxies with a power-law dust distribution have much less flux at mid-infrared (MIR) than the other. Then, we applied the model to a gravitationally lensed LBG \cb (cB58), a unique probe of the dust emission from LBGs. Our model framework well reproduced the hot dust temperature under a natural assumption for this galaxy. We also examined the detectability of LBGs at submm wavelengths in an eight-hour deep survey by ALMA. By integrating over their redshifted SEDs with the observed luminosity functions, we obtained the contribution of LBGs to the cosmic infrared background radiation (CIRB). Although they have non-negligible amount of dust, their contribution was found to be small, especially in the FIR $\sim 200 \mu$m. Thus, we need a strongly obscured population of galaxies which contains a large amount of star formation, at some epoch in the history of the universe.

Source confusion has been a long-standing problem in the astronomical history. In the previous formulation, sources are assumed to be distributed homogeneously on the sky. This fundamental assumption is not realistic in many applications. In this work, by making use of the point field theory, we derive general analytic formulae for the confusion problems with arbitrary distribution and correlation functions. As a typical example, we apply these new formulae to the source confusion of infrared galaxies. We first calculate the confusion statistics for power-law galaxy number counts as a test case. When the slope of differential number counts, \gamma, is steep, the confusion limits becomes much brighter and the probability distribution function (PDF) of the fluctuation field is strongly distorted. Then we estimate the PDF and confusion limits based on the realistic number count model for infrared galaxies. The gradual flattening of the slope of the source counts makes the clustering effect rather mild. Clustering effects result in an increase of the limiting flux density with \sim 10%. In this case, the peak probability of the PDF decreases up to \sim 15% and its tail becomes heavier.

The small grain sizes produced by Type II supernova (SN II) models in young, metal-poor galaxies make the appearance of their infrared (IR) spectral energy distribution (SED) quite different from that of nearby, older galaxies. To study this effect, we have developed a model for the evolution of dust content and the IR SED of low-metallicity, extremely young galaxies based on Hirashita et al. (2002). We find that, even in the intense ultraviolet (UV) radiation field of very young galaxies, small silicate grains are subject to stochastic heating resulting in a broad temperature distribution and substantial MIR continuum emission. Larger carbonaceous grains are in thermal equilibrium at T \simeq 50 - 100K, and they also contribute to the MIR. We present the evolution of SEDs and IR extinction of very young, low-metallicity galaxies. The IR extinction curve is also shown. In the first few Myrs, the emission peaks at \lambda \sim 30-50um at later times dust self-absorption decreases the apparent grain temperatures, shifting the bulk of the emission into the submillimetre band. We successfully apply the model to the IR SED of a low metallicity (1/41 Z_\odot) dwarf galaxy SBS0335-052. We find the SED, optical properties and extinction of the star forming region to be consistent with a very young and compact starburst. We also predict the SED of another extremely low-metallicity galaxy, I Zw 18, for future observational tests. Some prospects for future observations are discussed.

We study the outflow of dust particles on the surface layers of optically thick disks. At the surface of disks around young stars, small dust particles (size < 10 micron) experience stellar radiation pressure support and orbit more slowly than the surrounding gas. The resulting tail-wind imparts energy and angular momentum to the dust particles, moving them outward. This outflow occurs in the thin surface layer of the disk that is exposed to starlight, and the outward mass flux is carried primarily by particles of size ~0.1 micron. Beneath the irradiated surface layer, dust particles experience a head-wind, which drives them inward. For the specific case of a minimum-mass-solar-nebula, less than a thousandth of the dust mass experiences outward flow. If the stellar luminosity is 15 times brighter than the sun, however, or if the gas disk mass is as small as ~100 M_earth, then the surface outflow can dominate the inward flux in certain radial ranges, leading to the formation of rings or gaps in the dust disks.

We estimated the luminosity function (LF) of IRAS galaxies in the PSCz catalogue. The faint end of the PSCz LF is slightly steeper than that of the LF derived by Saunders et al. (1990; S90). Using an analytical form for the LF used by S90, we obtain the following parameters: \alpha = 1.23 \pm 0.04, L_*=(8.85 \pm 1.75) \times 10^8 h^{-2} L_\odot, \sigma =0.724 \pm 0.010, and \phi_* = (2.34 \pm 0.30) \times 10^{-2} h^3 Mpc^{-3}. We also examined the evolution in the sample by a simple assumption \phi_*(z) \propto (1+z)^P, and found P=3.40 \pm 0.70. It does not affect the three parameters, \alpha, L_*, and \sigma, but \phi_*(z=0) is overestimated up to \sim 15% if we ignore evolution. We estimated the temperature dependence of the LF. The LFs of warm and cool galaxies are quite different: the LF of warm galaxies has a very steep faint end with \alpha =1.37. We also discuss a lump found at the brightest end of the LF.

We present millimeter imaging observations in the 1200 micron continuum of the disk around beta Pictoris. With the 25 arcsec beam, the beta Pic disk is unresolved perpendicularly to the disk plane (< 10 arcsec), but slightly resolved in the northeast-southwest direction (26 arcsec). Peak emission is observed at the stellar position. A secondary maximum is found 1000 AU along the disk plane in the southwest, which does not positionally coincide with a similar feature reported earlier at 850 micron. Arguments are presented which could be seen in support of the reality of these features. The observed submm/mm emission is consistent with thermal emission from dust grains, which are significantly larger than those generally found in the interstellar medium, including mm-size particles, and thus more reminiscent of the dust observed in protostellar disks. Modelling the observed scattered light in the visible and the emission in the submm/mm provides evidence for the particles dominating the scattering in the visible/NIR and those primarily responsible for the thermal emission at longer wavelengths belonging to different populations.

We study the phenomenological implications of the classical limit of the "stringy" commutation relations [x_i,p_j]=i hbar[(1+beta p^2) delta_{ij} + beta' p_i p_j]. In particular, we investigate the "deformation" of Kepler's third law and apply our result to the rotation curves of gas and stars in spiral galaxies.

Before the end of 2002 will be launched the GALEX satellite (a NASA/SMEX project) which will observe all the sky in Ultraviolet (UV) through filters at 1500 and 2300 A down to m(AB) 21. In 2004 will be launched the ASTRO-F satellite which will perform an all sky survey at Far-Infrared (FIR) wavelengths. The cross-correlation of both surveys will lead to very large samples of galaxies for which FIR and UV fluxes will be available. Using the FIR to UV flux ratio as a quantitative tracer of the dust extinction we will be able to measure the extinction in the nearby universe (z<0.2) and to perform a statistically significant analysis of the extinction as a function of galactic properties. Of particular interest is the construction of pure FIR and UV selected samples for which the extinction will be measured as templates for the observation of high redshift galaxies.

We continue our investigation of the phenomenological implications of the "deformed" commutation relations [x_i,p_j]=i hbar[(1 + beta p^2) delta_{ij} + beta' p_i p_j]. These commutation relations are motivated by the fact that they lead to the minimal length uncertainty relation which appears in perturbative string theory. In this paper, we consider the effects of the deformation on the classical orbits of particles in a central force potential. Comparison with observation places severe constraints on the value of the minimum length.

We study the radial migration of dust particles in accreting protostellar disks analogous to the primordial solar nebula. This study takes account of the two dimensional (radial and normal) structure of the disk gas, including the effects of the variation in the gas velocity as a function of distance from the midplane. It is shown that the dust component of disks accretes slower than the gas component. At high altitude from the disk midplane, the gas rotates faster than particles because of the inward pressure gradient force, and its drag force causes particles to move outward in the radial direction. Viscous torque induces the gas within a scale height from the disk midplane to flow outward, carrying small (size < 100 micron at 10 AU) particles with it. Only particles at intermediate altitude or with sufficiently large sizes (> 1 mm at 10 AU) move inward. When the particles' radial velocities are averaged over the entire vertical direction, particles have a net inward flux. At large distances from the central star, particles migrate inward with a velocity much faster than the gas accretion velocity. However, their inward velocity is reduced below that of the gas in the inner regions of the disk. The rate of velocity decrease is a function of the particles' size. While larger particles retain fast accretion velocity until they approach closer to the star, 10 micron particles have slower velocity than the gas in the most part of the disk (r < 100 AU). This differential migration of particles causes the size fractionation. Dust disks composed mostly of small particles (size < 10 micron) accrete slower than gas disks, resulting in the increase in the dust-gas ratio during the gas accretion phase.

We propose a new empirical method to estimate the total far-infrared flux of galaxies from the spectral energy distribution (SED) at wavelengths lambda shorter than 100 microns. It is difficult to derive the total far-infrared luminosity from only the IRAS data, though it is one of the most important properties of galaxies. Observations by Infrared Telescope in Space (IRTS) indicate that the SED of the diffuse emission from the Galactic plane in this wavelength region can be derived from the 60 microns to 100 microns color. This empirical SED relation was improved in order to obtain a better fit to the Galactic plane data for I(60 microns)/I(100 microns) > 0.6, and applied to 96 IRAS galaxies for which ISOPHOT and KAO data are available at lambda > 100 microns. As a result, the empirical relation describes well the far-infrared (FIR) SED for a majority of the galaxies. Additionally, the total FIR flux for lambda longer than 40 microns is derived from the flux densities at 60 and 100 microns by using this model. For the 96 IRAS galaxies, the uncertainty in the total far-infrared flux of the present method is 26%. The present method is more accurate than the previous one widely used to derive the total infrared flux from the IRAS 60 and 100 microns data.

We give an explanation for the origin of various properties observed in local infrared galaxies, and make predictions for galaxy counts and cosmic background radiation (CBR), by a new model extended from that for optical/near-infrared galaxies. Important new characteristics of this study are that (1) mass scale dependence of dust extinction is introduced based on the size-luminosity relation of optical galaxies, and that (2) the big grain dust temperature T_dust is calculated based on a physical consideration for energy balance, rather than using the empirical relation between T_dust and total infrared luminosity L_IR found in local galaxies, which has been employed in most of previous works. Consequently, the local properties of infrared galaxies, i.e., optical/infrared luminosity ratios, L_IR-T_dust correlation, and infrared luminosity function are outputs predicted by the model. Our model indeed reproduces these local properties reasonably well. We then found considerably different results for MIR-submm counts and CBR from most of previous works based on the empirical L_IR-T_dust relation; especially, it is shown that the dust temperature of starbursting primordial elliptical galaxies is expected to be very high (40-80K). This indicates that intense starbursts of forming elliptical galaxies should have occurred at z~2-3, in contrast to the previous results that significant starbursts beyond z~1 tend to overproduce the far-infrared (FIR) CBR detected by COBE/FIRAS. On the other hand, our model predicts that the mid-infrared (MIR) flux from warm/nonequilibrium dust is relatively weak in such galaxies making FIR CBR, and this effect reconciles the prima facie conflict between the upper limit on MIR CBR from TeV gamma-ray observations and the COBE detections of FIR CBR. (abridged)

We investigate the observability of hydrogen molecules in absorption. The absorption efficiency of the hydrogen molecules become comparable with or larger than that of the dust grains in the metal-poor condition expected in the early Universe. If we can use bright infrared continuum sources behind the molecular gas clouds, the absorption measurement of the hydrogen molecules will be an important technique to explore the primordial gas clouds that are contracting into first-generation objects.

We re-examine the technique of determining the redshifts of galaxies from the ratio of the submillimeter-to-radio continuum flux densities based on a recently published catalog of 850~$\mu$m sources. We derived the expected variation of this ratio as a function of redshift incorporating the expected average luminosity and the spectral energy distribution (SED) of dust emission and the radio continuum. We find that the existing data for most of the high redshift ($z \gtsim 1$) sources correspond to our new calculation. Amongst the well-identified sources, there is none with an index significantly higher than predicted. Sources which have an index lower than predicted are either within the error zone or much lower, the latter presumably having AGN-dominated radio-continuum emission. We find the median redshift to be $\sim 2$, which is consistent with that deduced by previous work ($\sim 3$) within the error. We also discuss the various systematic effects that can affect the accuracy of the redshift estimate. We examine other methods of redshift estimation, like photometric ratio in the submillimeter and locating the peak of the SED in the rest system of the objects. We conclude that while the various methods are helpful in identifying high-redshift objects and making a crude estimate of the redshift, they are not, at present, accurate enough for a detailed study of redshift distribution of the submillimeter galaxies.

The virial mass ($M_{\rm vir}$)-metallicity relation among the Local Group dwarf spheroidal galaxies (dSphs) is examined. Hirashita, Takeuchi, & Tamura showed that the dSphs can be divided into two distinct classes with respect to the relation between their virial masses and luminosities: low-mass ($M_{\rm vir} \la 10^8 M_\odot$) and high-mass ($M_{\rm vir} \ga 10^8 M_\odot$) groups. We see that both the mass-metallicity and the mass-luminosity relations of the high-mass dSphs are understood as a low-mass extension of giant ellipticals. On the contrary, we find that the classical galactic-wind model is problematic to apply to the low-mass dSphs, whose low binding energy is comparable to that released by several supernova explosions. A strongly regulated star formation in their formation phase is required to reproduce their observed metallicity. Such regulation is naturally expected in a gas cloud with the primordial elemental abundance according to Nishi & Tashiro. A significant scatter in the mass-metallicity relation for the low-mass dSphs is also successfully explained along with the scenario of Hirashita and coworkers. We not only propose a new picture for a chemical enrichment of the dSphs, but also suggest that the mass-metallicity and the mass-luminosity relations be understood in a consistent context.

We investigate what we can learn about galaxy formation and evolution from the data which will be obtained by the forthcoming large submillimeter/millimeter facilities, mainly by the Atacama Submillimeter Telescope Experiment (ASTE) and the Atacama Large Millimeter Array/Large Millimeter and Submillimeter Array (ALMA/LMSA). We first calculate the source counts from 350 um to 3 mm using the empirical infrared galaxy number count model of Takeuchi et al. (2001). Based on the number counts, we evaluate the source confusion and determine the confusion limit at various wavebands as a function of the characteristic beam size. At submillimeter wavelengths, source confusion with the 10-15-m class facilities becomes severe at 0.1 to 1 mJy level, and astrometry and flux measurements will be difficult. However, we show that very a large-area survey of submillimeter sources brighter than 10-50 mJy can provide a unique constraint on infrared galaxy evolution at z = 1-2. We also find that the 5\sigma-confusion limit of the LMSA is fainter than 1 \mu Jy. The source counts at such faint flux levels give important information of the epoch of galaxy formation. We then show that multiband photometry from the infrared (by ASTRO-F) to the millimeter can be utilized as a redshift estimator. In addition, we compare the observed 1.4, 5, and 8-GHz source counts with our model counts to examine the contribution of star forming galaxies to the faint radio galaxies.

We model the star formation history (SFH) and the chemical evolution of the Galactic disk by combining an infall model and a limit-cycle model of the interstellar medium (ISM). Recent observations have shown that the SFH of the Galactic disk violently variates or oscillates. We model the oscillatory SFH based on the limit-cycle behavior of the fractional masses of three components of the ISM. The observed period of the oscillation ($\sim 1$ Gyr) is reproduced within the natural parameter range. This means that we can interpret the oscillatory SFH as the limit-cycle behavior of the ISM. We then test the chemical evolution of stars and gas in the framework of the limit-cycle model, since the oscillatory behavior of the SFH may cause an oscillatory evolution of the metallicity. We find however that the oscillatory behavior of metallicity is not prominent because the metallicity reflects the past integrated SFH. This indicates that the metallicity cannot be used to distinguish an oscillatory SFH from one without oscillations.

We analyze the dynamics of gas-dust coupling in the presence of stellar radiation pressure in circumstellar gas disks, which are in a transitional stage between the gas-dominated, optically thick, primordial nebulae, and the dust-dominated, optically thin Vega-type disks. Dust undergo radial migration, seeking a stable equilibrium orbit in corotation with gas. The migration of dust gives rise to radial fractionation of dust and creates a variety of possible observed disk morphologies, which we compute by considering the equilibrium between the dust production and the dust-dust collisions removing particles from their equilibrium orbits. Sand-sized and larger grains are distributed throughout most of the gas disk, with concentration near the gas pressure maximum in the inner disk. Smaller grains (typically in the range of 10 to 200 micron) concentrate in a prominent ring structure in the outer region of the gas disk (presumably at radius 100 AU), where gas density is rapidly declining with radius. The width and density, as well as density contrast of the dust ring with respect to the inner dust disk depend on the distribution of gas. Our results open the prospect for deducing the distribution of gas in circumstellar disks by observing their dust. We have qualitatively compared our models with two observed transitional disks around HR 4796A and HD 141569A. Dust migration can result in observation of a ring or a bimodal radial dust distribution, possibly very similar to the ones produced by gap-opening planet(s) embedded in the disk, or shepherding it from inside or outside. We conclude that a convincing planet detection via dust imaging should include specific non-axisymmetric structure following from the dynamical simulations of perturbed disks.

Recently reported infrared (IR) galaxy number counts and cosmic infrared background (CIRB) all suggest that galaxies have experienced a strong evolution sometime in their lifetime. We statistically estimate the galaxy evolution history from these data. We find that an order of magnitude increase of the far-infrared (FIR) luminosity at redshift z = 0.5 - 1.0 is necessary to reproduce the very high CIRB intensity at 140 um reported by Hauser et al. (1998). z \sim 0.75 and decreases to, even at most, a factor of 10 toward z \sim 5, though many variants are allowed within these constraints. This evolution history also satisfies the constraints from the galaxy number counts obtained by IRAS, ISO and, roughly, SCUBA. The rapid evolution of the comoving IR luminosity density required from the CIRB well reproduces the very steep slope of galaxy number counts obtained by ISO. We also estimate the cosmic star formation history (SFH) from the obtained FIR luminosity density, considering the effect of the metal enrichment in galaxies. The derived SFH increases steeply with redshift in 0 < z < 0.75, and becomes flat or even declines at z > 0.75. This is consistent with the SFH estimated from the reported ultraviolet luminosity density. In addition, we present the performance of the Japanese ASTRO-F FIR galaxy survey. We show the expected number counts in the survey. We also evaluate how large a sky area is necessary to derive a secure information of galaxy evolution up to z \sim 1 from the survey, and find that at least 50 - 300 deg^2 is required.

We analyze the star formation history (SFH) of the Galactic disk by using an infall model. Based on the observed SFH of the Galactic disk, we first determine the timescales of the gas infall into the Galactic disk ($t_{\rm in}$) and that of the gas consumption to form stars ($t_{\rm sf}$). Since each of the two timescales does not prove to be determined independently from the SFH, we first fix $t_{\rm sf}$. Then, $t_{\rm in}$ is determined so that we minimize $\chi^2$. Consequently, we choose three parameter sets: ($t_{\rm sf}$ [Gyr], $t_{\rm in}$ [Gyr]) = (6.0, 23), (11, 12), and (15, 9.0), where we set the Galactic age as 15 Gyr. All of the three cases predict almost identical star formation history. Next, we test the intermittency (or variability) of the star formation rate (SFR) along with the smooth SFH suggested from the infall model. The large value of the $\chi^2$ statistic supports the violent time-variation of the SFH. If we interpret the observed SFH with smooth and variable components, the amplitude of the variable component is comparable to the smooth component. Thus, intermittent SFH of the Galactic disk is strongly suggested. We also examined the metallicity distribution of G-dwarfs. We found that the true parameter set lies between ($t_{\rm sf}$ [Gyr], $t_{\rm in}$ [Gyr]) = (6, 23) and (11, 12), though we should need a more sophisticated model including the process of metal enrichment within the Galactic halo.

We studied the statistical methods for the estimation of the luminosity function (LF) of galaxies. We focused on four nonparametric estimators: $1/V_{\rm max}$ estimator, maximum-likelihood estimator of Efstathiou et al. (1988), Cho{\l}oniewski's estimator, and improved Lynden-Bell's estimator. The performance of the $1/V_{\rm max}$ estimator has been recently questioned, especially for the faint-end estimation of the LF. We improved these estimators for the studies of the distant Universe, and examined their performances for various classes of functional forms by Monte Carlo simulations. We also applied these estimation methods to the mock 2dF redshift survey catalog prepared by Cole et al. (1998). We found that $1/V_{\rm max}$ estimator yields a completely unbiased result if there is no inhomogeneity, but is not robust against clusters or voids. This is consistent with the well-known results, and we did not confirm the bias trend of $1/V_{\rm max}$ estimator claimed by Willmer (1997) in the case of homogeneous sample. We also found that the other three maximum-likelihood type estimators are quite robust and give consistent results with each other. In practice we recommend Cho{\l}oniewski's estimator for two reasons: 1. it simultaneously provides the shape and normalization of the LF; 2. it is the fastest among these four estimators, because of the algorithmic simplicity. Then, we analyzed the photometric redshift data of the Hubble Deep Field prepared by Fern\'{a}ndez-Soto et al. (1999) using the above four methods. We also derived luminosity density $\rho_{\rm L}$ at $B$- and $I$-band. Our $B$-band estimation is roughly consistent with that of Sawicki, Lin, & Yee (1997), but a few times lower at $2.0 < z < 3.0$. The evolution of $\rho_{\rm L}(I)$ is found to be less prominent.

We constructed a model of infrared and sub-mm (hereafter IR) galaxy number count and estimated history of the IR luminosity density. We treat the evolutionary change of galaxy luminosities as a stepwise nonparametric form, in order to explore the most suitable evolutionary history which reproduces the present observational results. We found the evolutionary patterns which satisfy both constraints required from Cosmic Infrared Background (CIRB) and IR galaxy number counts. One order of magnitude increase of luminosity at redshift $z=0.75 - 1.0$ was found in IR $60 \mu$m luminosity density evolution. We also found that a large number of galaxies ($\sim 10^7$ in the whole sky) will be detected in all-sky survey at far-infrared by Infrared Imaging Surveyor (IRIS); Japanese infrared satellite project Astro-F.

To determine the exact shape of the luminosity function (LF) of galaxies is one of the central problems in galactic astronomy and observational cosmology. The most popular method to estimate the LF is maximum likelihood, which is clearly understood with the concepts of the information theory. In the field of information theory and statistical inference, great advance has been made by the discovery of Akaike's Information Criterion (AIC). It enables us to perform a direct comparison among different types of models with different numbers of parameters, and provides us a common basis of the model adequacy. In this paper we applied AIC to the determination of the shape of the LF. We first treated the estimation using stepwise LF, and derived a formula to obtain the optimal bin number. In addition, we studied the method to compare the goodness-of-fit of the parametric form with stepwise LF.

We examine the correlation between physical quantities to explore the existence of dark matter (DM) in the Local Group dwarf spheroidal galaxies (dSphs). In order to clarify whether DM exists in the dSphs we compare two extreme models of their internal kinematics: [1] a tidal model (the mass of a dSph is estimated with the luminous mass, since the observed large velocity dispersions of the dSphs are attributed to the tidal force by the Galaxy), and [2] a DM model (the mass of a dSph is estimated with the virial mass, since the velocity dispersion is considered to reflect the DM potential). In both models, we find that the relation between the surface brightness of the dSphs and the tidal force by the Galaxy is consistently interpreted. This makes a critical remark about the previous studies that concluded that the tidal force is effective for the dSphs based only on model [1]. We also check the correlation between the Galactocentric distance and the tidal force. Consequently, both models are also supported in this test. Thus, we are unable to judge which of the two is more promising for the dSphs in correlation investigations. The physical process of tidal disruption is also considered. Since the timescale for the tidally perturbed stars to escape is much shorter than the orbital periods of the dSphs, it is difficult for the tidally perturbed dSphs to exist as an assembly. Thus, we suggest that the gravitational field of DM is necessary to bind the dSphs with a large velocity dispersion.

The cluster Zwicky 1615.8+3505 is considered to be a dynamically young poor cluster. To investigate the morphology and star-formation activity of galaxies under the environment of a dynamically young poor cluster, we have performed V, R, and I surface photometry, optical low-resolution spectroscopy, and 12CO (J=1-0) line observations for member galaxies. Our data show that more than 90% of the observed galaxies show regular morphologies and no star-formation activities, indicating that the environment does not affect these galaxy properties. Among sixteen galaxies observed, only NGC 6104 shows a significant star-formation activity, and shows a distorted morphology, indicating a tidal interaction. This galaxy contains double knots, and only one knot possesses Seyfert activity, though the sizes and luminosities are similar to each other; we also discuss this feature.

Infrared Imaging Surveyor (IRIS) is a satellite which will be launched at the beginning of 2003. One of the main purposes of the IRIS mission is an all-sky survey in the far-infrared region with a flux limit much deeper than that of IRAS. Detection of a large number of galaxies ($\sim several \times 10^6$ in the whole sky) is expected in this survey. We investigated the expected optical and near-infrared (NIR) number counts of galaxies detected by the far-infrared scanner (FIS) of IRIS (hereafter, IRIS galaxies) and possibility of optical and NIR follow-up of them. The spectral energy distribution and the luminosity function of the IRIS galaxies are modeled based on the properties of galaxies observed by IRAS. The IRIS galaxies are divided into two populations according to their infrared luminosities ($L_{IR}$); normal spirals ($L_{IR}<10^{10}\LO$) and starbursts ($L_{IR}>10^{10}\LO$). The expected number counts of IRIS galaxies for both of the populations are calculated in $B$ and $H$ bands. We show that about 60 normal galaxies and about 80 starburst galaxies are detected per square degree in both of the two bands, when galaxy evolution is not taken into account. All of the normal population of IRIS galaxies are located at the redshift $z\ltsim 0.1$. As for the starburst population, we also calculated number of galaxies with a simple model of evolution. The total number of starburst population predicted by the evolution model is larger by 20% than that expected from the non-evolution model. In the evolution model, the numbers of low-z ($z<1$), intermediate-z ($1<z<3$), and high-z ($z>3$) galaxies are 100, 20, and 0.2 per square degree, respectively.

Infrared Imaging Surveyor (IRIS) is a satellite which will be launched in the beginning of 2003. One of the main purposes of the IRIS mission is an all-sky survey at far-infrared (FIR) with a flux limit much deeper than that of IRAS. In order to examine the performance of the survey, we estimated the FIR galaxy counts in four (50, 70, 120, and 150 $\mu$m) bands based on some models. We adopted a multicomponent model which consists of cirrus and starburst components for galaxy spectra, and the nearby FIR luminosity function derived from that of IRAS galaxies. We derived the number counts, redshift distributions, and infrared diffuse background radiation spectra for i) no evolution, ii) pure luminosity evolution, iii) pure density evolution with $q_0 = 0.1$ and 0.5. We found that a large number of galaxies ($\sim {\rm a few} \times 10^6$ in the whole sky) will be detected in this survey. With the aid of a vast number of detection, we will detect the effect of galaxy evolution, and evaluate the amplitude of evolution at least in the nearby universe in the IRIS survey, though it will be still difficult to constrain which type of evolution takes place from the number count alone. We also studied the estimation of redshifts of detected galaxies by their infrared colors alone. Although significant contamination takes place among nearby faint galaxies and high-$z$ ones, we found that rough estimation of galaxy redshift can be practicable by jointly using present and future optical surveys.

We analyzed 10 dwarf spheroidal galaxies (dSphs) in the Local Group, and found two distinct sequences on the $M_{\rm vir} / L$ -- $M_{\rm vir}$ plane: $M_{\rm vir} / L \propto M_{\rm vir}^{1.6}$ for $M_{\rm vir} < 10^8 M_\odot$ whereas $M_{\rm vir} / L \simeq {\rm const.}$ for $M_{\rm vir} > 10^8 M_\odot$ ($M_{\rm vir}$ and $L$ are the virial mass and the total luminosity of a dSph, respectively). We interpret the disconuity as the threshold for the gas in dSphs to be blown away by successive supernovae. We succeeded to give a quantitative explanation of the discontinuity mass of $10^8M_\odot$ as blow-away condition. We further derived the above relation for the low-mass dSphs, assuming that the initial star formation rate of the dSphs is proportional to the inverse of the cooling time. The relation of high-mass dSphs is also explained along with the same consideration, with the condition that the gas cannot been blown away effectively.

We have performed a systematic study of several regions in the sky where the number of galaxies exhibiting star formation (SF) activity is greater than average. We used Kiso ultraviolet-excess galaxies (KUGs) as our SF-enhanced sample. By statistically comparing the KUG and non-KUG distributions, we discovered four KUG-rich regions with a size of $\sim 10^\circ \times 10^\circ$. One of these regions corresponds spatially to a filament of length $\sim 60 h^{-1}$ Mpc in the Lynx-Ursa Major region ($\alpha \sim 9^{\rm h} - 10^{\rm h}, \delta \sim 42^\circ - 48^\circ$). We call this ``the Lynx-Ursa Major (LUM) filament''. We obtained $V(RI)_{\rm C}$ surface photometry of 11 of the KUGs in the LUM filament and used these to investigate the integrated colors, distribution of SF regions, morphologies, and local environments. We found that these KUGs consist of distorted spiral galaxies and compact galaxies with blue colors. Their star formation occurs in the entire disk, and is not confined to just the central regions. The colors of the SF regions imply that active star formation in the spiral galaxies occurred $10^{7 - 8}$ yr ago, while that of the compact objects occurred $10^{6-7}$ yr ago. Though the photometric characteristics of these KUGs are similar to those of interacting galaxies or mergers, most of these KUGs do not show direct evidence of merger processes.

We discuss a physical interpretation of the relation between $M_{\rm vir}/L$ and $M_{\rm vir}$ of dwarf spheroidal galaxies (dSphs), where $M_{\rm vir}$ and $L$ are the virial mass and the total luminosity of a dSph, respectively. We used 11 dSphs in the Local Group as the sample. We find two distinct sequences on the $M_{\rm vir}/L$--$M_{\rm vir}$ plane: $M_{\rm vir}/L\propto M_{\rm vir}^{2.0}$ for dSphs with $M_{\rm vir}<10^8M_\odot$, whereas $M_{\rm vir}/L\simeq$ constant for $M_{\rm vir}>10^8M_\odot$. A ``discontinuity'' is seen at $M_{\rm vir}\simeq 10^8M_\odot$. We interpret the ``discontinuity'' as the threshold for the gas in dSphs to be blown away by successive supernovae. If a dSph has virial mass (most of which is dark mass) less than $10^8M_\odot$, the gas is blown away, while in a dSph of larger mass, the deep potential well prevents the blow-away mechanism from working effectively. Thus, large mass ratio of dark matter (DM) to baryonic matter (i.e., large $M_{\rm vir}/L$) is realized in a low-mass ($M_{\rm vir}<10^8M_\odot$) dSph through the gas depletion, whereas $M_{\rm vir}/L$ becomes lower in high-mass ($M_{\rm vir}>10^8M_\odot$) dSphs. We further make an attempt to explain the above relation for the low-mass dSphs, $M_{\rm vir}/L\propto M_{\rm vir}^{2.0}$, based on estimate of cooling time, using the scaling laws of virial temperature, virial mass and radius of a dSph and assuming that the heating by OB-star radiation terminates the star formation activity. We succeed in deriving the above relation for the mass-to-light ratio and luminosity.

We present H-alpha velocity fields of thirteen giant H II regions in four nearby dwarf irregular galaxies, NGC 2366, Holmberg II, IC 2574, and WLM. We classify the velocity features as well as the morphologies of the H II regions. The H-alpha velocity features are divided into three categories; three H II regions show chaotic feature with a typical scale of variation of a few 100 pc in size and a few 10 km s^{-1} in velocity, one shows expanding-bubble feature, and the remaining nine have calm velocity fields. There is a correlation between the H-alpha velocity feature and the morphology of H II regions. We measured bulk motion of the H II regions relative to the ambient H I velocity for the H II regions with calm velocity fields and found a typical velocity difference of about 5 km s^{-1}. We discuss a model for origin of star-forming regions based on the presence of the velocity difference between H-alpha and H I gas as well as the H I characteristics.

Sunspot proper motions and flares of a super active region NOAA 5395, which was the biggest and the most flare-active region in the 22nd sunspot cycle, were analyzed in details. We measured sunspot proper motions by using the H-alpha - 5.0 A images obtained with the 60-cm Domeless Solar Telescope (DST) at Hida Observatory, Kyoto University and found some peculiar vortex-like motions of small satellite spots successively emerged from the leading edge of this sunspot group. To explain these motions of small sunspots, we proposed a schematic model of the successive emergence of twisted and winding magnetic flux ropes coiling around a trunk of magnetic flux tube. The location of the strongest flare activity was found to coincide with very the site of the vortex-like motions of sunspots. We conclude that the flare-productive magnetic shear is produced by the emergence of the twisted magnetic flux bundle. Magnetic energy is stored in the twisted flux tubes which are originally formed in the convection zone and released as flares in the course of the emergence of the twisted flux tubes above the photosphere.

We examined the general characteristics of the Kiso Ultra-violet Excess Galaxies (KUGs). We present for the first time the quantitative expressions for the criteria of the KUGs; the boundary color separating the KUGs from the non-KUGs is (B-V)_{T} = 0.74 and the KUG degrees of UV strength are found to correlate with the mean (B-V)_{T} colors. We investigate the nature of the KUGs, a sample of blue galaxy population, and show that (1) about a half of the KUGs are spiral galaxies with Sb to Scd, (2) the KUGs are biased to late-type galaxies and include early-type galaxies with young star populations, and (3) the KUGs are preferably found among less luminous galaxies with log L(B) < 10. The KUGs also contain the post-starburst galaxies, many of which are found among the blue galaxy population at intermediate redshifts. The analysis of the far-infrared data shows that a typical present-to-past star formation rate for a KUG is 0.4.

The fraction of star-forming galaxies in rich clusters of galaxies increases rapidly with the redshift. This is interpreted as the result of a rapid evolution of cluster galaxies, though its mechanism is not yet clear. One hypothesis is that if galaxies run into the dense intracluster medium (ICM) regions, starbursts can be induced due to compression of molecular clouds in the galaxies by a raised "external" ICM pressure. In a merging cluster, there would be dense ICM regions and some galaxies may experience a rapid increase of external pressure. Thus, if the mechanism works, starbursts would occur in such galaxies. Such a situation is considered to be realized in Coma cluster, and galaxies showing a recent starburst have a characteristic spatial distribution; they populate an elongated region between two sub-clusters. We examined the above hypothesis by measuring the spatial distribution of blue galaxies, regarded as star-forming galaxies in Abell~168 to check whether a case like Coma is realized in this recently merging cluster. However, we could not find an enhanced blue-galaxy distribution between two sub-clusters. We discuss the causes of our result, including possibilities of surviving the hypothesis.