Source author record

Tsutomu T. Takeuchi

Tsutomu T. Takeuchi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

40works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

40 published item(s)

preprint2022arXiv

A new galaxy spectral energy distribution model consistent with the evolution of dust

The spectral energy distribution (SED) of galaxies provides fundamental information on the related physical processes. However, the SED is significantly affected by dust in its interstellar medium. Dust is mainly produced by asymptotic giant branch stars and Type II supernovae. In addition, the dust mass increases through the metal accretion, and the grain size changes by the collisions between the grains. The contribution of each process and the extinction depend on the size distribution. Therefore, the SED model should treat the evolution of the dust mass and size distribution. In spite of the importance of dust evolution, many previous SED models have not considered the evolution of the total mass and size distribution in a physically consistent manner. In this work, we constructed a new radiative transfer SED model, based on our dust evolution model consistent with the chemical evolution. To reduce the computational cost, we adopted the mega-grain and the one-dimensional plane parallel galaxy approximation. As a fiducial case, we calculated Milky Way-like galaxy SEDs at various ages under the closed-box model. We found that a galaxy at the age of 100~Myr does not produce small grains such as polycyclic aromatic hydrocarbons. After 1~Gyr, we observed a drastic increase of infrared emission and attenuation caused by a rapid increase of dust mass. This phenomenon can be treated appropriately for the first time by our new model. This model can be used for the SED fitting to a galaxy at any stage of evolution.

preprint2022arXiv

Dissecting Nearby Galaxies with piXedfit: II. Spatially Resolved Scaling Relations Among Stars, Dust, and Gas

We study spatially resolved scaling relations among stars, dust, and gas in ten nearby spiral galaxies. In a preceding paper Abdurro'uf et al. (2022), we have derived spatially resolved properties of the stellar population and dust by panchromatic spectral energy distribution (SED) fitting using piXedfit. Now, we investigate resolved star formation ($Σ_{\rm H_{2}}$--$Σ_{\rm SFR}$--$Σ_{*}$) and dust scaling relations. While the relations with all sub-galactic regions of the galaxies are reasonably tight ($σ\lesssim 0.3$ dex), we find that most of the scaling relations exhibit galaxy-to-galaxy variations in normalization and shape. Only two relations of $Σ_{\rm dust}$--$Σ_{\rm gas}$ and $Σ_{\rm dust}$--$Σ_{\rm H_{2}}$ do not show noticeable galaxy-to-galaxy variations among our sample galaxies. We further investigate correlations among the scaling relations. We find significant correlations among the normalization of the $Σ_{\rm H_{2}}$--$Σ_{\rm SFR}$--$Σ_{*}$ relations, which suggest that galaxies with higher levels of resolved $\text{H}_{2}$ fraction ($f_{\rm H_{2}}$) tend to have higher levels of resolved star formation efficiency (SFE) and specific star formation rate (sSFR). We also observe that galaxies with higher levels of resolved dust-to-stellar mass ratios tend to have higher levels of resolved sSFR, SFE, and $f_{\rm H_{2}}$. Moreover, we find that galaxies with higher global sSFR and less compact morphology tend to have higher levels of the resolved sSFR, SFE, and $f_{\rm H_{2}}$, which can explain the variations in the normalization of the $Σ_{\rm H_{2}}$--$Σ_{\rm SFR}$--$Σ_{*}$ relationships. Overall, we observe indications of the contributions of both global and local factors in governing the star formation process in galaxies.

preprint2022arXiv

Estimation of the Star Formation Rate of Galaxies with Radio Continuum Obtained with Murchison Widefield Array

We investigate the correlation between the integrated low-frequency and infrared (IR) emissions of star-forming galaxies extracted from the {\sl Herschel} Reference Survey. By taking advantage of the GaLactic Extragalactic All-sky MWA (GLEAM) survey operated by the Murchison Widefield Array (MWA) we examine how this correlation varies at a function of frequency across the 20 GLEAM narrow bands at $72\mbox{--}231\; [\mbox{MHz}]$. These examinations are important for ensuring the reliability of the radio luminosity as a SFR indicator. In this study, we focus on 18 star-forming galaxies whose radio emission is detected by the GLEAM survey. These galaxies show that a single power-law is sufficient to characterise the far-infrared-to-radio correlation across the GLEAM frequency bands and up to $1.5\; [\mbox{GHz}]$. Thus, the radio continuum in this wavelength range can serve as a good dust extinction-free SFR estimator. This is particularly important for future investigation of the cosmic SFR independently from other estimators, since the radio continuum can be detected from $z=0$ to high redshifts ($z \sim 5\mbox{--}10$) in a coherent manner.

preprint2022arXiv

The bimodality in the mass-metallicity relation in SDSS-MaNGA galaxy pairs

Aims: Interacting galaxies show a metallicity dilution compared to isolated galaxies of similar masses in the mass-metallicity space at the global scale. We investigate the spatially resolved mass-metallicity relation (MZR) of galaxy pairs in the SDSS-MaNGA survey to confirm that the local relation between the stellar mass surface density and the metallicity is consistent with the MZR at the global scale. Methods: We investigate the relationship between the stellar mass surface density and the metallicity abundance, 12+log (O/H), for star-forming spaxels belonging to 298 galaxy pairs identified using visual and kinematic indicators in the SDSS-MaNGA survey. We also investigate if a)the location of a spaxel relative to the galaxy centre and b)the galaxy pair separation have any effect on the local MZR. Results: We find that the correlation between mass and metallicity holds for interacting galaxies at the local level. However, we find two peaks in spaxel distribution, one peak with enriched metallicity and the other with diluted metallicity. We find that the spaxels belonging to the galaxy central regions (i.e. at lower R/Reff) are concentrated close to the two peaks. We also find that the metallicity-diluted spaxels belong to galaxy pairs with closer projected separations and that spaxels with enriched metallicity belong to galaxy pairs with greater projected separations. Conclusions: We find two discrete peaks in the spatially resolved MZR for star-forming spaxels that belong to galaxy pairs. The peaks are likely related to the galaxy projected separation or the stage of the interaction process of a galaxy pair.

preprint2022arXiv

Wavelets and sparsity for Faraday tomography

Faraday tomography through broadband polarimetry can provide crucial information on magnetized astronomical objects, such as quasars, galaxies, or galaxy clusters. However, the limited wavelength coverage of the instruments requires that we solve an ill-posed inverse problem when we want to obtain the Faraday dispersion function (FDF), a tomographic distribution of the magnetoionic media along the line of sight. This paper explores the use of wavelet transforms and the sparsity of the transformed FDFs in the form of wavelet shrinkage (WS) for finding better solutions to the inverse problem. We recently proposed the Constraining and Restoring iterative Algorithm for Faraday Tomography (CRAFT; Cooray et al. 2021), a new flexible algorithm that showed significant improvements over the popular methods such as Rotation Measure Synthesis. In this work, we introduce CRAFT+WS, a new version of CRAFT incorporating the ideas of wavelets and sparsity. CRAFT+WS exhibit significant improvements over the original CRAFT when tested for a complex FDF of realistic Galactic model. Reconstructions of FDFs demonstrate super-resolution in Faraday depth, uncovering previously unseen Faraday complexities in observations. The proposed approach will be necessary for effective cosmic magnetism studies using the Square Kilometre Array and its precursors.

preprint2021arXiv

Dissecting Nearby Galaxies with piXedfit: I. Spatially Resolved Properties of Stars, Dust, and Gas as Revealed by Panchromatic SED Fitting

We study spatially resolved properties (on spatial scales of $\sim 1-2$ kpc out to at least $3$ effective radii) of the stars, dust, and gas in ten nearby spiral galaxies. The properties of the stellar population and dust are derived by fitting the spatially resolved spectral energy distribution (SED) with more than 20 photometric bands ranging from far-ultraviolet to far-infrared. Our newly developed software piXedfit performs point spread function matching of images, pixel binning, and models the stellar light, dust attenuation, dust emission, and emission from a dusty torus heated by an active galactic nucleus simultaneously through the energy balance approach. With this self-consistent analysis, we present the spatially resolved version of the IRX--$β$ relation, finding that it is consistent with the relationship from the integrated photometry. We show that the old stellar populations contribute to the dust heating, which causes an overestimation of star formation rate (SFR) derived from the total ultraviolet and infrared luminosities on kpc scales. With archival high-resolution maps of atomic and molecular gas, we study the radial variation of the properties of the stellar populations (including stellar mass, age, metallicity, and SFR), dust (including dust mass, dust temperature, and abundance of polycyclic aromatic hydrocarbon), gas, as well as dust-to-stellar mass and dust-to-gas mass ratios. We observe a depletion of molecular gas mass fraction in the central region of the majority of the galaxies, suggesting that the lack of available fuel is an important factor in suppressing the specific SFR at the center.

preprint2020arXiv

A Method for Unmasking Incomplete Astronomical Signals: Application to CO Multi-line Imaging of Nearby Galaxies Project

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.

preprint2020arXiv

A Study on the Baryon Acoustic Oscillation with Topological Data Analysis

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

preprint2020arXiv

CO Multi-line Imaging of Nearby Galaxies (COMING). IX. $^{12}$CO($J$=2-1)/$^{12}$CO($J$=1-0) line ratio on kiloparsec scales

While molecular gas mass is usually derived from $^{12}$CO($J$=1-0) - the most fundamental line to explore molecular gas - it is often derived from $^{12}$CO($J$=2-1) assuming a constant $^{12}$CO($J$=2-1)/$^{12}$CO($J$=1-0) line ratio ($R_{2/1}$). We present variations of $R_{2/1}$ and effects of the assumption that $R_{2/1}$ is a constant in 24 nearby galaxies using $^{12}$CO data obtained with the Nobeyama 45-m radio telescope and IRAM 30-m telescope. The median of $R_{2/1}$ for all galaxies is 0.61, and the weighted mean of $R_{2/1}$ by $^{12}$CO($J$=1-0) integrated-intensity is 0.66 with a standard deviation of 0.19. The radial variation of $R_{2/1}$ shows that it is high (~0.8) in the inner ~1 kpc while its median in disks is nearly constant at 0.60 when all galaxies are compiled. In the case that the constant $R_{2/1}$ of 0.7 is adopted, we found that the total molecular gas mass derived from $^{12}$CO($J$=2-1) is underestimated/overestimated by ~20%, and at most by 35%. The scatter of a molecular gas surface density within each galaxy becomes larger by ~30%, and at most by 120%. Indices of the spatially resolved Kennicutt-Schmidt relation by $^{12}$CO($J$=2-1) are underestimated by 10-20%, at most 39% in 17 out of 24 galaxies. $R_{2/1}$ has good positive correlations with star-formation rate and infrared color, and a negative correlation with molecular gas depletion time. There is a clear tendency of increasing $R_{2/1}$ with increasing kinetic temperature ($T_{\rm kin}$). Further, we found that not only $T_{\rm kin}$ but also pressure of molecular gas is important to understand variations of $R_{2/1}$. Special considerations should be made when discussing molecular gas mass and molecular gas properties inferred from $^{12}$CO($J$=2-1) instead of $^{12}$CO($J$=1-0).

preprint2020arXiv

CO Multi-line Imaging of Nearby Galaxies (COMING). X. Physical conditions of molecular gas and the local SFR-Mass relation

We investigate the molecular gas properties of galaxies across the main sequence of star-forming (SF) galaxies in the local Universe using $^{12}$CO($J=1-0$) (hereafter $^{12}$CO) and $^{13}$CO($J=1-0$) ($^{13}$CO) mapping data of 147 nearby galaxies obtained in the COMING project, a legacy project of the Nobeyama Radio Observatory. In order to improve signal-to-noise ratios of both lines, we stack all the pixels where $^{12}$CO emission is detected after aligning the line center expected from the first-moment map of $^{12}$CO. As a result, $^{13}$CO emission is successfully detected in 80 galaxies with a signal-to-noise ratio larger than three. The error-weighted mean of integrated-intensity ratio of $^{12}$CO to $^{13}$CO lines ($R_{1213}$) of the 80 galaxies is $10.9$ with a standard deviation of $7.0$. We find that (1) $R_{1213}$ positively correlates to specific star-formation rate (sSFR) with a correlation coefficient of $0.46$, and (2) both flux ratio of IRAS 60~$μ$m to 100~$μ$m ($f_{60}/f_{100}$) and inclination-corrected linewidth of $^{12}$CO stacked spectra ($σ_{{\rm ^{12}CO},i}$) also correlate with sSFR for galaxies with the $R_{1213}$ measurement. Our results support the scenario where $R_{1213}$ variation is mainly caused by the changes in molecular-gas properties such as temperature and turbulence. The consequent variation of CO-to-H$_2$ conversion factor across the SF main sequence is not large enough to completely extinguish the known correlations between sSFR and $M_{\rm mol}/M_{\rm star}$ ($μ_{\rm mol}$) or star-formation efficiency (SFE) reported in previous studies, while this variation would strengthen (weaken) the sSFR-SFE (sSFR-$μ_{\rm mol}$) correlation.

preprint2020arXiv

Constructing a multivariate distribution function with a vine copula: toward multivariate luminosity and mass functions

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.

preprint2015arXiv

Can we use Weak Lensing to Measure Total Mass Profiles of Galaxies on 20 kiloparsec Scales?

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.

preprint2015arXiv

Evolution of mid-infrared galaxy luminosity functions from the entire AKARI NEP-Deep field with new CFHT photometry

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$μ$m) exists nowhere else, due to filter gaps of other space telescopes. It allows us to estimate restframe 8$μ$m and 12$μ$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$μ$m, 12$μ$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 ($Ω_{\mathrm{IR}}$), we found that the $Ω_{\mathrm{IR}}$ evolves as $\propto (1+z)^{4.2\pm 0.4}$.

preprint2014arXiv

Evolution of grain size distribution in high-redshift dusty quasars: Integrating large amounts of dust and unusual extinction curves

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.

preprint2014arXiv

Evolution of the extinction curves in galaxies

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\;μ$m). Then, grain growth becomes effective at small grain radii, forming a bump at $a \sim 10^{-3}${--}$10^{-2}\;μ$m on the grain size distribution. Consequently, the extinction curve at ultraviolet (UV) wavelengths becomes steep, and a bump at $1/λ\sim 4.5\;μ{\rm m}^{-1}\;(λ: \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.

preprint2014arXiv

Galaxy luminosity function and its cosmological evolution: Testing a new feedback model depending on galaxy-scale dust opacity

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.

preprint2013arXiv

CO Component Estimation Based on the Independent Component Analysis

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.

preprint2013arXiv

What determines the grain size distribution in galaxies?

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

preprint2012arXiv

Application of a Self-Organizing State Space Model to the Leonid Meteor Storm in 2001

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.

preprint2012arXiv

Dust formation history of galaxies: a critical role of metallicity for the dust mass growth by accreting materials in the interstellar medium

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.

preprint2012arXiv

Far-Ultraviolet and Far-Infrared Bivariate Luminosity Function of Galaxies: Complex Relation between Stellar and Dust Emission

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

preprint2012arXiv

Infrared Spectral Energy Distribution of Galaxies in the AKARI All Sky Survey: Correlations with Galaxy Properties, and Their Physical Origin

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.

preprint2012arXiv

Reexamination of the Infrared Excess-Ultraviolet Slope Relation of Local Galaxies

The relation between the ratio of infrared (IR) and ultraviolet (UV) flux densities (the infrared excess: IRX) and the slope of the UV spectrum (β) 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-βrelation depending on sample selection. We reexamined the IRX-β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 λ= 42--122\mum, using the data from AKARI which has wider wavelength coverage toward longer wavelengths, we obtained an appropriate IRX-βrelation with total dust emission (TIR): \log(L_{\rm TIR}/L_{\rm FUV}) = \log [10^{0.4(3.06+1.58β)}-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.

preprint2012arXiv

Test for anisotropy in the mean of the CMB temperature fluctuation in spherical harmonic space

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.

preprint2012arXiv

The North Ecliptic Pole Wide survey of AKARI: a near- and mid-infrared source catalog

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 μ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 μm, about 20,000 sources in the shorter parts of the mid-IR bands between 7 - 11 μm, and about 16,000 sources in the longer parts of the mid-IR bands, with \sim 4,000 sources at 24 μm. The estimated 5? detection limits are approximately 21 magnitude (mag) in the 2 - 5 μm bands, 19.5 - 19 mag in the 7 - 11 μm, and 18.8 - 18.5 mag in the 15 - 24 μ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 μm, 18.6 mag at 9 μm, and 18 mag at 18 μ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.

preprint2011arXiv

Infrared Luminosity Functions of AKARI-SDSS Galaxies

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

preprint2010arXiv

Constructing a bivariate distribution function with given marginals and correlation: application to the galaxy luminosity function

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

preprint2010arXiv

Copula Cosmology: Constructing a Likelihood Function

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.

preprint2010arXiv

Luminosity Functions of Local Infrared Galaxies with AKARI: Implications to the Cosmic Star Formation History and AGN Evolution

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.

preprint2010arXiv

Precise Estimation of Cosmological Parameters Using a More Accurate Likelihood Function

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.

preprint2010arXiv

Star-galaxy separation by far-infrared color-color diagrams for the AKARI FIS All-Sky Survey (Bright Source Catalogue Version beta-1)

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.

preprint2008arXiv

Deep 15um AKARI observations in the CDFS: estimating dust luminosities for a MIR-selected sample and for Lyman Break Galaxies and the evolution of L(dust)/L(UV) with the redshift

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.

preprint2007arXiv

IR and UV Galaxies at z=0.6 -- Evolution of Dust Attenuation and Stellar Mass as Revealed by SWIRE and GALEX

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 ($Ω= 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μ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.

preprint2006arXiv

Deep Extragalactic Surveys around the Ecliptic Poles with AKARI (ASTRO-F)

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.

preprint2006arXiv

UV and FIR selected star-forming galaxies at z=0: differences and overlaps

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.

preprint2002arXiv

A New Empirical Method of Estimation of the Far-Infrared Flux of Galaxies

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.

preprint2001arXiv

Absorption Measurement of Hydrogen Molecules in the Early Universe

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.

preprint2001arXiv

The Radio-to-Submillimeter Flux Density Ratio of Galaxies as a Measure of Redshift

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~$μ$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.

preprint2000arXiv

Exploring Galaxy Evolution from Infrared Number Counts and Cosmic Infrared Background

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.

preprint1998arXiv

Photometric Properties of Kiso Ultraviolet-Excess Galaxies in the Lynx-Ursa Major Region

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 ($α\sim 9^{\rm h} - 10^{\rm h}, δ\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.