Search

One of the fundamental goals of modern astrophysics is to estimate the physical parameters of galaxies. We present a hierarchical Bayesian model to compute age maps from images in the H α line (taken with Taurus tunable filter, TTF), ultraviolet band (GALEX far UV, FUV), and infrared bands (Spitzer 24, 70, and 160 |$\mu$| m). We present the burst ages for young stellar populations in a sample of nearby and nearly face-on galaxies. The H α to FUV flux ratio is a good relative indicator of the very recent star formation history (SFH). As a nascent star-forming region evolves, the H α line emission declines earlier than the UV continuum, leading to a decrease in the H α/FUV ratio. Using star-forming galaxy models, sampled with a probabilistic formalism, and allowing for a variable fraction of ionizing photons in the clusters, we obtain the corresponding theoretical ratio H α/FUV to compare with our observed flux ratios, and thus to estimate the ages of the observed regions. We take into account the mean uncertainties and the interrelationships between parameters when computing H α/FUV. We propose a Bayesian hierarchical model where a joint probability distribution is defined to determine the parameters (age, metallicity, IMF) from the observed data (the observed flux ratios H α/FUV). The joint distribution of the parameters is described through independent and identically distributed (i.i.d.) random variables generated through MCMC (Markov Chain Monte Carlo) techniques. [ABSTRACT FROM AUTHOR]

The H2O Southern Galactic Plane Survey (HOPS) has mapped 100 deg² of the Galactic plane for water masers and thermal molecular line emission using the 22 m Mopra telescope. We describe the automated spectral-line fitting pipelines used to determine the properties of emission detected in HOPS data cubes, and use these to derive the physical and kinematic properties of gas in the survey. A combination of the angular resolution, sensitivity, velocity resolution and high critical density of lines targeted make the HOPS data cubes ideally suited to finding precursor clouds to the most massive and dense stellar clusters in the Galaxy. We compile a list of the most massive HOPS ammonia regions and investigate whether any may be young massive cluster progenitor gas clouds. HOPS is also ideally suited to trace the flows of dense gas in the Galactic Centre. We find the kinematic structure of gas within the inner 500 pc of the Galaxy is consistent with recent predictions for the dynamical evolution of gas flows in the centre of the Milky Way. We confirm a recent finding that the dense gas in the inner 100 pc has an oscillatory kinematic structure with characteristic length-scale of 20 pc, and also identify similar oscillatory kinematic structure in the gas at radii larger than 100 pc. Finally, we make all of the above fits and the remaining HOPS data cubes across the 100 deg² of the survey available to the community. [ABSTRACT FROM AUTHOR]

We study the distribution and evolution of the stellar mass and the star formation rate (SFR) of the brightest group galaxies (BGGs) over 0.04 < z < 1.3 using a large sample of 407 X-ray galaxy groups selected from the COSMOS, AEGIS, and XMM-LSS fields. We compare our results with predictions from the semi-analytic models based on the Millennium simulation. In contrast to model predictions, we find that, as the Universe evolves, the stellar mass distribution evolves towards a normal distribution. This distribution tends to skew to low-mass BGGs at all redshifts implying the presence of a star-forming population of the BGGs with MS ∼ 1010.5 M⊙ which results in the shape of the stellar mass distribution deviating from a normal distribution. In agreement with the models and previous studies, we find that the mean stellar mass of BGGs grows with time by a factor of ∼2 between z = 1.3 and z = 0.1, however, the significant growth occurs above z = 0.4. The BGGs are not entirely a dormant population of galaxies, as low-mass BGGs in low-mass haloes are more active in forming stars than the BGGs in more massive haloes, over the same redshift range. We find that the average SFR of the BGGs evolves steeply with redshift and fraction of the passive BGGs increases as a function of increasing stellar mass and halo mass. Finally, we show that the specific SFR of the BGGs within haloes with M200 ≤ 1013.4 M⊙ decreases with increasing halo mass at z < 0.4. [ABSTRACT FROM AUTHOR]

This document is a correction to a paper titled "Determining star-formation rates in active galactic nuclei hosts via stellar population synthesis." The correction addresses an inconsistency in the value of |$H_0$| used in the analysis. The authors have now adopted a uniform value of |$H_0=73.0$| km s|$^{-1}$| Mpc|$^{-1}$| for both star-formation rates (SFR|$_\\star$| and SFR|$_{Gas}$|). As a result, some figures and the transformation equation have been updated. The conclusions of the paper remain unchanged. The document provides the corrected equations and figures for reference. [Extracted from the article]

This document is a correction to a previously published paper titled "Binary black hole mergers from population III stars: uncertainties from star formation and binary star properties." The authors discovered a minor bug in the code used to produce the results shown in various figures. The correction provides the updated figures and lists the changes made to the published version. The authors state that the bug did not significantly affect the results and conclusions of the paper. The main data presented in the paper are publicly available on Zenodo and GitLab. [Extracted from the article]

Tidal streams and stellar shells are naturally formed in galaxy interactions and mergers. The Giant Stellar Stream (GSS), the North-East (NE), and Western (W) stellar shelves observed in Andromeda Galaxy (M31) are examples of these structures and were formed through the merger of M31 and a satellite galaxy. Recent observational papers have provided strong evidence that the shells and GSS originate from a single progenitor. In this paper, we investigate the formation of these two stellar shelves and the detailed nature of their relationship to the GSS. We present numerical simulations of tidal disruption of a satellite galaxy assuming that it is a progenitor of the GSS and the shell system. We represent the progenitor as a dwarf spheroidal galaxy with the stellar mass of 109 M⊙ and evolve its merger with M31 for 3 Gyr to reproduce the chemodynamical properties of the NE and W shelves. We find that an initial metallicity of the progenitor with a negative radial gradient of ΔFeH  = −0.3 ± 0.2 successfully reproduces observed metallicities of the NE, W shelves, and the GSS, showing that all these structures can originate from the same merger event. [ABSTRACT FROM AUTHOR]

In our previous study (Tsukamoto et al. 2023b, PASJ, 75, 835), we investigated the formation and early evolution of protoplanetary disks with 3D non-ideal magnetohydrodynamics simulations considering dust growth, and found that the modified equations of the conventional steady accretion disk model that consider magnetic braking, dust growth, and ambipolar diffusion reproduce the disk structure (such as density and vertical magnetic field) obtained from simulations very well. In this paper, as a sequel to our previous study, we analytically investigate the structure and evolution of protoplanetary disks corresponding to Class 0/I young stellar objects using the modified steady accretion disk model combining an analytical model of envelope accretion. We estimate that the disk radius is several astronomical units at the disk formation epoch and increases to several hundred astronomical units at the end of the accretion phase. The disk mass is estimated to be |$0.01 \lesssim M_{\rm disk} \lesssim 0.1 \, M_\odot$| for a disk with a radius of several tens of astronomical units and a mass accretion rate of |$\dot{M}_{\rm disk} \sim 10^{-6} \, M_\odot \,\, {\rm yr^{-1}}$|⁠. These estimates seems to be consistent with recent observations. We also found that, with typical disk ionization rates (ζ ≳ 10−19 s−1) and a moderate mass accretion rate (⁠|$\dot{M}_{\rm disk}\gtrsim 10^{-8} \, M_\odot \,\, {\rm yr^{-1}}$|⁠), magnetorotational instability is suppressed in the disk because of low plasma β and efficient ambipolar diffusion. We argue that the radial profile of specific angular momentum (or rotational velocity) at the disk outer edge should be continuously connected to that of the envelope if the disk evolves by magnetic braking, and should be discontinuous if the disk evolves by an internal angular momentum transport process such as gravitational instability or magnetorotational instability. Future detailed observations of the specific angular momentum profile around the disk outer edge are important for understanding the angular momentum transport mechanism of protoplanetary disks. [ABSTRACT FROM AUTHOR]

Our ability to trace the star-forming molecular gas is important to our understanding of the Universe. We can trace this gas using CO emission, converting the observed CO intensity into the H |$_2$| gas mass of the region using the CO-to-H |$_2$| conversion factor (⁠|$X_{\rm{{\small CO}}}$|⁠). In this paper, we use simulations to study the conversion factor and the molecular gas within galaxies. We analysed a suite of simulations of isolated disc galaxies, ranging from dwarfs to Milky Way-mass galaxies, that were run using the fire-2 subgrid models coupled to the chimes non-equilibrium chemistry solver. We use the non-equilibrium abundances from the simulations, and we also compare to results using abundances assuming equilibrium, which we calculate from the simulation in post-processing. Our non-equilibrium simulations are able to reproduce the relation between CO and H |$_2$| column densities, and the relation between |$X_{\rm{{\small CO}}}$| and metallicity, seen within observations of the Milky Way. We also compare to the xCOLD GASS survey, and find agreement with their data to our predicted CO luminosities at fixed star formation rate. We also find the multivariate function used by xCOLD GASS overpredicts the H |$_2$| mass for our simulations, motivating us to suggest an alternative multivariate function of our fitting, though we caution that this fitting is uncertain due to the limited range of galaxy conditions covered by our simulations. We also find that the non-equilibrium chemistry has little effect on the conversion factor (<5 per cent) for our high-mass galaxies, though still affects the H |$_2$| mass and |$L_{\rm{{\small CO}}}$| by |$\approx$| 25  per cent. [ABSTRACT FROM AUTHOR]

It is well known that stars are subject to radial migration, i.e. over time, they move away from their birth location. This dynamical process tends to mix different stellar populations and hence hinders the determination of the true chemical evolution of a galaxy (e.g. metallicity gradients). One way to account for radial migration is to infer stellar birth radii for individual stars. Many attempts to do so have been performed over the last few years, but are limited to the Milky Way, as computing the birth position of stars requires precise measurements of stellar metallicity and age for individual stars that cover large Galactic radii. Fortunately, recent and future surveys will provide numerous opportunities for inferring birth radii for external galaxies such as the LMC. In this paper, we investigate the possibility of doing so using the NIHAO cosmological zoom-in simulations. We find that it is theoretically possible to infer birth radii with a ∼25 per cent median uncertainty for individual stars in galaxies with i) orderliness of the orbits, |$\langle v_\phi \rangle /\sigma _{v} > 2 $|⁠ , ii) a dark matter halo mass greater or equal to approximately the LMC mass (∼2 × 1011  |${\rm M}_\odot$|⁠), and iii) after the average azimuthal velocity of the stellar disc reaches ∼70 per cent of its maximum. From our analysis, we conclude that it is possible and useful to infer birth radii for the LMC and other external galaxies that satisfy the above criteria. [ABSTRACT FROM AUTHOR]

In order to develop a complete theory of star formation, one essentially needs to know two things: what collapses and how long it takes. This is the second paper in a series, where we query how long a parcel of gas takes to collapse and the process it undergoes. We embed pseudo-Lagrangian tracer particles in simulations of collapsing molecular clouds, identify the particles that end in dense knots, and then examine the collapse history of the gas. We find a nearly universal behaviour of cruise-then-collapse, wherein a core stays at intermediate densities for a significant fraction of its life before finally collapsing. We identify time immediately before each core collapses, |$t_{\rm {sing}}$|⁠ , and examine how it transitions to high density. We find that the time to collapse is uniformly distributed between |$0.25 t_{\rm {ff}}$| and the end of the simulation at |$\sim\!\! 1 t_{\rm {ff}}$|⁠ , and that the duration of collapse is universally short, |$\Delta t \sim 0.1 t_{\rm {ff}}$|⁠ , where |$t_{\rm {ff}}$| is the free-fall time at the mean density. We describe the collapse in three stages: collection, hardening, and singularity. Collection sweeps low-density gas into moderate density. Hardening brings kinetic and gravitational energies into quasi-equipartition. Singularity is the free-fall collapse, forming an envelope in rough energy balance and central overdensity in |$\sim\!\! 0.1 t_{\rm {ff}}$|⁠. [ABSTRACT FROM AUTHOR]

To advance our understanding of the evolution of the interstellar medium (ISM) of our Galaxy, numerical models of Milky Way (MW) type galaxies are widely used. However, most models only vaguely resemble the MW (e.g. in total mass), and often use imposed analytic potentials (which cannot evolve dynamically). This poses a problem in asserting their applicability for the interpretation of observations of our own Galaxy. The goal of this work is to identify a numerical model that is not only an MW-type galaxy, but one that can mimic some of the main observed structures of our Galaxy, using dynamically evolving potentials, so that it can be used as a base model to study the ISM cycle in a galaxy like our own. This paper introduces a suite of 15 MW-type galaxy models developed using the arepo numerical code, that are compared to Galactic observations of |$^{12}$| CO and H  i emission via longitude–velocity plots, from where we extract and compare the skeletons of major galactic features and the terminal gas velocities. We found that our best-fitting model to the overall structure, also reproduces some of the more specific observed features of the MW, including a bar with a pattern speed of |$30.0 \pm 0.2$| km s |$^{-1}$|  kpc |$^{-1}$|⁠ , and a bar half-length of |$3.2 \pm 0.8$|  kpc. Our model shows large streaming motions around spiral arms, and strong radial motions well beyond the inner bar. This model highlights the complex motions of a dynamic MW-type galaxy and has the potential to offer valuable insight into how our Galaxy regulates the ISM and star formation. [ABSTRACT FROM AUTHOR]

In recent years, the study of the Milky Way has significantly advanced due to extensive spectroscopic surveys of its stars, complemented by astroseismic and astrometric data. However, it remains disjoint from recent advancements in understanding the physics of the Galactic interstellar medium (ISM). This paper introduces a new model for the chemical evolution of the Milky Way that can be constrained on stellar data, because it combines a state-of-the-art ISM model with a Milky Way stellar disc model. Utilizing a data set of red clump stars from APOGEE, known for their precise ages and metallicities, we concentrate on the last 6 billion years – a period marked by Milky Way's secular evolution. We examine the oxygen abundance in the low- |$\alpha$| disc stars relative to their ages and birth radii, validating or constraining critical ISM parameters that remain largely unexplored in extragalactic observations. The models that successfully reproduce the radius–metallicity distribution and the age–metallicity distribution of stars without violating existing ISM observations indicate a need for modest differential oxygen enrichment in Galactic outflows, meaning that the oxygen abundance of outflows is higher than the local ISM abundance, irrespective of outflow mass loading. The models also suggest somewhat elevated ISM gas velocity dispersion levels over the past 6 billion years compared to galaxies of similar mass. The extra turbulence necessary could result from energy from gas accretion onto the Galaxy, supernovae clustering in the ISM, or increased star formation efficiency per freefall time. This work provides a novel approach to constraining the Galactic ISM and outflows, leveraging the detailed insights available from contemporary Milky Way surveys. [ABSTRACT FROM AUTHOR]

It is crucial for a deeper understanding of the formation and evolution of galaxies in the Universe to study stellar mass (M *) and star formation rate (SFR). Traditionally, astronomers infer the properties of galaxies from spectra, which are highly informative, but expensive and hard to be obtained. Fortunately, modern sky surveys obtained a vast amount of high-spatial-resolution photometric images. The photometric images are obtained relatively economically than spectra, and it is very helpful for related studies if M * and SFR can be estimated from photometric images. Therefore, this paper conducted some preliminary researches and explorations on this regard. We constructed a deep learning model named Galaxy Efficient Network (GalEffNet) for estimating integrated M * and specific star formation rate (sSFR) from Dark Energy Spectroscopic Instrument galaxy images. The GalEffNet primarily consists of a general feature extraction module and a parameter feature extractor. The research results indicate that the proposed GalEffNet exhibits good performance in estimating M * and sSFR, with σ reaching 0.218 and 0.410 dex. To further assess the robustness of the network, prediction uncertainty was performed. The results show that our model maintains good consistency within a reasonable bias range. We also compared the performance of various network architectures and further tested the proposed scheme using image sets with various resolutions and wavelength bands. Furthermore, we conducted applicability analysis on galaxies of various sizes, redshifts, and morphological types. The results indicate that our model performs well across galaxies with various characteristics and indicate its potentials of broad applicability. [ABSTRACT FROM AUTHOR]

Fast radio bursts (FRBs) are extremely powerful sources of radio waves observed at cosmological distances. We use a sophisticated model of FRB observations – presented in detail in a companion paper – to fit FRB population parameters using large samples of FRBs detected by ASKAP and Parkes, including seven sources with confirmed host galaxies. Our fitted parameters demonstrate that the FRB population evolves with redshift in a manner consistent with, or faster than, the star formation rate (SFR), ruling out a non-evolving population at better than 98 per cent CL (depending on modelling uncertainties). Our estimated maximum FRB energy is |$\log _{10} E_{\rm max} [{\rm erg}] = 41.70_{-0.06}^{+0.53}$| (68 per cent CL) assuming a 1 GHz emission bandwidth, with slope of the cumulative luminosity distribution |$\gamma =-1.09_{-0.10}^{+0.14}$|⁠. We find a log-mean host DM contribution of |$129_{-48}^{+66}$|  pc cm−3 on top of a typical local (interstellar medium and halo) contribution of ∼80 pc cm−3, which is higher than most literature values. These results are insensitive to assumptions of the FRB spectral index, and are consistent with the model of FRBs arising as the high-energy limit of magnetar bursts, but allow for FRB progenitors that evolve faster than the SFR. [ABSTRACT FROM AUTHOR]

In this paper, we have conducted an investigation focused on a segment of the Spitzer mid-infrared bubble N59, specifically referred to as R1 within our study. Situated in the inner Galactic plane, this region stands out for its hosting of five 6.7 GHz methanol masers, as well as numerous compact |$\mathrm{H}\, \rm {{\small II}}$| regions, massive clumps, filaments, and prominent bright rims. As 6.7 GHz masers are closely linked to the initial phases of high-mass star formation, exploring regions that exhibit a high abundance of these maser detections provides an opportunity to investigate relatively young massive star-forming sites. To characterize the R1 region comprehensively, we utilize multiwavelength (archival) data from optical to radio wavelengths, together with 13CO and C18O data. Utilizing the Gaia DR3 data, we estimate the distance towards the bubble to be 4.66 ± 0.70 kpc. By combining near-infrared (NIR) and mid-infrared (MIR) data, we identify 12 Class I and 8 Class II sources within R1. Furthermore, spectral energy distribution (SED) analysis of selected sources reveals the presence of four embedded high-mass sources with masses ranging from 8.70 to 14.20 M⊙. We also identified several O and B-type stars from radio continuum analysis. Our molecular study uncovers two distinct molecular clouds in the region, which, although spatially close, occupy different regions in velocity space. We also find indications of a potential hub-filament system fostering star formation within the confines of R1. Finally, we propose that the feedback from the |$\mathrm{H}\, \rm {{\small II}}$| regions has led to the formation of prominent Bright Rimmed Clouds (BRCs) within our region of interest. [ABSTRACT FROM AUTHOR]

Thermal dust polarization is a powerful tool to probe magnetic fields (⁠|$\boldsymbol{B}$|⁠) and grain properties. However, a systematic study of the dependence of dust polarization on grain properties in protostellar environments is not yet available. In this paper, we post-process a non-ideal MHD simulation of a collapsing protostellar core with our updated POLARIS code to study in detail the effects of iron inclusions and grain growth on thermal dust polarization. We found that superparamagnetic (SPM) grains can produce high polarization degree of |$p \sim 10\!-\!40~{{\ \rm per\ cent}}$| beyond ∼500 au from the protostar because of their efficient alignment by magnetically enhanced radiative torque mechanism. The magnetic field turbulence in the envelope causes the decrease in p with increasing emission intensity I as p ∝ I α with the slope α ∼ −0.3. But within 500 au, SPM grains tend to have inefficient internal alignment and be aligned with |$\boldsymbol{B}$| by RATs only, producing lower |$p \sim 1~{{\ \rm per\ cent}}$| and a steeper slope of α ∼ −0.6. For paramagnetic (PM) grains, the alignment loss of grains above |$1\, {\mu \rm {m}}$| in the inner ∼200 au produces |$p \lt \lt 1~{{\ \rm per\ cent}}$| and the polarization hole with α ∼ −0.9. Grain growth can increase p in the envelope for SPM grains, but cause stronger depolarization for SPM grains in the inner ∼500 au and for PM grains in the entire protostellar core. Finally, we found the increase of polarization angle dispersion function S with iron inclusions and grain growth, implying the dependence of B-field strength measured using the David–Chandrasekhar–Fermi technique on grain alignment and grain properties. [ABSTRACT FROM AUTHOR]

In this paper, we discuss measurements of the stellar population and star-forming properties for 43 spectroscopically confirmed publicly available high-redshift z > 7 JWST galaxies in the JADES and CEERS observational programs. We carry out a thorough study investigating the relationship between spectroscopic features and photometrically derived ones, including from spectral energy distribution (SED) fitting of models, as well as morphological and structural properties. We find that the star formation rates (SFRs) measured from H β line emission are higher than those estimated from Bayesian SED fitting and ultraviolet (UV) luminosity, with ratios SFRH β/SFRUV ranging from ∼2 to 13. This is a sign that the star formation history is consistently rising given the time-scales of H β versus UV star formation probes. In addition, we investigate how well equivalent widths (EWs) of H β λ4861, [O  iii ] λ4959, and [O  iii ] λ5007 can be measured from photometry, finding that, on average, the EW derived from photometric excesses in filters is 30 per cent smaller than the direct spectroscopic measurement. We also discover that a stack of the line emitting galaxies shows a distinct morphology after subtracting imaging that contains only the continuum. This gives us a first view of the line or ionized gas emission from z > 7 galaxies, demonstrating that this material has a similar distribution, statistically, as the continuum. We also compare the derived SFRs and stellar masses for both parametric and non-parametric star formation histories, where we find that 35 per cent of our sample formed at least 30 per cent of their stellar mass in recent (<10 Myr) starburst events. [ABSTRACT FROM AUTHOR]

The fragmentation of gas to form stars in molecular clouds is intrinsically linked to the turbulence within them. These internal motions are set at the birth of the cloud and may vary with galactic environment and as the cloud evolves. In this paper, we introduce a new suite of 15 high-resolution 3D molecular cloud simulations using the moving mesh code arepo to investigate the role of different decaying turbulent modes (mixed, compressive, and solenoidal) and virial ratios on the evolution of a |$10^4\, \mathrm{M}_{\odot }$| molecular cloud. We find that diffuse regions maintain a strong relic of the initial turbulent mode, whereas the initial gravitational potential dominates dense regions. Solenoidal seeded models thus give rise to a diffuse cloud with filament-like morphology, and an excess of brown dwarf mass fragments. Compressive seeded models have an early onset of star-formation, centrally condensed morphologies and a higher accretion rate, along with overbound clouds. 3D filaments identified using disperse and analysed through a new python toolkit we develop and make publicly available with this work called fiesta , show no clear trend in lengths, masses and densities between initial turbulent modes. Overbound clouds, however, produce more filaments and thus have more mass in filaments. The hubs formed by converging filaments are found to favour star-formation, with surprisingly similar mass distributions independent of the number of filaments connecting the hub. [ABSTRACT FROM AUTHOR]

The primordial clouds in the mini-haloes hatch the first generation stars of the Universe, which play a crucial role in cosmic evolution. In this paper, we investigate how turbulence impacts the structure of primordial star-forming clouds. Previous cosmological simulations of the first star formation predicted a typical mass of around |$\mathrm{ 100 \, M_\odot }$|⁠. This conflicts with recent observations of extremely metal-poor stars, suggesting a lower mass scale of about |$\mathrm{25 \, M_\odot }$|⁠. The discrepancy may arise from unresolved turbulence in the star-forming cloud, driven by primordial gas accretion during mini-halo formation in the previous simulations. To quantitatively examine the turbulence effect on the primordial cloud formation, we employ the adaptive mesh refinement code Enzo to model the gas cloud with primordial composition, including artificially driven turbulence on the cloud scale and relevant gas physics. This artificially driven turbulence utilizes a stochastic forcing model to mimic the unresolved turbulence inside mini-haloes. Our results show that the turbulence with high Mach number and compressional mode effectively fragments the cloud into several clumps, each with dense cores of |$\mathrm{22.7 - 174.9 \, M_\odot }$| that undergo Jeans instability to form stars. Fragmentation caused by intense and compressive turbulence prevents a runaway collapse of the cloud. The self-bound clumps with smaller masses in the turbulent primordial clouds suggest a possible pathway to decrease the theoretical mass scale of the first stars, further reconciling the mass discrepancy between simulations and observations. [ABSTRACT FROM AUTHOR]

Observations of star-forming regions provide snapshots in time of the star formation process, and can be compared with simulation data to constrain the initial conditions of star formation. In order to make robust inferences, different metrics must be used to quantify the spatial and kinematic distributions of stars. In this paper, we assess the suitability of the INdex to Define Inherent Clustering And TEndencies (INDICATE) method as a diagnostic to infer the initial conditions of star-forming regions that subsequently undergo dynamical evolution. We use INDICATE to measure the degree of clustering in N -body simulations of the evolution of star-forming regions with different initial conditions. We find that the clustering of individual stars, as measured by INDICATE, becomes significantly higher in simulations with higher initial stellar densities, and is higher in subvirial star-forming regions where significant amounts of dynamical mixing have occurred. We then combine INDICATE with other methods that measure the mass segregation (ΛMSR), relative stellar surface density ratio (ΣLDR), and the morphology (Q -parameter) of star-forming regions, and show that the diagnostic capability of INDICATE increases when combined with these other metrics. [ABSTRACT FROM AUTHOR]

Young stellar clusters are predominantly the hub of star formation and hence, ideal to perform comprehensive studies over the least explored substellar regime. Various unanswered questions like the mass distribution in brown dwarf regime and the effect of diverse cluster environment on brown dwarf formation efficiency still plague the scientific community. The nearby young cluster, IC 1396 with its feedback-driven environment, is ideal to conduct such study. In this paper, we adopt a multiwavelength approach, using deep Subaru HSC along with other data sets and machine learning techniques to identify the cluster members complete down to ∼ 0.03 M⊙ in the central 22 arcmin area of IC 1396. We identify 458 cluster members including 62 brown dwarfs which are used to determine mass distribution in the region. We obtain a star-to-brown dwarf ratio of ∼ 6 for a stellar mass range 0.03–1 M⊙ in the studied cluster. The brown dwarf fraction is observed to increase across the cluster as radial distance from the central OB-stars increases. This study also compiles 15 young stellar clusters to check the variation of star-to-brown dwarf ratio relative to stellar density and ultraviolet (UV) flux ranging within 4–2500 stars pc−2 and 0.7–7.3  G 0, respectively. The brown dwarf fraction is observed to increase with stellar density but the results about the influence of incident UV flux are inconclusive within this range. This is the deepest study of IC 1396 as of yet and it will pave the way to understand various aspects of brown dwarfs using spectroscopic observations in future. [ABSTRACT FROM AUTHOR]

Stellar feedback, the energetic interaction between young stars and their birthplace, plays an important role in the star formation history of the Universe and the evolution of the interstellar medium. Correctly interpreting the observations of star-forming regions is essential to understand stellar feedback, but it is a non-trivial task due to the complexity of the feedback processes and degeneracy in observations. In our recent paper, we introduced a conditional invertible neural network (cINN) that predicts seven physical properties of star-forming regions from the luminosity of 12 optical emission lines as a novel method to analyse degenerate observations. We demonstrated that our network, trained on synthetic star-forming region models produced by the warpfield -emission predictor (warpfield-emp), could predict physical properties accurately and precisely. In this paper, we present a new updated version of the cINN that takes into account the observational uncertainties during network training. Our new network named Noise-Net reflects the influence of the uncertainty on the parameter prediction by using both emission-line luminosity and corresponding uncertainties as the necessary input information of the network. We examine the performance of the Noise-Net as a function of the uncertainty and compare it with the previous version of the cINN, which does not learn uncertainties during the training. We confirm that the Noise-Net outperforms the previous network for the typical observational uncertainty range and maintains high accuracy even when subject to large uncertainties. [ABSTRACT FROM AUTHOR]

Westerlund 1 (Wd 1) is one of the most relevant star clusters in the Milky Way to study massive star formation, although it is still poorly known. Here, we used photometric and spectroscopic data to model the eclipsing binary W36, showing that its spectral type is O6.5 III  +  O9.5 IV, hotter and more luminous than thought before. Its distance d W36 = 4.03 ± 0.25 kpc agrees, within the errors, with three recent Gaia-EDR3-based distances reported in Paper I, Beasor & Davies, and by Negueruela's group. However, they follow different approaches to fix the zero-points for red sources such as those in Wd 1, and to select the best approach, we used an accurate modelling of W36. The weighted mean distance of our parallax (Paper I) and binary distances results in d wd1 = 4.05 ± 0.20 kpc, with an unprecedented accuracy of 5 per cent. We adopted isochrones based on the Geneva code with supersolar abundances to infer the age of W36B as 6.4 ± 0.7 Myr. This object seems to be part of the prolific star formation burst represented by OB giants and supergiants that occurred at 7.1 ± 0.5 Myr ago, which coincides with the recently published PMS isochrone with age 7.2 Myr. Other BA-type luminous evolved stars and yellow hypergiants spread in the age range of 8–11 Myr. The four red supergiants discussed in paper I represent the oldest population of the cluster with an age of 10.7 ± 1 Myr. The multiple episodes of star formation in Wd 1 are reminiscent of that reported for the R136/30 Dor LMC cluster. [ABSTRACT FROM AUTHOR]

There has been a recent resurgence in hydroxyl (OH) megamaser research driven by Square Kilometre Array (SKA) precursor/pathfinder telescopes. This will continue in the lead-up to the SKA mid-frequency array, which will greatly expand our view of OH megamasers and their cosmic evolution over ≳80 per cent of the age of the Universe. This is expected to yield large scientific returns as OH megamasers trace galaxy mergers, extreme star formation, high molecular gas densities, and potentially binary/dual supermassive black hole systems. In this paper, we predict the distortion to the OH luminosity function that a magnification bias will inflict, and in turn, predict the distortion on the OH megamaser number counts as a function of redshift. We identify spectral flux density thresholds that will enable efficient lensed OH megamaser selection in large spectral line surveys with MeerKAT and SKA. The surface density of lensed galaxies that could be discovered in this way is a strong function of the redshift evolution of the OH megamaser luminosity function, with predictions as high as ∼1 lensed OH source per square degree at high redshifts (z ≳ 1) for anticipated SKA spectral line survey designs. This could enable efficient selection of some of the most highly obscured galaxies in the Universe. This high-redshift selection efficiency, in combination with the large survey speed of the SKA at ≲1 GHz frequencies and the high magnifications possible with compact OH emission regions (μOH ≫ 10), will enable a transformational view of OH in the Universe. [ABSTRACT FROM AUTHOR]

Filamentary flows towards the centre of molecular clouds have been recognized as a crucial process in the formation and evolution of stellar clusters. In this paper, we present a comprehensive observational study that investigates the gas properties and kinematics of the Giant Molecular Cloud G148.24+00.41 using the observations of CO (1-0) isotopologues. We find that the cloud is massive (105 M⊙) and is one of the most massive clouds of the outer Galaxy. We identified six likely velocity coherent filaments in the cloud having length, width, and mass in the range of 14–38 pc, 2.5–4.2 pc, and (1.3–6.9) × 103 M⊙, respectively. We find that the filaments are converging towards the central area of the cloud, and the longitudinal accretion flows along the filaments are in the range of ∼ 26–264 M⊙ Myr−1. The cloud has fragmented into seven clumps having mass in the range of ∼ 260–2100 M⊙ and average size around ∼ 1.4 pc, out of which the most massive clump is located at the hub of the filamentary structures, near the geometric centre of the cloud. Three filaments are found to be directly connected to the massive clump and transferring matter at a rate of ∼ 675 M⊙ Myr−1. The clump hosts a near-infrared cluster. Our results show that large-scale filamentary accretion flows towards the central region of the collapsing cloud is an important mechanism for supplying the matter necessary to form the central high-mass clump and subsequent stellar cluster. [ABSTRACT FROM AUTHOR]

While dwarf galaxies observed in the field are overwhelmingly star forming, dwarf galaxies in environments as dense or denser than the Milky Way are overwhelmingly quenched. In this paper, we explore quenching in the lower density environment of the Small-Magellanic-Cloud-mass galaxy NGC 3109 (M |$_* \sim 10^8 \, \text{M}_\odot$|⁠), which hosts two known dwarf satellite galaxies (Antlia and Antlia B), both of which are |${\rm H}\, \rm{\small I}$| deficient compared to similar galaxies in the field and have recently stopped forming stars. Using a new semi-analytic model in concert with the measured star formation histories and gas masses of the two dwarf satellite galaxies, we show that they could not have been quenched solely by direct ram pressure stripping of their interstellar media, as is common in denser environments. Instead, we find that separation of the satellites from pristine gas inflows, coupled with stellar-feedback-driven outflows from the satellites (jointly referred to as the starvation quenching model), can quench the satellites on time-scales consistent with their likely infall times into NGC 3109's halo. It is currently believed that starvation is caused by 'weak' ram pressure that prevents low-density, weakly bound gas from being accreted on to the dwarf satellite, but cannot directly remove the denser interstellar medium. This suggests that star-formation-driven outflows serve two purposes in quenching satellites in low-mass environments: outflows from the host form a low-density circumgalactic medium that cannot directly strip the interstellar media from its satellites, but is sufficient to remove loosely bound gaseous outflows from the dwarf satellites driven by their own star formation. [ABSTRACT FROM AUTHOR]

Models derived in 2009 to fit mid-infrared (8–24 micron) source counts from the IRAS, ISO, and Spitzer missions, provide an excellent fit to deep counts with JWST , demonstrating that the evolution of dusty star-forming galaxies is well understood. The evolution of dust in galaxies at high redshifts is discussed and a simple prescription is proposed to model this. This allows more realistic models for source-counts at submillimetre wavelength. A reasonable fit to 250, 500, 850, and 1100 micron counts is obtained. This paper therefore draws together the IRAS, ISO, Spitzer, Akari, Herschel , submillimetre ground-based, and JWST surveys into a single picture. [ABSTRACT FROM AUTHOR]

Multiphase galactic outflows, generated by supernova (SN) feedback, are likely to be more metal rich than the interstellar media from which they are driven due to incomplete mixing between SN ejecta and the ambient interstellar medium. This enrichment is important for shaping galactic metallicities and metallicity gradients, but measuring it quantitatively from simulations requires resolution high enough to resolve mass, momentum and energy exchanges between the different phases of the outflows. In this context, we present QED, which are simulations of outflows, driven by SN feedback, conducted using Quokka , a new GPU-optimized adaptive mesh refinement radiation-hydrodynamics code. This code allows us to reach combinations of resolution, simulation volume, and simulation duration larger than those that have previously been possible, and to resolve all gas phases from cold neutral medium, T ∼ 100 K, to hot ionized gas, T ≳ 107 K. In this, a first of a series of papers exploring generation and evolution of multiphase outflows from a wide range of galactic environments and star formation rates, we quantify the extent of selective metal loading in solar neighbourhood-like environments. We explain the selective metal loading, we find as a result of the transport of metals within and between phases, a phenomenon we can study owing to the parsec-scale resolution that our simulations achieve. We also quantify the sensitivity of metal loading studies to numerical resolution, and present convergence criteria for future studies. [ABSTRACT FROM AUTHOR]

The fraction of ionizing photons that escape their host galaxies to ionize hydrogen in the intergalactic medium (IGM) is a critical parameter in analyses of the reionization era. In this paper, we use the meraxes semi-analytic galaxy formation model to infer the mean ionizing photon escape fraction and its dependence on galaxy properties through joint modelling of the observed high redshift galaxy population and existing constraints on the reionization history. Using a Bayesian framework, and under the assumption that escape fraction is primarily related to halo mass, we find that the joint constraints of the ultraviolet luminosity function, cosmic microwave background optical depth, and the Ly α forest require an escape fraction of |$(18\pm 5)$| per cent for galaxies within haloes of M ≲ 109 M⊙ and |$(5\pm 2)$| per cent for more massive haloes. In terms of galaxy properties, this transition in escape fraction occurs at stellar masses of M ⋆ ∼ 107 M⊙, nearly independent of redshift. As a function of redshift, reionization is dominated by the smaller M ⋆ ≲ 107 M⊙ galaxies with high escape fractions at z ≳ 6 and by the larger M ⋆ ≳ 107 M⊙ galaxies with lower escape fractions at z ≲ 6. Galaxies with star formation rates of 10−2.5 M⊙yr−1 to 10−1.5 M⊙yr−1 provide the dominant source of ionizing photons throughout reionization. Our results are consistent with recent direct measurements of a |$\sim 5~{{\ \rm per\ cent}}$| escape fraction from massive galaxies at the end of reionization and support the picture of low mass galaxies being the dominant sources of ionizing photons during reionization. [ABSTRACT FROM AUTHOR]

This paper presents the third data release of the INvestigating Stellar Population In RElics (INSPIRE) project, comprising 52 ultra-compact massive galaxies (UCMGs) observed with the X-Shooter spectrograph. We measure integrated stellar velocity dispersion, [Mg/Fe] abundances, ages, and metallicities for all the INSPIRE objects. We thus infer star formation histories and confirm the existence of a degree of relicness (DoR), defined in terms of the fraction of stellar mass formed by z  = 2, the time at which a galaxy has assembled 75 per cent of its mass, and the final assembly time. Objects with a high DoR assembled their stellar mass at early epochs, while low-DoR objects show a non-negligible fraction of later formed populations and hence a spread in ages and metallicities. A higher DoR correlates with larger [Mg/Fe], supersolar metallicity, and larger velocity dispersion values. The 52 UMCGs span a large range of DoR from 0.83 to 0.06, with 38 of them having formed more than 75  per cent of their mass by z  = 2. Of these, nine are extreme relics (DoR>0.7), since they formed the totality (⁠|$\gt 99~{{\ \rm per\ cent}}$|⁠) of their stellar mass by redshift z  = 2. The remaining 14 UCMGs cannot be considered relics, as they are characterized by more extended star formation histories. With INSPIRE  we built the first sizeable sample of relics outside the local Universe, up to z ∼ 0.4, increasing the number of confirmed relics by a factor of >10, and opening up an important window to explain the mass assembly of massive galaxies in the high- z Universe. [ABSTRACT FROM AUTHOR]

During star formation, dust plays the crucial role of coupling gas to stellar radiation fields, allowing radiation feedback to influence gas fragmentation and thus the stellar initial mass function (IMF). Variations in dust abundance therefore provide a potential avenue by which variation in galaxy metallicity might affect the IMF. In this paper, we present a series of radiation-magnetohydrodynamic simulations in which we vary the metallicity and thus the dust abundance from 1 per cent of solar to 3× solar, spanning the range from the lowest metallicity dwarfs to the most metal-rich early-type galaxies (ETGs) found in the local Universe. We design the simulations to keep all dimensionless parameters constant so that the interaction between feedback and star-forming environments of varying surface density and metallicity is the only factor capable of breaking the symmetry between the simulations and modifying the IMF, allowing us to isolate and understand the effects of each environmental parameter cleanly. We find that shifts in the IMF with varying metallicity at a fixed surface density are smaller than the shifts in varying surface density at a fixed surface metallicity. We also find that metallicity-induced IMF variations are too small to explain the mass-to-light ratio shifts seen in the ETGs. We therefore conclude that metallicity variations are much less important than variations in surface density in driving changes in the IMF and that the latter rather than the former are most likely responsible for the IMF variations found in ETGs. [ABSTRACT FROM AUTHOR]

With the launch of JWST and other scheduled missions aimed at probing the distant universe, we are entering a new promising era for high- z astronomy. One of our main goals is the detection of the first population of stars (Population III or Pop III stars), and models suggest that Pop III star formation is allowed well into the Epoch of Reionization (EoR), rendering this an attainable achievement. In this paper, we focus on our chance of detecting massive Pop IIIs at the moment of their death as Pair-Instability Supernovae (PISNe). We estimate the probability of discovering PISNe during the EoR in galaxies with different stellar masses (7.5 ≤ Log(M ⋆/M⊙) ≤ 10.5) from six dustyGadget simulations of 50  h −1 cMpc per side. We further assess the expected number of PISNe in surveys with JWST/NIRCam and Roman/WFI. On average, less than one PISN is expected in all examined JWST fields at z ≃ 8 with Δ z  = 1, and O(1) PISN may be found in a ∼1 deg2 Roman field in the best-case scenario, although different assumptions on the Pop III IMF and/or Pop III star formation efficiency can decrease this number substantially. Including the contribution from unresolved low-mass haloes holds the potential for increased discoveries. JWST/NIRCam and Roman/WFI allow the detection of massive-progenitor (∼250 M⊙) PISNe throughout all the optimal F200W-F356W, F277W-F444W, and F158-F213 colours. PISNe are also pre-dominantly located at the outskirts of their hosting haloes, facilitating the disentangling of underlying stellar emission thanks to the spatial-resolution capabilities of the instruments. [ABSTRACT FROM AUTHOR]

The estimation of galaxy stellar masses depends on the assumed prior of the star-formation history (SFH) and spatial scale of the analysis (spatially resolved versus integrated scales). In this paper, we connect the prescription of the SFH in the spectral energy distribution (SED) fitting to spatially resolved scales (∼kpc) to shed light on the systematics involved when estimating stellar masses. Specifically, we fit the integrated photometry of ∼970 massive (log (M⋆/M⊙) = 9.8–11.5), intermediate redshift (z  = 0.5–2.0) galaxies with PROSPECTOR , assuming both exponentially declining tau model and flexible SFHs. We complement these fits with the results of spatially resolved SFH estimates obtained by pixel-by-pixel SED fitting, which assume tau models for individual pixels. These spatially resolved SFHs show a large diversity in shapes, which can largely be accounted for by the flexible SFHs with PROSPECTOR. The differences in the stellar masses from those two approaches are overall in good agreement (average difference of ∼0.07 dex). Contrarily, the simpler tau model SFHs typically miss the oldest episode of star formation, leading to an underestimation of the stellar mass by ∼0.3 dex. We further compare the derived global specific star-formation rate (sSFR), the mass-weighted stellar age (t50), and the star-formation time-scale (τSF) obtained from the different SFH approaches. We conclude that the spatially resolved scales within galaxies motivate a flexible SFH on global scales to account for the diversity of SFHs and counteract the effects of outshining of older stellar populations by younger ones. [ABSTRACT FROM AUTHOR]

Protostellar jets and outflows are signposts of active star formation. In H  ii regions, molecular tracers like CO only reveal embedded portions of the outflow. Outside the natal cloud, outflows are dissociated, ionized, and eventually completely ablated, leaving behind only the high-density jet core. Before this process is complete, there should be a phase where the outflow is partially molecular and partially ionized. In this paper, we capture the HH 900 outflow while this process is in action. New observations from the Enhanced Resolution Imager and Spectrograph/SPIFFIER near-infrared (IR) integral field unit spectrograph using the K-middle filter (λ = 2.06–2.34 μ m) reveal H2 emission from the dissociating outflow and Br-γ tracing its ionized skin. Both lines trace the wide-angle outflow morphology but H2 only extends ∼5000 au into the H  ii region while Br-γ extends the full length of the outflow (∼12 650 au), indicating rapid dissociation of the molecules. H2 has higher velocities further from the driving source, consistent with a jet-driven outflow. Diagnostic line ratios indicate that photoexcitation, not just shocks, contributes to the excitation in the outflow. We argue that HH 900 is the first clear example of an evaporating molecular outflow and predict that a large column of neutral material that may be detectable with Atacama Large Millimeter Array accompanies the dissociating molecules. Results from this study will help guide the interpretation of near-IR images of externally irradiated jets and outflows such as those obtained with the JWST in high-mass star-forming regions where these conditions may be common. [ABSTRACT FROM AUTHOR]

Arkenstone is a new model for multiphase, stellar feedback-driven galactic winds designed for inclusion in coarse resolution cosmological simulations. In this first paper of a series, we describe the features that allow Arkenstone to properly treat high specific energy wind components and demonstrate them using idealized non-cosmological simulations of a galaxy with a realistic circumgalactic medium (CGM), using the arepo code. Hot, fast gas phases with low mass loadings are predicted to dominate the energy content of multiphase outflows. In order to treat the huge dynamic range of spatial scales involved in cosmological galaxy formation at feasible computational expense, cosmological volume simulations typically employ a Lagrangian code or else use adaptive mesh refinement with a quasi-Lagrangian refinement strategy. However, it is difficult to inject a high specific energy wind in a Lagrangian scheme without incurring artificial burstiness. Additionally, the low densities inherent to this type of flow result in poor spatial resolution. Arkenstone addresses these issues with a novel scheme for coupling energy into the transition region between the interstellar medium (ISM) and the CGM, while also providing refinement at the base of the wind. Without our improvements, we show that poor spatial resolution near the sonic point of a hot, fast outflow leads to an underestimation of gas acceleration as the wind propagates. We explore the different mechanisms by which low and high specific energy winds can regulate the star formation rate of galaxies. In future work, we will demonstrate other aspects of the Arkenstone model. [ABSTRACT FROM AUTHOR]

All modern galaxy formation models employ stochastic elements in their sub-grid prescriptions to discretize continuous equations across the time domain. In this paper, we investigate how the stochastic nature of these models, notably star formation, black hole accretion, and their associated feedback, that act on small (< kpc) scales, can back-react on macroscopic galaxy properties (e.g. stellar mass and size) across long (> Gyr) time-scales. We find that the scatter in scaling relations predicted by the EAGLE model implemented in the swift code can be significantly impacted by random variability between re-simulations of the same object, even when galaxies are resolved by tens of thousands of particles. We then illustrate how re-simulations of the same object can be used to better understand the underlying model, by showing how correlations between galaxy stellar mass and black hole mass disappear at the highest black hole masses (M BH > 108 M⊙), indicating that the feedback cycle may be interrupted by external processes. We find that although properties that are collected cumulatively over many objects are relatively robust against random variability (e.g. the median of a scaling relation), the properties of individual galaxies (such as galaxy stellar mass) can vary by up to 25 per cent, even far into the well-resolved regime, driven by bursty physics (black hole feedback), and mergers between galaxies. We suggest that studies of individual objects within cosmological simulations be treated with caution, and that any studies aiming to closely investigate such objects must account for random variability within their results. [ABSTRACT FROM AUTHOR]

The role of massive (≥ 8 |$\, {\rm M}_{\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* ∼ 3 × 103 to 1 × 105 L⊙) sample of 38 massive star-forming regions observed with ALMA at 1.3 mm 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 and 15 per field). The flux density budget is split evenly (53 per cent–47 per cent) 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, λ 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.7 GHz 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 generalized approach used to distinguishing sources which are potentially star-forming from those which are not, utilizing interferometric visibility properties. [ABSTRACT FROM AUTHOR]

This document is a correction to a previous article by Bellstedt et al. regarding the cosmic star formation histories of bulges, discs, and spheroids in galaxies. The correction addresses an error in the dust parameters presented in tables 2, 3, and 4 of the original paper. Specifically, the |$\\tau _{\\rm birth}$| and |$\\tau _{\\rm screen}$| values were mistakenly swapped for each component. The corrected values are provided in updated tables 1, 2, and 3. These corrections do not impact any other aspect of the original paper and reflect the settings used to generate the original results. [Extracted from the article]

Many galaxies display clear bulges and discs, and understanding how these components form is a vital step towards understanding how the galaxy has evolved into what we see today. The BUDDI-MaNGA project aims to study galaxy evolution and morphological transformations through the star-formation histories of the bulges and discs. We have applied our buddi software to galaxies from the MaNGA Survey in the SDSS DR15 in order to isolate their bulge and disc spectra, from which we derived their stellar populations. To date, this work provides the largest sample of clean bulge and disc spectra extracted from IFU datacubes using the galaxies light profile information, and will form the basis for a series of papers aiming to answer open questions on how galaxies have formed and evolved, and the role of their individual structures. This paper presents an introduction to the project, including an overview of these fits, a characterization of the sample, and a series of tests on the fits to ensure reliability. [ABSTRACT FROM AUTHOR]

This is the second paper of the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) Dynamics and stellar Population (DynPop) series, which analyses the global stellar population, radial gradients, and non-parametric star-formation history of ∼10K galaxies from the MaNGA Survey final data release 17 and relates them with dynamical properties of galaxies. We confirm the correlation between the stellar population properties and the stellar velocity dispersion σe, but also find that younger galaxies are more metal-poor at fixed σe. Stellar age, metallicity, and mass-to-light ratio (M */ L) all decrease with increasing galaxy rotation, while their radial gradients become more negative (i.e. lower value at the outskirts). The exception is the slow rotators, which also appear to have significantly negative metallicity gradients, confirming the mass–metallicity gradient correlation. Massive disc galaxies in the green valley, on the |$(\sigma _{\rm e},\rm age)$| plane, show the most negative age and metallicity gradients, consistent with their old central bulges surrounded by young star-forming discs and metal-poor gas accretion. Galaxies with high σe, steep total mass-density slope, low dark matter fraction, high M */ L , and high metallicity have the highest star-formation rate at earlier times, and are currently quenched. We also discover a population of low-mass star-forming galaxies with low rotation but physically distinct from the massive slow rotators. A catalogue of these stellar population properties is provided publicly. [ABSTRACT FROM AUTHOR]

The feedback loop between the galaxies producing the background radiation field for reionization and their growth is crucial, particularly for low-mass haloes. Despite this, the vast majority of galaxy formation studies employ a spatially uniform, time-varying reionizing background, with the majority of reionization studies employing galaxy formation models only required to work at high redshift. This paper uses the well-studied TNG galaxy formation model, calibrated at low redshift, coupled to the arepo-rt code, to self-consistently solve the coupled problems of galaxy evolution and reionization, evaluating the impact of patchy (and slow) reionization on early galaxies. thesan-hr is an extension of the thesan project to higher resolution (a factor of 50 increase, with a baryonic mass of m b ≈ 104 M⊙), to additionally enable the study of 'mini-haloes' with virial temperatures T vir < 104 K. Comparing the self-consistent model to a uniform UV background, we show that galaxies in thesan-hr are predicted to be larger in physical extent (by a factor ∼2), less metal enriched (by ∼0.2 dex), and less abundant (by a factor ∼10 at M 1500 =   − 10) by z  = 5. We show that differences in star formation and enrichment patterns lead to significantly different predictions for star formation in low mass haloes, low-metallicity star formation, and even the occupation fraction of haloes. We posit that cosmological galaxy formation simulations aiming to study early galaxy formation (z ≳ 3) must employ a spatially inhomogeneous UV background to accurately reproduce galaxy properties. [ABSTRACT FROM AUTHOR]

Radiation feedback from massive population III stars may have given rise to low mass star formation from primordial or nearly primordial material. If early Universe low mass stars did form, some should remain locally as white dwarfs, sub-giants, or main sequence stars. In this paper, we present model calculations for the evolution of single 0.8–3.0 M⊙ stars with primordial metallicity from pre-main sequence to the white dwarf cooling track, and calculations for the evolution of single 4.0–7.0 M⊙ stars which conclude in the giant phase. One goal of this work is to identify potential observable markers for potential observed progenitors of first or nearly first stars. We uncover a number of seemingly peculiar evolutionary differences between that of population III low mass stars compared with younger higher Z stars, as well as compared to other primordial evolution models. We also present an initial–final mass relationship and identify the minimum mass of a single white dwarf that could have had a population III progenitor. [ABSTRACT FROM AUTHOR]

The mere existence of multiple stellar generations in Milky Way globular clusters indicates that each generation was unable to stop star formation, that instead persisted unimpeded for several million years. This evidence argues for an extended stage of star formation within a forming globular cluster, during which stellar feedback was substantially ineffective and the nascent globular cluster was able to accrete processed gas from its surrounding, and efficiently convert it into successive stellar generations. It has been argued that such delayed feedback results from core collapse in most massive stars failing to trigger an energetic supernova explosion, but rather leading directly to black hole formation. Thus, globular clusters offer a concrete phenomenological example for the lack of feedback in young starbursts, an option that has been widely advocated to account for the unexpected abundance of ultraviolet-luminous galaxies at z = 9–16, as revealed by JWST observations. The paper is meant to attract attention to this opportunity for a synergic cooperation of globular cluster and high-redshift research. [ABSTRACT FROM AUTHOR]

The first stars, galaxies, star clusters, and direct-collapse black holes are expected to have formed in low-mass (∼105–109 M⊙) haloes at Cosmic Dawn (z ∼ 10–30) under conditions of efficient gas cooling, leading to gas collapse towards the centre of the halo. The halo mass cooling threshold has been analysed by several authors using both analytical models and numerical simulations, with differing results. Since the halo number density is a sensitive function of the halo mass, an accurate model of the cooling threshold is needed for (semi-)analytical models of star formation at Cosmic Dawn. In this paper, the cooling threshold mass is calculated (semi-)analytically, considering the effects of H2-cooling and formation (in the gas phase and on dust grains), cooling by atomic metals, Lyman-α cooling, photodissociation of H2 by Lyman–Werner photons (including self-shielding by H2), photodetachment of H− by infrared photons, photoevaporation by ionization fronts, and the effect of baryon streaming velocities. We compare the calculations to several high-resolution cosmological simulations, showing excellent agreement. We find that in regions of typical baryon streaming velocities, star formation is possible in haloes of mass ≳ 1–2 × 106 M⊙ for z ≳ 20. By z ∼ 8, the expected Lyman–Werner background suppresses star formation in all minihaloes below the atomic cooling threshold (T vir = 104 K). The halo mass cooling threshold increases by another factor of ∼4 following reionization, although this effect is slightly delayed (z ∼ 4–5) because of effective self-shielding. [ABSTRACT FROM AUTHOR]

In this paper, we extend the Tiered Radio Extragalactic Continuum Simulation (T-RECS) to include H  i emission. The H  i T-RECS model is based on the most recent H  i mass function estimates, combined with prescriptions to convert H  i mass to total integrated H  i flux. It further models source size, morphology, and kinematics, including rotational velocity and H  i line width. The continuum T-RECS model is updated to improve the agreement with deeper number counts available at 150 MHz. The model for star-forming galaxies (SFGs) is also modified according to the most recent indications of a star formation rate (SFR)–radio luminosity relation, which depends primarily on stellar mass rather than redshift. We further introduce prescriptions to associate an H  i mass with the T-RECS radio continuum SFG and active galactic nuclei (AGN) populations. This gives us a way to meaningfully associate counterparts between H  i and continuum catalogues, thus building H  i × continuum simulated observations. Clustering properties of the sources in both H  i and the continuum are reproduced by associating the galaxies with dark matter haloes in a cosmological simulation. We deliver a set of mock catalogues as well as the code to produce them, which can be used for simulating observations and predicting results from radio surveys with existing and forthcoming radio facilities, such as the Square Kilometre Array (SKA). [ABSTRACT FROM AUTHOR]

We study properties of the interstellar medium, an ingredient of stars, and star formation activity, in four nearby galaxy pairs in the early and mid stages of interaction for both a galaxy scale and a kpc scale. The galaxy-scale Kennicutt–Schmidt law shows that seven of eight interacting galaxies have a star formation rate within a factor of three compared with the best fit of the isolated galaxies, although we have shown that molecular hydrogen gas is efficiently produced from atomic hydrogen during the interaction in the previous paper (Kaneko et al. 2017 PASJ, 69, 66). The galaxy-scale specific star formation rate (sSFR) and star formation efficiency (SFE) in interacting galaxies are comparable to those in isolated galaxies. We also investigate SFE and the Kennicutt–Schmidt law on a kpc scale. The spatial distributions of SFE reveal that SFE is locally enhanced, and the enhanced regions take place asymmetrically or at off-centre regions. The local enhancement of SFE could be induced by shock. We find that the index of the Kennicutt–Schmidt law for the interacting galaxies in the early stage is 1.30 ± 0.04, which is consistent with that of the isolated galaxies. Since CO emission, which is used in the Kennicutt–Schmidt law, is a tracer of the amount of molecular gas, this fact suggests that dense gas, which is more directly connected to star formation, is not changed at the early stage of interaction. [ABSTRACT FROM AUTHOR]

Without active galactic nucleus (AGN) feedback, simulated massive, star-forming galaxies become too compact relative to observed galaxies at z ≲ 2. In this paper, we perform high-resolution re-simulations of a massive (⁠|$M_{\star }\sim 10^{11}\, \rm {{\rm M}_{\odot }}$|⁠) galaxy at z ∼ 2.3, drawn from the Feedback in Realistic Environments (FIRE) project. In the simulation without AGN feedback, the galaxy experiences a rapid starburst and shrinking of its half-mass radius. We experiment with driving mechanical AGN winds, using a state-of-the-art hyper-Lagrangian refinement technique to increase particle resolution. These winds reduce the gas surface density in the inner regions of the galaxy, suppressing the compact starburst and maintaining an approximately constant half-mass radius. Using radiative transfer, we study the impact of AGN feedback on the magnitude and extent of the multiwavelength continuum emission. When AGN winds are included, the suppression of the compact, dusty starburst results in lowered flux at FIR wavelengths (due to decreased star formation) but increased flux at optical-to-near-IR wavelengths (due to decreased dust attenuation, in spite of the lowered star formation rate), relative to the case without AGN winds. The FIR half-light radius decreases from ∼1 to |$\sim 0.1\, \rm {kpc}$| in |$\lesssim 40\, \rm {Myr}$| when AGN winds are not included, but increases to |$\sim 2\, \rm {kpc}$| when they are. Interestingly, the half-light radius at optical-NIR wavelengths remains approximately constant over |$35\, \rm {Myr}$|⁠ , for simulations with and without AGN winds. In the case without winds, this occurs despite the rapid compaction, and is due to heavy dust obscuration in the inner regions of the galaxy. This work highlights the importance of forward-modelling when comparing simulated and observed galaxy populations. [ABSTRACT FROM AUTHOR]

The distribution of baryons provides a significant way to understand the formation of galaxy clusters by revealing the details of its internal structure and changes over time. In this paper, we present theoretical studies on the scaled profiles of physical properties associated with the baryonic components, including gas density, temperature, metallicity, pressure and entropy as well as stellar mass, metallicity and satellite galaxy number density in galaxy clusters from z = 4 to z = 0 by tracking their progenitors. These mass-complete simulated galaxy clusters are coming from The Three Hundred with two runs: Gizmo-SIMBA and Gadget-X. Through comparisons between the two simulations, and with observed profiles that are generally available at low redshift, we find that (1) the agreements between the two runs and observations are mostly at outer radii r ≳ 0.3 r 500, in line with the self-similarity assumption. While Gadget-X shows better agreements with the observed gas profiles in the central regions compared to Gizmo-SIMBA ; (2) the evolution trends are generally consistent between the two simulations with slightly better consistency at outer radii. In detail, the gas density profile shows less discrepancy than the temperature and entropy profiles at high redshift. The differences in the cluster centre and gas properties imply different behaviours of the AGN models between Gadget-X and Gizmo-SIMBA , with the latter, maybe too strong for this cluster simulation. The high-redshift difference may be caused by the star formation and feedback models or hydrodynamics treatment, which requires observation constraints and understanding. [ABSTRACT FROM AUTHOR]

The EXperiment for Cryogenic Large-Aperture Intensity Mapping (EXCLAIM) is a balloon-borne cryogenic telescope that will survey the spectrum of diffuse emission from both the Milky Way and the cosmic web to probe star formation, the interstellar medium, and galaxy evolution across cosmic time. EXCLAIM's primary extragalactic science survey maps 305 deg2 along the celestial equator with an R  = 512 spectrometer over the frequency range ν = 420 − 540 GHz, targeting emission of the [C  ii ] line over redshifts 2.5 < z < 3.5 and several CO lines for z < 1. Cross-correlation with galaxy redshift catalogues isolates line emission from the large-scale structure at target redshifts. In this paper, we forecast the sensitivity for both the two-point and conditional one-point cross-correlation. We predict that EXCLAIM will detect both the [C  ii ]-QSO cross-power spectrum and the conditional voxel intensity distribution (CVID) at various redshifts under a broad range of [C  ii ] intensity models, allowing it to differentiate among these models in the literature. These forecasts for the power spectra include the effects of line interlopers and continuum foreground contamination. We then convert the joint [C  ii ] constraints from both the cross-power spectrum and the CVID into constraints on the [C  ii ] halo luminosity–mass relation |$L_\mathrm{[C\, \small {II}]}(M)$| model parameters and the star formation rate density (SFRD) from [C  ii ] emission. We also develop sensitivity estimates for CO, showing the ability to differentiate between models. [ABSTRACT FROM AUTHOR]

In this paper, we place the Milky Way (MW) in the context of similar-looking galaxies in terms of their star-formation and chemical evolution histories. We select a sample of 138 Milky Way analogues (MWAs) from the SDSS-IV/MaNGA survey based on their masses, Hubble types, and bulge-to-total ratios. To compare their chemical properties to the detailed spatially-resolved information available for the MW, we use a semi-analytic spectral fitting approach, which fits a self-consistent chemical-evolution and star-formation model directly to the MaNGA spectra. We model the galaxies' inner and outer regions assuming that some of the material lost in stellar winds falls inwards. We also incorporate chemical enrichment from type II and Ia supernovae to follow the alpha-element abundance at different metallicities and locations. We find some MWAs where the stellar properties closely reproduce the distribution of age, metallicity, and alpha enhancement at both small and large radii in the MW. In these systems, the match is driven by the longer time-scale for star formation in the outer parts, and the inflow of enriched material to the central parts. However, other MWAs have very different histories. These divide into two categories: self-similar galaxies where the inner and outer parts evolve identically; and centrally-quenched galaxies where there is very little evidence of late-time central star formation driven by material accreted from the outer regions. We find that, although selected to be comparable, there are subtle morphological differences between galaxies in these different classes, and that the centrally-quenched galaxies formed their stars systematically earlier. [ABSTRACT FROM AUTHOR]