ABSTRACT
Galaxy kinematics and gravitational lensing are two complementary ways to constrain the distribution of dark matter on galaxy scales. The typical dark matter density profiles adopted in ...dynamical studies cannot easily be adopted in lensing studies. Ideally, a mass model should be used that has the global characteristics of realistic dark matter distributions, and that allows for an analytical calculation of the magnifications and deflection angles. A simple model with these properties, the broken power-law (BPL) model, has very recently been introduced. We examine the dynamical structure of the family of BPL models. We derive simple closed expressions for basic dynamical properties, and study the distribution function under the assumption of velocity isotropy. We find that none of the BPL models with realistic parameters has an isotropic distribution function that is positive over the entire phase space, implying that the BPL models cannot be supported by an isotropic velocity distribution, or models with a more radially anisotropic orbital structure. This result limits the attractiveness of the BPL family as a tool for lensing studies to some degree. More generally, we find that not all members of the general family of double power-law or Zhao models, often used to model dark matter haloes, can be supported by an isotropic or radially anisotropic distribution function. In other words, the distribution function may become negative even for spherically symmetric models with a well-behaved density profile.
Computer simulations of photon transport through an absorbing and/or scattering medium form an important research tool in astrophysics. Nearly all software codes performing such simulations for ...three-dimensional geometries employ the Monte Carlo (MC) radiative transfer (RT) method, including various forms of biasing to accelerate the calculations. Because of the probabilistic nature of the MC technique, the outputs are inherently noisy, but it is often assumed that the average values provide the physically correct result. We show that this assumption is not always justified. Specifically, we study the intensity of radiation penetrating an infinite, uniform slab of material that absorbs and scatters the radiation with equal probability. The basic MCRT method, without any biasing mechanisms, starts to break down for transverse optical depths τ 20 because so few of the simulated photon packets reach the other side of the slab. When including biasing techniques such as absorption/scattering splitting and path length stretching, the simulated photon packets do reach the other side of the slab but the biased weights do not necessarily add up to the correct solution. While the noise levels seem to be acceptable, the average values sometimes severely underestimate the correct solution. Detecting these anomalies requires the judicious application of statistical tests, similar to those used in the field of nuclear particle transport, possibly in combination with convergence tests employing consecutively larger numbers of photon packets. In any case, for transverse optical depths τ 75 the MC methods used in our study fail to solve the one-dimensional slab problem, implying the need for approximations such as a modified random walk.
We revisit the differential energy distribution of steady-state dynamical models. It has been shown that the differential energy distribution of steady-state spherical models does not vary strongly ...with the anisotropy profile, and that it is hence mainly determined by the density distribution of the model. We explore this similarity in more detail. Through a worked example and a simple proof, we show that the mean binding energy per unit mass ⟨ℰ⟩, or equivalently the total integrated binding energy
B
tot
=
M
⟨ℰ⟩, is independent of the orbital structure, not only for spherical models but for any steady-state dynamical model. Only the higher-order moments of the differential energy distribution depend on the details of the orbital structure. We show that the standard deviation of the differential energy distribution of spherical dynamical models varies systematically with the anisotropy profile: radially anisotropic models tend to prefer more average binding energies, whereas models with a more tangential orbital distribution slightly favour more extreme binding energies. Finally, we find that the total integrated binding energy supplements the well-known trio consisting of total kinetic energy, total potential energy, and total energy on an equal footing. Knowledge of any one out of these four energies suffices to calculate the other three.
Interstellar dust in galaxies can be traced either through its extinction effects on the star light or through its thermal emission at infrared wavelengths. Recent radiative transfer studies of ...several nearby edge-on galaxies have found an apparent inconsistency in the dust energy balance: the radiative transfer models that successfully explain the optical extinction underestimate the observed fluxes by an average factor of 3. We investigate the dust energy balance for IC 4225 and NGC 5166, two edge-on spiral galaxies observed by the Herschel Space Observatory in the frame of the H-ATLAS survey. We start from models which were constrained from optical data and extend them to construct the entire spectral energy distribution of our galaxies. These predicted values are subsequently compared to the observed far-infrared fluxes. We find that including a young stellar population in the modelling is necessary as it plays a non-negligible part in the heating of the dust grains. While the modelling approach for both galaxies is nearly identical, we find two very different results. As is often seen in other edge-on spiral galaxies, the far-infrared emission of our radiative transfer model of IC 4225 underestimates the observed fluxes by a factor of about 3. For NGC 5166 on the other hand, we find that both the predicted spectral energy distribution as well as the simulated images match the observations particularly well. We explore possible reasons for this difference and conclude that it is unlikely that one single mechanism is the cause of the dust energy balance problem in spiral galaxies. We discuss the different approaches that can be considered in order to get a conclusive answer on the origin this discrepancy.
We present Herschel Spectral and Photometric Imaging Receiver (SPIRE) Fourier Transform Spectrometer (FTS) observations of the Antennae (NGC 4038/39), a well-studied, nearby (22Mpc), ongoing merger ...between two gas-rich spiral galaxies. We detect five CO transitions (J = 4-3 to J = 8-7), both CI transitions, and the NII 205 mu m transition across the entire system, which we supplement with ground-based observations of the CO J = 1-0, J = 2-1, and J = 3-2 transitions and Herschel Photodetecting Array Camera and Spectrometer (PACS) observations of CII and OI 63 mu m. We find that there are two components to the molecular gas: a cold (T sub(kin) ~ 10-30 K) and a warm (T sub(kin) > 100 K) component. Finally, we find that a combination of turbulent heating, due to the ongoing merger, and supernova and stellar winds are sufficient to heat the molecular gas.
We present the submillimeter spectra from 450 to 1550 GHz of 11 nearby active galaxies observed with the SPIRE Fourier Transform Spectrometer (SPIRE/FTS) on board Herschel. We detect CO transitions ...from J sub(up) = 4 to 12, as well as the two CI fine structure lines at 492 and 809 GHz and the NII1461 GHz line. We used radiative transfer models to analyze the observed CO spectral line energy distributions. The FTS CO data were complemented with ground-based observations of the low-J CO lines. We found that the warm molecular gas traced by the mid-J CO transitions has similar physical conditions (n sub(H2) ~ 10 super(3.2)-10 super(3.9) cm super(-3) and T sub(kin) ~ 300-800 K) in most of our galaxies. Furthermore, we found that this warm gas is likely producing the mid-IR rotational H sub(2) emission. We could not determine the specific heating mechanism of the warm gas, however, it is possibly related to the star formation activity in these galaxies. Our modeling of the CI emission suggests that it is produced in cold (T sub(kin) < 30 K) and dense (n sub(H2) > 10 super(3) cm super(-3)) molecular gas. Transitions of other molecules are often detected in our SPIRE/FTS spectra. The HF J = 1-0 transition at 1232 GHz is detected in absorption in UGC 05101 and in emission in NGC 7130. In the latter, near-infrared pumping, chemical pumping, or collisional excitation with electrons are plausible excitation mechanisms likely related to the active galactic nucleus of this galaxy. In some galaxies, few H sub(2)O emission lines are present. Additionally, three OH+ lines at 909, 971, and 1033 GHz are identified in NGC 7130.
Context.
The ultraviolet (UV) to sub-millimetre spectral energy distribution of galaxies can be roughly divided into two sections: the stellar emission (attenuated by dust) at UV to near-infrared ...wavelengths and dust emission at longer wavelengths. In Dobbels et al. (2020, A&A, 634, A57), we show that these two sections are strongly related, and we can predict the global dust properties from the integrated UV to mid-infrared emission with the help of machine learning techniques.
Aims.
We investigate if these machine learning techniques can also be extended to resolved scales. Our aim is to predict resolved maps of the specific dust luminosity, specific dust mass, and dust temperature starting from a set of surface brightness images from UV to mid-infrared wavelengths.
Methods.
We used a selection of nearby galaxies retrieved from the DustPedia sample, in addition to M31 and M33. These were convolved and resampled to a range of pixel sizes, ranging from 150 pc to 3 kpc. We trained a random forest model which considers each pixel individually.
Results.
We find that the predictions work well on resolved scales, with the dust mass and temperature having a similar root mean square error as on global scales (0.32 dex and 3.15 K on 18″ scales respectively), and the dust luminosity being noticeably better (0.11 dex). We find no significant dependence on the pixel scale. Predictions on individual galaxies can be biased, and we find that about two-thirds of the scatter can be attributed to scatter between galaxies (rather than within galaxies).
Conclusions.
A machine learning approach can be used to create dust maps, with its resolution being only limited to the input bands, thus achieving a higher resolution than
Herschel
. These dust maps can be used to improve global estimates of dust properties, they can lead to a better estimate of dust attenuation, and they can be used as a constraint on cosmological simulations that trace dust.
We expand our previous analytical and numerical studies of the family of Sérsic models, which are routinely used to describe early-type galaxies and the bulges of spiral galaxies. In particular, we ...focus on the total energy budget, an important dynamical property that has not been discussed in detail in previous works. We use two different methods to calculate the total energy for the Sérsic model family that result in two independent expressions that can be used along the entire sequence of Sérsic models. We use these expressions to investigate whether the Spitzer concentration index is a reliable measure for the intrinsic 3D concentration of galaxies, and we conclude that it is not a very useful measure for the central concentration. The popular Third Galaxy Concentration index, on the other hand, is shown to be a reliable measure for the intrinsic 3D concentration, even though it is based on the surface brightness distribution and not on the intrinsic 3D density.
The Sérsic or
R
1/
m
model has become the de facto standard model to describe the surface brightness profiles of early-type galaxies and the bulges of spiral galaxies. The photometric, intrinsic, and ...dynamical properties of this model have been investigated, but mainly for fairly large Sérsic indices
m
. For small values of
m
, appropriate for low-mass and dwarf ellipticals, a detailed investigation of these properties is still lacking. In this study, we used a combination of numerical and analytical techniques to investigate the Sérsic model over the entire range of Sérsic parameters, focussing on the small
m
regime, where a number of interesting and surprising properties are found. For all values
m
< 1, the model is characterised by a finite central luminosity density, and for
m
< 1/2, even a central depression in the luminosity density profile. This behaviour translates to the dynamical properties: we show that all Sérsic models with
m
⩾ 1/2 can be supported by an isotropic velocity dispersion tensor, and that these isotropic models are stable to both radial and non-radial perturbations. The models with
m
< 1/2, on the other hand, cannot be supported by an isotropic velocity dispersion tensor.
Dust attenuation is a crucial but highly uncertain parameter that hampers the determination of intrinsic galaxy properties, such as stellar masses, star formation rates and star formation histories. ...The shape of the dust attenuation law is not expected to be uniform between galaxies, nor within a galaxy. Our DustKING project was introduced at the first BINA workshop in 2016 and aims to study the variations of dust attenuation curves in nearby galaxies. At the second BINA workshop in 2018, I presented the results of our pilot study for the spiral galaxy NGC 628. We find that the average attenuation law of this galaxy is characterised by a MW-like bump and a steep UV slope. Furthermore, we observe intriguing variations within the galaxy, with regions of high AV exhibiting a shallower attenuation curve. Finally, we discuss how our work might benefit from data taken with the UVIT from the Indian AstroSat mission.