Low-mass End Research Articles

A Close Look at Lyα Emitters with JWST/NIRCam at z ≈ 3.1

We study 10 spectroscopically confirmed Lyα emitters (LAEs) at z ≈ 3.1 in the Ultra Deep Survey field, covered by the James Webb Space Telescope/NIRCam in the Public Release Imaging for Extragalactic Research program. All LAEs are detected in all NIRCam bands from F090W to F444W, corresponding to rest frame wavelengths of 2200 Å to 1.2 μm. Based on morphological analysis of the F200W images, three out of the 10 targets are resolved into pair-like systems with separations of <0.″9, and another three show asymmetric structures. We then construct the spectral energy distributions (SEDs) of these LAEs, which show little to no extinction. All sources, including the pairs, show similar SED shapes, with a prominent flux excess in the F200W band, corresponding to extremely strong [O iii]+Hβ emission lines (EWrest = 740–6500 Å). The median effective radii, stellar mass, and UV slope of our sample are 0.36 kpc, 3.8 × 107 M ⊙, and −2.48, respectively. The average burst age, estimated by stellar mass over star formation rate, is <40 Myr. These measurements reveal an intriguing starbursting dwarf galaxy population lying off the extrapolations of the z ∼ 3 scaling relations to the low-mass end: ∼0.7 dex above the star-forming main sequence, ∼0.35 dex below the mass–size relation, and bluer in the UV slope than typical high-z galaxies at similar UV luminosities. We speculate that these numbers may require a larger main-sequence scatter or tail in the dwarf galaxy regime toward the starburst outliers.

JADES: Insights into the low-mass end of the mass–metallicity–SFR relation at 3 &lt; z &lt; 10 from deep JWST/NIRSpec spectroscopy

We analysed the gas-phase metallicity properties of a sample of low-stellar-mass (log M⋆/M⊙ ≲ 9) galaxies at 3 &lt; z &lt; 10 observed with JWST/NIRSpec as part of the JADES programme in its deep GOODS-S tier. By combining this sample with more massive galaxies at similar redshifts from other programmes, we study the scaling relations between stellar mass (M⋆), oxygen abundance (O/H), and star-formation rate (SFR) for 146 galaxies spanning three orders of magnitude in stellar mass and out to the epoch of early galaxy assembly. We find evidence for a shallower slope at the low-mass end of the mass–metallicity relation (MZR), with 12+log(O/H) = (7.72 ± 0.02) + (0.17 ± 0.03) log(M⋆/108 M⊙), in good agreement with the MZR probed by local analogues of high-redshift systems, such as the ‘Green Pea’ and ‘Blueberry’ galaxies. The inferred slope is well matched by models including ‘momentum-driven’ supernova (SN) winds, suggesting that feedback mechanisms in dwarf galaxies (and at high z) might be different from those in place at higher masses. The evolution in the normalisation is observed to be relatively mild compared to previous determinations of the MZR at z ∼ 3 (∼0.1 − 0.2 dex across the explored mass regime). We observe a deviation from the local fundamental metallicity relation (FMR) for our sample at high redshift, especially at z &gt; 6, with galaxies significantly less enriched than predicted given their M⋆ and SFR (with a median offset in log(O/H) of ∼0.5 dex, significant at ∼5σ). These observations are consistent with an enhanced stochasticity in the gas accretion and star-formation history of high-redshift systems, prompting us to reconsider the nature of the relationship between M⋆, O/H, and SFR in the early Universe.

Recovery of the low- and high-mass end slopes of the IMF in massive early-type galaxies using detailed elemental abundances

ABSTRACT Star formation in the early Universe has left its imprint on the chemistry of observable stars in galaxies. We derive elemental abundances and the slope of the low-mass end of the initial mass function (IMF) for a sample of 25 very massive galaxies, separated into brightest cluster galaxies (BCGs) and their massive satellites. The elemental abundances of BGCs and their satellites are similar, but for some elements, satellite galaxies show a correlation with the global velocity dispersion. Using a subset of derived elemental abundances, we model the star formation histories of these galaxies with chemical evolution models, and predict the high-mass end slope of the IMF and star formation time-scales. The high-mass end IMF slope of the satellite galaxies correlates with the global velocity dispersion. The low- and the high-mass end IMF slopes are weakly correlated in a general sense that top heavy IMFs are paired with bottom heavy IMFs. Our results do not necessarily imply that the IMF was simultaneously bottom and top heavy. Instead, our findings can be considered consistent with a temporal variation in the IMF, where, for massive galaxies, the high-mass end IMF slope is representative of the very early age and the low-mass end slope of the later star formation. The small but noticeable differences between the BCGs and the satellites in terms of their elemental abundances and IMF slopes, together with their stellar kinematical properties, suggest somewhat different formation pathways, where BCGs experience more major, gas-free mergers.

The central black hole in the dwarf spheroidal galaxy Leo I: Not supermassive, at most an intermediate-mass candidate

It has recently been claimed that a surprisingly massive black hole (BH) is present in the core of the dwarf spheroidal galaxy (dSph) Leo I. This finding, based on integral field spectroscopy, challenges the typical expectation that dSphs host intermediate-mass BHs since such a BH would be classified as supermassive. Indeed, the analysis points toward Leo I harboring a BH with a lower mass limit exceeding a few $10^6 at $1 and the no-BH case excluded at 95&lt;!PCT!&gt; significance. Such a mass, which is comparable to the entire stellar mass of the galaxy, makes Leo I a unique system that warrants further investigation. Using equilibrium models based on distribution functions that depend on actions fJ coupled with the same integral field spectroscopy data and an extensive exploration of a very large parameter space, we demonstrate, within a comprehensive Bayesian model--data comparison framework, that the posterior on the BH mass is flat toward the low-mass end and, thus, that the kinematics of the central galaxy region only imposes an upper limit on the BH mass of few $10^5 (at $3 Such an upper limit indicates that the putative BH of Leo I is indeed an intermediate-mass BH, and it is also in line with formation scenarios and expectations from scaling relations at the mass regime of dwarf galaxies.

Closing the Gap between Observed Low-mass Galaxy H i Kinematics and Cold Dark Matter Predictions

Testing the standard cosmological model (ΛCDM) at small scales is challenging. Galaxies that inhabit low-mass dark matter halos provide an ideal test bed for dark matter models by linking observational properties of galaxies at small scales (low mass, low velocity) to low-mass dark matter halos. However, the observed kinematics of these galaxies do not align with the kinematics of the dark matter halos predicted to host them, obscuring our understanding of the low-mass end of the galaxy–halo connection. We use deep H i observations of low-mass galaxies at high spectral resolution in combination with cosmological simulations of dwarf galaxies to better understand the connection between dwarf galaxy kinematics and low-mass halos. Specifically, we use H i line widths to directly compare to the maximum velocities in a dark matter halo and find that each deeper measurement approaches the expected one-to-one relationship between the observed kinematics and the predicted kinematics in ΛCDM. We also measure baryonic masses and place these on the baryonic Tully–Fisher relation (BTFR). Again, our deepest measurements approach the theoretical predictions for the low-mass end of this relation, a significant improvement on similar measurements based on line widths measured at 50% and 20% of the peak. Our data also hint at the rollover in the BTFR predicted by hydrodynamical simulations of ΛCDM for low-mass galaxies.

Cloud properties across spatial scales in simulations of the interstellar medium

Molecular clouds (MCs) are structures of dense gas in the interstellar medium (ISM) that extend from ten to a few hundred parsecs and form the main gas reservoir available for star formation. Hydrodynamical simulations of a varying complexity are a promising way to investigate MCs evolution and their properties. However, each simulation typically has a limited range in resolution and different cloud extraction algorithms are used, which complicates the comparison between simulations. In this work, we aim to extract clouds from different simulations covering a wide range of spatial scales. We compare their properties, such as size, shape, mass, internal velocity dispersion, and virial state. We applied the Hop cloud detection algorithm on (M)HD numerical simulations of stratified ISM boxes and isolated galactic disk simulations that were produced using Flash Ramses and Arepo We find that the extracted clouds are complex in shape, ranging from round objects to complex filamentary networks in all setups. Despite the wide range of scales, resolution, and sub-grid physics, we observe surprisingly robust trends in the investigated metrics. The mass spectrum matches in the overlap between simulations without rescaling and with a high-mass power-law index of -1 for logarithmic bins of mass, in accordance with theoretical predictions. The internal velocity dispersion scales with the size of the cloud as $ $ for large clouds ($R pc $). For small clouds we find larger sigma compared to the power-law scaling, as seen in observations, which is due to supernova-driven turbulence. Almost all clouds are gravitationally unbound with the virial parameter scaling as $ vir $, which is slightly flatter compared to observed scaling but in agreement given the large scatter. We note that the cloud distribution towards the low-mass end is only complete if the more dilute gas is also refined, rather than only the collapsing regions.

Observational Properties of 155 β Cephei Pulsating Variable Stars

β Cephei pulsating variable (BCEP) stars are the most massive pulsating variable stars in the main sequence, exhibiting both p- and g-mode pulsations. In this study, we identified 155 BCEP stars or candidates using data from TESS and Gaia, of which 83 were first confirmed as BCEP stars. They have visual magnitudes ranging from 8 to 12 mag and effective temperatures between approximately 20,000 and 30,000 K, while the parallaxes of most targets are between 0.2 and 0.6 mas. The study indicates that these BCEP stars have pulsation periods ranging from 0.06 to 0.31 days, with amplitudes ranging from 0.1 to 55.8 mmag in the TESS band. Additionally, the number of BCEP stars increases as the pulsation amplitude decreases. These targets align with the distribution region of BCEP stars in the luminosity–period (L–P) and temperature–period diagrams. We have updated the L–P relation of BCEP stars. The Hertzsprung–Russell diagram indicates that these targets are in the main-sequence evolutionary phase, with masses ranging from 7 to 20 M ⊙ and luminosities between 2800 and 71,000 L ⊙. They are almost in the theoretical instability region of BCEP stars but as previously reported, this region at the low-mass end (red) is not filled. The distribution of the pulsation constant indicates that the dominant pulsation periods of BCEP stars consist mainly of low-order p-mode pulsations with a high proportion of radial fundamental modes. These BCEP stars are excellent objects for enhancing our understanding of the structure and evolution of massive stars through asteroseismology.

Influence of protostellar outflows on star and protoplanetary disk formation in a massive star-forming clump

Context. Due to the presence of magnetic fields, protostellar jets or outflows are a natural consequence of accretion onto protostars. They are expected to play an important role in star and protoplanetary disk formation. Aims. We aim to determine the influence of outflows on star and protoplanetary disk formation in star-forming clumps. Methods. Using RAMSES, we performed the first magnetohydrodynamics calculation of massive star-forming clumps with ambipolar diffusion; radiative transfer, including the radiative feedback of protostars; and protostellar outflows while systematically resolving the disk scales. We compared this simulation to a model without outflows. Results. We found that protostellar outflows have a significant impact on both star and disk formation. They provide a significant amount of additional kinetic energy to the clump, with typical velocities of around a few 10 km s−1; impact the clump and disk temperatures; reduce the accretion rate onto the protostars; and enhance fragmentation in the filaments. We found that they promote a more numerous stellar population. They do not impact the low-mass end of the IMF much, which is probably controlled by the mass of the first Larson core; however, they have an influence on its peak and high-mass end. Conclusions. Protostellar outflows appear to have a significant influence on both star and disk formation and should therefore be included in realistic simulations of star-forming environments.

Diffuse Light in Milky Way–like Haloes

We investigate the diffuse light (DL) content of dark matter haloes in the mass range 11.5≤logMhalo≤13 , a range that also includes the dark matter halo of the Milky Way, taking advantage of a state-of-the-art semianalytic model run on the merger trees extracted from a set of high-resolution cosmological simulations. The fraction of DL in such relatively small haloes is found to progressively decrease from the high- to the low-mass end, in good agreement with analytic and numerical results from simulations, and in good agreement also with the fraction of the DL observed in the Milky Way and M31. Haloes with different masses have different efficiencies in producing DL: logMhalo≃13 is found to be the characteristic halo mass where the production of DL is the most efficient, while the overall efficiency decreases at both larger and smaller scales (this work). The DL content in this range of halo mass is the result of stellar stripping due to tidal interaction between satellites and their hosts (95%) and mergers between satellites and the central galaxy (5%), with preprocessed material, a subchannel of mergers and stripping, and so already included in the 100%, that contributes no more than 8% on average. The halo concentration is the main driver of the DL formation: more concentrated haloes have higher DL fractions that come from the stripping of more massive satellites in the high-halo-mass end, while dwarfs contribute mostly in the low-halo-mass end.

Constraints on the Faint End of the Galaxy Stellar Mass Function at z ≃ 4–8 from Deep JWST Data

We analyze a sample of 3300 galaxies between redshifts z ≃ 3.5 and z ≃ 8.5 selected from James Webb Space Telescope (JWST) images in the Hubble Ultra Deep Field and UKIDSS Ultra Deep Survey field, including objects with stellar masses as low as ≃108 M ⊙ up to z ≃ 8. The depth and wavelength coverage of the JWST data allows us, for the first time, to derive robust stellar masses for such high-z, low stellar mass galaxies on an individual basis. We compute the galaxy stellar mass function, after complementing our sample with ancillary data from CANDELS to constrain the GMSF at high stellar masses ( M>M* ). Our results show a steepening of the low stellar mass end slope (α) with redshift, with α = −1.61 ± 0.05 at z ≃ 4 and α = −1.98 ± 0.14 at z ≃ 7. We also observe an evolution of the normalization ϕ * from z ≃ 7 to z ≃ 4, with ϕz≃4*/ϕz≃7*=130−50+210 . Our study incorporates a novel method for the estimation of the Eddington bias, which takes into account its possible dependence both on stellar mass and redshift, while allowing for skewness in the error distribution. We finally compute the resulting cosmic stellar mass density and find a flatter evolution with redshift than previous studies.

Identifying Quiescent Compact Objects in Massive Galactic Single-lined Spectroscopic Binaries

It is now clearly established that massive stars are predominantly found in multiple systems. While their period and eccentricity distributions are now well established across different metallicity regimes, the mass-ratio distribution has been mostly limited to double-lined spectroscopic binaries. The mass-ratio distribution therefore remains subject to significant uncertainties and open questions encompass the shape and extend of the companion mass-function towards its low-mass end and the nature of undetected companions in single-lined spectroscopic binaries, i.e. low-mass main-sequence stars, binary interaction products, compact objects, etc. We have conducted a large and systematic analysis of a sample of more than 80 single-line O-type spectroscopic binaries in the Milky Way and in the Large Magellanic Cloud using a sophisticated multi-technique, and multi-wavelength approach. We report on the developed methodology, the constraints obtained on the nature of SB1 companions, the distribution of O star mass-ratios at LMC metallicity and the occurrence of quiescent OB+black hole binaries.

Modelling post-reionization H i distributions in fuzzy dark matter cosmologies using conditional normalizing flows

ABSTRACT Upcoming 21-cm intensity mapping experiments like the Square Kilometer Array (SKA) hold significant potential to constrain the properties of dark matter. In this work, we model neutral hydrogen (H i) distributions using high-resolution hydrodynamical N-body simulations of both cold dark matter (CDM) and fuzzy dark matter (FDM) cosmologies in the post-reionization redshift range of z = 3.42−4.94. We show that the H i abundance decreases in FDM-like cosmologies. Extreme FDM models with m ∼ 10−22 eV are at odds with a range of measurements. Due to the increased halo bias, the H i bias increases, paralleled by the damped Lyman-α (DLA) bias which we infer from the cross-section of DLAs. The distribution of the latter in extreme FDM models has a high median at the low-mass end, which can be traced to the high column density of cosmic filaments. FDM models exhibit a very similar abundance of DLAs compared to CDM while sub-DLAs are already less abundant. We study the prospects of detecting the brightest H i peaks with SKA1-Low at z = 4.94, indicating moderate signal-to-noise ratios (SNRs) at angular resolution θA = 2 arcmin with a rapidly declining SNR for lower values of θA. After training the conditional normalizing flow network HIGlow on 2D H i maps, we interpolate its latent space of axion masses to predict the peak flux for a new, synthetic FDM cosmology, finding good agreement with expectations. This work thus underscores the potential of normalizing flows in capturing complex, non-linear structures within H i maps, offering a versatile tool for conditional sample generation and prediction tasks.

X-ray scaling relations of early-type galaxies in IllustrisTNG and a new way of identifying backsplash objects

ABSTRACT We investigate how feedback and environment shapes the X-ray scaling relations of early-type galaxies (ETGs), especially at the low-mass end. We select central-ETGs from the TNG100 box of IllustrisTNG that have stellar masses $\log _{10}(M_{\ast }/\mathrm{M_{\odot }})\in [10.7, 11.9]$. We derive mock X-ray luminosity (LX, 500) and spectroscopic-like temperature (Tsl, 500) of hot gas within R500 of the ETG haloes using the MOCK-X pipeline. The scaling between LX, 500 and the total mass within 5 effective radii ($M_{5R_{\rm e}}$) agrees well with observed ETGs from Chandra. IllustrisTNG reproduces the observed increase in scatter of LX, 500 towards lower masses, and we find that ETGs with $\log _{10} (M_{5R_{\rm e}}/\mathrm{M_{\odot }}) \leqslant 11.5$ with above-average LX, 500 experienced systematically lower cumulative kinetic AGN feedback energy historically (vice versa for below-average ETGs). This leads to larger gas mass fractions and younger stellar populations with stronger stellar feedback heating, concertedly resulting in the above-average LX, 500. The LX, 500–Tsl, 500 relation shows a similar slope to the observed ETGs but the simulation systematically underestimates the gas temperature. Three outliers that lie far below the LX–Tsl relation all interacted with larger galaxy clusters recently and demonstrate clear features of environmental heating. We propose that the distinct location of these backsplash ETGs in the LX–Tsl plane could provide a new way of identifying backsplash galaxies in future X-ray surveys.

Kinematical Fluctuations Vary with Galaxy Surface Mass Density

The Galaxy inner parts are generally considered to be optically symmetric, as well as kinematically symmetric for most massive early-type galaxies. At the lower-mass end, many galaxies contain lots of small patches in their velocity maps, causing their kinematics to be nonsmooth in small scales and far from symmetry. These small patches can easily be mistaken for measurement uncertainties and have not been well discussed. We used the comparison of observations and numerical simulations to demonstrate the small patches existence beyond uncertainties. For the first time we have found that the fluctuation degrees have an approximate inverse loglinear relation with the galaxy stellar surface mass densities. This tight relation among galaxies that do not show obvious optical asymmetry that traces environmental perturbations indicates that stellar motion in galaxies has inherent asymmetry besides external environment influences. The degree of the kinetic asymmetry is closely related to and constrained by the intrinsic properties of the host galaxy.

ELVES. IV. The Satellite Stellar-to-halo Mass Relation Beyond the Milky Way

Quantifying the connection between galaxies and their host dark matter halos has been key for testing cosmological models on various scales. Below M ⋆ ∼ 109 M ⊙, such studies have primarily relied on the satellite galaxy population orbiting the Milky Way (MW). Here we present new constraints on the connection between satellite galaxies and their host dark matter subhalos using the largest sample of satellite galaxies in the Local Volume (D ≲ 12 Mpc) to date. We use 250 confirmed and 71 candidate dwarf satellites around 27 MW-like hosts from the Exploration of Local VolumE Satellites (ELVES) Survey and use the semianalytical SatGen model for predicting the population of dark matter subhalos expected in the same volume. Through a Bayesian model comparison of the observed and the forward-modeled satellite stellar mass functions (SSMFs), we infer the satellite stellar-to-halo mass relation. We find that the observed SSMF is best reproduced when subhalos at the low-mass end are populated by a relation of the form , with a moderate slope of and a low scatter, constant as a function of the peak halo mass, of . A model with a steeper slope (α grow = 2.39 ± 0.06) and a scatter that grows with decreasing M peak is also consistent with the observed SSMF but is not required. Our new model for the satellite–subhalo connection, based on hundreds of Local Volume satellite galaxies, is in line with what was previously derived using only MW satellites.

Extending the unified subhalo model to warm dark matter

ABSTRACT Using a set of high-resolution N-body simulations, we extend the unified distribution model of cold dark matter (CDM) subhaloes to the warm dark matter (WDM) case. The same model framework combining the unevolved mass function, unevolved radial distribution, and tidal stripping can predict the mass function and spatial distribution of subhaloes in both CDM and WDM simulations. The dependence of the model on the DM particle property is universally parametrized through the half-mode mass of the initial power spectrum. Compared with the CDM model, the WDM model differs most notably in two aspects. (1) In contrast to the power-law form in CDM, the unevolved subhalo mass function for WDM is scale-dependent at the low mass end due to the cut-off in the initial power spectrum. (2) WDM subhaloes are more vulnerable to tidal stripping and disruption due to their lower concentrations at accretion time. Their survival rate is also found to depend on the infall mass. Accounting for these differences, the model predicts a final WDM subhalo mass function that is also proportional to the unevolved subhalo mass function. The radial distribution of WDM subhaloes is predicted to be mass-dependent. For low mass subhaloes, the radial distribution is flatter in the inner halo and steeper in the outer halo compared to the CDM counterpart, due to the scale-dependent unevolved mass function and the enhanced tidal stripping. The code for sampling subhaloes according to our generalized model is available at https://github.com/fhtouma/subgen2.

Shaping the unseen: the influence of baryons and environment on low-mass, high-redshift dark matter haloes in the SIEGE simulations

ABSTRACT We use zoom-in, hydrodynamical, cosmological N-body simulations tracing the formation of the first stellar clumps from the SImulating the Environments where Globular clusters Emerged project, to study key structural properties of dark matter haloes when the Universe was only $0.92\, {\rm Gyr}$ old. The very high resolution (maximum physical resolution $0.3\, {h}^{-1}\, {\rm pc}$ at z = 6.14, smallest dark matter particle mass $164\, {\rm M}_{\odot }$) allows us to reach the very low mass end of the stellar-to-halo mass relation ($M_{\rm vir}=10^{7.5{\!-\!}9.5}\, {\rm M}_{\odot }$) to study the processes that mould dark matter haloes during the first stages of structure formation. We investigate the role of baryonic cooling and stellar feedback, modelled from individual stars, in shaping haloes, and of environmental effects as accretion of dark matter along cosmic filaments and mergers. We find that the onset of star formation (typically for $\log M_{\rm vir}/\, {\rm M}_{\odot }\simeq 7.6$) causes the inner cusp in the haloes’ density profile to flatten into a core with constant density and size proportionally to the halo virial mass. Even at these mass scales, we confirm that baryons make haloes that have formed stars rounder in the central regions than haloes that have not formed stars yet, with median minor-to-major 〈q〉 and intermediate-to-major 〈s〉 axes 0.66 and 0.84, respectively. Our morphological analysis shows that, at z = 6.14, haloes are largely prolate in the outer parts, with the major axis aligned along filaments of the cosmic web or towards smaller sub-haloes, with the degree of elongation having no significant dependence on the halo mass.

The Mass–Metallicity Relation of Dwarf Galaxies at Cosmic Noon from JWST Observations

We present a study of the mass–metallicity relation (MZR) of 51 dwarf galaxies (M ⋆ ≈ 106.5–109.5 M ⊙) at z = 2–3 from the A2744 and SMACS J0723-3732 galaxy cluster fields. These dwarf galaxies are identified and confirmed by deep JWST/NIRISS imaging and slitless grism spectroscopic observations. By taking advantage of the superior performance of JWST and the gravitational lensing effect, we extend the previous MZR relation at z = 2–3 to a much lower-mass regime down by ≈2.5 orders of magnitude as compared with previous studies. We find that the MZR has a shallower slope at the low-mass end (M ⋆ < 109 M ⊙), with a slope turnover point of ≈109 M ⊙. This implies that the dominating feedback processes in dwarf galaxies may be different from that in massive galaxies. From z = 3, to z = 2, the metallicity of the dwarf galaxies is enhanced by ≈0.09 dex for a given stellar mass, consistent with the mild evolution found in galaxies with higher mass. Furthermore, we confirm the existence of a fundamental metallicity relation (FMR) between the gas-phase metallicity, stellar mass, and star formation rate in dwarf galaxies at z = 2–3. Our derived FMR, which has no significant redshift evolution, can be used as a benchmark to understand the origin of the anticorrelation between the star formation rate and metallicity of dwarf galaxies in the high-z Universe.

Early Planet Formation in Embedded Disks (eDisk). VI. Kinematic Structures around the Very-low-mass Protostar IRAS 16253-2429

Precise estimates of protostellar masses are crucial to characterize the formation of stars of low masses down to brown dwarfs (BDs; M * < 0.08 M ☉). The most accurate estimation of protostellar mass uses the Keplerian rotation in the circumstellar disk around the protostar. To apply the Keplerian rotation method to a protostar at the low-mass end, we have observed the Class 0 protostar IRAS 16253-2429 using the Atacama Large Millimeter/submillimeter Array (ALMA) in the 1.3 mm continuum at an angular resolution of 0.″07 (10 au), and in the 12CO, C18O, 13CO (J = 2–1), and SO (J N = 65−54) molecular lines, as part of the ALMA Large Program Early Planet Formation in Embedded Disks project. The continuum emission traces a nonaxisymmetric, disk-like structure perpendicular to the associated 12CO outflow. The position–velocity (PV) diagrams in the C18O and 13CO lines can be interpreted as infalling and rotating motions. In contrast, the PV diagram along the major axis of the disk-like structure in the 12CO line allows us to identify Keplerian rotation. The central stellar mass and the disk radius are estimated to be ∼0.12–0.17 M ☉ and ∼13–19 au, respectively. The SO line suggests the existence of an accretion shock at a ring (r ∼ 28 au) surrounding the disk and a streamer from the eastern side of the envelope. IRAS 16253-2429 is not a proto-BD but has a central stellar mass close to the BD mass regime, and our results provide a typical picture of such very-low-mass protostars.

Modeling of the nebular-phase spectral evolution of stripped-envelope supernovae

We present an extended grid of multi-epoch 1D nonlocal thermodynamic equilibrium radiative transfer calculations for nebular-phase Type Ibc supernovae (SNe) from He-star explosions. Compared to our previous work, which was focused on a post-explosion epoch of 200 days, here we study the spectral evolution from 100 to about 450 days. We also augment the model set with progenitors that evolved without wind mass loss. Models with the same final, pre-SN mass have similar yields and produce essentially the same emergent spectra. Hence, the uncertain progenitor mass loss history compromises the inference of the initial, main sequence mass. This shortcoming does not affect Type IIb SNe in which mass-loss has left a small residual H-rich envelope in the progenitor star at core collapse and, hence, an intact He core. However, our 1D models with a different pre-SN mass tend to yield widely different spectra, as seen through variations in the strong emission lines due to [N II] λλ 6548, 6583, [O I] λλ 6300, 6364, [Ca II] λλ 7291, 7323, [Ni II] λ 7378, and the forest of Fe II lines below 5500 Å. At the lower mass end, the ejecta are He-rich, and at 100 days, they cool through He I, N II, Ca II, and Fe II lines, with N II and Fe II dominating at 450 days. These models, associated with He giants, stand in conflict to observed SNe Ib, which typically lack strong N II emission. Instead, they may lead to SNe Ibn or, because of additional stripping by a companion star, ultra-stripped SNe Ic. In contrast, for higher pre-SN masses, the ejecta are progressively He poor and cool at 100 days through O I, Ca II, and Fe II lines, with O I and Ca II dominating at 450 days. Non-uniform, aspherical, large-scale mixing is more likely to determine the SN type at intermediate pre-SN masses, rather than any compositional differences. Variations in clumping and mixing, as well as departures from spherical symmetry would increase the spectral diversity, but also introduce additional degeneracies. More robust predictions from spectral modeling thus require that careful attention be paid to the initial conditions by incorporating the salient features of physically consistent 3D explosion models.