Dust Opacity Research Articles

Peculiar Dust Emission within the Orion Molecular Cloud

It is widely assumed that dust opacities in molecular clouds follow a power-law profile with an index, β. Recent studies of the Orion Molecular Cloud (OMC) 2/3 complex, however, show a flattening in the spectral energy distribution (SED) at λ > 2 mm, implying nonconstant indices on scales ≳0.08 pc. The origin of this flattening is not yet known, but it may be due to the intrinsic properties of the dust grains or contamination from other sources of emission. We investigate the SED slopes in OMC 2/3 further using observations of six protostellar cores with Northern Extended Millimeter Array (NOEMA) from 2.9–3.6 mm and Atacama Large Millimeter/submillimeter Array (ALMA) Atacama Compact Array in Band 4 (1.9–2.1 mm) and Band 5 (1.6–1.8 mm) on core and envelope scales of ∼0.02–0.08 pc. We confirm flattened opacity indices between 2.9 mm and 3.6 mm for the six cores with β ≈ −0.16 to 1.45, which are notably lower than the β-values of >1.3 measured for these sources on 0.08 pc scales from single-dish data. Four sources have consistent SED slopes between the ALMA data and the NOEMA data. We propose that these sources may have a significant fraction of emission coming from large dust grains in embedded disks, which biases the emission more at longer wavelengths. Two sources, however, had inconsistent slopes between the ALMA and NOEMA data, indicating different origins of emission. These results highlight how care is needed when combining multiscale observations or extrapolating single-band observations to other wavelengths.

Dusty substructures induced by planets in ALMA disks: How dust growth and dynamics changes the picture

Abstract Protoplanetary disks exhibit a rich variety of substructure in millimeter continuum emission, often attributed to unseen planets. As these planets carve gaps in the gas, dust particles can accumulate in the resulting pressure bumps, forming bright features in the dust continuum. We investigate the role of dust dynamics in the gap-opening process with 2D radiation hydrodynamics simulations of planet–disk interaction and a two-population dust component modeled as a pressureless fluid. We consider the opacity feedback and backreaction due to drag forces as mm grains accumulate in pressure bumps at different stages of dust growth. We find that dust dynamics can significantly affect the resulting substructure driven by the quasi-thermal-mass planet with Mp/M⋆ = 10−4. Opacity feedback causes nonaxisymmetric features to become more compact in azimuth, whereas the drag-induced backreaction tends to dissolve nonaxisymmetries. For our fiducial model, this results in multiple concentric rings of dust rather than the expected vortices and corotating dust clumps found in models without dust feedback. A higher coagulation fraction disproportionately enhances the effect of dust opacity feedback, favoring the formation of crescents rather than rings. Our results suggest that turbulent diffusion is not always necessary to explain the rarity of observed nonaxisymmetric features, and that incorporating dust dynamics is vital for interpreting the observed substructure in protoplanetary disks. We also describe and test the implementation of the publicly-available dust fluid module in the PLUTO code.

Uncovering the truth about M101, NGC 3938, and their significant others through radiative transfer

ABSTRACT Solving the inverse problem in spiral galaxies, that allows the derivation of the spatial distribution of dust, gas, and stars, together with their associated physical properties, directly from panchromatic imaging observations, is one of the main goals of this work. To this end, we used radiative transfer models to decode the spatial and spectral distributions of the nearby face-on galaxies M101 and NGC 3938. In both cases, we provide excellent fits to the surface-brightness distributions derived from GALEX, SDSS, 2MASS, Spitzer, and Herschel imaging observations. Together with previous results from M33, NGC 628, M51, and the Milky Way, we obtain a small statistical sample of modelled nearby galaxies that we analyse in this work. We find that in all cases Milky Way-type dust with Draine-like optical properties provide consistent and successful solutions. We do not find any ‘submm excess’, and no need for modified dust-grain properties. Intrinsic fundamental quantities like star-formation rates (SFR), specific SFR (sSFR), dust opacities, and attenuations are derived as a function of position in the galaxy and overall trends are discussed. In the SFR surface density versus stellar mass surface density space, we find a structurally resolved relation (SRR) for the morphological components of our galaxies, that is steeper than the main sequence (MS). Exception to this is for NGC 628, where the SRR is parallel to the MS.

Dust mass in protoplanetary disks with porous dust opacities

Atacama Large Millimeter/submillimeter Array surveys have suggested that protoplanetary disks are not massive enough to form the known exoplanet population, based on the assumption that the millimeter continuum emission is optically thin. In this work, we investigate how the mass determination is influenced when the porosity of dust grains is considered in radiative transfer models. The results show that disks with porous dust opacities yield similar dust temperatures, but systematically lower millimeter fluxes, as compared to disks that incorporate compact dust grains. Moreover, we have recalibrated the relation between dust temperature and stellar luminosity for a wide range of stellar parameters. We also calculated the dust masses of a large sample of disks using the traditionally analytic approach. The median dust mass from our calculation is about six times higher than the literature result, and this is mostly driven by the different opacities of porous and compact grains. A comparison of the cumulative distribution function between disk dust masses and exoplanet masses shows that the median exoplanet mass is about two times lower than the median dust mass when grains are assumed to be porous and there are no exoplanetary systems with masses higher than the most massive disks. Our analysis suggests that adopting porous dust opacities may alleviate the mass budget problem for planet formation. As an example illustrating the combined effects of optical depth and porous dust opacities on the mass estimation, we conducted new IRAM/NIKA-2 observations toward the IRAS 04370+2559 disk and performed a detailed radiative transfer modeling of the spectral energy distribution (SED). The best-fit dust mass is roughly 100 times higher than the value given by a traditionally analytic calculation. Future spatially resolved observations at various wavelengths are required to better constrain the dust mass.

A Parameter Study of 1D Atmospheric Models of Pulsating AGB Stars

Using the atmospheric pulsation code written by George Bowen, we have performed a parameter study examining the effects of modifying various parameters of models of oxygen-rich AGB atmospheres pulsating in the fundamental and first-overtone modes. For each pulsation mode, we have examined the effects of adjusting the dust condensation temperature, dust condensation temperature range, pulsation amplitude, dust opacity, and metallicity. Our model grids are generated with the constraint that their luminosities are chosen to span the range of observed mass loss rates at a chosen mass. The dust condensation temperature, pulsation amplitude, and dust opacity have strong effects on the ultimate location and shape of the final model grids in the mass luminosity plane. The mass loss rate evolution of the fundamental and first-overtone mode models show a significant difference in behavior. While the fundamental mode models exhibit the typically assumed power–law relation with mass and luminosity, the first-overtone mode models show significant non-power law behavior at observed mass loss rates. Effectively, models in the first-overtone mode require somewhat higher luminosities to reach the same mass loss rates seen in fundamental mode models of the same mass, consistent with observed AGB stars.

A public grid of radiative transfer simulations for Lyα and metal lines in idealised galactic outflows

The vast majority of star-forming galaxies are surrounded by large reservoirs of gas ejected from the interstellar medium. Ultraviolet absorption and emission lines represent powerful diagnostics to constrain the cool phase of these outflows, through resonant transitions of hydrogen and metal ions. The interpretation of these observations is often remarkably difficult as it requires detailed modelling of the propagation of the continuum and emission lines in the gas. To this aim, we present a large public grid of ≈20 000 simulated spectra that includes H I Lyα and five metal transitions associated with Mg II, C II, Si II, and Fe II which is accessible online. The spectra have been computed with the RASCAS Monte Carlo radiative transfer code for 5760 idealised spherically symmetric configurations surrounding a central point source emission, and characterised by their column density, Doppler parameter, dust opacity, wind velocity, as well as various density and velocity gradients. Designed to predict and interpret Lyα and metal line profiles, our grid exhibits a wide diversity of resonant absorption and emission features, as well as fluorescent lines. We illustrate how it can help better constrain the wind properties by performing a joint modelling of observed Lyα, C II, and Si II spectra. Using CLOUDY simulations and virial scaling relations, we also show that Lyα is expected to be a faithful tracer of the gas at T ≈ 104 − 105 K, even if the medium is highly ionised. While C II is found to probe the same range of temperatures as Lyα, other metal lines merely trace cooler phases (T ≈ 104 K). As their gas opacity strongly depends on gas temperature, incident radiation field, metallicity and dust depletion, we caution that optically thin metal lines do not necessarily originate from low H I column densities and may not accurately probe Lyman continuum leakage.

Thermodynamics of Giant Molecular Clouds: The Effects of Dust Grain Size

The dust grain size distribution (GSD) likely varies significantly across star-forming environments in the Universe, but its impact on star formation remains unclear. This ambiguity arises because the GSD interacts nonlinearly with processes like heating, cooling, radiation, and chemistry, which have competing effects and varying environmental dependencies. Processes such as grain coagulation, expected to be efficient in dense star-forming regions, reduce the abundance of small grains and increase that of larger grains. Motivated by this, we investigate the effects of similar GSD variations on the thermochemistry and evolution of giant molecular clouds (GMCs) using magnetohydrodynamic simulations spanning a range of cloud masses and grain sizes, which explicitly incorporate the dynamics of dust grains within the full-physics framework of the STARFORGE project. We find that grain size variations significantly alter GMC thermochemistry: the leading-order effect is that larger grains, under fixed dust mass, GSD dynamic range, and dust-to-gas ratio, result in lower dust opacities. This reduced opacity permits interstellar radiation field and internal radiation photons to penetrate more deeply. This leads to rapid gas heating and inhibited star formation. Star formation efficiency is highly sensitive to grain size, with an order-of-magnitude reduction when grain size dynamic range increases from 10−3–0.1 μm to 0.1–10 μm. Additionally, warmer gas suppresses low-mass star formation, and decreased opacities result in a greater proportion of gas in diffuse ionized structures.

Impact of ice growth on the physical and chemical properties of dense cloud cores

We investigated the effect of time-dependent ice growth on dust grains on the opacity and hence on the dust temperature in a collapsing molecular cloud core, with the aim of quantifying the effect of the dust temperature variations on ice abundances as well as the evolution of the collapse. To perform the simulations, we employed a one-dimensional collapse model that self-consistently and time-dependently combines hydrodynamics with chemical and radiative transfer simulations. The dust opacity was updated on the fly based on the ice growth as a function of the location in the core. The results of the fully dynamical model were compared against simulations run with different values of fixed ice thickness. We found that the ice thickness increases quickly and reaches a saturation value (as a result of a balance between adsorption and desorption) of approximately 90 monolayers in the central core (volume density ~104 cm−3), and several tens of monolayers at a volume density of ~103 cm−3, after only a few 105 yr of evolution. The results thus exclude the adoption of thin (approximately ten monolayer) ices in molecular cloud simulations except at very short timescales. However, the differences in abundances and the dust temperature between the fully dynamic simulation and those with a fixed dust opacity are small; abundances change between the solutions generally within a factor of two only. The assumptions on the dust opacity do have an effect on the collapse dynamics through the influence of the photoelectric effect on the gas temperature, and the simulations take a different time to reach a common central density. This effect is, however, small as well. In conclusion, carrying out chemical simulations using a dust temperature corresponding to a fixed opacity seems to be a good approximation. Still, although at least in the present case its effect on the overall results is limited – as long as the grains are monodisperse – ice growth should be considered to obtain the most accurate representation of the collapse dynamics. We have found in a previous work that considering a grain size distribution leads to a complicated ice composition that depends on the grain size nonlinearly. With this in mind, we will carry out a follow-up study where the influence of the grain size on the present simulation setup is investigated.

The Physical Origin of the Stellar Initial Mass Function

Stars are among the most fundamental structures of our Universe. They comprise most of the baryonic and luminous mass of galaxies; synthesize heavy elements; and inject mass, momentum, and energy into the interstellar medium. They are also home to the planets. Because stellar properties are primarily decided by their mass, the so-called stellar initial mass function (IMF) is critical to the structuring of our Universe. We review the various physical processes and theories that have been put forward as well as the numerical simulations that have been carried out to explain the origin of the stellar IMF. Key messages from this review include the following: ▪Gravity and turbulence most likely determine the power-law, high-mass part of the IMF.▪Depending of the Mach number and the density distribution, several regimes are possible, including ΓIMF ≃ 0, −0.8, −1, or −1.3, where dN/d log M ∝ M ΓIMF . These regimes are likely universal; however, the transition between these regimes is not.▪Protostellar jets can play a regulating influence on the IMF by injecting momentum into collapsing clumps and unbinding gas.▪The peak of the IMF may be a consequence of dust opacity and molecular hydrogen physics at the origin of the first hydrostatic core. This depends weakly on large-scale environmental conditions such as radiation, magnetic field, turbulence, or metallicity. This likely constitutes one reason for the relative universality of the IMF.

The eventful life of GS-z14-0, the most distant galaxy at redshift z=14.32

We developed a model for the star formation history (SFH) of super-early galaxies and applied it to GS-z14-0, the most distant galaxy known, located at $z=14.32$ (294 million years after the Big Bang). The SFH, starting at $z=26.7$, is complex. Initially ($z>18$), the galaxy experiences feedback-regulated phases that are bursty, relatively faint (reaching $M_ UV =-18.4$), and unattenuated. When dust shielding allows for a smooth star formation rate (SFR), the galaxy quickly becomes heavily obscured. During this obscured phase, which lasts for approximately 20<!PCT!> of the total star-forming time, 70<!PCT!> of the observed stars are formed. Super-early galaxies in this phase should be detectable by ALMA. Twenty-six million years before observation, as the galaxy becomes super-Eddington, a powerful radiation-driven outflow clears most of the dust and significantly reduces the SFR by a factor of seven, from $100 The galaxy transitions into a "blue monster" dominating the bright end of the UV luminosity function. When the outflow ceases due to decreased dust opacity, the galaxy relaxes into a post-starburst phase, in which it is currently observed. Our model accurately reproduces all the observed and inferred properties of the galaxy. The analysis of this extreme system opens exciting opportunities for studying the beginnings of the luminous Universe.

Relationships Between HCl, H2O, Aerosols, and Temperature in the Martian Atmosphere: 2. Quantitative Correlations

AbstractThe detection of hydrogen chloride (HCl) in the atmosphere of Mars was among the primary objectives of the ExoMars Trace Gas Orbiter (TGO) mission. Its discovery using the Atmospheric Chemistry Suite mid‐infrared channel (ACS MIR) showed a distinct seasonality and possible link to dust activity. This paper is part 2 of a study investigating the link between HCl and aerosols by comparing gas measurements made with TGO to dust and water ice opacities measured with the Mars Climate Sounder (MCS). In part 1, we showed, and compared, the seasonal evolution of vertical profiles of HCl, water vapor, temperature, dust opacity, and water ice opacity over the dusty periods around perihelion (solar longitudes 180°–360°) across Mars years 34–36. In part 2, we investigated the quantitative correlations in the vertical distribution between each quantity, as well as ozone. We show that there is a strong positive correlation between HCl and water vapor, which is expected due to fast photochemical production rates for HCl when reacting with water vapor photolysis products. We also show a strong positive correlation between water vapor and temperature, but are unable to show any correlation between temperature and HCl. There are weak correlations between the opacities of dust and water ice, and dust and water vapor, but only very low correlations between dust and HCl. We close with a discussion of possible sources and sinks and that interactions between HCl and water ice are the most likely for both, given the inter‐comparison.

Relationships Between HCl, H2O, Aerosols, and Temperature in the Martian Atmosphere: 1. Climatological Outlook

AbstractDetecting trace gases such as hydrogen chloride (HCl) in Mars' atmosphere is among the primary objectives of the ExoMars Trace Gas Orbiter (TGO) mission. Terrestrially, HCl is closely associated with active volcanic activity, so its detection on Mars was expected to point to some form of active magmatism/outgassing. However, after its discovery using the mid‐infrared channel of the TGO Atmospheric Chemistry Suite (ACS MIR), a clear seasonality was observed, beginning with a sudden increase in HCl abundance from below detection limits to 1–3 ppbv in both hemispheres coincident with the start of dust activity, followed by very sudden and rapid loss at the southern autumnal equinox. In this study, we have investigated the relationship between HCl and atmospheric dust by making comparisons in the vertical distribution of gases measured with ACS and aerosols measured co‐located with the Mars Climate Sounder (MCS). This study includes HCl, water vapor, and ozone measured using ACS MIR, water vapor and temperature measured with the near infrared channel of ACS, and temperature, dust opacity, and water ice opacity measured with MCS. In part 1, we show that dust loading has a strong impact in temperature, which controls the abundance of water ice and water vapor, and that HCl is very closely linked to water activity. In part 2, we investigate the quantitative correlations between each quantity and discuss the possible source and sinks of HCl, their likelihood given the correlations, and any issues arising from them.

On the Doublet Flux Ratio of Mg ii Resonance Lines in and Around Galaxies

Observations of metallic doublet emission lines, particularly Mg ii λ λ2796, 2803, provide crucial information for understanding galaxies and their circumgalactic medium. This study explores the effects of resonant scattering on the Mg ii doublet lines and the stellar continuum in spherical and cylindrical geometries. Our findings show that under certain circumstances, resonance scattering can cause an increase in the doublet flux ratio and the escaping flux of the lines beyond what is expected in optically thin spherical media. As expected, the doublet ratio is consistently lower than the intrinsic ratio when the scattering medium is spherically symmetric and dusty. However, if the scattering medium has a disk shape, such as face-on disk galaxies, and is viewed face-on, the doublet ratio is predicted to be higher than 2. It is also shown that doublet ratios as low as those observed in compact star-forming galaxies cannot be explained solely by pure dust attenuation of intrinsic Mg ii emission lines in spherical models unless dust opacity deviates markedly from that expected based on the dust-to-Mg+ gas ratio of our Galaxy. The importance of the continuum-pumped emission lines and expanding media is discussed to understand observational aspects, including doublet flux ratios, which can be lower than 1.5 or higher than 2, as well as symmetric or asymmetric line profiles. It is also discussed that the diffuse warm neutral medium may be an important source of Mg ii emission. These results provide insight into the complexity of the shape and orientation of distant, spatially unresolved galaxies.

SMA 200–400 GHz Survey for Dust Properties in the Icy Class II Disks in the Taurus Molecular Cloud

We present a new Submillimeter Array survey of 47 Class II sources in the Taurus–Auriga region. Our observations made 12 independent samples of flux densities over the 200–400 GHz frequency range. We tightly constrained the spectral indices of most sources to a narrow range of 2.0 ± 0.2; only a handful of spatially resolved (e.g., diameter >250 au) disks present larger spectral indices. The simplest interpretation for this result is that the (sub)millimeter luminosities of all of the observed target sources are dominated by very optically thick (e.g., τ ≳ 5) dust thermal emission. Some previous works that were based on the optically thin assumption thus might have underestimated optical depths by at least 1 order of magnitude. Assuming DSHARP dust opacities, this corresponds to underestimates of dust masses by a similar factor. For our specific selected sample, the lower limits of dust masses implied by the optically thick interpretation are 1–3 times higher than those previous estimates that were made based on the optically thin assumption. Moreover, some population synthesis models show that, to explain the observed, narrowly distributed spectral indices, the disks in our selected sample need to have very similar dust temperatures (T dust). Given a specific assumption of median T dust, the maximum grain sizes ( amax ) can also be constrained, which is a few times smaller than 0.1 mm for T dust ∼ 100 K and a few millimeters for T dust ∼ 24 K. The results may indicate that dust grain growth outside the water snow line is limited by the bouncing/fragmentation barriers. This is consistent with the recent laboratory experiments, which indicated that the coagulation of water-ice-coated dust is not efficient, and the water-ice-free dust is stickier and thus can coagulate more efficiently. In the Class II disks, the dust mass budget outside of the water snow line may be largely retained instead of being mostly consumed by planet formation. While Class II disks still possess sufficient dust masses to feed planet formation at a later time, it is unknown whether or not dust coagulation and planet formation can be efficient or natural outside of the water snow line.

Identification of Hot Gas around Low-mass Protostars

The low carbon content of Earth and primitive meteorites compared to the Sun and interstellar grains suggests that carbon-rich grains were destroyed in the inner few astronomical units of the young solar system. A promising mechanism to selectively destroy carbonaceous grains is thermal sublimation within the soot line at ≳300 K. To address whether such hot conditions are common among low-mass protostars, we observe CH3CN transitions at 1, 2, and 3 mm with the NOrthern Extended Millimeter Array toward seven low-mass and one intermediate-mass protostar (L bol ∼ 2–300L ⊙), as CH3CN is an excellent temperature tracer. We find >300 K gas toward all sources, indicating that hot gas may be prevalent. Moreover, the excitation temperature for CH3OH obtained with the same observations is always lower (∼135–250 K), suggesting that CH3CN and CH3OH have a different spatial distribution. A comparison of the column densities at 1 and 3 mm shows a stronger increase at 3 mm for CH3CN than for CH3OH. Since the dust opacity is lower at longer wavelengths, this indicates that CH3CN is enhanced in the hot gas compared to CH3OH. If this CH3CN enhancement is the result of carbon-grain sublimation, these results suggest that Earth’s initial formation conditions may not be rare.

CAI formation in the early Solar System

Calcium-aluminium-rich inclusions (CAIs) are the oldest dated solid materials in the Solar System, and are found as light-coloured crystalline ingredients in carbonaceous chondrite meteorites. Their formation time is commonly associated with age zero of the Solar System. Nevertheless, the physical and chemical processes that once led to the formation of these submillimetre- to centimetre-sized mineral particles in the early solar nebula are still a matter of debate. In this paper, we propose a pathway to form such inclusions during the earliest phases of disc evolution. We combine 1D viscous disc evolutionary models with 2D radiative transfer, equilibrium condensation, and new dust opacity calculations. We show that the viscous heating associated with the high accretion rates in the earliest evolutionary phases causes the midplane inside of about 0.5 au to heat up to limiting temperatures of about 1500–1700 K, but no further. These high temperatures force all refractory material components of the inherited interstellar dust grains to sublimate – except for a few Al-Ca-Ti oxides, such as Al2O3, Ca2Al2SiO7, and CaTiO3. This is a recurring and very stable result in all our simulations, because these minerals form a natural thermostat. Once the Mg-Fe silicates are gone, the dust becomes more transparent and the heat is more efficiently transported to the disc surface, which prevents further warming. This thermostat mechanism keeps these minerals above their annealing temperature for hundreds of thousands of years, allowing them to form large pure crystalline particles. These particles are dragged out by the viscously spreading disc, and once they reach a distance of about 0.5 au, the silicates recondense on the surface of the Ca-Al-rich particles, adding an amorphous silicate matrix. We estimate that this mechanism of CAI production works during the first 50 000 yr of disc evolution. These particles then continue to move outward and populate the entire disc up to radii of about 50 au, before the accretion rate eventually subsides, the disc cools, and the particles start to drift inwards.

Bayesian analysis of the molecular emission and dust continuum of protoplanetary disks

Context. The MIRI instrument on board the James Webb Space Telescope probes the chemistry and dust mineralogy of the inner regions of protoplanetary disks. The observed spectra are unprecedented in their detail and reveal a rich chemistry with strong diversity between objects. This complicates interpretations that are mainly based on manual continuum subtraction and 0D slab models. Aims. We investigate the physical conditions under which the gas emits in protoplanetary disks. Based on MIRI spectra, we apply a full Bayesian analysis that provides the posterior distributions of dust and molecular properties, such as column densities and emission temperatures. Methods. To do so, we introduced the Dust Continuum Kit with Line emission from Gas (DuCKLinG), a Python-based model simultaneously describing the molecular line emission and the dust continuum of protoplanetary disks without large computational cost. The model describes the dust continuum emission by dust models with precomputed dust opacities. The molecular emission is based on LTE slab models but from extended radial ranges with gradients in column densities and emission temperatures. We compare the model to observations using Bayesian analysis with linear regression techniques to reduce the dimension of the parameter space. We benchmarked this model to a complex thermo-chemical ProDiMo model of AATau and fit the MIRI spectrum of GW Lup. The latter allowed for a comparison to the previous results obtained with single slab models and hand-fitted continuum. Results. We successfully decrease the computational time of the fitting method by a factor of 80 by eliminating linear parameters, such as the emission areas, from the Bayesian run. This approach does not significantly change the retrieved molecular parameters, and only the calculated errors on the optically thin dust masses slightly decrease. For an AA Tau ProDiMo mock observation, we find that the retrieved molecular conditions from DuCKLinG (column densities from 3 × 1018 cm−2 to 4 × 1020 cm−2, radial range from 0.2 au to 1.2 au, and temperature range from about 200 K to 400 K) fall within the true values from ProDiMo (column densities between 4 × 1017 cm-2 to 5 × 1020 cm−2, radial extent 0.1 au to 6.6 au, and temperature range from about 120 to 1000 K). The smaller DuCKLinG ranges can be explained by the relative flux contributions of the different parts of ProDiMo. The parameter posterior of GW Lup reinforces previously found results. The previously determined column densities fall within the retrieved ranges in this study for all examined molecules (CO2, H2O, HCN, and C2H2). Similar overlap is found for the temperatures with only the temperature range of HCN (from 570−60+60 to 750−70+90 K) not including the previously found value (875 K). This discrepancy may be due to the simultaneous fitting of all molecules compared to the step-by-step fitting of the previous study. There is statistically significant evidence for radial temperature and column density gradients for H2O and CO2 compared to the constant temperature and column density assumed in the 0D slab models. Additionally, HCN and C2H2 emit from a small region with near constant conditions. Due to the small selected wavelength range 13.6–16.3 µm, the dust properties are not well constrained for GW Lup. DuCKL inG can become an important tool to analyse the molecular emission and dust mineralogy of large samples based on JWST /MIRI spectra in an automated way.

Planet Formation by Gas-assisted Accretion of Small Solids

We compute the accretion efficiency of small solids, with radii 1 cm ≤ R s ≤ 10 m, on planets embedded in gaseous disks. Planets have masses 3 ≤ M p ≤ 20 Earth masses (M ⊕) and orbit within 10 au of a solar mass star. Disk thermodynamics is modeled via 3D radiation-hydrodynamics calculations that typically resolve the planetary envelopes. Both icy and rocky solids are considered, explicitly modeling their thermodynamic evolution. The maximum efficiencies of 1 ≤ R s ≤ 100 cm particles are generally ≲10%, whereas 10 m solids tend to accrete efficiently or be segregated beyond the planet’s orbit. A simplified approach is applied to compute the accretion efficiency of small cores, with masses M p ≤ 1 M ⊕ and without envelopes, for which efficiencies are approximately proportional to Mp2/3 . The mass flux of solids, estimated from unperturbed drag-induced drift velocities, provides typical accretion rates dM p /dt ≲ 10−5 M ⊕ yr−1. In representative disk models with an initial gas-to-dust mass ratio of 70–100 and total mass of 0.05–0.06 M ⊙, the solids’ accretion falls below 10−6 M ⊕ yr−1 after 1–1.5 Myr. The derived accretion rates, as functions of time and planet mass, are applied to formation calculations that compute dust opacity self-consistently with the delivery of solids to the envelope. Assuming dust-to-solid coagulation times of ≈0.3 Myr and disk lifetimes of ≈3.5 Myr, heavy-element inventories in the range 3–7 M ⊕ require that ≈90–150 M ⊕ of solids cross the planet’s orbit. The formation calculations encompass a variety of outcomes, from planets a few times M ⊕, predominantly composed of heavy elements, to giant planets. The peak luminosities during the epoch of the solids’ accretion range from ≈10−7 to ≈10−6 L ⊙.

Spatial Variations of Dust Opacity and Grain Growth in Dark Clouds: L1689, L1709, and L1712

The far-infrared (FIR) opacity of dust in dark clouds within the Ophiuchus molecular cloud is investigated through multiwavelength infrared observations from UKIDSS, Spitzer, and Herschel. Employing the infrared color excess technique with both near-infrared and mid-infrared photometric data, a high-resolution extinction map in the K band (A K ) is constructed for three dark clouds: L1689, L1709, and L1712. The derived extinction map has a resolution of 1′ and reaches a depth of A K ∼ 3 mag. The FIR optical depths τ 250 at a reference wavelength of 250 μm are obtained by fitting the Herschel PACS and SPIRE continuum data at 100, 160, 250, 350, and 500 μm using a modified blackbody model. The average dust opacity per unit gas mass at 250 μm, r κ 250, is determined through a pixel-by-pixel correlation of τ 250 with A K , yielding a value of approximately 0.09 cm2 g−1, which is about 2–3 times higher than the typical value in the diffuse interstellar medium. Additionally, an independent analysis across 16 subregions within the Ophiuchus cloud indicates spatial variations in dust opacity, with values ranging from 0.07 to 0.12 cm2 g−1. Although the observed trend of increasing dust opacity with higher extinction implies grain growth, our findings indicate that rapid grain growth has clearly not yet occurred in the dark clouds studied in this work.

A deep search for large complex organic species toward IRAS16293-2422 B at 3 mm with ALMA

Complex organic molecules (COMs) have been detected ubiquitously in protostellar systems. However, at shorter wavelengths ($ sim 0.8$\,mm), it is generally more difficult to detect larger molecules than at longer wavelengths sim 3$\,mm) because of the increase in millimeter dust opacity, line confusion, and unfavorable partition function. We aim to search for large molecules (more than eight atoms) in the Atacama Large Millimeter/submillimeter Array (ALMA) Band 3 spectrum of IRAS 16293-2422 B. In particular, the goal is to quantify the usability of ALMA Band 3 for molecular line surveys in comparison to similar studies at shorter wavelengths. We used deep ALMA Band 3 observations of IRAS 16293-2422 B to search for more than 70 molecules and identified as many lines as possible in the spectrum. The spectral settings were set to specifically target three-carbon species such as i- and n-propanol and glycerol, the next step after glycolaldehyde and ethylene glycol in the hydrogenation of CO. We then derived the column densities and excitation temperatures of the detected species and compared the ratios with respect to methanol between Band 3 ($ sim 3$\,mm) and Band 7 sim 1$\,mm, Protostellar Interferometric Line Survey) observations of this source to examine the effect of the dust optical depth. We identified lines of 31 molecules including many oxygen-bearing COMs such as CH$_3$OH, CH$_2$OHCHO, CH$_3$CH$_2$OH, and c-C$_2$H$_4$O and a few nitrogen- and sulfur-bearing ones such as HOCH$_2$CN and CH$_3$SH. The largest detected molecules are gGg-(CH$_2$OH)$_2$ and CH$_3$COCH$_3$. We did not detect glycerol or i- and n-propanol, but we do provide upper limits for them which are in line with previous laboratory and observational studies. The line density in Band 3 is only $ sim 2.5$ times lower in frequency space than in Band 7. From the detected lines in Band 3 at a $ level sim 25-30<!PCT!>$ of them could not be identified indicating the need for more laboratory data of rotational spectra. We find similar column densities and column density ratios of COMs (within a factor $ sim 2$) between Band 3 and Band 7. The effect of the dust optical depth for IRAS 16293-2422 B at an off-source location on column densities and column density ratios is minimal. Moreover, for warm protostars, long wavelength spectra ($ sim 3$\,mm) are not only crowded and complex, but they also take significantly longer integration times than shorter wavelength observations sim 0.8$\,mm) to reach the same sensitivity limit. The 3\,mm search has not yet resulted in the detection of larger and more complex molecules in warm sources. A full deep ALMA Band $2-3$ (i.e. sim 3-4$\,mm wavelengths) survey is needed to assess whether low frequency data have the potential to reveal more complex molecules in warm sources.