The discovery of the most metal-poor stream, C-19, provides us with a fossil record of a stellar structure born very soon after the Big Bang. In this work, we search for new C-19 members over the whole sky by combining two complementary stream-searching algorithms, STREAMFINDER and StarGO,, and utilizing low-metallicity star samples from the Pristine survey as well as Gaia BP/RP spectro-photometric catalogues. We confirm twelve new members, spread over more than 100$^\circ$, using velocity and metallicity information from a set of spectroscopic follow-up programs that targeted a quasi-complete sample of our bright candidates ($G \lesssim 16.0$). From the updated set of stream members, we confirm that the stream is wide, with a stream width of $\sim200$ pc, and dynamically hot, with a derived velocity dispersion of $10.9^{+2.1}_{-1.5}$ km/s. The tension remains between these quantities and a purely baryonic scenario in which the relatively low-mass stream (even updated to a few $10^4M_{\odot}$) stems from a globular cluster progenitor, as suggested by its chemical abundances. Some heating mechanism, such as preheating of the cluster in its own dark matter halo or through interactions with halo sub-structures appears necessary to explain the tension. The impact of binaries on the measured dispersion also remains unknown. Detailed elemental abundances of more stream members as well as multi-epoch radial velocities from spectroscopic observations are therefore crucial to fully understand the nature and past history of the most metal-poor stream of the Milky Way.

We analyze the merger and assembly histories of Milky Way (MW) and Andromeda (M31)-like galaxies to quantify how, and how often, disk galaxies of this mass can survive recent major mergers (stellar mass ratio $\ge$ 1:4). For this, we use the cosmological magneto-hydrodynamical simulation TNG50 and identify 198 analog galaxies, selected based on their $z=0$ stellar mass ($10^{10.5-11.2} {\rm M_{\odot}}$), disky stellar morphology and local environment. Firstly, major mergers are common: 85 per cent (168) of MW/M31-like galaxies in TNG50 have undergone at least one major merger across their lifetime. In fact, 31 galaxies (16 per cent) have undergone a recent major merger, i.e. in the last 5 Gyr. The gas available during the merger suffices to either induce starbursts at pericentric passages or to sustain prolonged star formation after coalescence: in roughly half of the cases, the pre-existing stellar disk is destroyed because of the merger but reforms thanks to star formation. Moreover, higher merger mass ratios are more likely to destroy the stellar disks. In comparison to those with more ancient massive mergers, MW/M31-like galaxies with recent major mergers have, on average, somewhat thicker stellar disks, more massive and somewhat shallower stellar haloes, larger stellar ex-situ mass fractions, but similarly massive kinematically-defined bulges. All this is qualitatively consistent with the different observed properties of the Galaxy and Andromeda and with the constraints on their most recent major mergers, 8-11 and ~2 Gyr ago, respectively. According to contemporary cosmological simulations, a recent quiet merger history is not a pre-requisite for obtaining a relatively-thin stellar disk at $z=0$.

We present the full public release of all data from the Illustris simulation project. Illustris is a suite of large volume, cosmological hydrodynamical simulations run with the moving-mesh code Arepo and including a comprehensive set of physical models critical for following the formation and evolution of galaxies across cosmic time. Each simulates a volume of (106.5 Mpc)^3 and self-consistently evolves five different types of resolution elements from a starting redshift of z=127 to the present day, z=0. These components are: dark matter particles, gas cells, passive gas tracers, stars and stellar wind particles, and supermassive black holes. This data release includes the snapshots at all 136 available redshifts, halo and subhalo catalogs at each snapshot, and two distinct merger trees. Six primary realizations of the Illustris volume are released, including the flagship Illustris-1 run. These include three resolution levels with the fiducial "full" baryonic physics model, and a dark matter only analog for each. In addition, we provide four distinct, high time resolution, smaller volume "subboxes". The total data volume is ~265 TB, including ~800 full volume snapshots and ~30,000 subbox snapshots. We describe the released data products as well as tools we have developed for their analysis. All data may be directly downloaded in its native HDF5 format. Additionally, we release a comprehensive, web-based API which allows programmatic access to search and data processing tasks. In both cases we provide example scripts and a getting-started guide in several languages: currently, IDL, Python, and Matlab. This paper addresses scientific issues relevant for the interpretation of the simulations, serves as a pointer to published and on-line documentation of the project, describes planned future additional data releases, and discusses technical aspects of the release.

Magnetic fields are known to be dynamically important in the interstellar medium of our own Galaxy, and they are ubiquitously observed in diffuse gas in the halos of galaxies and galaxy clusters. Yet, magnetic fields have typically been neglected in studies of the formation of galaxies, leaving their global influence on galaxy formation largely unclear. We extend our MHD implementation in the moving-mesh code Arepo to cosmological problems which include radiative cooling and the formation of stars. In particular, we replace our previously employed divergence cleaning approach with a Powell 8-wave scheme, which turns out to be significantly more stable, even in very dynamic environments. We verify the improved accuracy through simulations of the MRI in accretion disks, that reproduce its correct linear growth rate. Using this new MHD code, we simulate the formation of isolated disk galaxies similar to the Milky Way using idealized initial conditions with and without magnetic fields. We find that the magnetic field is quickly amplified in the initial starburst and the differential rotation of the forming disk until it eventually saturates when it becomes comparable to the thermal pressure. The additional pressure component leads to a lower star formation rate at late times compared to simulations without magnetic fields, and induces changes in the spiral arm structures of the gas disk. In addition, we observe highly magnetized fountain-like outflows from the disk. These results are robust with numerical resolution and are largely independent of the initial magnetic seed field assumed in the initial conditions, as the amplification process is rapid and self-regulated. Our findings suggest an important influence of magnetic fields on galaxy formation and evolution, cautioning against their neglect in theoretical models of structure formation.

We implemented sink particles in the adaptive mesh refinement (AMR) hydrodynamics code FLASH. Sink particles are created in regions of local gravitational collapse, and their trajectories and accretion can be followed over many dynamical times. We perform a series of tests including the time integration of circular and elliptical orbits, the collapse of a Bonnor-Ebert sphere and a rotating, fragmenting cloud core. We compare the collapse of a highly unstable singular isothermal sphere to the theory by Shu (1977), and show that the sink particle accretion rate is in excellent agreement with the theoretical prediction. To model eccentric orbits and close encounters of sink particles accurately, we show that a very small timestep is often required, for which we implemented subcycling of the N-body system. We emphasize that a sole density threshold for sink particle creation is insufficient in supersonic flows, if the density threshold is below the opacity limit. In that case, the density can exceed the threshold in strong shocks that do not necessarily lead to local collapse. Additional checks for bound state, gravitational potential minimum, Jeans instability and converging flows are absolutely necessary for a meaningful creation of sink particles. We apply our new sink particle module for FLASH to the formation of a stellar cluster, and compare to a smoothed particle hydrodynamics (SPH) code with sink particles. Our comparison shows encouraging agreement of gas properties, indicated by column density distributions and radial profiles, and of sink particle formation times and positions. We find excellent agreement in the number of sink particles formed, and in their accretion and mass distributions.

We study the population of supermassive black holes (SMBHs) and their effects on massive central galaxies in the IllustrisTNG cosmological hydrodynamical simulations of galaxy formation. The employed model for SMBH growth and feedback assumes a two-mode scenario in which the feedback from active galactic nuclei occurs through a kinetic, comparatively efficient mode at low accretion rates relative to the Eddington limit, and in the form of a thermal, less efficient mode at high accretion rates. We show that the quenching of massive central galaxies happens coincidently with kinetic-mode feedback, consistent with the notion that active supermassive black cause the low specific star formation rates observed in massive galaxies. However, major galaxy mergers are not responsible for initiating most of the quenching events in our model. Up to black hole masses of about $10^{8.5}\,{\rm M}_\odot$, the dominant growth channel for SMBHs is in the thermal mode. Higher mass black holes stay mainly in the kinetic mode and gas accretion is self-regulated via their feedback, which causes their Eddington ratios to drop, with SMBH mergers becoming the main channel for residual mass growth. As a consequence, the quasar luminosity function is dominated by rapidly accreting, moderately massive black holes in the thermal mode. We show that the associated growth history of SMBHs produces a low-redshift quasar luminosity function and a redshift zero black hole mass-stellar bulge mass relation in good agreement with observations, whereas the simulation tends to over-predict the high-redshift quasar luminosity function.

Context. The dynamical evolution of open clusters in the tidal field of the Milky Way and the feeding of the disc field star population depend strongly on the initial conditions at the time of gas removal. Detailed dynamical models tailored to individual clusters help us understand the role of open clusters in the Galactic disc evolution. Aims. We present a detailed dynamical model of Praesepe, which reproduces the mass profile, the stellar mass function, and the mass segregation observed with the help of Gaia EDR3 data. Based on this model, we investigate the kinematic properties of the tidal tail stars in detail. Methods. We used direct N-body simulations along the eccentric orbit of Praesepe in the tidal field of the Milky Way, where each particle represents one star. The initial mass and size of the cluster, the dynamical state, and the initial mass function were adapted to reach the best-fitting model. Based on this model and a comparison model on a circular orbit, we analysed the stars in the tidal tails in terms of density, angular momentum, and orbit shapes. Results. Praesepe can be well reproduced by a cluster model with concentrated star formation in a supervirial state after instantaneous gas expulsion, adopting a global star formation efficiency of 17%. About 75% of the initially 7500 MSol are lost in the violent relaxation phase, and the observed mass segregation can be understood by two-body relaxation. We find that the self-gravity of the tail stars is the dominant force altering the angular momentum of the tail stars. For a typical star, the total change after escaping is about 1.6 kpc km/s. This corresponds to an offset in guiding radius of 7 pc, where tail stars contribute up to 70% to the alteration. The total radial shift of the orbit of the cluster in the Galactic plane can exceed 50 pc. This effect is not a result of the eccentricity of the orbit.

Knowledge of abundance ratios as functions of metallicity can lead to insights on the origin and evolution of our Galaxy and its stellar populations. We aim to trace the chemical evolution of the neutron-capture elements Sr, Zr, La, Ce, Nd, Sm, and Eu in the Milky Way stellar disk to constrain the formation sites of these elements as well as to probe the evolution of the Galactic thin and thick disks. Using spectra of high resolution and high signal-to-noise we determine Sr, Zr, La, Ce, Nd, Sm, and Eu abundances for a sample of 593 F and G dwarf stars in the Solar neighbourhood. We present abundance results for Sr, Zr, La, Ce, Nd, Sm and Eu. We find that Nd, Sm, and Eu show trends similar to what is observed for the alpha-elements when compared to [Fe/H]. [Sr/Fe] and [Zr/Fe] show decreasing abundance ratios for increasing metallicity, reaching sub-solar values at super-solar metallicities. [La/Fe] and [Ce/Fe] do not show any clear trend with metallicity. The rapid neutron-capture process is active early in the Galaxy, mainly in type II supernovae from stars in the mass range 8-10 M_sun. Eu is almost completely produced by r-process but Nd and Sm show similar trends to Eu even if their s-process component is higher. Sr and Zr show significant enrichment at low metallicity that requires extra r-process production, that probably is different from the classical r-process. Finally, La and Ce are mainly produced via s-process from AGB stars in mass range 2-4 M_sun. The trend of [X/Fe] with age found could be explained by considering that the decrease in [X/Fe] for the thick disk stars can be due to the decrease of type II supernovae with time meaning a reduced enrichment of r-process elements in the interstellar medium. In the thin disk the trends are flatter that probably is due to that the main production from s-process is balanced by Fe production from type Ia supernovae.

Fundamental questions about the physics responsible for fragmenting molecular parsec-scale clumps into cores of ~1000 au are still open, that only a statistically significant investigation with ALMA is able to address: what are the dominant agents that determine the core demographics, mass, and spatial distribution as a function of the physical properties of the hosting clumps, their evolutionary stage and the different Galactic environments in which they reside? To what extent extent is fragmentation driven by clumps dynamics or mass transport in filaments? With ALMAGAL we observed the 1.38 mm continuum and lines toward more than 1000 dense clumps in our Galaxy, with M>500M_sun, surface density > 0.1 g/cm2 and d<7.5 kpc. The ACA and two 12-m array setups were used to deliver a minimum resolution of ~1000 au over the entire sample distance range. The sample covers all evolutionary stages from infrared dark clouds (IRDCs) to HII regions from the tip of the Galactic bar to the outskirts of the Galaxy. The spectral setup includes several molecular lines to trace the multiscale physics and dynamics of gas, notably CH3CN, H2CO, SiO, CH3OH, DCN, HC3N, SO etc. We present an initial overview of the observations and the early science product and results, with a first characterization of the morphological properties of the continuum emission. We use "perimeter-versus-area" and convex hull-versus-area metrics to classify the different morphologies. More extended and morphologically complex shapes are found toward clumps that are relatively more evolved and have higher surface densities.

Feedback from Young Massive Clusters (YMCs) is an important driver of galaxy evolution. In the first few Myr, mechanical feedback is dominated by collective effects of the massive stellar winds in the YMC. The mass-loss rates and terminal wind velocities of these stars change by orders of magnitude over pre-SN timescales as the massive stars evolve, and mass-loss rates of Cool Supergiant (CSG) stars in particular are uncertain by a factor $\sim~20$ or more. In this work we perform a first study of the time evolution of average cluster wind velocity $\bar{V}_{\mathrm{cl}}$ as a function of stellar metallicity $Z$, assuming single star evolution. We also check the validity of assuming Wolf-Rayet stars dominate the feedback effects of a YMC, as often done when interpreting X-ray and $\gamma$-ray observations, and test how sensitive $\bar{V}_{\mathrm{cl}}$ is to current uncertainties in mass-loss rates. We use pySTARBURST99 to calculate integrated properties of YMCs for $Z$ in the range of $0.0004-0.02$, representing a range of environments from IZw18 to the Galactic Centre. We find that $\bar{V}_{\mathrm{cl}}$ drops off rapidly for sub-LMC $Z$, and we recommend a value of $500-1000\,~\textrm{km~s}^{-1}$ be used in this regime. We show accounting only for WR stars can overestimate $\bar{V}_{\mathrm{cl}}$ by $500-2000\,~\textrm{km~s}^{-1}$ at $Z \geq Z_\text{LMC}$. We also find that different RSG mass-loss assumptions can change the inferred $\bar{V}_{\mathrm{cl}}$ by $\sim1000\,~\textrm{km~s}^{-1}$, highlighting the need for improved observational constraints for RSGs in YMCs.

(Abridged) We study the galaxy mass-size relation in CARLA spectroscopically confirmed clusters at $1.410.5$ are systematically $\gtrsim 0.2-0.3~{\rm dex}$ larger than field ETGs. The passive ETG average size evolution is slower at $1

In the last decade, the photospheric solar metallicity as determined from spectroscopy experienced a remarkable downward revision. Part of this effect can be attributed to an improvement of atomic data and the inclusion of NLTE computations, but also the use of hydrodynamical model atmospheres seemed to play a role. This "decrease" with time of the metallicity of the solar photosphere increased the disagreement with the results from helioseismology. With a CO5BOLD 3D model of the solar atmosphere, the CIFIST team at the Paris Observatory re-determined the photospheric solar abundances of several elements, among them C, N, and O. The spectroscopic abundances are obtained by fitting the equivalent width and/or the profile of observed spectral lines with synthetic spectra computed from the 3D model atmosphere. We conclude that the effects of granular fluctuations depend on the characteristics of the individual lines, but are found to be relevant only in a few particular cases. 3D effects are not reponsible for the systematic lowering of the solar abundances in recent years. The solar metallicity resulting from this analysis is Z=0.0153, Z/X=0.0209.

Context. Teegarden's Star is the brightest and one of the nearest ultra-cool dwarfs in the solar neighbourhood. For its late spectral type (M7.0V), the star shows relatively little activity and is a prime target for near-infrared radial velocity surveys such as CARMENES. Aims. As part of the CARMENES search for exoplanets around M dwarfs, we obtained more than 200 radial-velocity measurements of Teegarden's Star and analysed them for planetary signals. Methods. We find periodic variability in the radial velocities of Teegarden's Star. We also studied photometric measurements to rule out stellar brightness variations mimicking planetary signals. Results. We find evidence for two planet candidates, each with $1.1M_\oplus$ minimum mass, orbiting at periods of 4.91 and 11.4 d, respectively. No evidence for planetary transits could be found in archival and follow-up photometry. Small photometric variability is suggestive of slow rotation and old age. Conclusions. The two planets are among the lowest-mass planets discovered so far, and they are the first Earth-mass planets around an ultra-cool dwarf for which the masses have been determined using radial velocities.

Context. B hypergiants (BHGs) are important for understanding high-mass stellar evolution. While they are in a similar parameter space of B supergiants (BSGs), some BHGs are known to be luminous blue variables (LBVs). Their spectra with absorption and emission features resemble Of/WNh stars. Yet, their wind physics and their evolutionary connections are not clear. Aims. In this study, we aim to understand (i) the atmospheric and wind structure, (ii) the wind-launching and -driving mechanisms, and (iii) the spectrum formation of early-type BHGs. As an observational prototype, we use zet1 Sco (B1.5Ia+). Methods. Using the atmosphere code PoWRhd, we calculated the first hydrodynamically consistent models at the BHG domain. They give insights into the radiative driving of the calculated wind regimes and enable us to study the influence of clumping and X-rays on the resulting wind structure. Results. Our consistent model reproduces the main spectral features of zet1 Sco. The obtained mass-loss rate is higher than that of similar spectral type BSGs. However, the wind optical depth of BHGs is way below unity, making them less of a transition type. To reproduce zet1 Sco's spectrum, we needed low clumping with subsonic onset. The wind has a shallow-gradient velocity profile, deviating from the beta law, and is mainly driven by Fe III opacity. Conclusions. Our study suggests that despite more mass loss, early-type Galactic BHGs have winds relatively similar to BSGs. Their winds are not thick enough to characterize them as "transition-type" stars, unlike Of/WNh, implying that emission features arise more easily in cooler than in hotter stars. The spectral BHG appearance is likely connected to atmospheric inhomogeneities below the sonic point. To reach an appearance similar to LBVs, BHGs need to be either closer to the Eddington limit or have higher wind clumping than inferred for zeta1 Sco.

We report the analysis of a planetary microlensing event AT2021uey. The event was observed outside the Galactic bulge and was alerted by both space- (Gaia) and ground-based (ZTF and ASAS-SN) surveys. From the observed data, we find that the lens system is located at a distance of 1 kpc and comprises an M-dwarf host star of about half a solar mass, orbited by a Jupiter-like planet beyond the snowline. The source star could be a metal-poor giant located in the halo according to the spectral analyses and modelling. Hence, AT2021uey is a unique example of the binary-lens event outside the bulge that is offered by a disc-halo lens-source combination.

Context. Galaxies fly inside galaxy clusters and ram pressure by the ICM can remove a large amount of the ISM from the galaxy, and deposit the gas in the ICM. The ISM decoupled from the host galaxy leaves a long trail following the moving galaxy. Such long trails track the galaxy motion and can be detected with sensitive data in Halpha. Aims. We study the Halpha tail trailing NGC 4569 in the Virgo cluster. Methods. The initial discovery was made with the deep Halpha imaging data with CFHT, from the VESTIGE project. The follow-up spectroscopic observations were made with APO/DIS, MMT/Binospec and CFHT/SITELLE. Results. Besides the known 80 kpc Halpha tail downstream of NGC 4569, the deep Halpha imaging data allow the Halpha tail detected to at least 230 kpc from the galaxy. More importantly, the Halpha clumps implied from the imaging data are confirmed with the spectroscopic data. The Halpha clumps show a smooth radial velocity gradient across about 1300 km/s, eventually reaching the velocity of the cluster. We build a simple model to explain the deceleration of stripped clumps and constrain the age to about 0.9 Gyr. Conclusions. This discovery, for the first time, demonstrates the full deceleration process of the stripped ISM. This discovery also showcases the potential with wide-field Halpha survey on galaxy clusters to discover intracluster optical emission-line clouds originated from cluster galaxies. These clouds provide kinematic tracers to the infall history of cluster galaxies and the turbulence in the ICM. They are also excellent multi-phase objects to study the relevant important physical processes.

JWST observations, when combined with HST data, promise to improve age estimates of star clusters in nearby spiral galaxies. However, feedback from young cluster stars pushes out the natal gas and dust, making cluster formation and evolution a challenge to model. Here, we use JWST + HST observations of the nearby spiral galaxy NGC 628 to produce spectral energy distribution (SED) templates of compact star clusters spanning 275 nm through 21 {\mu}m. These preliminary SEDs capture the cluster stars and associated gas and dust within radii of 0.12" to 0.67" (corresponding to 6 to 33 pc at the distance of NGC 628). One important finding is that the SEDs of 1, 2, 3, and 4 Myr clusters can be differentiated in the infrared. Another is that in 80-90% of the cases we study, the PAH and H_alpha emission track one another, with the dust responsible for the 3.3 {\mu}m PAH emission largely removed by 4 Myr, consistent with pre-supernova stellar feedback acting quickly on the surrounding gas and dust. Nearly-embedded cluster candidates have infrared SEDs which are quite similar to optically visible 1 to 3 Myr clusters. In nearly all cases we find there is a young star cluster within a few tenths of an arcsec (10 - 30 pc) of the nearly embedded cluster, suggesting the formation of the cluster was triggered by its presence. The resulting age estimates from the empirical templates are compatible both with dynamical estimates based on CO superbubble expansion velocities, and the TODDLERS models which track spherical evolution of homogeneous gas clouds around young stellar clusters.

Modern machine learning will allow for simulation-based inference from reionization-era 21cm observations at the Square Kilometre Array. Our framework combines a convolutional summary network and a conditional invertible network through a physics-inspired latent representation. It allows for an efficient and extremely fast determination of the posteriors of astrophysical and cosmological parameters, jointly with well-calibrated and on average unbiased summaries. The sensitivity to non-Gaussian information makes our method a promising alternative to the established power spectra.

We analyze the IllustrisTNG simulations to study the mass, volume fraction and phase distribution of gaseous baryons embedded in the knots, filaments, sheets and voids of the Cosmic Web from redshift $z=8$ to redshift $z=0$. We find that filaments host more star-forming gas than knots, and that filaments also have a higher relative mass fraction of gas in this phase than knots. We also show that the cool, diffuse Intergalactic Medium (IGM; T&lt;10^5 \, {\rm K}, n_{\rm H}&lt;10^{-4}(1+z) \, {\rm cm^{-3}}) and the Warm-Hot Intergalactic Medium (WHIM; $ 10^5 \, {\rm K}