Astrophysics of Galaxies (astro-ph.GA)

Abstract: We conduct a comparison of the merging galaxy populations detected by a sample of visual identification of tidal features around galaxies as well as spectroscopically-detected close pairs of galaxies to determine whether our method of selecting merging galaxies biases our understanding of galaxy interactions. Our volume-limited parent sample consists of 852 galaxies from the Galaxy And Mass Assembly (GAMA) survey in the redshift range $0.04 \leq z \leq 0.20$ and stellar mass range $9.50 \leq$ log$_{10}(M_{\star}/\rm{M}_{\odot})\leq 11.0$. We conduct our comparison using images from the Ultradeep layer of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) to visually-classify galaxies with tidal features and compare these to the galaxies in the GAMA spectroscopic close-pair sample. We identify 198 galaxies possessing tidal features, resulting in a tidal feature fraction $f_{\rm{tidal}}$ = 0.23 $\pm$ 0.02. We also identify 80 galaxies involved in close pairs, resulting in a close pair fraction $f_{\rm{pair}}$ = 0.09 $\pm$ 0.01. Upon comparison of our tidal feature and close pair samples we identify 42 galaxies that are present in both samples, yielding a fraction $f_{\rm{both}}$ = 0.05 $\pm$ 0.01. We find evidence to suggest that the sample of close pairs of galaxies is more likely to detect early-stage mergers, where two separate galaxies are still visible, and the tidal feature sample detects later-stage mergers, where only one galaxy nucleus remains visible. The overlap of the close pair and tidal feature samples likely detect intermediate-stage mergers. Our results are in good agreement with the predictions of cosmological hydrodynamical simulations regarding the populations of merging galaxies detected by close pair and tidal feature samples.

Abstract: Using a semi-analytic galaxy-formation model, we study analogues of 8 recently discovered JWST galaxies at $z>{\sim}12$. We select analogues from a cosmological simulation with a $(311{\rm cMpc})^3$ volume and an effective particle number of $10^{12}$ enabling resolution of every atomic-cooling galaxy at $z{\le}20$. We vary model parameters to reproduce the observed UV luminosity function of $5{<}z{<}13$, aiming for a statistically representative high-redshift galaxy mock catalogue. Using the forward-modelled JWST photometry, we identify analogues from this catalogue and study their properties as well as possible evolutionary paths and local environments. We find faint JWST galaxies ($M_{\rm UV}>{\sim}-19.5$) remain consistent with standard galaxy-formation model and that our fiducial catalogue includes large samples of their analogues. The properties of these analogues broadly agree with conventional SED fitting results, except for having systematically lower redshifts due to the evolving UV luminosity function, and for having higher specific star formation rates as a result of burstier histories in our model. On the other hand, only a handful of bright galaxy analogues can be identified for observed $z{\sim}12$ galaxies. Moreover, in order to reproduce $z>{\sim}16$ JWST galaxy candidates, boosted star-forming efficiencies and reduced feedback regulation are necessary relative to models of lower-redshift populations. This suggests star formation in the first galaxies could differ significantly from their lower-redshift counterparts. We also find that these candidates are subject to low-redshift contamination, which is present in our fiducial results as both the dusty or quiescent galaxies at $z{\sim}5$.

Abstract: The putative host galaxy of FRB 20171020A was first identified as ESO 601-G036 in 2018, but as no repeat bursts have been detected, direct confirmation of the host remains elusive. In light of recent developments in the field, we re-examine this host and determine a new association confidence level of 98%. At 37 Mpc, this makes ESO 601-G036 the third closest FRB host galaxy to be identified to date and the closest to host an apparently non-repeating FRB (with an estimated repetition rate limit of < 0.011 bursts per day above 10 erg). Due to its close distance, we are able to perform detailed multi-wavelength analysis on the ESO 601-G036 system. Follow-up observations confirm ESO 601-G036 to be a typical star-forming galaxy with HI and stellar masses of log(M_HI/M_sol) ~ 9.2 and log(M_*/M_sol) = 8.64, and a star formation rate of SFR = 0.09 +/- 0.01 M_sol/yr. We detect, for the first time, a diffuse gaseous tail (log(M_HI/M_sol) ~ 8.3) extending to the south-west that suggests recent interactions, likely with the confirmed nearby companion ESO 601-G037. ESO 601-G037 is a stellar shred located to the south of ESO 601-G036 that has an arc-like morphology, is about an order of magnitude less massive, and has a lower gas metallicity that is indicative of a younger stellar population. The properties of the ESO 601-G036 system indicate an ongoing minor merger event, which is affecting the overall gaseous component of the system and the stars within ESO 601-G037. Such activity is consistent with current FRB progenitor models involving magnetars and the signs of recent interactions in other nearby FRB host galaxies.

Abstract: The upcoming ByCycle project on the VISTA/4MOST multi-object spectrograph will offer new prospects of using a massive sample of $\sim 1$ million high spectral resolution ($R$ = 20,000) background quasars to map the circumgalactic metal content of foreground galaxies (observed at $R$ = 4000 - 7000), as traced by metal absorption. Such large surveys require specialized analysis methodologies. In the absence of early data, we instead produce synthetic 4MOST high-resolution fibre quasar spectra. To do so, we use the TNG50 cosmological magnetohydrodynamical simulation, combining photo-ionization post-processing and ray tracing, to capture MgII ($\lambda2796$, $\lambda2803$) absorbers. We then use this sample to train a Convolutional Neural Network (CNN) which searches for, and estimates the redshift of, MgII absorbers within these spectra. For a test sample of quasar spectra with uniformly distributed properties ($\lambda_{\rm{MgII,2796}}$, $\rm{EW}_{\rm{MgII,2796}}^{\rm{rest}} = 0.05 - 5.15$ \AA, $\rm{SNR} = 3 - 50$), the algorithm has a robust classification accuracy of 98.6 per cent and a mean wavelength accuracy of 6.9 \AA. For high signal-to-noise spectra ($\rm{SNR > 20}$), the algorithm robustly detects and localizes MgII absorbers down to equivalent widths of $\rm{EW}_{\rm{MgII,2796}}^{\rm{rest}} = 0.05$ \AA. For the lowest SNR spectra ($\rm{SNR=3}$), the CNN reliably recovers and localizes EW$_{\rm{MgII,2796}}^{\rm{rest}}$ $\geq$ 0.75 \AA\, absorbers. This is more than sufficient for subsequent Voigt profile fitting to characterize the detected MgII absorbers. We make the code publicly available through GitHub. Our work provides a proof-of-concept for future analyses of quasar spectra datasets numbering in the millions, soon to be delivered by the next generation of surveys.

Abstract: The observed star formation rate in galaxies is well below what it should be if only gravitational collapse was at play. There is still no consensus on what is the main process responsible for the regulation the star formation rate. It has recently been shown that one candidate to regulate star formation, the feedback from massive stars, is suitable only if the mean column density at the kiloparsec scale is not too high, under $\approx 20 \mathrm{M}_\odot\cdot\mathrm{pc}^{-2}$. On the other-hand, intense large scale turbulent driving could possibly slow down star formation in high density environment to values compatible with observations. In this work we explore the effect of the nature and strength of the turbulent driving, as well as the effect of the magnetic field. We perform a large series of feedback regulated numerical simulations of the interstellar medium (ISM) in which bidimensional large scale turbulent driving is also applied. We determine the driving intensity needed to reproduce the Schmidt-Kennicutt (SK) relation for several gas column densities, magnetization and driving compressibility. We confirm that in the absence of turbulent forcing, and even with substantial magnetic field, the SFR is too high, particularly at high column density, compared to the SK relation. We find that the SFR outcome depends strongly on the initial magnetic field and on the compressibility of the turbulent driving. As a consequence, higher magnetic field in high column density environment may lower the energy necessary to sustain a turbulence sufficiently intense to regulate star formation.

Abstract: In this paper we present the distribution of molecular gas in the Milky Way Galactic plane from $l$ = [59.75, 74.75]$^{\circ}$ and $b$ = [${-}$5.25, +5.25]$^{\circ}$, using the MWISP $^{12}$CO/$^{13}$CO/$\rm {C}^{18}{O}$ emission line data. The molecular gas in this region can be mainly attributed to the Local spur, Local arm, Perseus arm, and Outer arm. Statistics of the physical properties of the molecular gas in each arm, such as excitation temperature, optical depth, and column density, are presented. Using the DBSCAN algorithm, we identified 15 extremely distant molecular clouds with kinematic distances of 14.72$-$17.77 kpc and masses of 363$-$520 M$_{\odot}$, which we find could be part of the Outer Scutum-Centaurus (OSC) arm identified by \cite{2011ApJ...734L..24D} and \cite{2015ApJ...798L..27S}. It is also possible that, 12 of these 15 extremely distant molecular clouds constitute an independent structure between the Outer and the OSC arms or a spur. There exist two Gaussian components in the vertical distribution of the molecular gas in the Perseus spiral arm. These two Gaussian components correspond to two giant filaments parallel to the Galactic plane. We find an upward warping of the molecular gas in the Outer spiral arm with a displacement of around 270 pc with respect to the Galactic mid-plane.

Abstract: We study the underlying physics of cosmic-ray (CR) driven instabilities that play a crucial role for CR transport across a wide range of scales, from interstellar to galaxy cluster environments. By examining the linear dispersion relation of CR-driven instabilities in a magnetised electron-ion background plasma, we establish that both, the intermediate and gyroscale instabilities have a resonant origin and show that these resonances can be understood via a simple graphical interpretation. These instabilities destabilise wave modes parallel to the large-scale background magnetic field at significantly distinct scales and with very different phase speeds. Furthermore, we show that approximating the electron-ion background plasma with either magnetohydrodynamics (MHD) or Hall-MHD fails to capture the fastest growing instability in the linear regime, namely the intermediate-scale instability. This finding highlights the importance of accurately characterising the background plasma for resolving the most unstable wave modes. Finally, we discuss the implications of the different phase speeds of unstable modes on particle-wave scattering. Further work is needed to investigate the relative importance of these two instabilities in the non-linear, saturated regime and to develop a physical understanding of the effective CR transport coefficients in large-scale CR hydrodynamics theories.

Abstract: Radio galaxies are a subclass of active galactic nuclei the supermassive black hole releases energy into the environment via relativistic jets. The jets are not constantly active throughout the life of the host galaxy and alternate between active and quiescent phases. Remnant radio galaxies are detected during a quiescent phase and define a class of unique sources to constrain the AGN duty cycle. We present, a spatially resolved radio analysis of the radio galaxy associated with NGC 6086 and constraints on the spectral age of the diffuse emission to investigate the duty cycle and evolution of the source. We use three new low-frequency, high-sensitivity observations, performed with the Low Frequency Array at 144 MHz and with the upgraded Giant Metrewave Radio Telescope at 400 MHz and 675 MHz. To these, we add two Very Large Array archival observations at 1400 and 4700 MHz. In the new observations, we detect a second pair of larger lobes and three regions with a filamentary morphology. We analyse the spectral index trend in the inner remnant lobes and see systematic steeper values at the lower frequencies compared to the GHz ones. Steeper spectral indices are found in the newly detected outer lobes (up to 2.1), as expected if they trace a previous phase of activity of the AGN. However, the differences between the spectra suggest different dynamical evolution within the intragroup medium during their expansion and/or different magnetic field values. We place constraints on the age of the inner and outer lobes and derive the duty cycle of the source. This results in a total active time of $\sim$39%. The filamentary structures have a steep spectral index ($\sim$1) without any spectral index trend and only one of them shows a steepening in the spectrum. Their origin is not yet clear, but they may have formed due to the compression of the plasma or due to magnetic field substructures.

Abstract: Chemical elements in the hot medium permeating early-type galaxies, groups, and clusters make them an excellent laboratory for studying metal enrichment and cycling processes in the largest scales of the Universe. Here, we report the XMM-Newton RGS analysis of 14 early-type galaxies, including the well-known brightest cluster galaxies of Perseus, for instance. The spatial distribution of the O/Fe, Ne/Fe, and Mg/Fe ratios is generally flat at the central 60 arcsecond regions of each object, irrespective of whether or not a central Fe abundance drop has been reported. Common profiles between noble gas and normal metal suggest that the dust depletion process does not work predominantly in these systems. Therefore, observed abundance drops are possibly attributed to other origins, like systematics in the atomic codes. Giant systems of high gas mass-to-luminosity ratio tend to hold a hot gas ($\sim$ 2 keV) yielding the solar N/Fe, O/Fe, Ne/Fe, Mg/Fe, and Ni/Fe ratios. Contrarily, light systems at a subkiloelectronvolt temperature regime, including isolated or group-centered galaxies, generally exhibit super-solar N/Fe, Ni/Fe, Ne/O, and Mg/O ratios. We find that the latest supernova nucleosynthesis models fail to reproduce such a super-solar abundance pattern. Possible systematic uncertainties contributing to these high abundance ratios of cool objects are also discussed in tandem with the crucial role of future X-ray missions.

Abstract: The history of reionisation is highly dependent on the ionising properties of high-redshift galaxies. It is therefore important to have a solid understanding of how the ionising properties of galaxies are linked to physical and observable quantities. In this paper, we use the First Light and Reionisation Epoch Simulations (FLARES) to study the Lyman-continuum (LyC, i.e. hydrogen-ionising) emission of massive ($M_*>10^8\,\mathrm{M_\odot}$) galaxies at redshifts $z=5-10$. We find that the specific ionising emissivity (i.e. intrinsic ionising emissivity per unit stellar mass) decreases as stellar mass increases, due to the combined effects of increasing age and metallicity. FLARES predicts a median ionising photon production efficiency (i.e. intrinsic ionising emissivity per unit intrinsic far-UV luminosity) of $\log_{10}(\xi_{\rm ion}\rm{/erg^{-1}Hz})=25.40^{+0.16}_{-0.17}$, with values spanning the range $\log_{10}(\xi_{\rm ion}\rm{/erg^{-1}Hz})=25-25.75$. This is within the range of many observational estimates, but below some of the extremes observed. We compare the production efficiency with observable properties, and find a weak negative correlation with the UV-continuum slope, and a positive correlation with the OIII equivalent width. We also consider the dust-attenuated production efficiency (i.e. intrinsic ionising emissivity per unit dust-attenuated far-UV luminosity), and find a median of $\log_{10}(\xi_{\rm ion}\rm{/erg^{-1}Hz})\sim25.5$. Within our sample of $M_*>10^8\,\mathrm{M_\odot}$ galaxies, it is the stellar populations in low mass galaxies that contribute the most to the total ionising emissivity. Active galactic nuclei (AGN) emission accounts for $10-20$ % of the total emissivity at a given redshift, and extends the LyC luminosity function by $\sim0.5$ dex.

Abstract: With the successful launch and commissioning of JWST we are now able to routinely spectroscopically probe the rest-frame optical emission of galaxies at $z>6$ for the first time. Amongst the most useful spectral diagnostics used in the optical is the Balmer/4000~\AA\ break; this is, in principle, a diagnostic of the mean ages of composite stellar populations. However, the Balmer break is also sensitive to the shape of the star formation history, the stellar (and gas) metallicity, the presence of nebular continuum emission, and dust attenuation. In this work we explore the origin of the Balmer/4000~\AA\ break using the SYNTHESIZER synthetic observations package. We then make predictions of the Balmer/4000~\AA\ break using the First Light and Reionisation Epoch Simulations (FLARES) at $5<z<10$. We find that the average break strength weakly correlates with stellar mass and rest-frame far-UV luminosity, but that this is predominantly driven by dust attenuation. We also find that break strength provides a weak diagnostic of the age but performs better as a means to constrain star formation and stellar mass, alongside the UV and optical luminosity, respectively.

Abstract: Context. Observations of $\rm ^{14}N/^{15}N$ in the interstellar medium are becoming more frequent thanks to the increased telescope capabilities. However, interpreting these data is still puzzling. In particular, measurements of $\rm ^{14}N/^{15}N$ in diazenylium revealed high levels of anti-fractionation in cold cores. Aims. Furuya & Aikawa (2018), using astrophysical simulations coupled with a gas-grain chemical code, concluded that the $^{15}$N-depletion in prestellar cores could be inherited from the initial stages, when $\rm ^{14}N^{15}N$ is selectively photodissociated and 15N atoms deplete onto the dust grain, forming ammonia ices. We aim to test this hypothesis. Methods. We targeted three sources (the prestellar core L1544, the protostellar envelope IRAS4A, and the shocked region L1157-B1) with distinct degrees of desorption or sputtering of the ammonia ices. We observed the NH3 isotopologues with the GBT, and we inferred the $\rm ^{14}N/^{15}N$ via a spectral fitting of the observed inversion transitions. Results. $^{15}$NH3(1,1) is detected in L1544 and IRAS4A, whilst only upper limits are deduced in L1157-B1. The NH3 isotopic ratio is significantly lower towards the protostar than at the centre of L1544, where it is consistent with the elemental value. We also present the first spatially resolved map of NH3 nitrogen isotopic ratio towards L1544. Conclusions. Our results are in agreement with the hypothesis that ammonia ices are enriched in $^{15}$N, leading to a decrease of the $\rm ^{14}N/^{15}N$ ratio when the ices are sublimated into the gas phase for instance due to the temperature rise in protostellar envelopes. The ammonia $\rm ^{14}N/^{15}N$ value at the centre of L1544 is a factor of 2 lower than that of N2H+, suggesting that the dominant formation pathway is hydrogenation of N atoms on dust grains, followed by non-thermal desorption.

Abstract: We present a spatially resolved study of stellar populations in 6 galaxies with stellar masses $M_*\sim10^{10}M_\odot$ at $z\sim3.7$ using 14-filter JWST/NIRCam imaging from the JADES and JEMS surveys. The 6 galaxies are visually selected to have clumpy substructures with distinct colors over rest-frame $3600-4100$\r{A}, including a bright dominant stellar core that is close to their stellar-light centroids. With 23-filter photometry from HST to JWST, we measure the stellar-population properties of individual structural components via SED fitting using Prospector. We find that the central stellar cores are $\gtrsim2$ times more massive than the Toomre mass, indicating they may not form via in-situ fragmentation. The stellar cores have stellar ages of $0.4-0.7$ Gyr that are similar to the timescale of clump inward migration due to dynamical friction, suggesting that they likely instead formed through the coalescence of giant stellar clumps. While they have not yet quenched, the 6 galaxies are below the star-forming main sequence by $0.2-0.7$ dex. Within each galaxy, we find that the specific star formation rate is lower in the central stellar core, and the stellar-mass surface density of the core is already similar to quenched galaxies of the same masses and redshifts. Meanwhile, the stellar ages of the cores are either comparable to or younger than the extended, smooth parts of the galaxies. Our findings are consistent with model predictions of the gas-rich compaction scenario for the buildup of galaxies' central regions at high redshifts. We are likely witnessing the coeval formation of dense central cores, along with the onset of galaxy-wide quenching at $z>3$.