Search

Context. Multiple systems are common in field stars, and the frequency is found to be higher in early evolutionary stages. Thus, the study of young multiple systems during the embedded stages is key to have a comprehensive understanding of star formation. In particular, how material accretes from the large-scale envelope into the inner region and how this flow interacts with the system physically and chemically has not yet been well characterized observationally. Aims. We aim to provide a snapshot of the forming protobinary system SVS13A, consisting of VLA4A and VLA4B. This includes clear pictures of its kinematic structures, physical conditions, and chemical properties. Methods. We conducted ALMA observations toward SVS13A targeting CH3CN and CH3-13CN J=12-11 K-ladder line emission with a high spatial resolution of ~30 au at a spectral resolution of ~0.08 km s-1 Results. We perform LTE radiative transfer models to fit the spectral features of the line emission. We find the two-layer LTE radiative model including dust absorption is essential to interpret the CH3CN and CH3-13-CN line emission. We identify two major and four small kinematic components, and derive their physical and chemical properties. Conclusions. We find a possible infalling signature toward the bursting secondary source VLA4A, which may be fed by an infalling streamer from the large-scale envelope. The mechanical heating in the binary system, as well as the infalling shocked gas, likely play a role in the thermal structure of the protobinary system. By accumulating mass from the streamer, the system might have experienced a gravitationally unstable phase before the accretion outburst. Finally, the derived CH3CN/CH3-13-CN ratio is lower than the canonical ratio in the ISM and is different between VLA4A and VLA4B., Comment: 14 pages, 13 figures

Pre-stellar cores are the first steps in the process of star and planet formation. However, the dynamical and chemical evolution of pre-stellar cores is still not well understood. We aim at estimating the central density of the pre-stellar core IRAS16293E and at carrying out an inventory of molecular species towards the density peak of the core. We observed high-$J$ rotational transitions of N$_2$H$^+$ and N$_2$D$^+$, and several other molecular lines towards the dust emission peak using the Atacama Pathfinder EXperiment (APEX) telescope, and derived the density and temperature profiles of the core using far-infrared surface brightness maps from $Herschel$. The N$_2$H$^+$ and N$_2$D$^+$ lines were analysed by non-LTE radiative transfer modelling. Our best-fit core model consists in a static inner region, embedded in an infalling envelope with an inner radius of approximately 3000 au (21" at 141 pc). The observed high-J lines of N$_2$H$^+$ and N$_2$D$^+$ (with critical densities greater than 10$^6$ cm$^{-3}$) turn out to be very sensitive to depletion; the present single-dish observations are best explained with no depletion of N$_2$H$^+$ and N$_2$D$^+$ in the inner core. The N$_2$D$^+$/N$_2$H$^+$ ratio that best reproduces our observations is 0.44, one of the largest observed to date in pre-stellar cores. Additionally, half of the molecules that we observed are deuterated isotopologues, confirming the high-level of deuteration towards this source. Non-LTE radiative transfer modelling of N$_2$H$^+$ and N$_2$D$^+$ lines proved to be an excellent diagnostic of the chemical structure and dynamics of a pre-stellar core. Probing the physical conditions immediately before the protostellar collapse is a necessary reference for theoretical studies and simulations with the aim of understanding the earliest stages of star and planet formation and the time scale of this process., Comment: Accepted for publication in A&A

We present 5-1x10^7 Angstrom spectral energy distributions (SEDs) for twelve M dwarf stars covering spectral types M0-M8. Our SEDs are provided for community use as a sequel to the Measurements of the Ultraviolet Spectral Characteristics of Low-mass Exoplanetary Systems (MUSCLES) survey. The twelve stars include eight known exoplanet hosts and four stars chosen to fill out key parameter space in spectral type and rotation period. The SEDs are constructed from Hubble Space Telescope ultraviolet spectroscopy and XMM Newton, Chandra and/or Swift X-ray observations and completed with various model data, including Lyman alpha reconstructions, PHOENIX optical models, APEC coronal models and Differential Emission Measure models in the currently-unobservable Extreme Ultraviolet. We provide a complete overview of the Mega-MUSCLES program, including a description of the observations, models, and SED construction. The SEDs are available as MAST High-Level Science Products and we describe the various data products here. We also present ensemble measurements from our sample that are of particular relevance to exoplanet science, including the high-energy fluxes in the habitable zone and the FUV/NUV ratio. Combined with MUSCLES, Mega-MUSCLES provides SEDs covering a wide range of M\,dwarf spectral types and ages such that suitable proxies for any M dwarf planet host of interest may be found in our sample. However, we find that ultraviolet and X-ray fluxes can vary even between stars with similar parameters, such that observations of each exoplanet host star will remain the gold standard for interpreting exoplanet atmosphere observations., Comment: Updated Mega-MUSCLES SEDs available at https://doi.org/10.5281/zenodo.4556130

Context. The formation of stars has been subject to extensive studies in the past decades from molecular cloud to protoplanetary disk scales. It is still not fully understood how the surrounding material in a protostellar system, that often shows asymmetric structures with complex kinematic properties, feeds the central protostar(s) and their disk(s). Aims. We study the spatial morphology and kinematic properties of the molecular gas surrounding the IRS3A and IRS3B protostellar systems in the L1448N region located in the Perseus molecular cloud. Methods. We present 1 mm NOEMA observations of the PRODIGE large program and analyze the kinematic properties of molecular lines. Given the complexity of the spectral profiles, the lines are fitted with up to three Gaussian velocity components. The clustering algorithm DBSCAN is used to disentangle the velocity components into the underlying physical structure. Results. We discover an extended gas bridge (~3000 au) surrounding both the IRS3A and IRS3B systems in six molecular line tracers (C18O, SO, DCN, H2CO, HC3N, and CH3OH). This gas bridge is oriented along the northeast-southwest direction and shows clear velocity gradients on the order of 100 km/s/pc towards the IRS3A system. We find that the observed velocity profile is consistent with analytical streamline models of gravitational infall towards IRS3A. The high-velocity C18O (2-1) emission towards IRS3A indicates a protostellar mass of ~1.2 Msun. Conclusions. While high angular resolution continuum data often show IRS3A and IRS3B in isolation, molecular gas observations reveal that these systems are still embedded within a large-scale mass reservoir with a complex spatial morphology as well as velocity profiles. The kinematic properties of the extended gas bridge are consistent with gravitational infall toward the IRS3A protostar., Comment: 16 pages, 10 figures; accepted for publication in A&A

(Abridged) We explore the chemistry of the most abundant C, O, S, and N bearing species in molecular clouds, in the context of the IRAM 30 m Large Programme Gas phase Elemental abundances in Molecular Clouds (GEMS). In this work, we aim to assess the limitations introduced in the observational works when a uniform density is assumed along the line of sight for fitting the observations, developing a very simple numerical model of a turbulent box. We perform a MHD simulation in order to reproduce the turbulent steady-state of a turbulent box with properties typical of a molecular filament before collapse. We post-process the results of the MHD simulation with a chemical code to predict molecular abundances, and then post-process this cube with a radiative transfer code to create synthetic emission maps for a series of rotational transitions observed during the GEMS project. From the chemical point of view, we find that turbulence produces variations on the predicted abundances, but they are more or less critical depending on the chosen transition and the chemical age. When compared to real observations, the results from the turbulent simulation provides a better fit than when assuming a uniform gas distribution along the line of sight. In the view of our results, we conclude that taking into account turbulence when fitting observations might significantly improve the agreement with model predictions. This is especially important for sulfur bearing species that are very sensitive to the variations of density produced by turbulence at early times (0.1 Myr). The abundance of CO is also quite sensitive to turbulence when considering the evolution beyond a few 0.1 Myr., Comment: Published in Astronomy and Astrophysics

The characterization of young planets (< 300 Myr) is pivotal for understanding planet formation and evolution. We present the 3-5$\mu$m transmission spectrum of the 17 Myr, Jupiter-size ($R$ $\sim$10$R_{\oplus}$) planet, HIP 67522 b, observed with JWST/NIRSpec/G395H. To check for spot contamination, we obtain a simultaneous $g$-band transit with SOAR. The spectrum exhibits absorption features 30-50% deeper than the overall depth, far larger than expected from an equivalent mature planet, and suggests that HIP 67522 b's mass is $<$20 $M_{\oplus}$ irrespective of cloud cover and stellar contamination. A Bayesian retrieval analysis returns a mass constraint of $13.8\pm1.0M_{\oplus}$. This challenges the previous classification of HIP 67522 b as a hot Jupiter and instead, positions it as a precursor to the more common sub-Neptunes. With a density of $<$0.10g/cm$^{3}$, HIP 67522 b is one of the lowest density planets known. We find strong absorption from H$_{2}$O and CO$_{2}$ ($\ge7\sigma$), a modest detection of CO (3.5$\sigma$), and weak detections of H$_2$S and SO$_2$ ($\simeq2\sigma$). Comparisons with radiative-convective equilibrium models suggest supersolar atmospheric metallicities and solar-to-subsolar C/O ratios, with photochemistry further constraining the inferred atmospheric metallicity to 3$\times$10 Solar due to the amplitude of the SO$_2$ feature. These results point to the formation of HIP 67522 b beyond the water snowline, where its envelope was polluted by icy pebbles and planetesimals. The planet is likely experiencing substantial mass loss (0.01-0.03 M$_{\oplus}$ Myr$^{-1}$), sufficient for envelope destruction within a Gyr. This highlights the dramatic evolution occurring within the first 100 Myr of its existence., Comment: Accepted for publication in The Astronomical Journal; 32 pages, 18 figures, 7 tables

Radio detections of stellar systems provide a window onto stellar magnetic activity and the space weather conditions of extrasolar planets, information that is difficult to attain at other wavelengths. There have been recent advances observing auroral emissions from radio-bright low-mass stars and exoplanets largely due to the maturation of low-frequency radio instruments and the plethora of wide-field radio surveys. To guide us in placing these recent results in context, we introduce the foremost local analogues for the field: Solar bursts and the aurorae found on Jupiter. We detail how radio bursts associated with stellar flares are foundational to the study of stellar coronae, and time-resolved radio dynamic spectra offers one of the best prospects of detecting and characterising coronal mass ejections from other stars. We highlight the prospects of directly detecting coherent radio emission from exoplanetary magnetospheres, and early tentative results. We bridge this discussion to the field of brown dwarf radio emission, in which their larger and stronger magnetospheres are amenable to detailed study with current instruments. Bright, coherent radio emission is also predicted from magnetic interactions between stars and close-in planets. We discuss the underlying physics of these interactions and implications of recent provisional detections for exoplanet characterisation. We conclude with an overview of outstanding questions in theory of stellar, star-planet interaction, and exoplanet radio emission, and the prospects of future facilities in answering them., Comment: Accepted to Nature Astronomy. The manuscript is designed to be a primer for new doctoral students and scholars to the field of radio stars and exoplanets. 36 pages, 3 figures

The Milky Way's Central Molecular Zone (CMZ) differs dramatically from our local solar neighbourhood, both in the extreme interstellar medium conditions it exhibits (e.g. high gas, stellar, and feedback density) and in the strong dynamics at play (e.g. due to shear and gas influx along the bar). Consequently, it is likely that there are large-scale physical structures within the CMZ that cannot form elsewhere in the Milky Way. In this paper, we present new results from the Atacama Large Millimeter/submillimeter Array (ALMA) large programme ACES (ALMA CMZ Exploration Survey) and conduct a multi-wavelength and kinematic analysis to determine the origin of the M0.8$-$0.2 ring, a molecular cloud with a distinct ring-like morphology. We estimate the projected inner and outer radii of the M0.8$-$0.2 ring to be 79" and 154", respectively (3.1 pc and 6.1 pc at an assumed Galactic Centre distance of 8.2 kpc) and calculate a mean gas density $> 10^{4}$ cm$^{-3}$, a mass of $\sim$ $10^6$ M$_\odot$, and an expansion speed of $\sim$ 20 km s$^{-1}$, resulting in a high estimated kinetic energy ($> 10^{51}$ erg) and momentum ($> 10^7$ M$_\odot$ km s$^{-1}$). We discuss several possible causes for the existence and expansion of the structure, including stellar feedback and large-scale dynamics. We propose that the most likely cause of the M0.8$-$0.2 ring is a single high-energy hypernova explosion. To viably explain the observed morphology and kinematics, such an explosion would need to have taken place inside a dense, very massive molecular cloud, the remnants of which we now see as the M0.8$-$0.2 ring. In this case, the structure provides an extreme example of how supernovae can affect molecular clouds., Comment: 17 pages, 13 figures, and 2 tables. Published in Astronomy & Astrophysics

When the planet formation process begins in the disks surrounding young stars is still an open question. Annular substructures such as rings and gaps in disks are intertwined with planet formation, and thus their presence or absence is commonly used to investigate the onset of this process. Current observations show a limited number of disks surrounding protostars exhibiting annular substructures, all of them in the Class I stage. The lack of observed features in most of these sources may indicate a late emergence of substructures, but it could also be an artifact of these disks being optically thick. To mitigate the problem of optical depth, we investigate substructures within a very young Class 0 disk characterized by a low inclination using observations at longer wavelengths. We use 3 mm ALMA observations tracing dust emission at a resolution of 7 au to search for evidence of annular substructures in the disk around the deeply embedded Class 0 protostar Oph A SM1. The observations reveal a nearly face-on disk (i$\sim$16$^{\circ}$) extending up to 40 au. The radial intensity profile shows a clear deviation from a smooth profile near 30 au, which we interpret as the presence of either a gap at 28 au or a ring at 34 au with Gaussian widths of $\sigma=1.4^{+2.3}_{-1.2}$ au and $\sigma=3.9^{+2.0}_{-1.9}$ au, respectively. The 3 mm emission at the location of the possible gap or ring is determined to be optically thin, precluding the possibility that this feature in the intensity profile is due to the emission being optically thick. Annular substructures resembling those in the more evolved Class I and II disks could indeed be present in the Class 0 stage, earlier than previous observations suggested. Similar observations of embedded disks in which the high optical depth problem can be mitigated are clearly needed to better constrain the onset of substructures in the embedded stages., Comment: 12 pages, 7 figures, 2 tables. Accepted for publication in A&A Letters

Recent observations suggest that planets formation starts early, in protostellar disks of $\le10^5$ yrs, which are characterized by strong interactions with the environment, e.g., through accretion streamers and molecular outflows. To investigate the impact of such phenomena on disk physical and chemical properties it is key to understand what chemistry planets inherit from their natal environment. In the context of the ALMA Large Program Fifty AU STudy of the chemistry in the disk/envelope system of Solar-like protostars (FAUST), we present observations on scales from ~1500 au to ~60 au of H$_2$CO, HDCO, and D$_2$CO towards the young planet-forming disk IRS~63. H$_2$CO probes the gas in the disk as well as in a large scale streamer (~1500 au) impacting onto the South-East (SE) disk side. We detect for the first time deuterated formaldehyde, HDCO and D$_2$CO, in a planet-forming disk, and HDCO in the streamer that is feeding it. This allows us to estimate the deuterium fractionation of H$_2$CO in the disk: [HDCO]/[H$_2$CO]$\sim0.1-0.3$ and [D$_2$CO]/[H$_2$CO]$\sim0.1$. Interestingly, while HDCO follows the H$_2$CO distribution in the disk and in the streamer, the distribution of D$_2$CO is highly asymmetric, with a peak of the emission (and [D]/[H] ratio) in the SE disk side, where the streamer crashes onto the disk. In addition, D$_2$CO is detected in two spots along the blue- and red-shifted outflow. This suggests that: (i) in the disk, HDCO formation is dominated by gas-phase reactions similarly to H$_2$CO, while (ii) D$_2$CO was mainly formed on the grain mantles during the prestellar phase and/or in the disk itself, and is at present released in the gas-phase in the shocks driven by the streamer and the outflow. These findings testify on the key role of streamers in the build-up of the disk both concerning the final mass available for planet formation and its chemical composition., Comment: 12 pages, 10 figures, accepted for publication on A&A

The high-energy radiative output, from the X-ray to the ultraviolet, of exoplanet host stars drives photochemical reactions and mass loss in the upper regions of planetary atmospheres. In order to place constraints on the atmospheric properties of the three closest terrestrial exoplanets transiting M dwarfs, we observe the high-energy spectra of the host stars LTT1445A and GJ486 in the X-ray with XMM-Newton and Chandra and in the ultraviolet with HST/COS and STIS. We combine these observations with estimates of extreme ultraviolet flux, reconstructions of the Ly-a lines, and stellar models at optical and infrared wavelengths to produce panchromatic spectra from 1A--20um for each star. While LTT1445Ab, LTT1445Ac, and GJ486b do not possess primordial hydrogen-dominated atmospheres, we calculate that they are able to retain pure CO2 atmospheres if starting with 10, 15, and 50% of Earth's total CO2 budget, respectively, in the presence of their host stars' stellar wind. We use age-activity relationships to place lower limits of 2.2 and 6.6 Gyr on the ages of the host stars LTT1445A and GJ486. Despite both LTT1445A and GJ486 appearing inactive at optical wavelengths, we detect flares at ultraviolet and X-ray wavelengths for both stars. In particular, GJ486 exhibits two flares with absolute energies of 10^29.5 and 10^30.1 erg (equivalent durations of 4357+/-96 and 19724+/-169 s) occurring three hours apart, captured with HST/COS G130M. Based on the timing of the observations, we suggest that these high-energy flares are related and indicative of heightened flaring activity that lasts for a period of days, but our interpretations are limited by sparse time-sampling. Consistent high-energy monitoring is needed to determine the duration and extent of high-energy activity on individual M dwarfs, as well as the population as a whole., Comment: 21 pages, published in A&A

We aim to constrain the gas density and temperature distributions as well as gas masses in several T Tauri protoplanetary disks located in Taurus. We use the 12CO, 13CO, and C18O (2-1) isotopologue emission observed at 0.9 with the IRAM NOrthern Extended Millimeter Array (NOEMA) as part of the MPG-IRAM Observatory Program PRODIGE (PROtostars and DIsks: Global Evolution PIs: P. Caselli & Th. Henning). Our sample consists of Class II disks with no evidence of strong radial substructures. We use thesedata to constrain the thermal and chemical structure of these disks through theoretical models for gas emission. To fit the combined optically thick and thin CO line data in Fourier space, we developed the DiskCheF code, which includes the parameterized disk physical structure, machine-learning (ML) accelerated chemistry, and the RADMC-3D line radiative transfer module. A key novelty of DiskCheF is the fast and feasible ML-based chemistry trained on the extended grid of the disk physical-chemical models precomputed with the ANDES2 code. This ML approach allows complex chemical kinetics models to be included in a time-consuming disk fitting without the need to run a chemical code. We present a novel approach to incorporate chemistry into disk modeling without the need to explicitly calculate a chemical network every time. Using this new disk modeling tool, we successfully fit the 12CO, 13CO, and C18O (2-1) data from the CI, CY, DL, DM, DN, and IQ Tau disks. The combination of optically thin and optically thick CO lines allows us to simultaneously constrain the disk temperature and mass distribution, and derive the CO-based gas masses. These values are in reasonable agreement with the disk dust masses rescaled by a factor of 100 as well as with other indirect gas measurements., Comment: 25 pages, 19 figures

It has been consensus that star-disk systems accrete most of their mass and angular momentum during the collapse of a prestellar core, such that the rotational direction of a system is equivalent to the net rotation of the core. Recent results, however, indicate that stars experience post-collapse or late infall, during which the star and its disk is refreshed with material from the protostellar environment through accretion streamers. Apart from adding mass to the star-disk system, infall potentially supplies a substantial amount of angular momentum as the infalling material is initially not bound to the collapsing prestellar core. We investigate the orientation of infall on star-disk systems by analyzing the properties of accreting tracer particles in 3D magnetohydrodynamical simulations of a molecular cloud that is (4 pc)$^3$ in volume. In contrast to the traditional picture, where the rotational axis is inherited from the collapse of a coherent pre-stellar core, the orientation of star-disk systems can change substantially during the accretion process. In agreement with previous results that show larger contributions of late infall for increasing stellar masses, infall is more likely to lead to a prolonged change in orientation for stars of higher final mass. On average, brown dwarfs and very low mass stars are more likely to form and accrete all of their mass as part of a multiple system, while stars with final masses above a few 0.1 M$_{\odot}$ are more likely to accrete part of their mass as single stars. Finally, we find an overall trend: the post-collapse accretion phase is more anisotropic than the early collapse phase. This result is consistent with a scenario, where mass accretion from infall occurs via infalling streamers along a preferred direction, while the initial collapse is less anisotropic albeit the fact that material is funneled through accretion channels., Comment: 15 pages including 5 pages of appendix, submitted to A&A on April 16; the paper is a follow-up on the 2023 publication "Rejuvenating infall: a crucial yet overlooked source of mass and angular momentum" (https://ui.adsabs.harvard.edu/abs/2023EPJP..138..272K/abstract), comments are welcome

We present ALMA 12-m, 7-m & Total Power (TP) Array observations of the FU Orionis outbursting system, covering spatial scales ranging from 160 to 25,000 au. The high-resolution interferometric data reveals an elongated $^{12}$CO(2-1) feature previously observed at lower resolution in $^{12}$CO(3-2). Kinematic modeling indicates that this feature can be interpreted as an accretion streamer feeding the binary system. The mass infall rate provided by the streamer is significantly lower than the typical stellar accretion rates (even in quiescent states), suggesting that this streamer alone is not massive enough to sustain the enhanced accretion rates characteristic of the outbursting class prototype. The observed streamer may not be directly linked to the current outburst but rather a remnant of a previous, more massive streamer that may have contributed enough to the disk mass to render it unstable and trigger FU Ori's outburst. The new data detects, for the first time, a vast, slow-moving carbon monoxide molecular outflow emerging from this object. To accurately assess the outflow properties (mass, momentum, kinetic energy), we employed $^{13}$CO(2-1) data to correct for optical depth effects. The analysis indicates that the outflow corresponds to swept-up material not associated with the current outburst, similar to slow-molecular outflows observed around other FUor and Class I protostellar objects.

The origin of the chemical diversity observed around low-mass protostars probably resides in the earliest history of these systems. We aim to investigate the impact of protostellar feedback on the chemistry and grain growth in the circumstellar medium of multiple stellar systems. In the context of the ALMA Large Program FAUST, we present high-resolution (50 au) observations of CH$_3$OH, H$_2$CO, and SiO and continuum emission at 1.3 mm and 3 mm towards the Corona Australis star cluster. Methanol emission reveals an arc-like structure at $\sim$1800 au from the protostellar system IRS7B along the direction perpendicular to the major axis of the disc. The arc is located at the edge of two elongated continuum structures that define a cone emerging from IRS7B. The region inside the cone is probed by H$_2$CO, while the eastern wall of the arc shows bright emission in SiO, a typical shock tracer. Taking into account the association with a previously detected radio jet imaged with JVLA at 6 cm, the molecular arc reveals for the first time a bow shock driven by IRS7B and a two-sided dust cavity opened by the mass-loss process. For each cavity wall, we derive an average H$_2$ column density of $\sim$7$\times$10$^{21}$ cm$^{-2}$, a mass of $\sim$9$\times$10$^{-3}$ M$_\odot$, and a lower limit on the dust spectral index of $1.4$. These observations provide the first evidence of a shock and a conical dust cavity opened by the jet driven by IRS7B, with important implications for the chemical enrichment and grain growth in the envelope of Solar System analogues., Comment: 12 pages, 8 figures, 3 tables. Accepted Letter in Astronomy & Astrophysics

Complex Organic Molecules (COMs) have been found toward low-mass protostars but the origins of the COM emission are still unclear. It can be associated with, for example, hot corinos, outflows, and/or accretion shock/disk atmosphere. We have conducted NOEMA observations toward SVS13A from the PROtostars & DIsks: Global Evolution (PRODIGE) program. Our previous \ce{DCN} observations reveal a possible infalling streamer, which may affect the chemistry of the central protobinary by inducing accretion outbursts and/or shocked gas. Here, we further analyze six O-bearing COMs: CH3OH, aGg'-(CH2OH)2, C2H5OH, CH2(OH)CHO, CH3CHO, and CH3OCHO. Although the COM emission is not spatially resolved, we constrain the source sizes to $\lesssim0.3-0.4$ arcsec (90$-$120 au) by conducting uv-domain Gaussian fitting. Interestingly, the high-spectral resolution data reveal complex line profiles with multiple peaks showing differences between these six O-bearing COMs. The LTE fitting unveils differences in excitation temperatures and emitting areas among these COMs. We further conduct multiple-velocity-component LTE fitting to decompose the line emission into different kinematic components. Up to 6 velocity components are found from the LTE modeling. The temperature, column density, and source size of these components from each COM are obtained. We find a variety in excitation temperatures ($100-500$ K) and source sizes (D$\sim10-70$ au) from these kinematic components from different COMs. The emission of each COM can trace several components and different COMs most likely trace different regions. Given this complex structure, we suggest that the central region is inhomogeneous and unlikely to be heated by only protostellar radiation. We conclude that accretion shocks induced by the large-scale infalling streamer likely exist and contribute to the complexity of the COM emission., Comment: 29 pages, 18 figures, accepted to A&A

Despite a burgeoning set of ultracool dwarf ($\leq$M7) radio detections, their radio emissions remain enigmatic. Open questions include the plasma source and acceleration mechanisms for the non-auroral "quiescent" component of these objects' radio emissions, which can trace Jovian synchrotron radiation belt analogs. Ultracool dwarf binary systems can provide test beds for examining the underlying physics for these plasma processes. We extend a recently developed occurrence rate calculation framework to compare the quiescent radio occurrence rate of binary systems to single objects. This generalized and semi-analytical framework can be applied to any set of astrophysical objects conceptualized as unresolved binary systems with approximately steady-state emission or absorption. We combine data available in the literature to create samples of 179 single ultracool dwarfs (82 M dwarfs, 74 L dwarfs, and 23 T/Y dwarfs) and 27 binary ultracool dwarf systems. Using these samples, we show that quiescent radio emissions occur in $54^{+11}_{-11}$ - $65^{+11}_{-12}$ per cent of binaries where both components are ultracool dwarfs, depending on priors. We also show that binarity enhances the ultracool dwarf quiescent radio occurrence rate relative to their single counterparts. Finally, we discuss potential implications for the underlying drivers of ultracool dwarf quiescent radio emissions, including possible plasma sources.

Recent observations have demonstrated that very-low mass stars and brown dwarfs are capable of sustaining strong magnetic fields despite their cool and neutral atmospheres. These kG field strengths are inferred based on strong highly circularly polarized GHz radio emission, a consequence of the electron cyclotron maser instability. Crucially, these observations imply the existence of energetic non-thermal electron populations, associated with strong current systems, as are found in the auroral regions of the magnetized planets of the Solar System. Intense auroral electron precipitation will lead to electron collisions with the H$_{2}$ gas that should ultimately generate the ion H$_{3}^{+}$. With this motivation, we targeted a sample of ultracool dwarfs, known to exhibit signatures associated with aurorae, in search of the K-band emission features of H$_{3}^{+}$ using the Keck telescopes on Mauna Kea. From our sample of 9 objects, we found no clear indication of H$_{3}^{+}$ emission features in our low-medium resolution spectra (R$\sim$3600). We also modeled the impact of an auroral electron beam on a brown dwarf atmosphere, determining the depth at which energetic beams deposit their energy and drive particle impact ionization. We find that the H$_{3}^{+}$ non-detections can be explained by electron beams of typical energies $\gtrsim$2-10~keV, which penetrate deeply enough that any H$_{3}^{+}$ produced is chemically destroyed before radiating energy through its infrared transitions. Strong electron beams could further explain the lack of UV detections, and suggest that most or nearly all of the precipitating auroral energy must ultimately emerge as thermal emissions deep in brown dwarf atmospheres., Comment: accepted to ApJ

We are performing a line survey of 8 planet-forming Class II disks in Taurus with the IRAM NOrthern Extended Millimeter Array (NOEMA), as a part of the MPG-IRAM Observatory Program PRODIGE (PROtostars and DIsks: Global Evolution; PIs: P. Caselli and Th. Henning). Compact and extended disks around T Tauri stars CI, CY, DG, DL, DM, DN, IQ Tau, and UZ Tau E are observed in ~80 lines from >20 C-, O,- N-, and S-bearing species. The observations in four spectral settings at 210-280 GHz with $1\sigma$ rms sensitivity of ~ 8-12 mJy/beam at 0.9" and 0.3 km/s resolution will be completed in 2024. The uv-visibilities are fitted with the DiskFit model to obtain key stellar and disk properties. In this paper, the combined $^{12}$CO, $^{13}$CO and C$^{18}$O $J = 2-1$ data are presented. We find that the CO fluxes and disk masses inferred from dust continuum tentatively correlate with the CO emission sizes. We constrain dynamical stellar masses, geometries, temperatures, the CO column densities and gas masses for each disk. The best-fit temperatures at 100 au are ~ 17-37 K, and decrease radially with the power-law exponent q ~ 0.05-0.76. The inferred CO column densities decrease radially with the power-law exponent p ~ 0.2-3.1. The gas masses estimated from $^{13}$CO (2-1) are ~ $0.001-0.2 M_\textrm{Sun}$. The best-fit CO column densities point to severe CO freeze-out in the disks. The DL Tau disk is an outlier, and has either stronger CO depletion or lower gas mass than the rest of the sample. The CO isotopologue ratios are roughly consistent with the observed values in disks and the low-mass star-forming regions., Comment: 47 pages, 10 figures, 6 tables, accepted for publication by A&A on 22/02/2024

The ALMA interferometer has played a key role in revealing a new component of the Sun-like star forming process: the molecular streamers, i.e. structures up to thousands of au long funneling material non-axisymmetrically to disks. In the context of the FAUST ALMA LP, the archetypical VLA1623-2417 protostellar cluster has been imaged at 1.3 mm in the SO(5$_6$--4$_5$), SO(6$_6$--5$_5$), and SiO(5--4) line emission at the spatial resolution of 50 au. We detect extended SO emission, peaking towards the A and B protostars. Emission blue-shifted down to 6.6 km s$^{-1}$ reveals for the first time a long ($\sim$ 2000 au) accelerating streamer plausibly feeding the VLA1623 B protostar. Using SO, we derive for the first time an estimate of the excitation temperature of an accreting streamer: 33$\pm$9 K. The SO column density is $\sim$ 10$^{14}$ cm$^{-2}$, and the SO/H$_2$ abundance ratio is $\sim$ 10$^{-8}$. The total mass of the streamer is 3 $\times$ 10$^{-3}$ $Msun$, while its accretion rate is 3--5 $\times$ 10$^{-7}$ Msun yr$^{-1}$. This is close to the mass accretion rate of VLA1623 B, in the 0.6--3 $\times$ 10$^{-7}$ Msun yr$^{-1}$ range, showing the importance of the streamer in contributing to the mass of protostellar disks. The highest blue- and red-shifted SO velocities behave as the SiO(5--4) emission, the latter species detected for the first time in VLA1623-2417: the emission is compact (100-200 au), and associated only with the B protostar. The SO excitation temperature is $\sim$ 100 K, supporting the occurrence of shocks associated with the jet, traced by SiO., Comment: accepted by MNRAS

We aim at advancing our understanding of magnetic activity and the underlying dynamo mechanism at the end of the main-sequence. To this end, we collected simultaneous X-ray and radio observations for a sample of M7-L0 dwarfs using XMM-Newton jointly with the JVLA and the ATCA. We also included data from the all-sky surveys of eROSITA on board the Russian Spektrum-Roentgen-Gamma mission (SRG) and rotation periods from TESS together with an archival compilation of X-ray and radio data for UCDs from the literature. We limit the sample to objects with rotation period <1d, focusing on the study of a transition in magnetic behavior suggested by a drastic change in the radio detection rate at vsini \approx 38 km/s. We compiled the most up-to-date radio/X-ray luminosity relation for 26 UCDs with rotation periods lower than 1d, finding that rapid rotators lie the furthest away from the G\"udel-Benz relation previously studied for earlier-type stars. Radio bursts are mainly experienced by very fast UCDs, while X-ray flares are seen along the whole range of rotation. We examined the L_{\rm x}/L_{\rm bol} vs P_{\rm rot} relation, finding no evident relation between the X-ray emission and rotation, reinforcing previous speculations on a bimodal dynamo across late-type dwarfs. One radio-detected object has a rotation period consistent with the range of auroral bursting sources; while it displays moderately circularly polarized emission. A radio flare from this object is interpreted as gyrosynchrotron emission, and it displays X-ray and optical flares. We also found a slowly rotating apparent auroral emitter, that is also one of the X-ray brightest radio-detected UCDs. We speculate that this UCD is experiencing a transition of its magnetic behavior since it produces signatures expected from higher mass M dwarfs along with emerging evidence of auroral emission., Comment: A&A submitted

In this work, we study deuterium fractionation in four starless cores in the low-mass star-forming region L1688 in the Ophiuchus molecular cloud. We study how the deuterium fraction ($R_D$) changes with environment, compare deuteration of ions and neutrals, core centre and its envelope, and attempt to reproduce the observed results with a gas-grain chemical model. We chose high and low gas density tracers to study both core centre and the envelope. With the IRAM 30m antenna, we mapped N$_2$H$^+$(1-0), N$_2$D$^+$(1-0), H$^{13}$CO$^+$ (1-0) and (2-1), DCO$^+$(2-1), and $p$-NH$_2$D(1$_{11}$-1$_{01}$) towards the chosen cores. The missing $p$-NH$_3$ and N$_2$H$^+$(1-0) data were taken from the literature. To measure the molecular hydrogen column density, dust and gas temperature within the cores, we used the Herschel/SPIRE dust continuum emission data, the GAS survey data (ammonia), and the COMPLETE survey data to estimate the upper limit on CO depletion. We present the deuterium fraction maps for three species towards four starless cores. Deuterium fraction of the core envelopes traced by DCO$^+$/H$^{13}$CO$^+$ is one order of magnitude lower ($\sim$0.08) than that of the core central parts traced by the nitrogen-bearing species ($\sim$0.5). Deuterium fraction increases with the gas density as indicated by high deuterium fraction of high gas density tracers and low deuterium fraction of lower gas density tracers and by the decrease of $R_D$ with core radii, consistent with the predictions of the chemical model. Our model results show a good agreement with observations for $R_D$(N$_2$D$^+$/N$_2$H$^+$) and R$_D$(DCO$^+$/HCO$^+$) and underestimate the $R_D$(NH$_2$D/NH$_3$)., Comment: 30 pages with 22 figures, accepted for publication in MNRAS

We aim to investigate the origin of the HI Br$\gamma$ emission in young stars by using GRAVITY to image the innermost region of circumstellar disks, where important physical processes such as accretion and winds occur. With high spectral and angular resolution, we focus on studying the continuum and the HI Br$\gamma$-emitting area of the Herbig star HD58647. Using VLTI-GRAVITY, we conducted observations of HD58647 with both high spectral and high angular resolution. Thanks to the extensive $uv$ coverage, we were able to obtain detailed images of the circumstellar environment at a sub-au scale, specifically capturing the continuum and the Br$\gamma$-emitting region. Through the analysis of velocity-dispersed images and photocentre shifts, we were able to investigate the kinematics of the HI Br$\gamma$-emitting region. The recovered continuum images show extended emission where the disk major axis is oriented along a position angle of 14\degr. The size of the continuum emission at 5-sigma levels is $\sim$ 1.5 times more extended than the sizes reported from geometrical fitting (3.69 mas $\pm$ 0.02 mas). This result supports the existence of dust particles close to the stellar surface, screened from the stellar radiation by an optically thick gaseous disk. Moreover, for the first time with GRAVITY, the hot gas component of HD58647 traced by the Br$\gamma$ ,has been imaged. This allowed us to constrain the size of the Br$\gamma$-emitting region and study the kinematics of the hot gas; we find its velocity field to be roughly consistent with gas that obeys Keplerian motion. The velocity-dispersed images show that the size of the hot gas emission is from a more compact region than the continuum (2.3 mas $\pm$ 0.2 mas). Finally, the line phases show that the emission is not entirely consistent with Keplerian rotation, hinting at a more complex structure in the hot gaseous disk., Comment: accepted in A&A on 20/11/2023

V1298 Tau is a young pre-main sequence star hosting four known exoplanets that are prime targets for transmission spectroscopy with current-generation instruments. This work pieces together observations from the NICER X-ray telescope, the Space Telescope Imaging Spectrograph and Cosmic Origins Spectrograph instruments aboard Hubble Space Telescope, and empirically informed models to create a panchromatic spectral energy distribution for V1298 Tau spanning 1 to 100000 Angstroms. We describe the methods and assumptions used to assemble the panchromatic spectrum and show that despite this star's brightness, its high-energy spectrum is near the limit of present X-ray and ultraviolet observatories' abilities to characterize. We conclude by using the V1298 Tau spectrum as a benchmark for the activity saturation stage of high-energy radiation from solar-mass stars to compare the lifetime cumulative high-energy irradiation of the V1298 Tau planets to other planets orbiting similarly massive stars., Comment: The V1298 Tau panchromatic SED will be available via the Data Behind the Figure feature for Figure 1 of the final journal article, will be uploaded to the Exoplanet Archive as an IPAC table, and can be provided upon request

The role of massive ($\geq$ 8M$_{\odot}$) stars in defining the energy budget and chemical enrichment of the interstellar medium in their host galaxy is significant. In this first paper from the Tracing Evolution in Massive Protostellar Objects (TEMPO) project we introduce a colour-luminosity selected (L$_*$ $\sim$ 3$\times10^3$ to 1$\times10^5$ L$_{\odot}$) sample of 38 massive star forming regions observed with ALMA at 1.3mm and explore the fragmentation, clustering and flux density properties of the sample. The TEMPO sample fields are each found to contain multiple fragments (between 2-15 per field). The flux density budget is split evenly (53%-47%) between fields where emission is dominated by a single high flux density fragment and those in which the combined flux density of fainter objects dominates. The fragmentation scales observed in most fields are not comparable with the thermal Jeans length, $\lambda_J$, being larger in the majority of cases, suggestive of some non-thermal mechanism. A tentative evolutionary trend is seen between luminosity of the clump and the `spectral line richness' of the TEMPO fields; with 6.7GHz maser associated fields found to be lower luminosity and more line rich. This work also describes a method of line-free continuum channel selection within ALMA data and a generalised approach used to distinguishing sources which are potentially star-forming from those which are not, utilising interferometric visibility properties., Comment: 42 pages, 17 figures, 4 Tables including 3 Appendices (with additional tables and figures). Accepted for publication in MNRAS

We present the serendipitous discovery of 8 distant ($>$ 50 pc) late M dwarfs with plausible associated radio emission at 144 MHz. The M dwarf nature of our sources has been confirmed with optical spectroscopy performed using HET/LRS2 and Subaru/FOCAS, and their radio flux densities are within the range of 0.5-1.0 mJy at 144 MHz. Considering the radio-optical source separation and source densities of the parent catalogues, we suggest that it is statistically probable the M dwarfs are associated with the radio emission. However, it remains plausible that for some of the sources the radio emission originates from an optically faint and red galaxy hiding behind the M dwarf. The isotropic radio luminosities ($\sim10^{17-18}$ erg s$^{-1}$ Hz$^{-1}$) of the M dwarfs suggest that if the association is real, the radio emission is likely driven by a coherent emission process produced via plasma or electron-cyclotron maser instability processes, which is potentially caused by binary interaction. Long term monitoring in the radio and high-resolution radio follow-up observations are necessary to search for any variability and pinpoint the radio emission to determine whether our tentative conclusion that these ultracool dwarfs are radio emitting is correct. If the low-frequency radio emission is conclusively associated with the M dwarfs, this would reveal a new population of optically faint and distant ($>$ 50 pc) radio emitting M dwarfs., Comment: 10 pages, 5 figures, accepted for publication in A&A

The infall of gas from outside natal cores has proven to feed protostars after the main accretion phase (Class 0). This changes our view of star formation to a picture that includes asymmetric accretion (streamers), and a larger role of the environment. However, the connection between streamers and the filaments that prevail in star-forming regions is unknown. We investigate the flow of material toward the filaments within Barnard 5 (B5) and the infall from the envelope to the protostellar disk of the embedded protostar B5-IRS1. Our goal is to follow the flow of material from the larger, dense core scale, to the protostellar disk scale. We present new HC$_3$N line data from the NOEMA and 30m telescopes covering the coherence zone of B5, together with ALMA H$_2$CO and C$^{18}$O maps toward the protostellar envelope. We fit multiple Gaussian components to the lines so as to decompose their individual physical components. We investigate the HC$_3$N velocity gradients to determine the direction of chemically-fresh gas flow. At envelope scales, we use a clustering algorithm to disentangle the different kinematic components within H$_2$CO emission. At dense core scales, HC$_3$N traces the infall from the B5 region toward the filaments. HC$_3$N velocity gradients are consistent with accretion toward the filament spines plus flow along them. We found a $\sim2800$ au streamer in H$_2$CO emission which is blueshifted with respect to the protostar and deposits gas at outer disk scales. The strongest velocity gradients at large scales curve toward the position of the streamer at small scales, suggesting a connection between both flows. Our analysis suggests that the gas can flow from the dense core to the protostar. This implies that the mass available for a protostar is not limited to its envelope, and can receiving chemically-unprocessed gas after the main accretion phase., Comment: 25 pages, 27 figures, accepted for publication on Astronomy and Astrophysics. The scripts used for analysis can be seen at https://github.com/tere-valdivia/Barnard_5_infall

Studies of dense core morphologies and their orientations with respect to gas flows and the local magnetic field have been limited to only a small sample of cores with spectroscopic data. Leveraging the Green Bank Ammonia Survey alongside existing sub-millimeter continuum observations and Planck dust polarization, we produce a cross-matched catalogue of 399 dense cores with estimates of core morphology, size, mass, specific angular momentum, and magnetic field orientation. Of the 399 cores, 329 exhibit 2D $\mathrm{v}_\mathrm{LSR}$ maps that are well fit with a linear gradient, consistent with rotation projected on the sky. We find a best-fit specific angular momentum and core size relationship of $J/M \propto R^{1.82 \pm 0.10}$, suggesting that core velocity gradients originate from a combination of solid body rotation and turbulent motions. Most cores have no preferred orientation between the axis of core elongation, velocity gradient direction, and the ambient magnetic field orientation, favouring a triaxial and weakly magnetized origin. We find, however, strong evidence for a preferred anti-alignment between the core elongation axis and magnetic field for protostellar cores, revealing a change in orientation from starless and prestellar populations that may result from gravitational contraction in a magnetically-regulated (but not dominant) environment. We also find marginal evidence for anti-alignment between the core velocity gradient and magnetic field orientation in the L1228 and L1251 regions of Cepheus, suggesting a preferred orientation with respect to magnetic fields may be more prevalent in regions with locally ordered fields., Comment: 33 pages, 28 figures, accepted to MNRAS

Characterizing rocky exoplanet atmospheres is a key goal of exoplanet science, but interpreting such observations will require understanding the stellar UV irradiation incident on the planet from its host star. Stellar UV mediates atmospheric escape, photochemistry, and planetary habitability, and observations of rocky exoplanets can only be understood in the context of the UV SED of their host stars. Particularly important are SEDs from observationally favorable but poorly understood low-mass M-dwarf stars, which are the only plausible targets for rocky planet atmospheric characterization for the next 1-2 decades. In this work, we explore the utility of AstroSat UVIT for the characterization of the UV SEDs of low-mass stars. We present observations of the nearby M0 star HIP 23309 in the FUV and NUV gratings of UVIT. Our FUV spectra are consistent with contemporaneous HST data and our NUV spectra are stable between orbits, suggesting UVIT is a viable tool for the characterization of the SEDs of low-mass stars. We apply our measured spectra to simulations of photochemistry and habitability for a hypothetical rocky planet orbiting HIP 23309 and elucidate the utility and limitations of UVIT in deriving UV SEDs of M-dwarf exoplanet hosts. Our work validates UVIT as a tool to complement HST in the characterization of exoplanet host stars and carries implications for its successor missions like INSIST., Comment: Accepted to AJ

In the centre of pre-stellar cores, deuterium fractionation is enhanced due to the low temperatures and high densities. Therefore, the chemistry of deuterated molecules can be used to study the earliest stages of star formation. We analyse the deuterium fractionation of simple molecules, comparing the level of deuteration in the envelopes of the pre-stellar core L1544 in Taurus and the protostellar core HH211 in Perseus. We used single-dish observations of CCH, HCN, HNC, HCO$^+$, and their $^{13}$C-, $^{18}$O- and D-bearing isotopologues, detected with the Onsala 20m telescope. We derived the column densities and the deuterium fractions of the molecules. Additionally, we used radiative transfer simulations and results from chemical modelling to reproduce the observed molecular lines. We used new collisional rate coefficients for HNC, HN$^{13}$C, DNC, and DCN that consider the hyperfine structure of these molecules. We find high levels of deuteration for CCH (10%) in both sources, consistent with other carbon chains, and moderate levels for HCN (5-7%) and HNC (8%). The deuterium fraction of HCO$^+$ is enhanced towards HH211, most likely caused by isotope-selective photodissociation of C$^{18}$O. Similar levels of deuteration show that the process is likely equally efficient towards both cores, suggesting that the protostellar envelope still retains the chemical composition of the original pre-stellar core. The fact that the two cores are embedded in different molecular clouds also suggests that environmental conditions do not have a significant effect on the deuteration within dense cores. Radiative transfer modelling shows that it is necessary to include the outer layers of the cores to consider the effects of extended structures. Besides HCO$^+$ observations, HCN observations towards L1544 also require the presence of an outer diffuse layer where the molecules are relatively abundant., Comment: 27 pages, 17 figures, accepted for publication in A&A

White dwarfs stars are known to be polluted by their active planetary systems, but little attention has been paid to the accretion of wind from low-mass companions. The capture of stellar or substellar wind by white dwarfs is one of few methods available to astronomers which can assess mass-loss rates from unevolved stars and brown dwarfs, and the only known method to extract their chemical compositions. In this work, four white dwarfs with closely-orbiting, L-type brown dwarf companions are studied to place limits on the accretion of a substellar wind, with one case of a detection, and at an extremely non-solar abundance $m_{\rm Na}/m_{\rm Ca}>900$. The mass-loss rates and upper limits are tied to accretion in the white dwarfs, based on limiting cases for how the wind is captured, and compared with known cases of wind pollution from close M dwarf companions, which manifest in solar proportions between all elements detected. For wind captured in a Bondi-Hoyle flow, mass-loss limits $\dot M\lesssim 5\times10^{-17}$ M$_\odot$ yr$^{-1}$ are established for three L dwarfs, while for M dwarfs polluting their hosts, winds in the range $10^{-13} - 10^{-16}$ M$_\odot$ yr$^{-1}$ are found. The latter compares well with the $\dot M\sim 10^{-13} - 10^{-15}$ M$_\odot$ yr$^{-1}$ estimates obtained for nearby, isolated M dwarfs using Ly$\alpha$ to probe their astropsheres. These results demonstrate that white dwarfs are highly-sensitive stellar and substellar wind detectors, where further work on the actual captured wind flow is needed., Comment: 13 pages, 2 figures, 1 table, accepted to MNRAS

X-ray through infrared spectral energy distributions (SEDs) are essential for understanding a star's effect on exoplanet atmospheric composition and evolution. We present a catalog of panchromatic SEDs, hosted on the Barbara A. Mikulski Archive for Space Telescopes (MAST), for 11 exoplanet hosting stars which have guaranteed JWST observation time as part of the ERS or GTO programs but have no previous UV characterization. The stars in this survey range from spectral type F4-M6 (0.14-1.57 M$_\odot$), rotation periods of ~4-132 days, and ages of approximately 0.5-11.4 Gyr. The SEDs are composite spectra using data from the Chandra X-ray Observatory and XMM-Newton, the Hubble Space Telescope, BT-Settl stellar atmosphere models, and scaled spectra of proxy stars of similar spectral type and activity. From our observations, we have measured a set of UV and X-ray fluxes as indicators of stellar activity level. We compare the chromospheric and coronal activity indicators of our exoplanet-hosting stars to the broader population of field stars and find that a majority of our targets have activity levels lower than the average population of cool stars in the solar neighborhood. This suggests that using SEDs of stars selected from exoplanet surveys to compute generic exoplanet atmosphere models may underestimate the typical host star's UV flux by an order of magnitude or more, and consequently, that the observed population of exoplanetary atmospheres receive lower high-energy flux levels than the typical planet in the solar neighborhood., Comment: 28 pages, 10 figures, 7 tables, accepted for publication in AJ

Early dust grain growth in protostellar envelopes infalling on young discs has been suggested in recent studies, supporting the hypothesis that dust particles start to agglomerate already during the Class 0/I phase of young stellar objects (YSOs). If this early evolution were confirmed, it would impact the usually assumed initial conditions of planet formation, where only particles with sizes $\lesssim 0.25 \mu$m are usually considered for protostellar envelopes. We aim to determine the maximum grain size of the dust population in the envelope of the Class 0/I protostar L1527 IRS, located in the Taurus star-forming region (140 pc). We use Atacama Large millimetre/sub-millimetre Array (ALMA) and Atacama Compact Array (ACA) archival data and present new observations, in an effort to both enhance the signal-to-noise ratio of the faint extended continuum emission and properly account for the compact emission from the inner disc. Using observations performed in four wavelength bands and extending the spatial range of previous studies, we aim to place tight constraints on the spectral ($\alpha$) and dust emissivity ($\beta$) indices in the envelope of L1527 IRS. We find a rather flat $\alpha \sim$ 3.0 profile in the range 50-2000 au. Accounting for the envelope temperature profile, we derive values for the dust emissivity index, 0.9 < $\beta$ < 1.6, and reveal a tentative, positive outward gradient. This could be interpreted as a distribution of mainly ISM-like grains at 2000 au, gradually progressing to (sub-)millimetre-sized dust grains in the inner envelope, where at R=300 au, $\beta$ = 1.1 +/- 0.1. Our study supports a variation of the dust properties in the envelope of L1527 IRS. We discuss how this can be the result of in-situ grain growth, dust differential collapse from the parent core, or upward transport of disc large grains., Comment: Accepted for publication in A&A. Contains 18 pages, 21 figures, 5 tables Replacement on Nov 22 to change title number of FAUST series from "X" to "IX."

Observing magnetic star-planet interactions (SPI) offers promise for determining magnetic fields of exoplanets. Models of sub-Alfv\'enic SPI predict that terrestrial planets in close-in orbits around M~dwarfs can induce detectable stellar radio emission, manifesting as bursts of strongly polarized coherent radiation observable at specific planet orbital positions. We present 2-4 GHz detections of coherent radio bursts on the slowly-rotating M dwarf YZ Ceti, which hosts a compact system of terrestrial planets, the innermost orbiting with a 2-day period. Two coherent bursts occur at similar orbital phases of YZ Cet b, suggestive of an enhanced probability of bursts near that orbital phase. We model the system's magnetospheric environment in the context of sub-Alfv\'enic SPI and determine that YZ Ceti b can plausibly power the observed flux densities of the radio detections. However, we cannot rule out stellar magnetic activity, without a well characterized rate of non-planet-induced coherent radio bursts on slow rotators. YZ Ceti is therefore a candidate radio SPI system, with unique promise as a target for the long-term monitoring., Comment: Authors' version of article published in Nature Astronomy, see their website for official version of scientific record

More than 50% of solar-mass stars form in multiple systems. It is therefore crucial to investigate how multiplicity affects the star and planet formation processes at the protostellar stage. We report continuum and C$^{18}$O (2-1) observations of the VLA 1623-2417 protostellar system at 50 au angular resolution as part of the ALMA Large Program FAUST. The 1.3 mm continuum probes the disks of VLA 1623A, B, and W, and the circumbinary disk of the A1+A2 binary. The C$^{18}$O emission reveals, for the first time, the gas in the disk-envelope of VLA 1623W. We estimate the dynamical mass of VLA 1623W, $M_{\rm dyn}=0.45\pm0.08$ M$_{\odot}$, and the mass of its disk, $M_{\rm disk}\sim6\times10^{-3}$ M$_{\odot}$. C$^{18}$O also reveals streamers that extend up to 1000 au, spatially and kinematically connecting the envelope and outflow cavities of the A1+A2+B system with the disk of VLA 1623W. The presence of the streamers, as well as the spatial ($\sim$1300 au) and velocity ($\sim$2.2 km/s) offset of VLA 1623W suggest that either sources W and A+B formed in different cores, interacting between them, or that source W has been ejected from the VLA 1623 multiple system during its formation. In the latter case, the streamers may funnel material from the envelope and cavities of VLA 1623AB onto VLA 1623W, thus concurring to set its final mass and chemical content., Comment: 10 pages, 5 figures, accepted by MNRAS

Context: Core accretion models of massive star formation require the existence of massive, starless cores within molecular clouds. Yet, only a small number of candidates for such truly massive, monolithic cores are currently known. Aims: Here we analyse a massive core in the well-studied infrared-dark cloud (IRDC) called the 'dragon cloud' (also known as G028.37+00.07 or 'Cloud C'). This core (C2c1) sits at the end of a chain of a roughly equally spaced actively star-forming cores near the centre of the IRDC. Methods: We present new high-angular resolution 1 mm ALMA dust continuum and molecular line observations of the massive core. Results: The high-angular resolution observations show that this region fragments into two cores C2c1a and C2c1b, which retain significant background-subtracted masses of 23 Msun and 2 Msun (31 Msun and 6 Msun without background subtraction), respectively. The cores do not appear to fragment further on the scales of our highest angular resolution images (0.200 arcsec, 0.005 pc ~ 1000 AU). We find that these cores are very dense (nH2 > 10^6 cm-3) and have only trans-sonic non-thermal motions (Ms ~ 1). Together the mass, density and internal motions imply a virial parameter of < 1, which suggests the cores are gravitationally unstable, unless supported by strong magnetic fields with strengths of ~ 1 - 10 mG. From CO line observations, we find that there is tentative evidence for a weak molecular outflow towards the lower-mass core, and yet the more massive core remains devoid of any star formation indicators. Conclusions: We present evidence for the existence of a massive, pre-stellar core, which has implications for theories of massive star formation. This source warrants follow-up higher-angular-resolution observations to further assess its monolithic and pre-stellar nature., Comment: 8+4 pages, 4+2 Figures, 2 Tables, accepted for publication in Astronomy & Astrophysics

The high energy X-ray and ultraviolet (UV) radiation fields of exoplanet host stars play a crucial role in controlling the atmospheric conditions and the potential habitability of exoplanets. Major surveys of the X-ray/UV emissions from late-type (K and M spectral type) exoplanet hosts have been conducted by the MUSCLES and Mega-MUSCLES Hubble Space Telescope (HST) Treasury programs. These samples primarily consist of relatively old, ``inactive'', low mass stars. In this paper we present results from X-ray observations of the coronal emission from these stars obtained using the Chandra X-ray Observatory, the XMM-Newton Observatory, and the Neil Gehrels Swift Observatory. The stars effectively sample the coronal activity of low-mass stars at a wide range of masses and ages. The vast majority (21 of 23) of the stars are detected and their X-ray luminosities measured. Short-term flaring variability is detected for most of the fully-convective (M $\leq$ 0.35 M$_{\odot}$) stars but not for the more massive M dwarfs during these observations. Despite this difference, the mean X-ray luminosities for these two sets of M dwarfs are similar with more massive (0.35 M$_{\odot}$ $\leq$ M $\leq$ 0.6 M$_{\odot}$) M dwarfs at $\sim$5 $\times$ 10$^{26}$ erg s$^{-1}$ compared to $\sim$2 $\times$ 10$^{26}$ erg s$^{-1}$ for fully-convective stars older than 1 Gyr. Younger, fully-convective M dwarfs have X-ray luminosities between 3 and 6 $\times$ 10$^{27}$ erg s$^{-1}$.The coronal X-ray spectra have been characterized and provide important information that is vital for the modeling of the stellar EUV spectra., Comment: 39 pages, 15 figures. Accepted for publication in The Astronomical Journal

We address the problem of optimally identifying all kilonovae detected via gravitational wave emission in the upcoming LIGO/Virgo/KAGRA Collaboration observing run, O4, which is expected to be sensitive to a factor of $\sim 7$ more Binary Neutron Stars alerts than previously. Electromagnetic follow-up of all but the brightest of these new events will require $>1$ meter telescopes, for which limited time is available. We present an optimized observing strategy for the Dark Energy Camera during O4. We base our study on simulations of gravitational wave events expected for O4 and wide-prior kilonova simulations. We derive the detectabilities of events for realistic observing conditions. We optimize our strategy for confirming a kilonova while minimizing telescope time. For a wide range of kilonova parameters, corresponding to a fainter kilonova compared to GW170817/AT2017gfo we find that, with this optimal strategy, the discovery probability for electromagnetic counterparts with the Dark Energy Camera is $\sim 80\%$ at the nominal binary neutron star gravitational wave detection limit for the next LVK observing run (190 Mpc), which corresponds to a $\sim 30\%$ improvement compared to the strategy adopted during the previous observing run. For more distant events ($\sim 330$ Mpc), we reach a $\sim 60\%$ probability of detection, a factor of $\sim 2$ increase. For a brighter kilonova model dominated by the blue component that reproduces the observations of GW170817/AT2017gfo, we find that we can reach $\sim 90\%$ probability of detection out to 330 Mpc, representing an increase of $\sim 20 \%$, while also reducing the total telescope time required to follow-up events by $\sim 20\%$., Comment: 26 pages, 11 figures, accepted by ApJ

Characterizing the physical conditions at disk scales in Class 0 sources is crucial for constraining the protostellar accretion process and the initial conditions for planet formation. We use ALMA 1.3 mm and 3 mm observations to investigate the physical conditions of the dust around the Class 0 binary IRAS 16293-2422 A (sep <100 au) down to ~10 au scales. The circumbinary material's spectral index, alpha, has a median of 3.1 and a dispersion of ~0.2, providing no firm evidence of mm-sizes grains therein. Continuum substructures with brightness temperature peaks of T_b~60-80 K at 1.3 mm are observed near the disks at both wavelengths. These peaks do not overlap with strong variations of alpha, indicating they trace high-temperature spots instead of regions with significant optical depth variations. The lower limits to the inferred dust temperature in the hot spots are 122, 87 and 49 K. Depending on the assumed dust opacity index, these values can be several times higher. They overlap with high gas temperatures and enhanced complex organic molecular (COM) emission. This newly resolved dust temperature distribution is in better agreement with the expectations from mechanical instead of the most commonly assumed radiative heating. In particular, we find that the temperatures agree with shock heating predictions. This evidence and recent studies highlighting accretion heating in Class 0 disks suggest that mechanical heating (shocks, dissipation powered by accretion, etc.) is important during the early stages and should be considered when modeling and measuring properties of deeply embedded protostars and disks., Comment: 18 pages, 7 figures. Published in ApJL

Aims. We present high sensitivity and high-spectral resolution NOEMA observations of the Class 0/I binary system SVS13A, composed of the low-mass protostars VLA4A and VLA4B with a separation of ~90 au. VLA4A is undergoing an accretion burst that enriches the chemistry of the surrounding gas. This gives us an excellent opportunity to probe the chemical and physical conditions as well as the accretion process. Methods. We observe the (12K-11K) lines of CH3CN and CH313CN, the DCN (3-2) line, and the C18O (2-1) line toward SVS13A using NOEMA. Results. We find complex line profiles at disk scales which cannot be explained by a single component or pure Keplerian motion. By adopting two velocity components to model the complex line profiles, we find that the temperatures and densities are significantly different between these two components. This suggests that the physical conditions of the emitting gas traced via CH3CN can change dramatically within the circumbinary disk. In addition, combining our observations of DCN (3-2) with previous ALMA high-angular-resolution observations, we find that the binary system (or VLA4A) might be fed by an infalling streamer from envelope scales (~700 au). If this is the case, this streamer contributes to the accretion of material onto the system with a rate of at least 1.4x10-6 Msun yr-1. Conclusions. We conclude that the CH3CN emission in SVS13A traces hot gas from a complex structure. This complexity might be affected by a streamer that is possibly infalling and funneling material into the central region., Comment: 20 pages, 19 figures, accepted to A&A

Prestellar cores, the birthplace of Sun-like stars, form from the fragmentation of the filamentary structure that composes molecular clouds, from which they must inherit at least partially the kinematics. Furthermore, when they are on the verge of gravitational collapse, they show signs of subsonic infall motions. How extended these motions are, which depends on how the collapse occurs, remains largely unknown. We want to investigate the kinematics of the envelope that surrounds the prototypical prestellar core L1544, studying the cloud-core connection. To our aims, we observed the $\rm HCO^+$(1-0) transition in a large map. \hcop is expected to be abundant in the envelope, making it an ideal probe of the large-scale kinematics in the source. We modelled the spectrum at the dust peak by means of a non local-thermodynamical-equilibrium radiative transfer. In order to reproduce the spectrum at the dust peak, a large ($\sim 1\, \rm pc$) envelope is needed, with low density (tens of $\rm cm^{-3}$ at most) and contraction motions, with an inward velocity of $\approx 0.05\,\rm km \, s^{-1}$. We fitted the data cube using the Hill5 model, which implements a simple model {for the optical depth and excitation temperature profiles along the line-of-sight,} in order to obtain a map of the infall velocity. This shows that the infall motions are extended, with typical values in the range $0.1-0.2\,\rm km \, s^{-1}$. Our results suggest that the contraction motions extend in the diffuse envelope surrounding the core, which is consistent with recent magnetic field measurements in the source, which showed that the envelope is magnetically supercritical., Comment: Accepted for publication on ApJ, 24 Oct. 2022

We obtained ultraviolet and optical spectra for 9 M~dwarfs across a range of rotation periods to determine whether they showed stochastic intrinsic variability distinguishable from flares. The ultraviolet spectra were observed during the Far Ultraviolet M~Dwarf Evolution Survey \emph{Hubble~Space~Telescope} program using the Space Telescope Imaging Spectrograph. The optical observations were taken from the Apache Point Observatory 3.5-meter telescope using the Dual Imaging Spectrograph and from the Gemini South Observatory using the Gemini Multi-Object Spectrograph. We used the optical spectra to measure multiple chromospheric lines: the Balmer series from H$\alpha$ to H$10$ and the \ion{Ca}{2}~H and K lines. We find that after excising flares, these lines vary on the order of $1-20\%$ at minute-cadence over the course of an hour. The absolute amplitude of variability was greater for the faster rotating M~dwarfs in our sample. Among the 5 stars for which measured the weaker Balmer lines, we note a tentative trend that the fractional amplitude of the variability increases for higher order Balmer lines. We measured the integrated flux of multiple ultraviolet emission features formed in the transition region: the \ion{N}{5}, \ion{Si}{4} and \ion{C}{4} resonance line doublets, and the \ion{C}{2} and \ion{He}{2} multiplets. The signal-to-noise (S/N) ratio of the UV data was too low for us to detect non-flare variability at the same scale and time cadence as the optical. We consider multiple mechanisms for the observed stochastic variability and propose both observational and theoretical avenues of investigation to determine the physical causes of intrinsic variability in the chromospheres of M~dwarfs.

The near-ultraviolet (NUV) spectral region is a useful diagnostic for stellar flare physics and assessing the energy environment of young exoplanets, especially as relates to prebiotic chemistry. We conducted a pilot NUV spectroscopic flare survey of the young M dwarf AU Mic with the Neil Gehrels Swift Observatory's UltraViolet and Optical Telescope. We detected four flares and three other epochs of significantly elevated count rates during the 9.6 hours of total exposure time, consistent with a NUV flare rate of $\sim$0.5 hour$^{-1}$. The largest flare we observed released a minimum energy of 6$\times$10$^{33}$ erg between 1730-5000 \r{A}. All flares had durations longer than the $\sim$14-17 minute duration of each Swift visit, making measuring total flare energy and duration infeasible., Comment: Published in Research Notes of the AAS (RNAAS)

Context. In the past few years, there has been a rise in the detection of streamers, asymmetric flows of material directed toward the protostellar disk with material from outside the star's natal core. It is unclear how they affect the process of mass accretion, in particular beyond the Class 0 phase. Aims. We investigate the gas kinematics around Per-emb-50, a Class I source in the crowded star-forming region NGC 1333. Our goal is to study how the mass infall proceeds from envelope to disk scales in this source. Results. We discover a streamer delivering material toward Per-emb-50 in H$_2$CO and C$^{18}$O emission. The streamer's emission can be well described by the analytic solutions for an infalling parcel of gas along a streamline with conserved angular momentum, both in the image plane and along the line of sight velocities. The streamer has a mean infall rate of $1.3 \times 10^{ -6}$ M$_{ \odot}$ yr$^{ -1}$, $5 -10$ times higher than the current accretion rate of the protostar. SO and SO$_2$ emission reveal asymmetric infall motions in the inner envelope, additional to the streamer around Per-emb-50. Furthermore, the presence of SO$_2$ could mark the impact zone of the infalling material. Conclusions. The streamer delivers sufficient mass to sustain the protostellar accretion rate and might produce an accretion burst, which would explain the protostar's high luminosity with respect to other Class I sources. Our results highlight the importance of late infall for protostellar evolution: streamers might provide a significant amount of mass for stellar accretion after the Class 0 phase., Comment: 20 pages, 14 figures, accepted for publication in A&A

Methanol is a ubiquitous species commonly found in the molecular interstellar medium. It is also a crucial seed species for the building-up of the chemical complexity in star forming regions. Thus, understanding how its abundance evolves during the star formation process and whether it enriches the emerging planetary system is of paramount importance. We used new data from the ALMA Large Program FAUST (Fifty AU STudy of the chemistry in the disk/envelope system of Solar-like protostars) to study the methanol line emission towards the [BHB2007] 11 protobinary system (sources A and B), where a complex structure of filaments connecting the two sources with a larger circumbinary disk has been previously detected. Twelve methanol lines have been detected with upper energies in the range [45-537] K along with one 13CH3OH transition. The methanol emission is compact and encompasses both protostars, separated by only 28 au and presents three velocity components, not spatially resolved by our observations, associated with three different spatial regions, with two of them close to 11B and the third one associated with 11A. A non-LTE radiative transfer analysis of the methanol lines concludes that the gas is hot and dense and highly enriched in methanol with an abundance as high as 1e-5. Using previous continuum data, we show that dust opacity can potentially completely absorb the methanol line emission from the two binary objects. Although we cannot firmly exclude other possibilities, we suggest that the detected hot methanol is resulting from the shocked gas from the incoming filaments streaming towards [BHB2007] 11 A and B, respectively. Higher spatial resolution observations are necessary to confirm this hypothesis., Comment: 17 pages, accepted in A&A

Well-characterized binary systems will provide valuable opportunities to study the conditions that are necessary for the onset of both auroral and non-auroral magnetospheric radio emission in the ultracool dwarf regime. We present new detections of non-auroral "quiescent" radio emission at 4-8 GHz of the three ultracool dwarf binary systems GJ 564 BC, LP 415-20, and 2MASS J21402931+1625183. We also tentatively detect a highly circularly polarized pulse at 4-6 GHz that may indicate aurorae from GJ 564 BC. Finally, we show that the brightest binary ultracool dwarf systems may be more luminous than predictions from single-object systems.

Protostellar jets are an important agent of star formation feedback, tightly connected with the mass-accretion process. The history of jet formation and mass-ejection provides constraints on the mass accretion history and the nature of the driving source. We want to characterize the time-variability of the mass-ejection phenomena at work in the Class 0 protostellar phase, in order to better understand the dynamics of the outflowing gas and bring more constraints on the origin of the jet chemical composition and the mass-accretion history. We have observed the emission of the CO 2-1 and SO N_J=5_4-4_3 rotational transitions with NOEMA, towards the intermediate-mass Class 0 protostellar system Cep E. The CO high-velocity jet emission reveals a central component associated with high-velocity molecular knots, also detected in SO, surrounded by a collimated layer of entrained gas. The gas layer appears to accelerate along the main axis over a length scale delta_0 ~700 au, while its diameter gradually increases up to several 1000au at 2000au from the protostar. The jet is fragmented into 18 knots of mass ~10^-3 Msun, unevenly distributed between the northern and southern lobes, with velocity variations up to 15 km/s close to the protostar, well below the jet terminal velocities. The knot interval distribution is approximately bimodal with a scale of ~50-80yr close to the protostar and ~150-200yr at larger distances >12". The mass-loss rates derived from knot masses are overall steady, with values of 2.7x10^-5 Msun/yr (8.9x10^-6 Msun/yr) in the northern (southern) lobe. The interaction of the ambient protostellar material with high-velocity knots drives the formation of a molecular layer around the jet, which accounts for the higher mass-loss rate in the north. The jet dynamics are well accounted for by a simple precession model with a period of 2000yr and a mass-ejection period of 55yr., Comment: 13 pages, 9 figures, 3 table. Accepted in A&A

Observations of H I Lyman alpha, the brightest UV emission line of late-type stars, are critical for understanding stellar chromospheres and transition regions, modeling photochemistry in exoplanet atmospheres, and measuring the abundances of neutral hydrogen and deuterium in the interstellar medium. Yet, Lyman alpha observations are notoriously challenging due to severe attenuation from interstellar gas, hindering our understanding of this important emission line's basic morphology. We present high-resolution far- and near-UV spectroscopy of five G, K, and M dwarfs with radial velocities large enough to Doppler shift the stellar Lyman alpha emission line away from much of the interstellar attenuation, allowing the line core to be directly observed. We detect self-reversal in the Lyman alpha emission line core for all targets, and we show that the self-reversal depth decreases with increasing surface gravity. Mg II self-reversed emission line profiles provide some useful information to constrain the Lyman alpha line core, but the differences are significant enough that Mg II cannot be used directly as an intrinsic Lyman alpha template during reconstructions. We show that reconstructions that neglect self-reversal could overestimate intrinsic Lyman alpha fluxes by as much as 60%-100% for G and K dwarfs and 40%-170% for M dwarfs. The five stars of our sample have low magnetic activity and sub-solar metallicity; a larger sample size is needed to determine how sensitive these results are to these factors., Comment: 17 pages, 6 figures, accepted to ApJ

The Orion molecular cloud complex harbours the nearest GMCs and site of high-mass star formation. Its YSO populations are thoroughly characterized. The region is therefore a prime target for the study of star formation. Here, we verify the performance of the SuperCAM 64 pixel heterodyne array on APEX. We give a descriptive overview of a set of wide-field CO(3-2) spectral cubes obtained towards the Orion GMC complex, aimed at characterizing the dynamics and structure of the extended molecular gas in diverse regions of the clouds, ranging from very active sites of clustered star formation in Orion B to comparatively quiet regions in southern Orion A. We present a 2.7 square degree (130pc$^2$) mapping survey in the CO(3-2) transition, obtained using SuperCAM on APEX at an angular resolution of 19'' (7600AU or 0.037pc at a distance of 400pc), covering L1622, NGC2071, NGC2068, OriB9, NGC2024, and NGC2023 in Orion B, and the southern part of the L1641 cloud in Orion A. We describe CO integrated emission and line moment maps and position-velocity diagrams and discuss a few sub-regions in some detail. Evidence for expanding bubbles is seen with lines splitting into double components, most prominently in NGC2024, where we argue that the bulk of the molecular gas is in the foreground of the HII region. High CO(3-2)/CO(1-0) line ratios reveal warm CO along the western edge of Orion B in the NGC2023/NGC2024 region facing the IC434 HII region. Multiple, well separated radial velocity components seen in L1641-S suggest that it consists of a sequence of clouds at increasingly larger distances. We find a small, spherical cloud - the 'Cow Nebula' globule - north of NGC2071. We trace high velocity line wings for the NGC2071-IR outflow and the NGC2024 CO jet. The protostellar dust core FIR4 (rather than FIR5) is the true driving source of the NGC2024 monopolar outflow., Comment: Accepted for publication in Astronomy and Astrophysics