Is the Galactic submillimeter dust emissivity underestimated?

We present detailed modeling of the spectral energy distribution (SED) of the spiral galaxies NGC 891, NGC 4013, and NGC 5907 in the far-infrared and submm wavelengths. The model takes into account the emission produced by the diffuse dust and the star forming HII complexes. The dust mass is constrained by radiative transfer simulations in the optical (Xilouris et al. 1999). We find that the submm emission predicted by our model cannot account for the observed fluxes. Two scenarios may account for the "missing" submm flux. In the first scenario (Popescu et al. 2000), additional dust (to that derived from the optical, and associated with young stars) is embedded in the galaxy in the form of a thin disk and gives rise to additional submm emission. The other scenario investigates whether the average submm emissivity of the dust grains is higher than the values widely used in Galactic environments. In this case, the dust mass is equal to that derived from the optical observations, and the submm emissivity is treated as a free parameter calculated by fitting our model to the observed SED. We find the submm emissivity value to be ~3 times that often used for our Galaxy. While both scenarios reproduce the observed 850 micron surface brightness, the extra embedded dust model is not supported by the near infrared observations. We, thus, find that the enhanced dust submm emissivity scenario is the most plausible. [abridged]Comment: 12 pages, 10 figures, accepted for publication in Astronomy and Astrophysic

Modelling the dust content of spiral galaxies: More dust mass vs. enhanced dust grain emissivity

We present detailed modelling of the spectral energy distribution (SED) of the spiral galaxies NGC 891, NGC 4013, and NGC 5907 in the far-infrared (FIR) and sub-millimeter (submm) wavelengths. The model takes into account the emission of the diffuse dust component, which is heated by the UV and optical radiation fields produced by the stars, as well as the emission produced locally in star forming HII complexes. The radiative transfer simulations of Xilouris et al. (1999) in the optical bands are used to constrain the stellar and dust geometrical parameters, as well as the total amount of dust. We find that the submm emission predicted by our model can not account for the observed fluxes at these wavelengths. We examine two cases, one having more dust embedded in a second thin disk and another allowing for an enhanced submillimeter emissivity of the dust grains. We argue that both cases can equally well reproduce the observed SED. The case of having more dust embedded in a second disk though, is not supported by the near-infrared observations and thus more realistic distributions of the dust (i.e., in spiral arms and clumps) have to be examined in order to better fit the surface brightness of each galaxy.Comment: To appear in the proceedings of: "The Spectral Energy Distribution of Gas-Rich Galaxies: Confronting Models with Data", Heidelberg, 4-8 Oct. 2004, eds. C.C. Popescu and R.J. Tuffs, AIP Conf. Ser., in pres

The influence of clumping on surface brightness fits of edge-on spiral galaxies

We have used a Monte Carlo radiative transfer code to produce edge-on images of dusty galactic disks, allowing a fraction of the dust to be distributed in clumps. Synthetic images of edge-on galaxies have been constructed for different amounts of dust, distributions of clumps and fractions of dust in clumps, following the formalism of Bianchi et al. (1999). We have also considered models with stellar emission embedded in the clumps. The synthetic images have been fitted with analytical models made with smooth distributions of dust, adopting the procedure developed by Xilouris et al. (1999) to fit optical images of real edge-on galaxies. We have compared the parameters determined by the fit with the input parameters of the models. For the clumping distributions adopted in this paper, the neglect of clumping results in underestimating the amount of dust in a galaxy. However, the underestimation is never larger than 40%.Comment: 6 pages, 4 pictures, A&A accepte

Extinction Effects in Spiral Galaxy Rotation Curves

Observations show that the slope of the inner part of the H$\alpha$/[NII] rotation curves of disk galaxies is depressed by extinction: at fixed luminosity, the observed slope is in fact seen to depend on the disk inclination to the line of sight. Using a simple extinction model, we are able to reproduce well the observed trends. The model assumes an exponential distribution, both in the radial and $z$ directions, identical for star--forming regions and dust. Fits to the data are optimized by varying the scale height and scale length of absorbers, as well as the disk's central optical depth $\tau_\circ$, as seen face--on. The observations indicate that disk central opacity increases with total luminosity. Model fits for the most luminous galaxies (brighter than $M-5\log h=-21.4$ in the I band) yield $\tau_\circ \simeq 3.5^{+4.0}_{-2.0}$, near $\lambda=0.66$ $\mu$m. The large uncertainty on the inferred value of $\tau_\circ$ is due to the poorly known characteristics of the distribution of absorbers in the disk, as well as to the likelihood of strong departures from an exponential radial distribution near the galaxy centers.Comment: 19 pages, 5 figures, accepted for publication in A

An evolutionary disc model of the edge-on galaxy NGC 5907

We present a physical model that explains the two disparate observational facts: 1) the exponential vertical disc structure in the optical and NIR of the non-obscured part of the stellar disc and 2) the enhanced FIR/submm luminosity by about a factor of four near the obscured mid-plane, which requires additional dust and also stellar light to heat the dust component. We use multi-band photometry in U, B, V, R, and I- band combined with radiative transfer through a dust component to fit simultaneously the vertical surface-brightness and colour index profiles in all bands adopting a reasonable star formation history and dynamical heating function. The final disc model reproduces the surface-brightness profiles in all bands with a moderately declining star formation rate and a slowly starting heating function for young stars. The total dust mass is 57 million solar masses as required from the FIR/submm measurements. Without a recent star burst we find in the midplane an excess of 5.2-, 4.0-, and 3.0-times more stellar light in the U-, B-, and V-band, respectively. The corresponding stellar mass-to-light ratios are 0.91 in V- and 1.0 in R-band. The central face-on optical depth in V-band is 0.81 and the radial scale length of the dust is 40% larger than that of the stellar disc. Evolutionary disc models are a powerful method to understand the vertical structure of edge-on galaxies. Insights to the star formation history and the dynamical evolution of stellar discs can be gained. FIR/submm observations are necessary to restrict the parameter space for the models.Comment: 17 pages, 12 figures (24 files), A&A in pres

Dust in a merging galaxy sequence: the SCUBA view

We investigate the cold and warm dust properties during galaxy interactions using a merging galaxy sample ordered into a chronological sequence from pre- to post-mergers. Our sample comprises a total of 29 merging systems selected to have far-infrared and sub-millimeter observations. We use the 100-to-850 micron flux density ratio, f100/f850, as a proxy to the mass fraction of the warm and the cold dust in these systems. We find evidence for an increase in f100/f850 along the merging sequence from early to advanced mergers and interpret this trend as an increase of the warm relative to the cold dust mass. We argue that the two key parameters affecting the f100/f850 flux ratio is the star-formation rate and the dust content of individual systems relative to the stars. Using a sophisticated model for the absorption and re-emission of the stellar UV radiation by dust we show that these parameters can indeed explain both the increase and the observed scatter in the f100/f850 along the merging galaxy sequence.Comment: To appear in the proceedings of: "The Spectral Energy Distribution of Gas-Rich Galaxies: Confronting Models with Data", Heidelberg, 4-8 Oct. 2004, eds. C.C. Popescu and R.J. Tuffs, AIP Conf. Ser., in pres

Modelling the spectral energy distribution of galaxies. I. Radiation fields and grain heating in the edge-on spiral NGC891

We describe a new tool for the analysis of the UV to the sub-millimeter (sub-mm) spectral energy distribution (SED) of spiral galaxies. We use a consistent treatment of grain heating and emission, solve the radiation transfer problem for a finite disk and bulge, and self-consistently calculate the stochastic heating of grains placed in the resulting radiation field. We use this tool to analyse the well-studied nearby edge- on spiral galaxy NGC 891. First we investigate whether the old stellar population in NGC 891, along with a reasonable assumption about the young stellar population, can account for the heating of the dust and the observed far-infrared and sub- mm emission. The dust distribution is taken from the model of Xilouris et al. (1999), who used only optical and near-infrared observations to determine it. We have found that such a simple model cannot reproduce the SED of NGC 891, especially in the sub-mm range. It underestimates by a factor of 2–4 the observed sub-mm flux. A number of possible explanations exist for the missing sub-mm flux. We investigate a few of them and demonstrate that one can reproduce the observed SED in the far-infrared and the sub-mm quite well, as well as the observed radial profile at 850 µm. For the models calculated we give the relative proportion of the dust radiation powered by the old and young stellar popula- tions as a function of FIR/sub-mm wavelength. In all models we find that the dust is predominantly heated by the young stellar population

The distribution of the ISM in the Milky Way A three-dimensional large-scale model

We use the COBE/DIRBE (1.2, 2.2, 60, 100, 140, and 240 $\mu$m) maps and the COBE/FIRAS spectra (for the wavelength range 100 - 1000 $\mu$m) to constrain a model for the spatial distribution of the dust, the stars, and the gas in the Milky Way. By assuming exponential axisymmetric distributions for the dust and the stars and by performing the corresponding radiative transfer calculations we closely (given the simple geometry of the model) reproduce the FIR and NIR maps of the Milky Way. Similar distributions for the atomic and molecular hydrogen in the disk are used (with an inner cut-off radius for the atomic hydrogen) to fit the gas data. The star formation rate as a function of the Galactic radius is derived from the FIR emission and is well in agreement with existing estimates from various star formation tracers. The gas surface density is plotted against the star formation rate density and an ``intrinsic'' Galactic Schmidt law is derived with excellent agreement with the ``external'' Schmidt law found for spiral galaxies. The Milky Way is found to consume $\sim 1$ and $\sim 10$ of its gas in the outer and inner regions respectively (for a period of 0.1 Gyr) to make stars. The dust-induced B-V color excess observed in various directions and distances (up to $\sim 6.5$ kpc) with well-studied Cepheid stars is compared with the model predictions showing a good agreement. The simple assumption of exponential distributions of stars and dust in the Galaxy is found to be quite instructive and adequate in modeling all the available data sets from 0.45 $\mu$m (B-band) to 1000 $\mu$m.Comment: 14 Pages, 10 figure

Cold and warm dust along a merging galaxy sequence

We investigate the cold and warm dust properties during galaxy interactions using a merging galaxy sample ordered into a chronological sequence from pre- to post-mergers. Our sample comprises a total of 29 merging systems selected to have far-infrared and sub-millimeter observations. The sub-millimeter data are mainly culled from the literature while for 5 galaxies (NGC 3597, NGC 3690, NGC 6090, NGC 6670 and NGC 7252) the sub-millimeter observations are presented here for the first time. We use the 100-to-850 micron flux density ratio, f_{100}/f_{850}, as a proxy to the mass fraction of the warm and the cold dust in these systems. We find evidence for an increase in f_{100}/f_{850} along the merging sequence from early to advanced mergers and interpret this trend as an increase of the warm relative to the cold dust mass. We argue that the two key parameters affecting the f_{100}/f_{850} flux ratio is the star-formation rate and the dust content of individual systems relative to the stars. Using a sophisticated model for the absorption and re-emission of the stellar UV radiation by dust we show that these parameters can indeed explain both the increase and the observed scatter in the f_{100}/f_{850} along the merging galaxy sequence. We also discuss our results under the hypothesis that elliptical galaxies are formed via disc galaxy mergers.Comment: Accepted for publication in MNRA