Astrophysics of Galaxies (astro-ph.GA)

Abstract: This work presents the magnetic field geometry in a complex of three cometary (with head-tail morphology) globules, namely LDN 323, LDN 328, and LDN 331, using R-band polarization measurements of background stars. These observations were combined with a Planck sky survey to study the large-scale morphology of the magnetic fields in the region. The distances of the target stars were adopted from the Gaia catalog. The variation of degree of polarization and polarization position angle with distances of stars is analyzed. The field geometry is mostly found to follow the cometary shape of the cloud, with some randomness at certain locations. For studying the correlation between cloud morphology and magnetic field orientations, a modified version of the Histogram of Relative Orientation analysis was employed.

Abstract: We conducted 3D-MHD simulations to investigate the feedback processes in the central 1kpc scale of galaxies hosting both active star formation (SF) and an AGN wind. Our simulations naturally generated a turbulent and clumpy interstellar medium driven by SF evolution. We found that the AGN wind duty cycle plays a crucial role in shaping the evolution of the outflows. This cycle consists of an active, a remnant and an inactive phase, lasting up to 1.5 Myr. The duration of the cycle increases with larger star formation rate (SFR) and smaller AGN wind power (tested for luminosities log L = 42-44 ergs per second and SFR=1-1000 solar masses per year. The feedback on SF, whether positive or negative, depends on various factors, including the AGN outflow opening angle, power, and phase of activity, as well as the initial SFR. The passage of the AGN wind enhances SF in a ring around it, resembling the structures observed in ULIRGs, and is stronger for larger AGN power or SFR. Also, a higher SFR enhances the mixing of interstellar matter with the AGN wind, resulting in a greater number of colder, denser structures with volume filling factors ~ 0.02 to 0.12 and velocities comparable to those observed in Seyferts and LINERs, but smaller than those observed in ULIRGs. The efficiency of the AGN wind in transporting mass to kiloparsec distances diminishes with increasing SFR. The mass loss rates range from 50 to 250 solar masses per year within the initial 2 Myr of evolution, which aligns with observed rates in nearby Seyferts and ULIRGs.

Abstract: We present a new dust extinction technique with which we are able to retrieve parsec-scale gas surface density maps for entire nearby galaxies. The method measures the dust attenuation in optical bands on a pixel-by-pixel basis against a smoothed, reconstructed stellar distribution. The contribution of foreground light along the line-of-sight is calibrated using dust emission observations, assuming that the dust sits in a layer close to the mid-plane of the face-on galaxy. Here, we apply this technique to M51 (NGC 5194) as a proof-of-concept, obtaining a resolution of 0.14" (5 pc). Our dust (and gas) surface density map is consistent with independent dust- and CO-based studies at lower resolution. We find that discrepancies between our estimates of surface density and other studies stem primarily from the choice of dust model (i.e. different dust absorption coefficients). When assuming the same dust opacity law, our technique produces surface densities that are consistent with independent studies. This dust extinction technique provides us with gas surface density maps at an unprecedented resolution for full disc coverage studies of nearby galaxies. The resulting well-resolved spatial information opens the possibility for more in-depth examination of the influence of large-scale dynamics (and also stellar feedback mechanisms) on the interstellar medium at parsec-scales, and consequently star formation in nearby galaxies.

Abstract: Diagnostic diagrams of emission-line ratios have been used extensively to categorize extragalactic emission regions; however, these diagnostics are occasionally at odds with each other due to differing definitions. In this work, we study the applicability of supervised machine-learning techniques to systematically classify emission-line regions from the ratios of certain emission lines. Using the Million Mexican Model database, which contains information from grids of photoionization models using \texttt{cloudy}, and from shock models, we develop training and test sets of emission line fluxes for three key diagnostic ratios. The sets are created for three classifications: classic \hii{} regions, planetary nebulae, and supernova remnants. We train a neural network to classify a region as one of the three classes defined above given three key line ratios that are present both in the SITELLE and MUSE instruments' band-passes: [{\sc O\,iii}]$\lambda5007$/H$\beta$, [{\sc N\,ii}]$\lambda6583$/H$\alpha$, ([{\sc S\,ii}]$\lambda6717$+[{\sc S\,ii}]$\lambda6731$)/H$\alpha$. We also tested the impact of the addition of the [{\sc O\,ii}]$\lambda3726,3729$/[{\sc O\,iii}]$\lambda5007$ line ratio when available for the classification. A maximum luminosity limit is introduced to improve the classification of the planetary nebulae. Furthermore, the network is applied to SITELLE observations of a prominent field of M33. We discuss where the network succeeds and why it fails in certain cases. Our results provide a framework for the use of machine learning as a tool for the classification of extragalactic emission regions. Further work is needed to build more comprehensive training sets and adapt the method to additional observational constraints.

Abstract: We perform numerical simulations of supersonic magnetohydrodynamic (MHD) turbulence and calculate Fourier power spectra of E and B modes arising from dust polarization. We pay close attention to the ratio of E-mode to B-mode spectra (a.k.a. E/B power asymmetry) on small spatial scales. We find that the ratio depends on the strength of the mean magnetic field: the stronger the mean magnetic field is, the smaller the ratio is. More precisely speaking, the ratio scales with the Alfv\'en Mach number $M_A$, the root-mean-square velocity divided by the Alfv\'en speed of the mean magnetic field, when it lies in the range $1\lesssim M_A \lesssim 30$. This result implies that we can use the E/B power asymmetry to constrain the strength of the mean magnetic field in supersonic and super-Alfv\'enic MHD turbulence.

Abstract: A chemodynamical model of our galaxy is fitted to data from DR17 of the APOGEE survey supplemented with data from the StarHorse catalogue and gaia DR3. Dynamically, the model is defined by action-based distribution functions for dark matter and six stellar components plus a gas disc. The gravitational potential jointly generated by the model's components is used to examine the galaxy's chemical composition within action space. The observational data probably cover all parts of action space that are populated by stars. The overwhelming majority of stars have angular momentum J_\phi>0 implying that they were born in the Galactic disc. High-alpha stars dominate in a region that is sharply bounded by J_\phi \la J_\phi(solar). Chemically the model is defined by giving each stellar component a Gaussian distribution in ([Fe/H],[Mg/Fe]) space about a mean that is a linear function of the actions. The model's 47 dynamical and 70 chemical parameters are chosen to maximise the likelihood of the data given the model in 72 three-dimensional velocity spaces and 30 two-dimensional chemical spaces. The circular speed falls steadily from 237\kms at R=4\kpc to 218\kms at R=20\kpc. Dark matter contributes half the radial force on the Sun and has local density 0.011\msun\pc^{-3}, there being 24.5\msun\pc^{-2} in dark matter and 26.5\msun\pc^{-2} in stars within 1.1\kpc of the plane.

Abstract: The first bright objects to form in the Universe at redshift $z \sim 10-20$ might have been Dark Stars, made primarily of hydrogen and helium but powered by dark matter. In this study, we investigate the detectability of Supermassive Dark Stars (SMDS) by the Roman Space Telescope. RST will be able to detect SMDSs at redshifts as high as $z\simeq 14$. In cases with gravitational lensing factors of $\mu\sim 100$, RST will be able to find SMDS as small as $\sim10^4 M_{\odot}$ at $z\sim 12$ with $\sim 10^6$ s of exposure. To differentiate SMDS from early galaxies containing zero metallicity stars at similar redshifts, we compare their spectra, photometry in RST bands, color indexes and image morphology. With RST alone, the differentiation is possible only for limited cases: SMDS formed via "adiabatic contraction" (DM pulled into the star via gravity alone) with $M\gtrsim 10^5M_{\odot}$ and lensed by $\mu\gtrsim 30$ have distinct photometric signatures from those of the first galaxies. For SMDSs formed via "dark matter capture," their spectra are degenerate to those of many galaxies with little to no nebular emission. Thus with RST alone, the only way to tell them apart from first galaxies would be via image morphology: i.e. point object (SMDSs) vs. extended object (sufficiently magnified galaxies). However, if the same objects are further examined by JWST spectroscopy, a "smoking gun" for detection of SMDS is the HeII $\lambda$1640 absorption line. While RST does not cover the wavelength band required to find this line (for $z_{\rm emi}\gtrsim 10$), JWST does. Hence the two detectors can be used together in identifying SMDS. The confirmed detection of any SMDSs will provide evidence for a new type of star, powered by dark matter. Moreover, such massive stars can also be natural progenitors of the supermassive black holes powering the extremely bright quasars observed at $z\gtrsim 6$.

Abstract: The bulk kinematics and thermodynamics of hot supernovae-driven galactic winds is critically dependent on both the amount of swept up cool clouds and non-spherical collimated flow geometry. However, accurately parameterizing these physics is difficult because their functional forms are often unknown, and because the coupled non-linear flow equations contain singularities. We show that deep neural networks embedded as individual terms in the governing coupled ordinary differential equations (ODEs) can robustly discover both of these physics, without any prior knowledge of the true function structure, as a supervised learning task. We optimize a loss function based on the Mach number, rather than the explicitly solved-for 3 conserved variables, and apply a penalty term towards near-diverging solutions. The same neural network architecture is used for learning both the hidden mass-loading and surface area expansion rates. This work further highlights the feasibility of neural ODEs as a promising discovery tool with mechanistic interpretability for non-linear inverse problems.

Abstract: As the nearest accessible massive early-type galaxy, NGC 5128 presents an exceptional opportunity to measure dark matter halo parameters for a representative elliptical galaxy. Here we take advantage of rich new observational datasets of large-radius tracers to perform dynamical modeling of NGC 5128, using a discrete axisymmetric anisotropic Jeans approach with a total tracer population of nearly 1800 planetary nebulae, globular clusters, and dwarf satellite galaxies extending to a projected distance of $\sim250$ kpc from the galaxy center. We find that a standard NFW halo provides an excellent fit to nearly all the data, excepting a subset of the planetary nebulae that appear to be out of virial equilibrium. The best-fit dark matter halo has a virial mass of ${\rm M}_{vir}=4.4^{+2.4}_{-1.4}\times10^{12} {\rm M}_{\odot}$, and NGC 5128 appears to sit below the mean stellar mass--halo mass and globular cluster mass--halo mass relations, which both predict a halo virial mass closer to ${\rm M}_{vir} \sim 10^{13} {\rm M}_{\odot}$. The inferred NFW virial concentration is $c_{vir}=5.6^{+2.4}_{-1.6}$, nominally lower than $c_{vir} \sim 9$ predicted from published $c_{vir}$--${\rm M}_{vir}$ relations, but within the $\sim 30\%$ scatter found in simulations. The best-fit dark matter halo constitutes only $\sim10\%$ of the total mass at 1 effective radius but $\sim50\%$ at 5 effective radii. The derived halo parameters are relatively insensitive to reasonable variations in the tracer population considered, tracer anisotropies, and system inclination. Our analysis highlights the value of comprehensive dynamical modeling of nearby galaxies, and the importance of using multiple tracers to allow cross-checks for model robustness.

Abstract: Star-forming galaxies like the Milky Way are surrounded by a hot gaseous halo at the virial temperature - the so-called galactic corona - that plays a fundamental role in their evolution. The interaction between the disc and the corona has been shown to have a direct impact on accretion of coronal gas onto the disc with major implications for galaxy evolution. In this work, we study the gas circulation between the disc and the corona of star-forming galaxies like the Milky Way. We use high-resolution hydrodynamical N-body simulations of a Milky Way-like galaxy with the inclusion of an observationally-motivated galactic corona. In doing so, we use SMUGGLE, an explicit interstellar medium (ISM) and stellar feedback model coupled with the moving-mesh code Arepo. We find that the reservoir of gas in the galactic corona is sustaining star formation: the gas accreted from the corona is the primary fuel for the formation of new stars, helping in maintaining a nearly constant level of cold gas mass in the galactic disc. Stellar feedback generates a gas circulation between the disc and the corona (the so-called galactic fountain) by ejecting different gas phases that are eventually re-accreted onto the disc. The accretion of coronal gas is promoted by its mixing with the galactic fountains at the disc-corona interface, causing the formation of intermediate temperature gas that enhance the cooling of the hot corona. We find that this process acts as a positive feedback mechanism, increasing the accretion rate of coronal gas onto the galaxy.

Abstract: Radio-mode feedback is a key ingredient in galaxy formation and evolution models, required to reproduce the observed properties of massive galaxies in the local Universe. We study the cosmic evolution of radio-AGN feedback out to $z\sim2.5$ using a sample of 9485 radio-excess AGN. We combine the evolving radio luminosity functions with a radio luminosity scaling relationship to estimate AGN jet kinetic powers and derive the cosmic evolution of the kinetic luminosity density, $\Omega_{\rm{kin}}$ (i.e. the volume-averaged heating output). Compared to all radio-AGN, low-excitation radio galaxies (LERGs) dominate the feedback activity out to $z\sim2.5$, with both these populations showing a constant heating output of $\Omega_{\rm{kin}} \approx 4-5 \times 10^{32}\,\rm{W\,Mpc^{-3}}$ across $0.5 < z < 2.5$. We compare our observations to predictions from semi-analytical and hydrodynamical simulations, which broadly match the observed evolution in $\Omega_{\rm{kin}}$, although their absolute normalisation varies. Comparison to the Semi-Analytic Galaxy Evolution (SAGE) model suggests that radio-AGN may provide sufficient heating to offset radiative cooling losses, providing evidence for a self-regulated AGN feedback cycle. We integrate the kinetic luminosity density across cosmic time to obtain the kinetic energy density output from AGN jets throughout cosmic history to be $\sim 10^{50}\,\rm{J\,Mpc^{-3}}$. Compared to AGN winds, the kinetic energy density from AGN jets dominates the energy budget at $z \lesssim 2$; this suggests that AGN jets play an important role in AGN feedback across most of cosmic history.

Abstract: We investigate the incidence and properties of ionized gas outflows in a sample of 52 galaxies with stellar mass between $10^7$ M$_{\odot}$ and $10^9$ M$_{\odot}$ observed with ultra-deep JWST/NIRSpec MSA spectroscopy as part of the JWST Advanced Deep Extragalactic Survey (JADES). The high-spectral resolution (R2700) NIRSpec observations allowed us to identify for the first time the signature of outflows in the rest-frame optical nebular lines in low-mass galaxies at $z>3$. The incidence fraction of ionized outflows, traced by broad components, is about 25-40$\%$ depending on the intensity of the emission lines. The low incidence fraction might be due to both the sensitivity limit and the fact that outflows are not isotropic but have a limited opening angle which results in a detection only when this is directed toward our line of sight. Evidence for outflows increases slightly with stellar mass and star-formation rate. The median velocity and mass loading factor (i.e., the ratio between mass outflow rate and star formation rate) of the outflowing ionized gas are 500 km s$^{-1}$ and $\eta=2.1^{+2.5}_{-1.6}$, respectively. These are two and 100 times higher, respectively than the typical values observed in local dwarf galaxies. These outflows are able to escape the gravitational potential of the galaxy and enrich the circum-galactic medium and, potentially, the inter-galactic medium. Our results indicate that outflows can significantly impact the star formation activity in low-mass galaxies within the first 2 Gyr of the Universe.

Abstract: Measurements of the star formation activity on cloud scales are fundamental to uncovering the physics of the molecular cloud, star formation, and stellar feedback cycle in galaxies. Infrared (IR) emission from small dust grains and polycyclic aromatic hydrocarbons (PAHs) are widely used to trace the obscured component of star formation. However, the relation between these emission features and dust attenuation is complicated by the combined effects of dust heating from old stellar populations and an uncertain dust geometry with respect to heating sources. We use images obtained with NIRCam and MIRI as part of the PHANGS--JWST survey to calibrate dust emission at 21$\rm \mu m$, and the emission in the PAH-tracing bands at 3.3, 7.7, 10, and 11.3$\rm \mu m$ as tracers of obscured star formation. We analyse $\sim$ 16000 optically selected HII regions across 16 nearby star-forming galaxies, and benchmark their IR emission against dust attenuation measured from the Balmer decrement. We model the extinction-corrected H$\alpha$ flux as the sum of the observed H$\alpha$ emission and a term proportional to the IR emission, with $a_{IR}$ as the proportionality coefficient. A constant $a_{IR}$ leads to extinction-corrected H$\alpha$ estimates which agree with those obtained with the Balmer decrement with a scatter of $\sim$ 0.1 dex for all bands considered. Among these bands, 21$\rm \mu m$ emission is demonstrated to be the best tracer of dust attenuation. The PAH-tracing bands underestimate the correction for bright HII regions, since in these environments the ratio of PAH-tracing bands to 21$\rm \mu m$ decreases, signalling destruction of the PAH molecules. For fainter HII regions all bands suffer from an increasing contamination from the diffuse infrared background.

Abstract: JWST observations of the strongly lensed galaxy The Sparkler have revealed a population of gravitationally bound globular cluster (GC) candidates. Different analyses have resulted in broadly similar ages but significantly different metallicities, questioning the assembly history that has led to the formation of such a population. In this letter, we re-analyse the two sets of photometry available in the literature with the code MCMAME especially tailored to fit physical properties of GCs. We find the ages and metallicities from both datasets are consistent within 1 $\sigma$ uncertainties. A significant group of GCs is consistent with being old and metal poor ([Fe/H] $\sim -1.7$). For this group, the ages do not converge, hence, we conclude that they are definitively older than 1 Gyr and can be as old as the age of the Universe. The remaining GCs have younger ages and a metallicity spread. The ages and metallicities distribution of GCs in the Sparkler are consistent with those observed in Local Group's galaxies at similar lookback times. Comparing with predictions from E-MOSAICS simulations we confirm that the Sparkler GC population traces the self-enrichment history of a galaxy which might become a few times $10^9$ M$_{\odot}$ massive system at redshift $z = 0$