Star formation in galaxies is a complex phenomenon occurring on a very wide range of scales, and molecular clouds are at the heart of this process. The formation of these structures and the subsequent collapse of the gas within them to form new stars remain unresolved scientific questions. In particular, the role and importance of gravity at between the disk scale height and prestellar cores (100 to 0.01 pc) are still topics of debate. In this work, we conduct a case study examining the mass assembly and evolution of a giant molecular cloud complex in a numerical stratified-box simulation of the interstellar medium with photo-ionizing and supernova driving and resolving down to scales $\gtrsim 1$ pc and densities up to $10^3$ cm$^{-3}$. By introducing tracer particles to precisely track the forces acting on the gas during its evolution towards and within the clouds, we are able to quantify how much of the mass inflow is driven by the self-gravity of the gas and the gravity from the stellar disk. We find that up to 20% of the gas is gravity-driven at a scale of 100 pc, contributing 10% of the inflow from the warm to the cold phase and 20% from the cold phase to the individual molecular clouds, reaching up to 45% inside the molecular gas, at densities $\gtrsim 400$ cm$^{-3}$. However, at the 100 pc scale, the contribution of gravity-driven gas on the linewidth is negligible. We conclude that the bulk of the gas motions assembling the clouds in our simulation are caused by the supernova-driven supersonic turbulence, which results in locally convergent flows, with a small contribution from the stellar gravitational potential.

Our understanding of strong gravity near supermassive compact objects has recently improved thanks to the measurements made by the Event Horizon Telescope (EHT). We use here the M87* shadow size to infer constraints on the physical charges of a large variety of nonrotating or rotating black holes. For example, we show that the quality of the measurements is already sufficient to rule out that M87* is a highly charged dilaton black hole. Similarly, when considering black holes with two physical and independent charges, we are able to exclude considerable regions of the space of parameters for the doubly-charged dilaton and the Sen black holes.

The 2017 observing campaign of the Event Horizon Telescope (EHT) delivered the first very long baseline interferometry (VLBI) images at the observing frequency of 230 GHz, leading to a number of unique studies on black holes and relativistic jets from active galactic nuclei (AGN). In total, eighteen sources were observed: the main science targets, Sgr A* and M87 along with various calibrators. We investigated the morphology of the sixteen AGN in the EHT 2017 data set, focusing on the properties of the VLBI cores: size, flux density, and brightness temperature. We studied their dependence on the observing frequency in order to compare it with the Blandford-Königl (BK) jet model. We modeled the source structure of seven AGN in the EHT 2017 data set using linearly polarized circular Gaussian components and collected results for the other nine AGN from dedicated EHT publications, complemented by lower frequency data in the 2-86 GHz range. Then, we studied the dependences of the VLBI core flux density, size, and brightness temperature on the frequency measured in the AGN host frame. We compared the observations with the BK jet model and estimated the magnetic field strength dependence on the distance from the central black hole. Our results indicate a deviation from the standard BK model, particularly in the decrease of the brightness temperature with the observing frequency. Either bulk acceleration of the jet material, energy transfer from the magnetic field to the particles, or both are required to explain the observations.

Low- to intermediate-mass ($\sim$0.8$-$8 M$_\odot$) evolved stars contribute significantly to the chemical enrichment of the interstellar medium in the local Universe, making accurate mass-return estimates in their final stages crucial. The Nearby Evolved Stars Survey (NESS) is a large multi-telescope project targeting a volume-limited sample of $\sim$850 stars within 3 kpc in order to derive the dust and gas return rates in the Solar Neighbourhood, and to constrain the physics underlying these processes. We present an initial analysis of the CO-line observations, including detection statistics, carbon isotopic ratios, initial mass-loss rates, and gas-to-dust ratios. We describe a new data reduction pipeline to analyse the available NESS CO data from the JCMT, measuring line parameters and calculating empirical gas mass-loss rates. We present the first release of the available data on 485 sources, one of the largest homogeneous samples of CO data to date. Comparison with a large literature sample finds that high mass-loss rate and especially carbon-rich sources are over-represented in literature, while NESS is probing significantly more sources at low mass-loss rates, detecting 59 sources in CO for the first time and providing useful upper limits. CO line detection rates are 81% for the CO (2--1) line and 75% for CO (3--2). The majority (82%) of detected lines conform to the expected soft parabola shape, while eleven sources show a double wind. Calculated mass-loss rates show power-law relations with both the dust-production rates and expansion velocities up to $\sim 5 \times 10^{-6}$~\msunyr. Median gas-to-dust ratios of 250 and 680 are found for oxygen-rich and carbon-rich sources, respectively. Our analysis of CO observations in this first data release highlights the importance of our volume-limited approach in characterizing the local AGB population as a whole.

The study of hot corinos in Solar-like protostars has been so far mostly limited to the Class 0 phase, hampering our understanding of their origin and evolution. In addition, recent evidence suggests that planet formation starts already during Class I phase, which, therefore, represents a crucial step in the future planetary system chemical composition. Hence, the study of hot corinos in Class I protostars has become of paramount importance. Here we report the discovery of a hot corino towards the prototypical Class I protostar L1551 IRS5, obtained within the ALMA Large Program FAUST. We detected several lines from methanol and its isopotologues ($^{13}$CH$_{\rm 3}$OH and CH$_{\rm 2}$DOH), methyl formate and ethanol. Lines are bright toward the north component of the IRS5 binary system, and a possible second hot corino may be associated with the south component. The methanol lines non-LTE analysis constrains the gas temperature ($\sim$100 K), density ($\geq$1.5$\times$10$^{8}$ cm$^{-3}$), and emitting size ($\sim$10 au in radius). All CH$_{\rm 3}$OH and $^{13}$CH$_{\rm 3}$OH lines are optically thick, preventing a reliable measure of the deuteration. The methyl formate and ethanol relative abundances are compatible with those measured in Class 0 hot corinos. Thus, based on the present work, little chemical evolution from Class 0 to I hot corinos occurs.

The rate of long-duration gamma-ray bursts (GRBs) from isolated Pop III stars is not well known, as it depends on our poor understanding of their initial mass function (IMF), rotation rates, stellar evolution, and mass loss. A sub-population of massive ($M_{\rm ZAMS}\gtrsim20M_\odot$) Pop III stars is expected to suffer core-collapse and launch a relativistic jet that would power a GRB. In the collapsar scenario, a key requirement is that the pre-supernova star imparts sufficient angular momentum to the remnant black hole to form an accretion disc and launch a relativistic jet, which demands rapid initial rotation of the progenitor star and suppression of line-driven mass loss during its chemically homogeneous evolution. Here we explore a grid of stellar evolution models of Pop III stars with masses $20\leq M_{\rm ZAMS}/M_\odot \leq 100$, which are initially rotating with surface angular velocities $0.6\leq \Omega_0/\Omega_{\rm crit}\leq 0.9$, where centrifugally-driven mass loss ensues for $\Omega>\Omega_{\rm crit}$. Realistic accretion and jet propagation models are used to derive the initial black hole masses and spins, and jet breakout times for these stars. The GRB production efficiency is obtained over a phase space comprising progenitor initial mass, rotation, and wind efficiency. For modest wind efficiency of $\eta_{\rm wind}=0.45-0.35$, the Pop III GRB production efficiency is $\eta_{\rm GRB}\sim10^{-5}-3\times10^{-4}\,M_\odot^{-1}$, respectively, for a top-heavy IMF. This yields an observable all-sky equivalent rate of $\sim2-40\,{\rm yr}^{-1}$ by \textit{Swift}, with 75\% of the GRBs located at $z\lesssim8$. If the actual observed rate is much lower, then this would imply $\eta_{\rm wind}>0.45$, which leads to significant loss of mass and angular momentum that renders isolated Pop III stars incapable of producing GRBs and favours a binary scenario instead.

A fundamental prediction of the LambdaCDM cosmology is the hierarchical build-up of structure and therefore the successive merging of galaxies into more massive ones. As one can only observe galaxies at one specific time in cosmic history, this merger history remains in principle unobservable. By using the TNG100 simulation of the IllustrisTNG project, we show that it is possible to infer the unobservable stellar assembly and merger history of central galaxies from their observable properties by using machine learning techniques. In particular, in this first paper of ERGO-ML (Extracting Reality from Galaxy Observables with Machine Learning), we choose a set of 7 observable integral properties of galaxies (i.e. total stellar mass, redshift, color, stellar size, morphology, metallicity, and age) to infer, from those, the stellar ex-situ fraction, the average merger lookback times and mass ratios, and the lookback time and stellar mass of the last major merger. To do so, we use and compare a Multilayer Perceptron Neural Network and a conditional Invertible Neural Network (cINN): thanks to the latter we are also able to infer the posterior distribution for these parameters and hence estimate the uncertainties in the predictions. We find that the stellar ex-situ fraction and the time of the last major merger are well determined by the selected set of observables, that the mass-weighted merger mass ratio is unconstrained, and that, beyond stellar mass, stellar morphology and stellar age are the most informative properties. Finally, we show that the cINN recovers the remaining unexplained scatter and secondary cross-correlations. Our tools can be applied to large galaxy surveys in order to infer unobservable properties of galaxies' past, enabling empirical studies of galaxy evolution enriched by cosmological simulations.

The Very Large Array Sky Survey (VLASS) is a synoptic, all-sky radio sky survey with a unique combination of high angular resolution ($\approx$2.5"), sensitivity (a 1$\sigma$ goal of 70 $\mu$Jy/beam in the coadded data), full linear Stokes polarimetry, time domain coverage, and wide bandwidth (2-4 GHz). The first observations began in September 2017, and observing for the survey will finish in 2024. VLASS will use approximately 5500 hours of time on the Karl G. Jansky Very Large Array (VLA) to cover the whole sky visible to the VLA (Declination $>-40^{\circ}$), a total of 33,885 deg$^2$. The data will be taken in three epochs to allow the discovery of variable and transient radio sources. The survey is designed to engage radio astronomy experts, multi-wavelength astronomers, and citizen scientists alike. By utilizing an "on the fly" interferometry mode, the observing overheads are much reduced compared to a conventional pointed survey. In this paper, we present the science case and observational strategy for the survey, and also results from early survey observations.

Spectroscopic observations have shown for decades that the Wolf-Rayet (WR) phenomenon is ubiquitous among stars with different initial masses. Although much effort to understand the winds from massive WR stars has been presented in the literature, not much has been done for such type of stars in the low-mass range. Here we present an attempt to understand the winds from [WR]-type stars using results from spectral analyses with the full non-LTE stellar atmosphere code PoWR. These results are put into context with the properties of massive WR stars. We found that WC+[WC] stars and WO+[WO] stars create independent sequences in the mass-loss rate ($\dot{M}$) and modified wind momentum ($D_\mathrm{mom}$) versus luminosity ($L$) diagrams. Our analysis indicates that even when the winds of WR and [WR] stars become optically thin, there is no breakdown of the general mass-loss trend, contrary to the observed ``weak wind phenomenon'' in OB stars. We report that all WR-type stars studied here broadly define single sequences in the wind efficiency ($\eta$) versus transformed mass-loss rate ($\dot{M}_\mathrm{t}$), the $\dot{M}_\mathrm{t}$-$T_\mathrm{eff}$ diagram, and the $(L, T_\mathrm{eff}, \dot{M})$ space, which suggest these to be fundamental properties of the WR phenomenon (regardless of the mass range), at least for WR-type stars of the O and C sequences. Our analytical estimations could drive computations of future stellar evolution models for WR-type stars.

Radio used as a star formation rate (SFR) tracer presents enormous advantages by being unaffected by dust and radio sources being pinpointed at the sub-arc-second level. The interpretation of the low frequency 1.4 GHz luminosity is hampered by the difficulty in modeling the cosmic ray paths in the interstellar medium, and their interactions with the magnetic field. In this work, we compare the SFR derived from radio observations, and the ones derived from spectral energy distribution (SED) modeling. We aim at better understand the behavior of the SFR radio tracer, with a specific emphasis on the link with star-formation histories. We used the SED modeling code Code Investigating GALaxy Emission, CIGALE, with a non-parametric star formation history model (SFH) and fit the data over the wavelength range from the ultraviolet (UV) up to the mid-infrared (mid-IR). We interpret the difference between radio and SED-based SFR tracers in the light of recent gradients in the derived SFH. To validate the robustness of the results, we checked for any remaining active galaxy nuclei (AGN) contribution and tested the impact of our SFH modeling approach. Approximately 27% our galaxies present a radio SFR (SFR$_{\rm radio}$) at least ten times larger than the instantaneous SFR from SED-fitting (SFR$_{\rm SED}$). This trend affects primarily the galaxies that show a declining SFH activity over the last 300 Myr. Both SFR indicators converge toward a consistent value, when the SFHs are averaged over a period larger than 150 Myr to derive SFR$_{\rm SED}$. Although the radio at low frequency 1.4 GHz is a good tracer of the star formation activity of galaxies with constant or increasing SFH, our results indicate that this is not the case for galaxies that are quenching. Our analysis suggests that the star formation time sensitivity of the radio low frequency could be longer than 150 Myr.

Two-dimensional hydrodynamical simulations are presented from the formation up to the late evolution of planetary nebula, for 6 different stellar models from 1 to 5 Mo. The hydrodynamical models use stellar evolution calculations as inner boundary conditions and updated values for the number of ionizing photons. Special emphasis is placed on the formation of neutral spikes, as recently observed by the James Webb Space Telescope. The results indicate that neutral spikes can be detected either at the formation of planetary nebulae or in their decline. In the first case, the temporal window decreases with the mass of the model, ranging from 3,000 years in the 1 Mo case to 0 for 5 Mo. In the second case, only the 1.5, 2.0, and 2.5 Mo cases allow us to detect the neutral spikes for most of the remaining time.

Brown dwarfs are the bridge between low-mass stars and giant planets. One way of shedding light on their dominant formation mechanism is to study them at the earliest stages of their evolution, when they are deeply embedded in their parental clouds. Several works have identified pre- and proto-brown dwarfs candidates using different observational approaches. The aim of this work is to create a database with all the objects classified as very young substellar candidates in the litearature in order to study them in an homogeneous way. We have gathered all the information about very young substellar candidates available in the literature until 2020. We have retrieved their published photometry from the optical to the centimeter regime, and we have written our own codes to derive their bolometric temperatures and luminosities, and their internal luminosities. We have also populated the database with other parameters extracted from the literature, like e.g. the envelope masses, their detection in some molecular species, and presence of outflows. The result of our search is the SUCANES database, containing 174 objects classified as potential very young substellar candidates in the literature. We present an analysis of the main properties of the retrieved objects. Since we have updated the distances to several star forming regions, this has allowed us to reject some candidates based on their internal luminosities. We have also discussed the derived physical parameters and envelope masses for the best substellar candidates isolated in SUCANES. As an example of a scientific exploitation of this database, we present a feasibility study for the detection of radiojets with upcoming facilities: the ngVLA and the SKA interferometers. The SUCANES database is accessible through a Graphical User Interface and it is open to any potential user.

At distances from the active galaxy nucleus (AGN) where the ambient temperature falls below ~1500-1800 K, dust is able to survive. It is thus possible to have a large dusty structure present which surrounds the AGN. This is the first of two papers aiming at comparing six dusty torus models with available SEDs, namely Fritz et al. (2006), Nenkova et al. (2008B), Hoenig & Kishimoto (2010), Siebenmorgen et al. (2015), Stalevski et al. (2016), and Hoenig & Kishimoto (2017). In this first paper we use synthetic spectra to explore the discrimination between these models and under which circumstances they allow to restrict the torus parameters, while our second paper analyzes the best model to describe the mid-infrared spectroscopic data. We have produced synthetic spectra from current instruments: GTC/CanariCam and Spitzer /IRS and future JWST/MIRI and JWST/NIRSpec instruments. We find that for a reasonable brightness (F12um > 100mJy) we can actually distinguish among models except for the two pair of parent models. We show that these models can be distinguished based on the continuum slopes and the strength of the silicate features. Moreover, their parameters can be constrained within 15% of error, irrespective of the instrument used, for all the models but Hoenig & Kishimoto (2017). However, the parameter estimates are ruined when more than 50% of circumnuclear contributors are included. Therefore, future high spatial resolution spectra as those expected from JWST will provide enough coverage and spatial resolution to tackle this topic.

The spatial properties of small star-clusters suggest that they may originate from a fragmentation cascade of the cloud for which there might be traces up to a few dozen of kAU. Our goal is to investigate the multi-scale spatial structure of gas clumps, to probe the existence of a hierarchical cascade and to evaluate its possible link with star production in terms of multiplicity. From the Herschel emission maps of NGC 2264, clumps are extracted using getsf software at each of their associated spatial resolution, respectively [8.4, 13.5, 18.2, 24.9, 36.3]". Using the spatial distribution of these clumps and the class 0/I Young Stellar Object (YSO) from Spitzer data, we develop a graph-theoretic analysis to represent the multi-scale structure of the cloud as a connected network. From this network, we derive three classes of multi-scale structure in NGC 2264 depending on the number of nodes produced at the deepest level: hierarchical, linear and isolated. The structure class is strongly correlated with the column density $N_{\rm H_2}$ since the hierarchical ones dominate the regions whose N$_{\rm H_2} > 6 \times 10^{22}$cm$^{-2}$. Although the latter are in minority, they contain half of the class 0/I YSOs proving that they are highly efficient in producing stars. We define a novel statistical metric, the fractality coefficient F that measure the fractal index describing the scale-free process of the cascade. For NGC 2264, we estimate F = 1.45$\pm$0.12. However, a single fractal index fails to fully describe a scale-free process since the hierarchical cascade starts at a 13 kAU characteristic spatial scale. Our novel methodology allows us to correlate YSOs with their multi-scale gaseous environment. This hierarchical cascade that drives efficient star formation is suspected to be both hierarchical and rooted by the larger-scale gas environment up to 13 kAU.

We present an analysis of the optical properties of three Ultra Luminous X-ray (ULX) sources identified in NGC 925. We use Integral field unit data from the George Mitchel spectrograph in the context of the Metal-THINGS survey. The optical properties for ULX-1 and ULX-3 are presented, while the spaxel associated with ULX-2 had a low S/N, which prevented its analysis. We also report the kinematics and dimensions of the optical nebula associated with each ULX using ancillary data from the PUMA Fabry-Perot spectrograph. A BPT analysis demonstrates that most spaxels in NGC 925 are dominated by star-forming regions, including those associated with ULX-1 and ULX-3. Using the resolved gas-phase metallicities, a negative metallicity gradient is found, consistent with previous results for spiral galaxies, while the ionization parameter tends to increase radially throughout the galaxy. Interestingly, ULX-1 shows a very low gas metallicity for its galactocentric distance, identified by two independent methods, while exhibiting a typical ionization. We find that such low gas metallicity is best explained in the context of the high-mass X-ray binary population, where the low-metallicity environment favours active Roche lobe overflows that can drive much higher accretion rates. An alternative scenario invoking accretion of a low-mass galaxy is not supported by the data in this region. Finally, ULX-3 shows both a high metallicity and ionization parameter, which is consistent with the progenitor being a highly-accreting neutron star within an evolved stellar population region.

It is well known that galaxies falling into clusters can experience gas stripping due to ram-pressure by the intra-cluster medium (ICM). The most spectacular examples are galaxies with extended tails of optically-bright stripped material known as "jellyfish". We use the first large homogeneous compilation of jellyfish galaxies in clusters from the WINGS and OmegaWINGS surveys, and follow-up MUSE observations from the GASP MUSE programme to investigate the orbital histories of jellyfish galaxies in clusters and reconstruct their stripping history through position vs. velocity phase- space diagrams. We construct analytic models to define the regions in phase-space where ram-pressure stripping is at play. We then study the distribution of cluster galaxies in phase-space and find that jellyfish galaxies have on average higher peculiar velocities (and higher cluster velocity dispersion) than the overall population of cluster galaxies at all clustercentric radii, which is indicative of recent infall into the cluster and radial orbits. In particular, the jellyfish galaxies with the longest gas tails reside very near the cluster cores (in projection) and are moving at very high speeds, which coincides with the conditions of the most intense ram-pressure. We conclude that many of the jellyfish galaxies seen in clusters likely formed via fast (~1- 2 Gyr), incremental, outside-in ram-pressure stripping during first infall into the cluster in highly radial orbits.

The ALMA interferometer, with its unprecedented combination of high-sensitivity and high-angular resolution, allows for (sub-)mm wavelength mapping of protostellar systems at Solar System scales. Astrochemistry has benefited from imaging interstellar complex organic molecules in these jet-disk systems. Here we report the first detection of methanol (CH3OH) and methyl formate (HCOOCH3) emission towards the triple protostellar system VLA1623-2417 A1+A2+B, obtained in the context of the ALMA Large Program FAUST. Compact methanol emission is detected in lines from Eu = 45 K up to 61 K and 537 K towards components A1 and B, respectively. LVG analysis of the CH3OH lines towards VLA1623-2417 B indicates a size of 0.11-0.34 arcsec (14-45 au), a column density N(CH3OH) = 10^16-10^17 cm-2, kinetic temperature > 170 K, and volume density > 10^8 cm-3. An LTE approach is used for VLA1623-2417 A1, given the limited Eu range, and yields Trot < 135 K. The methanol emission around both VLA1623-2417 A1 and B shows velocity gradients along the main axis of each disk. Although the axial geometry of the two disks is similar, the observed velocity gradients are reversed. The CH3OH spectra from B shows two broad (4-5 km s-1) peaks, which are red- and blue-shifted by about 6-7 km s-1 from the systemic velocity. Assuming a chemically enriched ring within the accretion disk, close to the centrifugal barrier, its radius is calculated to be 33 au. The methanol spectra towards A1 are somewhat narrower (about 4 km s-1), implying a radius of 12-24 au.

Using observations with e-MERLIN and the VLA, together with archival data from ALMA, we obtain high-resolution radio images of two binary YSOs: L1551 IRS 5 and L1551 NE, covering a wide range of frequencies from 5 - 336 GHz, and resolving emission from the radio jet on scales of only ~15 au. By comparing these observations to those from a previous epoch, it is shown that there is a high degree of variability in the free-free emission from the jets of these sources. In particular, the northern component of L1551 IRS 5 shows a remarkable decline in flux density of a factor of ~5, suggesting that the free-free emission of this source has almost disappeared. By fitting the spectra of the sources, the ionised mass-loss rates of the jets are derived and it is shown that there is significant variability of up to a factor of ~6 on timescales of ~20 years. Using radiative transfer modelling, we also obtained a model image for the jet of the southern component of L1551 IRS 5 to help study the inner region of the ionised high-density jet. The findings favour the X-wind model launched from a very small innermost region.

In this work we present chemical abundances of heavy elements (Z$>$28) for a homogeneous sample of 1059 stars from HARPS planet search program. We also derive ages using parallaxes from Hipparcos and Gaia DR1 to compare the results. We study the [X/Fe] ratios for different populations and compare them with models of Galactic chemical evolution. We find that thick disk stars are chemically disjunt for Zn and Eu. Moreover, the high-alpha metal-rich population presents an interesting behaviour, with clear overabundances of Cu and Zn and lower abundances of Y and Ba with respect to thin disk stars. Several abundance ratios present a significant correlation with age for chemically separated thin disk stars (regardless of their metallicity) but thick disk stars do not present that behaviour. Moreover, at supersolar metallicities the trends with age tend to be weaker for several elements.