Mass Surface Density Maps Research Articles

Denoising Diffusion Probabilistic Models to Predict the Density of Molecular Clouds

We introduce the state-of-the-art deep-learning denoising diffusion probabilistic model as a method to infer the volume or number density of giant molecular clouds (GMCs) from projected mass surface density maps. We adopt magnetohydrodynamic simulations with different global magnetic field strengths and large-scale dynamics, i.e., noncolliding and colliding GMCs. We train a diffusion model on both mass surface density maps and their corresponding mass-weighted number density maps from different viewing angles for all the simulations. We compare the diffusion model performance with a more traditional empirical two-component and three-component power-law fitting method and with a more traditional neural network machine-learning approach. We conclude that the diffusion model achieves an order-of-magnitude improvement on the accuracy of predicting number density compared to that by other methods. We apply the diffusion method to some example astronomical column density maps of Taurus and the infrared dark clouds G28.37+0.07 and G35.39-0.33 to produce maps of their mean volume densities.

GMC collisions as triggers of star formation – VIII. The core mass function

ABSTRACTCompression in giant molecular cloud (GMC) collisions is a promising mechanism to trigger the formation of massive star clusters and OB associations. We simulate colliding and non-colliding magnetized GMCs and examine the properties of pre-stellar cores, selected from projected mass surface density maps, including after synthetic ALMA observations. We then examine core properties, including mass, size, density, velocity, velocity dispersion, temperature, and magnetic field strength. After 4 Myr, ∼1000 cores have formed in the GMC collision, and the high-mass end of the core mass function (CMF) can be fit by a power-law dN/dlogM ∝ M−α with α ≃ 0.7, i.e. relatively top heavy compared to a Salpeter mass function. Depending on how cores are identified, a break in the power law can appear around a few $\times 10\, \mathrm{M}_\odot$. The non-colliding GMCs form fewer cores with a CMF with α ≃ 0.8–1.2, i.e. closer to the Salpeter index. We compare the properties of these CMFs to those of several observed samples of cores. Considering other properties, cores formed from colliding clouds are typically warmer, have more disturbed internal kinematics, and are more likely to be gravitational unbound, than cores formed from non-colliding GMCs. The dynamical state of the protocluster of cores formed in the GMC–GMC collision is intrinsically subvirial but can appear to be supervirial if the total mass measurement is affected by observations that miss mass on large scales or at low densities.

NIHAO XIII: Clumpy discs or clumpy light in high-redshift galaxies?

Many massive star forming disc galaxies in the redshift range 3 to 0.5 are observed to have a clumpy morphology showing giant clumps of size $\sim$1 kpc and masses of about $10^7M_{\odot}$ to $10^{10} M_{\odot}$. The nature and fate of these giant clumps is still under debate. In this work we use 19 high-resolution simulations of disc galaxies from the NIHAO sample to study the formation and the evolution of clumps in the discs of high redshift galaxies. We use mock HST - CANDELS observations created with the radiative transfer code GRASIL-3D to carry out, for the first time, a quantitative comparison of the observed fraction of clumpy galaxies and its evolution with redshift with simulations. We find a good agreement between the observed clumpy fraction and the one of the NIHAO galaxies. We find that dust attenuation can suppress intrinsically bright clumps and enhance less luminous ones. In our galaxy sample we only find clumps in light (u-band) from young stars but not in stellar mass surface density maps. This means that the NIHAO sample does not show clumpy stellar discs but rather a clumpy light distribution originating from clumpy star formation events. The clumps found in the NIHAO sample match observed age/color gradients as a function of distance from the galaxy center but they show no sign of inward migration. Clumps in our simulations disperse on timescales of a about a hundred Myr and their contribution to bulge growth is negligible.

High-resolution mapping of dust via extinction in the M31 bulge

We map the dust distribution in the central 180" (~680 pc) region of the M31 bulge, based on HST/WFC3 and ACS observations in ten bands from near-ultraviolet (2700 A) to near-infrared (1.5 micron). This large wavelength coverage gives us great leverage to detect not only dense dusty clumps, but also diffuse dusty molecular gas. We fit a pixel-by-pixel spectral energy distributions to construct a high-dynamic-range extinction map with unparalleled angular resolution (~0.5" , i.e., ~2 pc) and sensitivity (the extinction uncertainty, \delta A_V~0.05). In particular, the data allow to directly fit the fractions of starlight obscured by individual dusty clumps, and hence their radial distances in the bulge. Most of these clumps seem to be located in a thin plane, which is tilted with respect to the M31 disk and appears face-on. We convert the extinction map into a dust mass surface density map and compare it with that derived from the dust emission as observed by Herschel . The dust masses in these two maps are consistent with each other, except in the low-extinction regions, where the mass inferred from the extinction tends to be underestimated. Further, we use simulations to show that our method can be used to measure the masses of dusty clumps in Virgo cluster early-type galaxies to an accuracy within a factor of ~2.

Star formation laws in the Andromeda galaxy: gas, stars, metals and the surface density of star formation

We use hierarchical Bayesian regression analysis to investigate star formation laws in the Andromeda galaxy (M31) in both local (30, 155, and 750pc) and global cases. We study and compare the well-known Kennicutt-Schmidt law, the extended Schmidt law and the metallicity/star formation correlation. Using a combination of H$\alpha$ and 24 $\mu$m emission, a combination of far-ultraviolet and 24$\mu$m, and the total infrared emission, we estimate the total star formation rate (SFR) in M31 to be between $0.35\pm 0.04$M$_{\odot}$yr$^{-1}$ and $0.4\pm 0.04$M$_{\odot}$yr$^{-1}$. We produce a stellar mass surface density map using IRAC 3.6$\mu$m emission and measured the total stellar mass to be $6.9 \times 10^{10}$M$_{\odot}$. For the Kennicutt-Schmidt law in M31, we find the power-law index $N$ to be between 0.49 and 1.18, for all the laws, the power-law index varies more with changing gas tracer than with SFR tracer. The power-law index also changes with distance from the centre of the galaxy. We also applied the commonly-used ordinary least squares fitting method and showed that using different fitting methods leads to different power-law indices. There is a correlation between the surface density of SFR and the stellar mass surface density, which confirms that the Kennicutt-Schmidt law needs to be extended to consider the other physical properties of galaxies. We found a weak correlation between metallicity, the SFR and the stellar mass surface density.

Galaxy gas flows inferred from a detailed, spatially resolved metal budget

We use the most extensive integral field spectroscopic map of a local galaxy, NGC 628, combined with gas and stellar mass surface density maps, to study the distribution of metals in this galaxy out to 3 effective radii ($\rm R_e$). At each galactocentric distance, we compute the metal budget and thus constrain the mass of metals lost. We find that in the disc about 50% of the metals have been lost throughout the lifetime of the galaxy. The fraction of metals lost is higher in the bulge ($\sim$70%) and decreases towards the outer disc ($\rm \sim 3 \ R_e$). In contrast to studies based on the gas kinematics, which are only sensitive to ongoing outflow events, our metal budget analysis enables us to infer the average outflow rate during the galaxy lifetime. By using simple physically motivated models of chemical evolution we can fit the observed metal budget at most radii with an average outflow loading factor of order unity, thus clearly demonstrating the importance of outflows in the evolution of disc galaxies of this mass range ($\rm log(M_\star/M_\odot) \sim 10)$. The observed gas phase metallicity is higher than expected from the metal budget and suggests late-time accretion of enriched gas, likely raining onto the disc from the metal-enriched halo.

Dynamical signatures of a ΛCDM-halo and the distribution of the baryons in M 33

We map the neutral atomic gas content of M33 using high resolution VLA and GBT observations and fit a tilted ring model to determine the orientation of the extended gaseous disk and its rotation curve. The disk of M33 warps from 8 kpc outwards without substantial change of its inclination with respect to the line of sight. Rotational velocities rise steeply with radius in the inner disk, reaching 100 km/s in 4 kpc, then the rotation curve becomes more perturbed and flatter with velocities as high as 120-130 km/s out to 23 kpc. We derive the stellar mass surface density map of M33's optical disk, via pixel -SED fitting methods based on population synthesis models, which highlights variations in the mass-to-light ratio. The stellar mass surface further out is estimated from deep images of outer disk fields. Stellar and gas maps are then used in the dynamical analysis of the rotation curve to constrain the dark matter distribution which is relevant at all radii. A dark matter halo with a Navarro-Frenk-White density profile in a LCDM cosmology, provides the best fit to the rotation curve for a dark halo concentration C=10 and a total halo mass of 4.3 10^{11}Msun. This imples a baryonic fraction of order 0.02 and the evolutionary history of this galaxy should account for loss of a large fraction of its original baryonic content.

Denosumab rapidly increases cortical bone in key locations of the femur: a 3D bone mapping study in women with osteoporosis.

Women with osteoporosis treated for 36 months with twice-yearly injections of denosumab sustained fewer hip fractures compared with placebo. Treatment might improve femoral bone at locations where fractures typically occur. To test this hypothesis, we used 3D cortical bone mapping of postmenopausal women with osteoporosis to investigate the timing and precise location of denosumab versus placebo effects in the hips. We analyzed clinical computed tomography scans from 80 female participants in FREEDOM, a randomized trial, wherein half of the study participants received subcutaneous denosumab 60 mg twice yearly and the others received placebo. Cortical 3D bone thickness maps of both hips were created from scans at baseline, 12, 24, and 36 months. Cortical mass surface density maps were also created for each visit. After registration of each bone to an average femur shape model followed by statistical parametric mapping, we visualized and quantified statistically significant treatment effects. The technique allowed us to pinpoint systematic differences between denosumab and control and to display the results on a 3D average femur model. Denosumab treatment led to an increase in femoral cortical mass surface density and thickness, already evident by the third injection (12 months). Overall, treatment with denosumab increased femoral cortical mass surface density by 5.4% over 3 years. One-third of the increase came from increasing cortical density, and two-thirds from increasing cortical thickness, relative to placebo. After 36 months, cortical mass surface density and thickness had increased by up to 12% at key locations such as the lateral femoral trochanter versus placebo. Most of the femoral cortex displayed a statistically significant relative difference by 36 months. Osteoporotic cortical bone responds rapidly to denosumab therapy, particularly in the hip trochanteric region. This mechanism may be involved in the robust decrease in hip fractures observed in denosumab-treated women at increased risk of fracture.

SPATIALLY RESOLVED STELLAR, DUST, AND GAS PROPERTIES OF THE POST-INTERACTING WHIRLPOOL GALAXY SYSTEM

Using infrared imaging from the Herschel Space Observatory, observed as part of the VNGS, we investigate the spatially resolved dust properties of the interacting Whirlpool galaxy system (NGC 5194 and NGC 5195), on physical scales of 1 kpc. Spectral energy distribution modelling of the new infrared images in combination with archival optical, near- through mid-infrared images confirms that both galaxies underwent a burst of star formation ~370-480 Myr ago and provides spatially resolved maps of the stellar and dust mass surface densities. The resulting average dust-to-stellar mass ratios are comparable to other spiral and spheroidal galaxies studied with Herschel, with NGC 5194 at log M(dust)/M(star)= -2.5+/-0.2 and NGC 5195 at log M(dust)/M(star)= -3.5+/-0.3. The dust-to-stellar mass ratio is constant across NGC 5194 suggesting the stellar and dust components are coupled. In contrast, the mass ratio increases with radius in NGC 5195 with decreasing stellar mass density. Archival mass surface density maps of the neutral and molecular hydrogen gas are also folded into our analysis. The gas-to-dust mass ratio, 94+/-17, is relatively constant across NGC 5194. Somewhat surprisingly, we find the dust in NGC 5195 is heated by a strong interstellar radiation field, over 20 times that of the ISRF in the Milky Way, resulting in relatively high characteristic dust temperatures (~30 K). This post-starburst galaxy contains a substantial amount of low-density molecular gas and displays a gas-to-dust ratio (73+/-35) similar to spiral galaxies. It is unclear why the dust in NGC 5195 is heated to such high temperatures as there is no star formation in the galaxy and its active galactic nucleus is 5-10 times less luminous than the one in NGC 5194, which exhibits only a modest enhancement in the amplitude of its ISRF.

The dust and gas properties of M83

We examine the dust and gas properties of the nearby, barred galaxy M83, which is part of the Very Nearby Galaxy Survey. Using images from the PACS and SPIRE instruments of Herschel, we examine the dust temperature and dust mass surface density distribution. We find that the nuclear, bar and spiral arm regions exhibit higher dust temperatures and masses compared to interarm regions. However, the distribution of dust temperature and mass are not spatially coincident. Assuming a trailing spiral structure, the dust temperature peaks in the spiral arms lie ahead of the dust surface density peaks. The dust mass surface density correlates well with the distribution of molecular gas as traced by CO (J=3-2) images (JCMT) and the star formation rate as traced by H?2 with a correction for obscured star formation using 24 micron emission. Using HI images from THINGS to trace the atomic gas component, we make total gas mass surface density maps and calculate the gas-to-dust ratio. We find a mean gas-to-dust ratio of 84 \pm 4 with higher values in the inner region assuming a constant CO-to-H2 conversion factor. We also examine the gas-to-dust ratio using CO-to-H2 conversion factor that varies with metallicity.

Spiral structure in nearby galaxies - I. Sample, data analysis and overview of results

This paper, the first of two, introduces an observational study of spiral structure in galaxies chosen from the SINGS survey. Near infrared (NIR) and optical data are used to produce mass surface density maps, and from these the morphology of the disc is examined. The aim of this work is to characterise the prevalence of spiral structure in this sample and, in the cases where a clear spiral pattern is found, include the findings in a comparative study (reported in paper II). A two-armed (`grand design') spiral pattern is found in approximately half the galaxies studied, including all those that are designated as grand design in the optical, but also including some, but not all, optically flocculent galaxies. It is found that the level of non-axisymmetric structure in the galaxies' mass distributions is only modestly higher in those galaxies that are classified as `grand design' compared with those that are not, implying that non-grand design galaxies possess significant power in higher order modes. There is no evidence that bars preferentially trigger the spirals, but they do appear to stir up non-axisymmetric structure in the disc. In contrast, there is evidence that strong/close tidal interactions with companion galaxies are associated with strong two-armed spiral structure in the infrared, though there are a number of galaxies with relatively weak infrared spiral structure that do not possess such companions.

Tracing spiral density waves in M81

We use SPITZER IRAC 3.6 and 4.5micron near infrared data from the Spitzer Infrared Nearby Galaxies Survey (SINGS), optical B, V and I and 2MASS Ks band data to produce mass surface density maps of M81. The IRAC 3.6 and 4.5micron data, whilst dominated by emission from old stellar populations, is corrected for small-scale contamination by young stars and PAH emission. The I band data are used to produce a mass surface density map by a B-V colour-correction, following the method of Bell and de Jong. We fit a bulge and exponential disc to each mass map, and subtract these components to reveal the non-axisymmetric mass surface density. From the residual mass maps we are able to extract the amplitude and phase of the density wave, using azimuthal profiles. The response of the gas is observed via dust emission in the 8micron IRAC band, allowing a comparison between the phase of the stellar density wave and gas shock. The relationship between this angular offset and radius suggests that the spiral structure is reasonably long lived and allows the position of corotation to be determined.

Molecular gas in the galaxy M 83

We present the kinematics of the molecular gas in the barred spiral galaxy M 83 (NGC 5236). The study is based on 12CO( 1–0 and 2–1) observations with the Swedish-ESO Submillimetre Telescope (SEST). Iso-velocity maps of the entire optical disk, 10′ 10′ or 13 13 kpc, are produced. They show the pattern of an inclined, rotating disk, but also the effects of streaming motions along the spiral arms. A dynamical mass of about 6 10 is estimated by fitting the rotation curve of an exponential disk model to these data. The gas constitutes about 13% of the disk mass. The pattern speed is determined from the residual velocity pattern. The locations of various resonances are discussed. The molecular gas velocity dispersion is determined, and a trend of decreasing dispersion with increasing galactocentric radius is found. A total gas (H2 + + He) mass surface density map is presented, and compared to the critical density for star formation of an isothermal gaseous disk. The critical density is exceeded in the spiral arms, but not in the interarm regions. The locations of Giant Molecular Associations (GMAs) and regions are consistent with this scenario of dynamically induced star formation.