Higher Dust Masses Research Articles

Gas conditions of a star-formation selected sample in the first billion years

ABSTRACT We present Atacama Large Millimetre/submillimetre Array (ALMA) observations of the [O iii] 88 $\mu {\rm m}$ emission of a sample of thirteen galaxies at $z = 6$ to 7.6 selected as [C ii]-emitting companion sources of quasars. To disentangle the origins of the luminous Oxygen line in the $z\, \gt\, 6$ Universe, we looked at emission-line galaxies that are selected through an excellent star-formation tracer [C ii] with star-formation rates between 9 and 162 $\rm {\rm M}_{\odot }\,yr^{-1}$. Direct observations reveal [O iii] emission in just a single galaxy ($L_\mathrm{[O\, {\small III}]}/L_\mathrm{[C\, {\small II}]}$$\, = 2.3$), and a stacked image shows no [O iii] detection, providing deep upper limits on the $L_\mathrm{[O\, {\small III}]}/L_\mathrm{[C\, {\small II}]}$ ratios in the $z \gt 6$ Universe ($L_\mathrm{[O\, {\small III}]}/L_\mathrm{[C\, {\small II}]}$$\, \lt 1.2$ at $3 \sigma$). While the fidelity of this sample is high, no obvious optical/near-infrared counterpart is seen in the JWST imaging available for four galaxies. Additionally accounting for low-z CO emitters, line stacking shows that our sample-wide result remains robust: The enhanced $L_\mathrm{[O\, {\small III}]}/L_\mathrm{[C\, {\small II}]}$ reported in the first billion years of the Universe is likely due to the selection towards bright, blue Lyman-break galaxies with high surface star-formation rates or young stellar populations. The deep upper limit on the rest-frame 90 μm continuum emission ($\lt 141 \mu$Jy at $3 \sigma$), implies a low average dust temperature ($T_\mathrm{dust} \lesssim 30\,$ K) and high-dust mass ($M_\mathrm{dust} \sim 10^8\, \mathrm{M}_\odot$). As more normal galaxies are explored in the early Universe, synergy between JWST and ALMA is fundamental to further investigate the ISM properties of the a broad range of samples of high-z galaxies.

A dusty streamer infalling onto the disk of a class I protostar

Context. Observations of interstellar material infalling onto star- and planet-forming systems have become increasingly common thanks to recent advancements in radio interferometry. These structures have the potential to alter the dynamics of protoplanetary disks significantly by triggering the formation of substructures, inducing shocks, and modifying their physical and chemical properties. Moreover, the protoplanetary disks are replenished with new material, which increases the overall mass budget for planet formation. Aims. In this study, we combine new ALMA band 3 and archival band 6 observations to characterize the dust content and infall rate of a 4000 au arc-like structure that is infalling onto [MGM2012] 512 (hereafter M512), a class I young stellar object located in the Lynds 1641 region of the Orion A molecular cloud. Methods. We detected the extended dust emission from this structure in both ALMA bands. We tested whether the velocity pattern of the streamer is consistent with infalling trajectories by means of analytical streamline models. We measured spectral index maps for the first time and derived a dust opacity-index profile along a streamer. We constrained its grain properties and mass. Results. We find that the arc structure is consistent with infalling motions. We measure a spectral index α ~ 3.2 across the entire structure and a dust opacity index β ~ 1.6. Considering grain properties consistent with the measured β, the structure can host up to 245 M⊕ of dust, which exceeds or is comparable to the mass of the inner unresolved 600 au, which contains the protoplanetary disk of M512. Assuming a typical dust-to-gas ratio of 1% for the streamer, the free-fall timescales (50 kyr) imply total mass-infall rates up to 1.5 × 10−6 M⊙ yr−1. M512 has been classified as an outbursting source with multi-epoch WISE photometry. It is thus an interesting case study for exploring the possible connection between infalling streamers and accretion outbursts. Conclusions. M512 is a unique source for which dust continuum emission of an arc-like streamer extending out to 4000 au can be characterized in a dual-band analysis. The dust properties are similar to those in the interstellar medium and imply a high dust mass. A massive streamer like this can affect the evolution of the star- and planet-forming inner system strongly.

The –M disk Relationship for Herbig Ae/Be Stars: A Lifetime Problem for Disks with Low Masses?

The accretion of material from protoplanetary disks onto their central stars is a fundamental process in the evolution of these systems and a key diagnostic in constraining the disk lifetime. We analyze the relationship between the stellar accretion rate and the disk mass in 32 intermediate-mass Herbig Ae/Be systems and compare them to their lower-mass counterparts, T Tauri stars. We find that the –M disk relationship for Herbig Ae/Be stars is largely flat at ∼10−7 M ☉ yr−1 over 3 orders of magnitude in dust mass. While most of the sample follows the T Tauri trend, a subset of objects with high accretion rates and low dust masses are identified. These outliers (12 out of 32 sources) have an inferred disk lifetime of less than 0.01 Myr and are dominated by objects with low infrared excess. This outlier sample is likely identified in part by the bias in classifying Herbig Ae/Be stars, which requires evidence of accretion that can only be reliably measured above a rate of ∼10−9 M ☉ yr−1 for these spectral types. If the disk masses are not underestimated and the accretion rates are not overestimated, this implies that these disks may be on the verge of dispersal, which may be due to efficient radial drift of material or outer disk depletion by photoevaporation and/or truncation by companions. This outlier sample likely represents a small subset of the larger young, intermediate-mass stellar population, the majority of which would have already stopped accreting and cleared their disks.

Can global warming bring more dust?

In the late twentieth century, global mean surface air temperature especially on land is continuously warming. Our analyses show that the global mean of dust increased since 1980, using the Modern-Era Retrospective Analysis version 2 for Research and Applications (MERRA-2) reanalysis data. This variation of global dust is mainly contributed by the dust increase outside of dust core areas (i.e. high dust mass concentration region). The causes to result in global dust variations are explored. In dust core areas, surface wind is the primary driving factor for surface dust, both of which show no remarkable trends of increase or decrease since 1980. In areas outside of the core areas, especially in arid and semi-arid areas in North and Middle Asia, surface air temperature warming is the primary impact factor causing the dust increase. An increase in surface air temperature is accompanied by enhancement of atmospheric instability which can trigger more upward motion and bring more dust. All 9 Earth System Models (ESMs) for the Aerosol Chemistry Model Intercomparison Project (AerChemMIP) reproduce the reasonable spatial distribution and seasonal cycle of dust in the present day. But only a few models such as BCC-ESM1 and GFDL-ESM4 simulate the increasing trend of dust similar to MERRA-2. While the primary impact of wind in dust core areas, and surface temperature outside of the core areas, especially in middle to high latitudes in Eurasian continent, are presented in most ESMs.

A Multifluid Dust Module in Athena++: Algorithms and Numerical Tests

We describe the algorithm, implementation, and numerical tests of a multifluid dust module in the Athena++ magnetohydrodynamic code. The module can accommodate an arbitrary number of dust species interacting with the gas via aerodynamic drag (characterized by the stopping time), with a number of numerical solvers. In particular, we describe two second-order accurate, two-stage, fully implicit solvers that are stable in stiff regimes, including short stopping times and high dust mass loading, and they are paired with the second-order explicit van Leer and Runge–Kutta gas dynamics solvers in Athena++, respectively. Moreover, we formulate a consistent treatment of dust concentration diffusion with dust back-reaction, which incorporates momentum diffusion and ensures Galilean invariance. The new formulation and stiff drag solvers are implemented to be compatible with most of the existing features of Athena++, including different coordinate systems, mesh refinement, and shearing box and orbital advection. We present a large suite of test problems, including the streaming instability in linear and nonlinear regimes, as well as local and global settings, which demonstrate that the code achieves the desired performance. This module will be particularly useful for studies of dust dynamics and planet formation in protoplanetary disks.

Dual constraints with ALMA: new [O iii] 88-μm and dust-continuum observations reveal the ISM conditions of luminous LBGs at z ∼ 7

ABSTRACT We present new [${\rm O\, {\small III}}$] 88-$\mu \mathrm{{m}}$ observations of five bright z ∼ 7 Lyman-break galaxies spectroscopically confirmed by ALMA through [${\rm C\, {\small II}}$] 158 $\mu \mathrm{{m}}$, unlike recent [${\rm O\, {\small III}}$] detections where Lyman α was used. This nearly doubles the sample of Epoch of Reionization galaxies with robust (5σ) [${\rm C\, {\small II}}$] and [${\rm O\, {\small III}}$] detections. We perform a multiwavelength comparison with new deep HST images of the rest-frame UV, whose compact morphology aligns well with [${\rm O\, {\small III}}$] tracing ionized gas. In contrast, we find more spatially extended [${\rm C\, {\small II}}$] emission likely produced in neutral gas, as indicated by an [${\rm N\, {\small II}}$] 205-$\mu \mathrm{{m}}$ non-detection in one source. We find a correlation between the optical ${[{\rm O\, {\small III}}]}+ {\mathrm{H\,\beta }}$ equivalent width and [${\rm O\, {\small III}}$]/[${\rm C\, {\small II}}$], as seen in local metal-poor dwarf galaxies. cloudy models of a nebula of typical density harbouring a young stellar population with a high-ionization parameter adequately reproduce the observed lines. Surprisingly, however, our models fail to reproduce the strength of [${\rm O\, {\small III}}$] 88-$\mu \mathrm{{m}}$, unless we assume an α/Fe enhancement and near-solar nebular oxygen abundance. On spatially resolved scales, we find [${\rm O\, {\small III}}$]/[${\rm C\, {\small II}}$] shows a tentative anticorrelation with infrared excess, LIR/LUV, also seen on global scales in the local Universe. Finally, we introduce the far-infrared spectral energy distribution fitting code mercurius to show that dust-continuum measurements of one source appear to favour a low dust temperature and correspondingly high dust mass. This implies a high stellar metallicity yield and may point towards the need of dust production or grain-growth mechanisms beyond supernovae.

The mass and size of Herbig disks as seen by ALMA

Context. Many population studies have been performed over the past decade with the Atacama Large millimeter/submillimeter Array (ALMA) to understand the bulk properties of protoplanetary disks around young stars. The studied populations have mostly consisted of late spectral type (i.e., G, K & M) stars, with relatively few more massive Herbig stars (spectral types B, A & F). With Gaia-updated distances, now is a good time to use ALMA archival data for a Herbig disk population study and take an important step forward in our understanding of planet formation. Aims. The aim of this work is to determine the masses and sizes of all Herbig dust disks observed with ALMA to date in a volume-limited sample out to 450 pc. These masses and sizes are put in the context of the Lupus and Upper Sco T Tauri disk populations. Methods. ALMA Band 6 and Band 7 archival data of 36 Herbig stars are used, making this work 64% complete out to 225 pc, and 38% complete out to 450 pc also including Orion. Using stellar parameters and distances, the dust masses and sizes of the disks are determined via a curve-of-growth method. Survival analysis is used to obtain cumulative distributions of the dust masses and radii. Results. Herbig disks have a higher dust mass than the T Tauri disk populations of Lupus and Upper Sco by factors of ~3 and ~7 respectively. In addition, Herbig disks are often larger than the typical T Tauri disk. Although the masses and sizes of Herbig disks extend over a similar range to those of T Tauri disks, the distributions of masses and sizes of Herbig disks are significantly skewed toward higher values. Lastly, group I disks are more massive than group II disks. An insufficient number of group II disks are observed at sufficient angular resolution to determine whether or not they are also small in size compared to group I disks. Conclusions. Herbig disks are skewed towards more massive and larger dust disks compared to T Tauri disks. Based on this we speculate that these differences find their origin in an initial disk mass that scales with the stellar mass, and that subsequent disk evolution enlarges the observable differences, especially if (sub)millimeter continuum optical depth plays a role. Moreover, the larger disk masses and sizes of Herbig stars could be linked to the increasing prevalence of giant planets with host star mass.

On the role of aerosol radiative effect in the wet season onset timing over the Congo rainforest during boreal autumn

Abstract. The boreal summer dry season length is reported to have been increasing in the last 3 decades over the Congo rainforest, which is the second-largest rainforest in the world. In some years, the wet season in boreal autumn starts early, while in others it arrives late. The mechanism behind such a change in the wet season onset date has not been investigated yet. Using multi-satellite data sets, we discover that the variation in aerosols in the dry season plays a major role in determining the subsequent wet season onset. Dry season aerosol optical depth (AOD) influences the strength of the southern African easterly jet (AEJ-S) and, thus, the onset of the wet season. Higher AOD associated with a higher dust mass flux reduces the net downward shortwave radiation and decreases the surface temperature over the Congo rainforest region, leading to a stronger meridional temperature gradient between the rainforest and the Kalahari Desert as early as in June. The latter, in turn, strengthens the AEJ-S, sets in an early and a stronger easterly flow, and leads to a stronger equatorward convergence and an early onset of the wet season in late August to early September. The mean AOD in the dry season over the region is strongly correlated (r=0.7) with the timing of the subsequent wet season onset. Conversely, in low AOD years, the onset of the wet season over the Congo basin is delayed to mid-October.

A Stellar Mass Dependence of Structured Disks: A Possible Link with Exoplanet Demographics

Abstract Gaps in protoplanetary disks have long been hailed as signposts of planet formation. However, a direct link between exoplanets and disks remains hard to identify. We present a large sample study of ALMA disk surveys of nearby star-forming regions to disentangle this connection. All disks are classified as either structured (transition, ring, extended) or nonstructured (compact) disks. Although low-resolution observations may not identify large-scale substructure, we assume that an extended disk must contain substructure from a dust evolution argument. A comparison across ages reveals that structured disks retain high dust masses up to at least 10 Myr, whereas the dust mass of compact, nonstructured disks decreases over time. This can be understood if the dust mass evolves primarily by radial drift, unless drift is prevented by pressure bumps. We identify a stellar mass dependence of the fraction of structured disks. We propose a scenario linking this dependence with that of giant exoplanet occurrence rates. We show that there are enough exoplanets to account for the observed disk structures if transitional disks are created by exoplanets more massive than Jupiter and ring disks by exoplanets more massive than Neptune, under the assumption that most of those planets eventually migrate inwards. On the other hand, the known anticorrelation between transiting super-Earths and stellar mass implies those planets must form in the disks without observed structure, consistent with formation through pebble accretion in drift-dominated disks. These findings support an evolutionary scenario where the early formation of giant planets determines the disk’s dust evolution and its observational appearance.

Impact of various air mass types on cloud condensation nuclei concentrations along coastal southeast Florida

Coastal southeast Florida experiences a wide range of aerosol conditions, including African dust, biomass burning (BB) aerosols, as well as sea salt and other locally-emitted aerosols. These aerosols are important sources of cloud condensation nuclei (CCN), which play an essential role in governing cloud radiative properties. As marine environments dominate the surface of Earth, CCN characteristics in coastal southeast Florida have broad implications for other regions with the added feature that this site is perturbed by both natural and anthropogenic emissions. This study investigates the influence of different air mass types on CCN concentrations at 0.2% (CCN0.2%) and 1.0% (CCN1.0%) supersaturation (SS) based on ground site measurements during selected months in 2013, 2017, and 2018. Average CCN0.2% and CCN1.0% concentrations were 373 ± 200 cm−3 and 584 ± 323 cm−3, respectively, for four selected days with minimal presence of African dust and BB (i.e., background days). CCN concentrations were not elevated on the four days with highest influence of African dust (289 ± 104 cm−3 [0.2% SS] and 591 ± 302 cm−3 [1.0% SS]), consistent with high dust mass concentrations comprised of coarse particles that are few in number. In contrast, CCN concentrations were substantially enhanced on the five days with the greatest impact from BB (1408 ± 976 cm−3 [0.2% SS] and 3337 ± 1252 cm−3 [1.0% SS]). Ratios of CCN0.2%:CCN1.0% were used to compare the hygroscopicity of the aerosols associated with African dust, BB, and background days. Average ratios were similar for days impacted by African dust and BB (0.54 ± 0.17 and 0.55 ± 0.17, respectively). A 29% higher average ratio was observed on background days (0.71 ± 0.14), owing in part to a strong presence of sea salt and reduced presence of more hydrophobic species such as those of a carbonaceous or mineral-dust nature. Finally, periods of heavy rainfall were shown to effectively decrease both CCN0.2% and CCN1.0% concentrations. However, the rate varied at which such concentrations increased after the rain. This work contributes knowledge on the nucleating ability of African dust and BB in a marine environment after varying periods of atmospheric transport (days to weeks). The results can be used to understand the hygroscopicity of these air mass types, predict how they may influence cloud properties, and provide a valuable model constraint when predicting CCN concentrations in comparable situations.

High Resolution Fine Dust Mass Concentration Calculation Using Two-wavelength Scanning Lidar System

High Resolution Fine Dust Mass Concentration Calculation Using Two-wavelength Scanning Lidar System

Reproducing submillimetre galaxy number counts with cosmological hydrodynamic simulations

ABSTRACT Matching the number counts of high-z submillimetre-selected galaxies (SMGs) has been a long-standing problem for galaxy formation models. In this paper, we use 3D dust radiative transfer to model the submm emission from galaxies in the simba cosmological hydrodynamic simulations, and compare predictions to the latest single-dish observational constraints on the abundance of 850 μm-selected sources. We find good agreement with the shape of the integrated 850 μm luminosity function, and the normalization is within 0.25 dex at >3 mJy, unprecedented for a fully cosmological hydrodynamic simulation, along with good agreement in the redshift distribution of bright SMGs. The agreement is driven primarily by simba’s good match to infrared measures of the star formation rate (SFR) function between z = 2 and 4 at high SFRs. Also important is the self-consistent on-the-fly dust model in simba, which predicts, on average, higher dust masses (by up to a factor of 2.5) compared to using a fixed dust-to-metals ratio of 0.3. We construct a light-cone to investigate the effect of far-field blending, and find that 52 per cent of sources are blends of multiple components, which makes a small contribution to the normalization of the bright end of the number counts. We provide new fits to the 850 μm luminosity as a function of SFR and dust mass. Our results demonstrate that solutions to the discrepancy between submm counts in simulations and observations, such as a top-heavy initial mass function, are unnecessary, and that submillimetre-bright phases are a natural consequence of massive galaxy evolution.

Massive dust clumps in the envelope of the red supergiant VY Canis Majoris

The envelope of the red supergiant VY CMa has long been considered an extreme example of episodic mass loss that is possibly taking place in other cool and massive evolved stars. Recent (sub-)millimeter observations of the envelope revealed the presence of massive dusty clumps within 800 mas of the star and reinforce the picture of drastic mass-loss phenomena in VY CMa. We present new ALMA observations at an angular resolution of 0.1″ and at an unprecedented sensitivity that reveal further details about the dusty clumps. We resolve more discrete features and identify a submillimeter counterpart of a more distant Southwest (SW) Clump known from visual observations. The brightest clump, named C, is marginally resolved in the observations. Gas seen against the resolved continuum emission of clump C produces a molecular spectrum in absorption, in lines of mainly sulfur-bearing species. Except for SW Clump, no molecular emission is found to be associated with the dusty clumps and we propose that the dusty structures have an atypically low gas content. We attempt to reproduce the properties of the dusty clumps through 3D radiative transfer modeling. Although a clump configuration explaining the observations is found, it is not unique. A very high optical depth of all clumps to the stellar radiation makes the modeling very challenging, and requires an unrealistically high dust mass for one of the clumps. We suggest that the dusty features have substructures (e.g., porosity) that allow deeper penetration of stellar photons within the clumps than in a homogeneous configuration. A comparison of the estimated clump ages to variations in the stellar visual flux for over a century suggests that the mechanism responsible for their formation is not uniquely manifested by enhanced or strongly diminished visual light. The study demonstrates that the dusty mass-loss episodes of VY CMa are indeed unparalleled among all known non-explosive stars. The origin of these episodes remains an unsolved problem.

The nature of submillimetre and highly star-forming galaxies in the EAGLE simulation

ABSTRACT We exploit EAGLE, a cosmological hydrodynamical simulation, to reproduce the selection of the observed submillimetre (submm) galaxy population by selecting the model galaxies at z ≥ 1 with mock submm fluxes $S_{850\, \mu \mathrm{m}}$ ≥ 1 mJy. We find a reasonable agreement between the model galaxies within this sample and the properties of the observed submm population, such as their star formation rates (SFRs) at z < 3, redshift distribution, and many integrated galaxy properties. We find that the median redshift of the $S_{850\, \mu \mathrm{m}}$ ≥ 1 mJy model population is z ≈ 2.5, and that they are massive galaxies (M* ∼ 1011 M⊙) with high dust masses (Mdust ∼ 108 M⊙), gas fractions (fgas ≈ 50 per cent), and SFRs ($\dot{M}_* \approx 100$ M⊙ yr−1). In addition, we find that they have major and minor merger fractions similar to the general population, suggesting that mergers are not the sole driver of the high SFRs in the model submm galaxies. Instead, the $S_{850\, \mu \mathrm{m}}$ ≥ 1 mJy model galaxies yield high SFRs primarily because they maintain a significant gas reservoir as a result of hosting an undermassive black hole relative to comparably massive galaxies. Not all ‘highly star-forming’ ($\dot{M}_* \ge 80$ M⊙ yr−1) eagle galaxies have submm fluxes $S_{850\, \mu \mathrm{m}}$ ≥ 1 mJy. We investigate the nature of these highly star-forming ‘Submm-Faint’ galaxies (i.e. $\dot{M}_* \ge 80$ M⊙ yr−1 but $S_{850\, \mu \mathrm{m}}$ < 1 mJy) and find that they are similar to the model submm galaxies, being gas rich and hosting undermassive black holes. However, they are also typically at higher redshifts (z > 4) and are lower mass (M* ∼ 1010 M⊙). These typically higher redshift galaxies show stronger evidence for having been triggered by major mergers, and critically, they are likely missed by most current submm surveys due to their higher dust temperatures and lower dust masses.

Stellar masses and disk properties of Lupus young stellar objects traced by velocity-aligned stacked ALMA 13CO and C18O spectra

Aims. Large samples of protoplanetary disks have been observed in recent ALMA surveys. The gas distributions and velocity structures of most of the disks can still not be imaged at high signal to noise ratios (S/Ns) because of the short integration time per source in these surveys. In this work, we apply the velocity-aligned stacking method to extract more information from molecular-line data of these ALMA surveys and to study the kinematics and disk properties traced by molecular lines. Methods. We re-analyzed the ALMA 13CO (3–2) and C18O (3–2) data of 88 young stellar objects (YSOs) in Lupus with the velocity-aligned stacking method. This method aligns spectra at different positions in a disk based on the projected Keplerian velocities at their positions and then stacks them. This method enhances the S/Ns of molecular-line data and allows us to obtain better detections and to constrain dynamical stellar masses and disk orientations. Results. We obtain 13CO detections in 41 disks and C18O detections in 18 disks with 11 new detections in 13CO and 9 new detections in C18O after applying the method. We estimate the disk orientations and the dynamical masses of the central YSOs from the 13CO data. Our estimated dynamical stellar masses correlate with the spectroscopic stellar masses, and in a subsample of 16 sources, where the inclination angles are better constrained, the two masses are in good agreement within the uncertainties and with a mean difference of 0.15 M⊙. With more detections of fainter disks, our results show that high gas masses derived from the 13CO and C18O lines tend to be associated with high dust masses estimated from the continuum emission. Nevertheless, the scatter is large and is estimated to be 0.9 dex, implying large uncertainties in deriving the disk gas mass from the line fluxes. We find that with such large uncertainties it is expected that there is no correlation between the disk gas mass and the mass accretion rate with the current data. Deeper observations to detect disks with gas masses <10−5 M⊙ in molecular lines are needed to investigate the correlation between the disk gas mass and the mass accretion rate.

GAMA/H-ATLAS: the local dust mass function and cosmic density as a function of galaxy type – a benchmark for models of galaxy evolution

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

Homogeneous Analysis of the Dust Morphology of Transition Disks Observed with ALMA: Investigating Dust Trapping and the Origin of the Cavities

Abstract We analyze the dust morphology of 29 transition disks (TDs) observed with Atacama Large (sub-)Millimeter Array (ALMA) at (sub-)millimeter emission. We perform the analysis in the visibility plane to characterize the total flux, cavity size, and shape of the ring-like structure. First, we found that the M dust–M ⋆ relation is much flatter for TDs than the observed trends from samples of class II sources in different star-forming regions. This relation demonstrates that cavities open in high (dust) mass disks, independent of the stellar mass. The flatness of this relation contradicts the idea that TDs are a more evolved set of disks. Two potential reasons (not mutually exclusive) may explain this flat relation: the emission is optically thick or/and millimeter-sized particles are trapped in a pressure bump. Second, we discuss our results of the cavity size and ring width in the context of different physical processes for cavity formation. Photoevaporation is an unlikely leading mechanism for the origin of the cavity of any of the targets in the sample. Embedded giant planets or dead zones remain as potential explanations. Although both models predict correlations between the cavity size and the ring shape for different stellar and disk properties, we demonstrate that with the current resolution of the observations, it is difficult to obtain these correlations. Future observations with higher angular resolution observations of TDs with ALMA will help discern between different potential origins of cavities in TDs.

Copious Amounts of Dust and Gas in a z = 7.5 Quasar Host Galaxy

Abstract We present IRAM/NOEMA and JVLA observations of the quasar J1342+0928 at z = 7.54 and report detections of copious amounts of dust and [C ii] emission in the interstellar medium (ISM) of its host galaxy. At this redshift, the age of the universe is 690 Myr, about 10% younger than the redshift of the previous quasar record holder. Yet, the ISM of this new quasar host galaxy is significantly enriched by metals, as evidenced by the detection of the [C ii] 158 μm cooling line and the underlying far-infrared (FIR) dust continuum emission. To the first order, the FIR properties of this quasar host are similar to those found at a slightly lower redshift ( ), making this source by far the FIR-brightest galaxy known at . The [C ii] emission is spatially unresolved, with an upper limit on the diameter of 7 kpc. Together with the measured FWHM of the [C ii] line, this yields a dynamical mass of the host of . Using standard assumptions about the dust temperature and emissivity, the NOEMA measurements give a dust mass of . The brightness of the [C ii] luminosity, together with the high dust mass, imply active ongoing star formation in the quasar host. Using [C ii]–SFR scaling relations, we derive star formation rates of 85–545 yr−1 in the host, consistent with the values derived from the dust continuum. Indeed, an episode of such past high star formation is needed to explain the presence of ∼108 M ☉ of dust implied by the observations.

Evidence of sub-surface energy storage in comet 67P from the outburst of 2016 July 03

On 3 July 2016, several instruments on board ESA's Rosetta spacecraft detected signs of an outburst event on comet 67P, at a heliocentric distance of 3.32 AU from the sun, outbound from perihelion. We here report on the inferred properties of the ejected dust and the surface change at the site of the outburst. The activity coincided with the local sunrise and continued over a time interval of 14 - 68 minutes. It left a 10m-sized icy patch on the surface. The ejected material comprised refractory grains of several hundred microns in size, and sub-micron-sized water ice grains. The high dust mass production rate is incompatible with the free sublimation of crystalline water ice under solar illumination as the only acceleration process. Additional energy stored near the surface must have increased the gas density. We suggest a pressurized sub-surface gas reservoir, or the crystallization of amorphous water ice as possible causes.

A Controlled Study of Cold Dust Content in Galaxies from z = 0–2

Abstract At , the formation of new stars is dominated by dusty galaxies whose far-IR emission indicates they contain colder dust than local galaxies of a similar luminosity. We explore the reasons for the evolving IR emission of similar galaxies over cosmic time using (1) local galaxies from GOALS ( ), (2) galaxies at from 5MUSES ( ), and (3) IR luminous galaxies spanning from GOODS and Spitzer xFLS ( ). All samples have Spitzer mid-IR spectra, and Herschel and ground-based submillimeter imaging covering the full IR spectral energy distribution, allowing us to robustly measure , , and for every galaxy. Despite similar infrared luminosities, dusty star-forming galaxies (DSFG) have a factor of 5 higher dust masses and 5 K colder temperatures. The increase in dust mass is linked to an increase in the gas fractions with redshift, and we do not observe a similar increase in stellar mass or star formation efficiency. , a proxy for , is strongly correlated with independently of redshift. We measure merger classification and galaxy size for a subsample, and there is no obvious correlation between these parameters and or . In DSFG, the change in can fully account for the observed colder dust temperatures, suggesting that any change in the spatial extent of the interstellar medium is a second-order effect.