Search

We introduce a novel orbit superposition method designed to reconstruct the stellar density structure, kinematics, and chemical abundance distribution of the entire Milky Way by leveraging 6D phase-space information from its resolved stellar populations, limited by the spatial coverage of APOGEE DR17., Comment: 17 pages, 13 figures, submitted to A&A

The Magellanic Stream has long been known to contain multiple HI strands and corresponding stellar populations are beginning to be discovered. Combining an H3-selected sample with stars drawn from the Gaia catalog, we trace stars along a sub-dominant strand of the Magellanic Stream, as defined by gas content, across 30$^\circ$ on the sky. We find that the dominant strand is devoid of stars with Galactocentric distance $\lesssim 55$ kpc while the subdominant strand shows a close correspondence to such stars. We conclude that (1) the two Stream strands have different origins, (2) they are likely only close in projection, (3) the subdominant strand is tidal in origin, and (4) the subdominant strand is composed of disk material, likely drawn from the disk of the Small Magellanic Cloud., Comment: 8 pages, 8 figures

The interplay between atomic gas, the star-formation history of a galaxy and its environment are intrinsically linked, and we need to decouple these dependencies to understand their role in galaxy formation and evolution. In this paper, we analyse the star formation histories (SFHs) of 187 galaxies from the MIGHTEE-HI Survey Early Science Release data, focusing on the relationships between HI properties and star formation. A strong correlation emerges between a galaxy's HI-to-stellar mass ratio and the time of formation, alongside an inverse correlation between stellar mass and time of formation, regardless of the inferred SFH. Additionally, galaxies with lower stellar masses and higher HI-to-stellar mass ratios exhibit longer gas depletion times compared to more massive galaxies, which appear to have depleted their gas and formed stars more efficiently. This suggests that smaller, gas-rich galaxies have higher depletion times due to shallower potential wells and less efficient star formation. Furthermore, we explore the connection between spin-filament alignment and HI content. We find no significant correlation between peak star formation activity and proximity to filaments. However, we do find that the two galaxies in our sample within 1 Mpc of a filament have very low gas-depletion timescales and have their spin axis misaligned with the filament, suggestive of a link between the galaxy properties and proximity to a filament., Comment: submitted to MNRAS; 14 pages, 14 figures

We investigate the alignment of non-red early-type galaxies (ETGs) with blue or green colors within large-scale filaments and compare this alignment pattern with that of red ETGs. Our analysis reveals a significant alignment of the major axes of red ETGs with the orientations of their host cosmic filaments, consistent with prior research. In contrast, non-red ETGs show no significant alignment signal. This divergence in alignment behavior between non-red and red ETGs implies a distinct evolutionary path for non-red ETGs, suggesting a formation process that may be independent of galaxy mergers or that recent mergers experienced by non-red ETGs may not follow the direction of the filament but rather be more random or even perpendicular to the filament orientation., Comment: Submitted to ApJ

Fast radio bursts (FRBs) are high-energy, short-duration phenomena in radio astronomy. Identifying their host galaxies can provide insights into their mysterious origins. In this paper, we introduce a novel approach to identifying potential host galaxies in three-dimensional space. We use FRB 20190425A and GW190425 as an example to illustrate our method. Recently, due to spatial and temporal proximity, the potential association of GW190425 with FRB 20190425A has drawn attention, leading to the identification of a likely host galaxy, UGC 10667, albeit without confirmed kilonova emissions. We search for the host galaxy of FRB 20190425A with a full CHIME localization map. Regardless of the validity of the association between GW190425 and FRB 20190425A, we identify an additional potential host galaxy (SDSS J171046.84+212732.9) from the updated GLADE galaxy catalog, supplementing the importance of exploring the new volume. We employed various methodologies to determine the most probable host galaxy of GW190424 and FRB 20190425A, including a comparison of galaxy properties and constraints on their reported observation limits using various Kilonova models. Our analysis suggests that current observations do not definitively identify the true host galaxy. Additionally, the Kilonova models characterized by a gradual approach to their peak are contradicted by the observational upper limits of both galaxies. Although the absence of optical emission detection raises doubts, it does not definitively disprove the connection between GW and FRB., Comment: 19 pages, 10 figures, 3 tables

We perform the first direct cosmological and astrophysical parameter inference from the combination of galaxy luminosity functions and colours using a simulation based inference approach. Using the Synthesizer code we simulate the dust attenuated ultraviolet--near infrared stellar emission from galaxies in thousands of cosmological hydrodynamic simulations from the CAMELS suite, including the Swift-EAGLE, Illustris-TNG, Simba & Astrid galaxy formation models. For each galaxy we calculate the rest-frame luminosity in a number of photometric bands, including the SDSS $\textit{ugriz}$ and GALEX FUV & NUV filters; this dataset represents the largest catalogue of synthetic photometry based on hydrodynamic galaxy formation simulations produced to date, totalling >200 million sources. From these we compile luminosity functions and colour distributions, and find clear dependencies on both cosmology and feedback. We then perform simulation based (likelihood-free) inference using these distributions, and obtain constraints on both cosmological and astrophysical parameters. Both colour distributions and luminosity functions provide complementary information on certain parameters when performing inference. Most interestingly we achieve constraints on $\sigma_8$, describing the clustering of matter. This is attributable to the fact that the photometry encodes the star formation--metal enrichment history of each galaxy; galaxies in a universe with a higher $\sigma_8$ tend to form earlier and have higher metallicities, which leads to redder colours. We find that a model trained on one galaxy formation simulation generalises poorly when applied to another, and attribute this to differences in the subgrid prescriptions, and lack of flexibility in our emission modelling. The photometric catalogues are publicly available at: https://camels.readthedocs.io/ ., Comment: 28 pages, 20 figures, submitted to MNRAS. Comments and feedback welcome!

We present the application of the image coaddition algorithm, Up-sampling and PSF Deconvolution Coaddition (UPDC), for stacking multiple exposure images captured by the James Webb Space Telescope (JWST) Near-Infrared Camera (NIRCam). By addressing the point spread function (PSF) effect, UPDC provides visually enhanced and sharper images. Furthermore, the anti-aliasing and super-resolution capabilities of UPDC make it easier to deblend sources overlapped on images, yielding a higher accuracy of aperture photometry. We apply this algorithm to the SMACS J0723 imaging data. Comparative analysis with the Drizzle algorithm demonstrates significant improvements in detecting faint sources, achieving accurate photometry, and effectively deblending (super-resolution) closely packed sources. {As a result, we have newly detected a pair of close binary stars that were previously unresolvable in the original exposures or the Drizzled image.} These improvements significantly benefit various scientific projects conducted by JWST. The resulting dataset, named "UPdec-Webb", can be accessible through the official website of the Chinese Virtual Observatory (ChinaVO)., Comment: 21 pages, 18 figures, accepted for publication in ApJS

Hub-filament systems (HFSs) are potential sites of massive star formation (MSF). To understand the role of filaments in MSF and the origin of HFSs, we conducted a multi-scale and multi-wavelength observational investigation of the molecular cloud G321.93-0.01. The $^{13}$CO($J$ = 2-1) data reveal multiple HFSs, namely, HFS-1, HFS-2, and a candidate HFS (C-HFS). HFS-1 and HFS-2 exhibit significant mass accretion rates ($\dot{M}_{||}$ $> 10^{-3}$ $M_{\odot}$ yr$^{-1}$) to their hubs (i.e., Hub-1 and Hub-2, respectively). Hub-1 is comparatively massive, having higher $\dot{M}_{||}$ than Hub-2, allowing to derive a relationship $\dot{M}_{||} \propto M^{\beta}_{\rm{hub}}$, with $\beta \sim1.28$. Detection of three compact HII regions within Hub-1 using MeerKAT 1.28 GHz radio continuum data and the presence of a clump (ATL-3), which meets Kauffmann & Pillai's criteria for MSF, confirm the massive star-forming activity in HFS-1. We find several low-mass ALMA cores (1-9 $M_{\odot}$) inside ATL-3. The presence of a compact HII region at the hub of C-HFS confirms that it is active in MSF. Therefore, HFS-1 and C-HFS are in relatively evolved stages of MSF, where massive stars have begun ionizing their surroundings. Conversely, despite a high $\dot{M}_{||}$, the non-detection of radio continuum emission toward Hub-2 suggests it is in the relatively early stages of MSF. Analysis of $^{13}$CO($J$ = 2-1) data reveals that the formation of HFS-1 was likely triggered by the collision of a filamentary cloud about 1 Myr ago. In contrast, the relative velocities ($\gtrsim 1$ km s$^{-1}$) among the filaments of HFS-2 and C-HFS indicate their formation through the merging of filaments., Comment: Accepted for publication in The Astronomical Journal, the manuscript consists of 27 pages, 11 figures, and 5 tables

We utilize a hybrid approach that integrates the traditional cross-correlation function (CCF) and machine learning to detect spectroscopic multi-systems, specifically focusing on double-line spectroscopic binary (SB2). Based on the ninth data release (DR9) of the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST), which includes a medium-resolution survey (MRS) containing 29,920,588 spectra, we identify 27,164 double-line and 3124 triple-line spectra, corresponding to 7096 SB2 candidates and 1903 triple-line spectroscopic binary (SB3) candidates, respectively, representing about 1% of the selection dataset from LAMOST-MRS DR9. Notably, 70.1% of the SB2 candidates and 89.6% of the SB3 candidates are newly identified. Compared to using only the traditional CCF technique, our method significantly improves the efficiency of detecting SB2, saves time on visual inspections by a factor of four., Comment: 18 pages, 11 figures, accepted by ApJS, Data available via China-VO PaperData repository

The progenitors of Type II-P supernovae (SN) have been confirmed to be red supergiants. However, the upper mass limit of the directly probed progenitors is much lower than that predicted by current theories, and the accurate determination of the progenitor masses is key to understand the final fate of massive stars. Located at a distance of only 6.72 Mpc, the Type II-P SN 2024ggi is one of the closest SN in the last decade. Previous studies have analyzed its progenitor by direct detection, but the derived progenitor mass may be influenced by the very uncertain circumstellar extinction and pulsational brightness variability. In this work, we try to constrain the progenitor mass with an environmental analysis based on images from the Hubble Space Telescope. We found that stars in the progenitor environment have a uniform spatial distribution without significant clumpiness, and we derived the star formation history of the environment with a hierarchical Bayesian method. The progenitor is associated with the youngest population in the SN environment with an age of log($t$/yr) = 7.41 (i.e. 25.7 Myr), which corresponds to an initial mass of $10.2^{+0.06}_{-0.09}$ $M_\odot$. Our work provides an independent measurement of the progenitor mass, which is not affected by circumstellar extinction and pulsational brightness variability., Comment: 8 pages, 3 figures; submitted to ApJL

Reproducing the physical characteristics of ultra-faint dwarf galaxies (UFDs) in cosmological simulations is challenging, particularly with respect to stellar metallicity and galaxy size. To investigate these difficulties in detail, we conduct high-resolution simulations ($M_{\rm gas} \sim 60 \, M_{\odot}$, $M_{\rm DM} \sim 370 \, M_{\odot}$ ) on six UFD analogs ($M_{\rm vir} \sim 10^8 - 10^9 \, M_{\odot}$, $M_{\rm \star} \sim 10^3 - 2.1 \times 10^4 \, M_{\odot}$). Our findings reveal that the stellar properties of UFD analogs are shaped by diverse star-forming environments from multiple progenitor halos in the early Universe. Notably, our UFD analogs exhibit a better match to the observed mass-metallicity relation (MZR), showing higher average metallicity compared to other theoretical models. The metallicity distribution functions (MDFs) of our simulated UFDs lack high-metallicity stars ($[\rm Fe/H] > -2.0$) while containing low-metallicity stars ($[\rm Fe/H] < -4.0$). Excluding these low-metallicity stars, our results align well with the MDFs of observed UFDs. However, forming stars with higher metallicity ($-2.0 \leq [\rm Fe/H]_{\rm max} \leq -1.5$) remains a challenge due to the difficulty of sustaining metal enrichment during their brief star formation period before cosmic reionization. Additionally, our simulations show extended outer structures in UFDs, resulting from dry mergers between progenitor halos. To ensure consistency, we adopt the same fitting method commonly used in observations to derive the half-light radius. We find that this method tends to produce lower values compared to direct calculations and struggles to accurately describe the extended outer structures. To address this, we employ a two-component density profile to obtain structural parameters, finding that it better describes the galaxy shape, including both inner and outer structures., Comment: 28 pages, 21 figures

We examine the quiescent fractions of massive galaxies in six $z\gtrsim3$ spectroscopically-confirmed protoclusters in the COSMOS field, one of which is newly confirmed and presented here. We report the spectroscopic confirmation of MAGAZ3NE~J100143+023021 at $z=3.122^{+0.007}_{-0.004}$ by the Massive Ancient Galaxies At $z>3$ NEar-infrared (MAGAZ3NE) survey. MAGAZ3NE~J100143+023021 contains a total of 79 protocluster members (28 spectroscopic and 51 photometric). Three spectroscopically-confirmed members are star-forming ultra-massive galaxies ($\log(M_{\star}/{\rm M}_\odot)>11$; UMGs), the most massive of which has $\log(M_{\star}/{\rm M}_\odot)=11.15^{+0.05}_{-0.06}$. Combining Keck/MOSFIRE spectroscopy and the COSMOS2020 photometric catalog, we use a weighted Gaussian kernel density estimator to map the protocluster and measure its total mass $2.25^{+1.55}_{-0.65}\times10^{14}~{\rm M}_{\odot}$ in the dense ``core'' region. For each of the six COSMOS protoclusters, we compare the quiescent fraction to the status of the central UMG as star-forming or quiescent. We observe that galaxies in these protoclusters appear to obey galactic conformity: elevated quiescent fractions are found in protoclusters with $UVJ$ quiescent UMGs and low quiescent fractions are found in protoclusters containing $UVJ$ star-forming UMGs. This correlation of star-formation/quiescence in UMGs and the massive galaxies nearby in these protoclusters is the first evidence for the existence of galactic conformity at $z>3$. Despite disagreements over mechanisms behind conformity at low redshifts, its presence at these early cosmic times would provide strong constraints on the physics proposed to drive galactic conformity., Comment: 6 figures, 5 tables. Accepted for publication in the Astrophysical Journal

Developing an efficient code for large, multiscale astrophysical simulations is crucial in preparing the upcoming era of exascale computing. RAMSES is an astrophysical simulation code that employs parallel processing based on the Message Passing Interface (MPI). However, it has limitations in computational and memory efficiency when using a large number of CPU cores. The problem stems from inefficiencies in workload distribution and memory allocation that inevitably occur when a volume is simply decomposed into domains equal to the number of working processors. We present RAMSES-yOMP, which is a modified version of RAMSES designed to improve parallel scalability. Major updates include the incorporation of Open Multi-Processing into the MPI parallelization to take advantage of both the shared and distributed memory models. Utilizing this hybrid parallelism in high-resolution benchmark simulations with full prescriptions for baryonic physics, we achieved a performance increase of a factor of 2 in the total run-time, while using 75% less memory and 30% less storage compared to the original code, when using the same number of processors. These improvements allow us to perform larger or higher-resolution simulations than what was feasible previously., Comment: 13 pages, 10 figures, submitted to ApJ

Upcoming surveys are likely to discover a new sample of interstellar objects (ISOs) within the Solar System, but questions remain about the origin and distribution of this population within the Galaxy. ISOs are ejected from their host systems with a range of velocities, spreading out into tidal streams - analogous to the stellar streams routinely observed from the disruption of star clusters and dwarf galaxies. We create a simulation of ISO streams orbiting in the Galaxy, deriving a simple model for their density distribution over time. We then construct a population model to predict the properties of the streams in which the Sun is currently embedded. We find that the number of streams encountered by the Sun is quite large, ~ 10^6 or more. However, the wide range of stream properties means that for reasonable future samples of ISOs observed in the Solar System, we may see ISOs from the same star ("siblings"), and we are likely to see ISOs from the same star cluster ("cousins"). We also find that ISOs are typically not traceable to their parent star, though this may be possible for ISO siblings. Any ISOs observed with a common origin will come from younger, dynamically colder streams., Comment: Submitted to AAS Journals, comments welcome

We present a systematic search for radio active galactic nuclei (AGNs) in dwarf galaxies using recent observations taken by the Very Large Array Sky Survey (VLASS). To select these objects, we first establish a criterion to identify radio-excess AGNs using the infrared-radio correlation (IRRC) parameter, $q$, that describes the tight relation between radio and IR emission in star forming (SF) galaxies. We find a $2\sigma$ threshold of $q < 1.94$ to select radio-excess AGNs, which is derived from a sample of $\sim 7,000$ galaxies across the full mass range in the NASA-Sloan Atlas (NSA) that have radio and IR detections from VLASS and the Wide-Field Infrared Survey Explorer, respectively. We create catalogs of radio-excess AGNs and SF galaxies and make these available to the community. Applying our criterion to dwarf galaxies with stellar masses $M_\star \lesssim 3 \times 10^9 M_\odot$ and redshifts $z \le 0.15$, and carefully removing interlopers, we find 10 radio-excess AGNs with radio-optical positional offsets between $\sim$ 0 and 2.3 arcseconds (0 - 2.7 kpc). Based on statistical arguments and emission line diagnostics, we expect the majority of these radio-excess AGNs to be associated with the dwarf host galaxies rather than background AGNs. Five of the objects have evidence for hosting AGNs at other wavelengths, and 5 objects are identified as AGNs in dwarf galaxies for the first time. We also identify 8 variable radio sources in dwarf galaxies by comparing the VLASS epoch 1 and epoch 2 observations to FIRST detections presented in arXiv:1909.04670., Comment: 27 pages, 21 figures, 3 tables, accepted for publication in The Astrophysical Journal

The UV continuum slope of galaxies, $\beta$, is a powerful diagnostic. Understanding the redshift evolution of $\beta$ and its dependence on key galaxy properties can shed light on the evolution of galaxy physical properties over cosmic time. In this study, we present $\beta$ measurements for 295 spectroscopically confirmed galaxies at $5.5 9.5$ to begin to redden, deviating from the trend observed at $z < 9.5$. By producing stacked spectra in bins of redshift and $\beta$, we derive trends between $\beta$ and dust attenuation, metallicity, ionization parameter, and stellar age indicators directly from spectra, finding a lack of dust attenuation to be the dominant driver of bluer $\beta$ values. We further report six galaxies with $\beta<-3.0$, which show a range of spectroscopic properties and signs of significant LyC photon leakage. Finally, we show that the redder $\beta$ values at $z > 9.5$ may require rapid build-up of dust reservoirs in the very early Universe or a significant contribution from the nebular continuum emission to the observed UV spectra, with the nebular continuum fraction depending on the gas temperatures and densities. Our modeling shows that in the absence of dust, nebular emission at $T > 15,000$ K can reproduce the range of $\beta$ that we see in our sample. Higher gas temperatures driven by hot, massive stars can boost the fraction of nebular continuum emission, potentially explaining the observed $\beta$ values as well as bright UV magnitudes seen across galaxies at $z > 10$., Comment: submitted to MNRAS, 22 pages including 1 appendix, 12 Figures (including 1 in appendix), comments welcome!

Deuterium fractionation is a well-established evolutionary tracer in low-mass star formation, but its applicability to the high-mass regime remains an open question. The abundances and ratios of deuterated species have often been proposed as reliable evolutionary indicators for different stages of the high-mass star formation. We investigate the role of N$_2$H$^+$ and key deuterated molecules as tracers of the different stages of the high-mass star formation, and test whether their abundance ratios can serve as reliable evolutionary indicators. We conducted APEX observations of o-H$_2$D$^+$ (1$_{10}$-1$_{11}$), N$_2$H$^+$ (4-3), and N$_2$d$^+$ (3-2) in 40 high-mass clumps at different evolutionary stages, selected from the ATLASGAL survey. Molecular column densities ($N$) and abundances ($X$), were derived through spectral line modelling, both under local thermodynamic equilibrium (LTE) and non-LTE conditions. The $N$(o-H$_2$D$^+$) show the smallest deviation from LTE results when derived under non-LTE assumptions. In contrast, N$_2$D$^+$ shows the largest discrepancy between the $N$ derived from LTE and non-LTE. In all the cases discussed, we found that $X$(o-H$_2$D$^+$) decreases more significantly with time than in the case of $X$(N$_2$D$^+$); whereas $X$(N$_2$H$^+$) increases slightly. Therefore, the validity of the recently proposed $X$(o-H$_2$D$^+$)/$X$(N$_2$D$^+$) ratio as a reliable evolutionary indicator was not observed for this sample. While the deuteration fraction derived from N$_2$D$^+$ and N$_2$H$^+$ clearly decreases with clump evolution, the interpretation of this trend is complex, given the different distribution of the two tracers. Our results suggest that a careful consideration of the observational biases and beam-dilution effects are crucial for an accurate interpretation of the evolution of the deuteration process during the high-mass star formation process., Comment: 13 pages, 5 figures, 3 tables. Accepted for publication in A&A

We present a catalog of 35 new candidate (13 high confidence) isolated, young stellar systems within the Virgo galaxy cluster identified through a citizen science search of public optical and ultraviolet imaging. "Blue blobs" are a class of blue, faint, isolated, extremely low stellar mass, and metal-rich star-forming clouds embedded in the hot intracluster medium of the Virgo cluster. Only six blue blobs were known previously and here we confirm an additional six of our candidates through velocity and metallicity measurements from follow-up optical spectroscopy on the Hobby-Eberly Telescope (HET). Our 13 high confidence candidates (including the six confirmed) have properties consistent with prior known blue blobs and are inconsistent with being low-mass galaxies. Most candidates are concentrated in relatively dense regions, roughly following filamentary structures within the cluster, but avoiding its center. Three of our candidates are likely the stellar counterparts of known 'optically dark' clouds of neutral hydrogen in the cluster, while a further four are widely separated extensions to previously known blue blobs. The properties of our new candidates are consistent with previous conclusions that blue blobs likely originated from ram pressure stripping events, however, their locations in velocity--projected cluster-centric radius phase-space imply that their parent galaxies are not on their first infall into the cluster. Through our ongoing follow-up program with HET we aim to confirm additional candidates, however, detailed understanding of the stellar populations and star formation histories of blue blobs will require JWST observations., Comment: Submitted to ApJ

We analyze 99 photometrically selected Little Red Dots (LRDs) at z ~ 4-8 in the GOODS fields, leveraging ultra-deep JADES NIRCam short-wavelength (SW) data. We examine the morphology of 30 LRDs; the remaining 69 are predominantly compact, i.e. are strongly dominated by sources < 400 pc in diameter and lack extended components even in stacked SW band images. Among the LRDs selected for morphological analysis, 50% show at least two distinct, associated sources or galaxy components, while the others appear as single sources with highly asymmetric structures. We find median stellar masses of log10(M*/Msun) = 9.07(-0.08)(+0.11) for pure stellar models with Av ~ 1.16(+0.11)(-0.21) mag, and log10(M*/Msun) = 9.67(+0.17)(-0.27) for models including AGNs, where Av ~ 2.74(+0.55)(-0.71) mag, consistent with recent results showing LRDs tend to have high stellar masses and dust content when fitted with AGN models. NIRSpec spectra are available for 15 sources, 6 of which fall within the morphological analysis sample and show multiple components. Among these 15, broad H-alpha emission is detected in 40%, with full-width half-maximum (FWHM) ranging from 1200 to 2900 km/s. One source exhibits broad H-beta emission with FWHM = 2000 +/- 500 km/s. Analysis of line ratios probing the interstellar medium (ISM) reveals a composite nature, indicating AGN activity combined with stellar processes. These findings suggest LRDs have a mixed nature, with AGN signatures in some cases linked to disturbed morphologies observed at rest-frame UV wavelengths., Comment: 28 pages, 9 figures, 3 tables. Submitted to ApJ. Comments are welcome!

The Hydrogen Lyman-alpha (Lya) line shows a large variety of shapes which is caused by factors at different scales, from the interstellar medium to the intergalactic medium. This work aims to provide a systematic inventory and classification of the spectral shapes of Lya emission lines to understand the general population of high-redshift Lya emitting galaxies (LAEs). Using the data from the MUSE eXtremely Deep Field, we select 477 galaxies at z=2.8-6.6. We develop a method to classify Lya emission lines in four spectral and three spatial categories, by combining a spectral analysis with a narrow-band image analysis. We measure spectral properties, such as the peak separation and the blue-to-total flux ratio. To ensure a robust sample for statistical analysis, we define a final unbiased sample of 206 galaxies by applying thresholds for signal-to-noise ratio, peak separation, and Lya luminosity. Our analysis reveals that between 32% and 51% of the galaxies exhibit double-peaked profiles. This fraction seems to evolve dependently with the Lya luminosity, while we don't notice a severe decrease of this fraction with redshift. A large amount of these double-peaked profiles shows blue-dominated spectra, suggesting unique gas dynamics and inflow characteristics in some high-redshift galaxies. Among the double-peaked galaxies, 4% are spurious detections. Around 20% out of the 477 sources of the parent sample lie in a complex environment, meaning there are other clumps or galaxies at the same redshift within a distance of 30kpc. Our results suggest that the Lya double-peak fraction may trace the evolution of IGM attenuation, but faintest galaxies are needed to be observed at high redshift. In addition, it is crucial to obtain secure systemic redshifts for LAEs to better constrain the nature of the double-peaks., Comment: 36 pages, 29 figures, 9 tables, accepted for publication in A&A

We present the large-scale distribution and kinematics of cold molecular gas across the compact galaxy group Stephan's Quintet, based on CO(2-1) observations performed with the Atacama Compact Array (ACA) and CO(1-0) data from the Combined Array for Research in Millimeter-wave Astronomy (CARMA). We find coherent structures of molecular gas associated with the galaxies and intra-group medium, which follow the distribution of warm H$_{2}$ previously seen with the James Webb Space Telescope (JWST). CO is associated with a ridge of shocked gas that crosses the galaxy group, and with a spiral arm of the intruding galaxy NGC7318b, which interacts with the intra-group medium along the ridge. Although the ridge contains widespread shocks, turbulent gas, and warm H$_{2}$, the CO lines are narrower than elsewhere in Stephan's Quintet (FWHM~25-65 km/s), indicative of settled cold gas. At a distinctly different velocity, CO is found in the active galaxy NGC7319 and Northern star-forming region SQ-A. A bridge of turbulent molecular gas connects NGC7319 with the ridge, covering a gap of ~700 km/s between these structures. The gas excitation ranges from $L'_{\rm CO(2-1)}$/$L'_{\rm CO(1-0)}$ ~ 0.3 in the bridge and SQ-A, to ~0.5 along the ridge, to near unity in the center of NGC7319. We also detect either a molecular outflow or turbulent molecular gas associated with the radio source in NGC7319. These ACA data are part of a program with the Atacama Large Millimeter/submillimeter Array (ALMA) and JWST to study molecular gas physics from the largest to the smallest scales across the intra-group medium of Stephan's Quintet., Comment: Accepted for publication in ApJ

The molecular gas in the interstellar medium (ISM) of star-forming galaxy populations exhibits diverse physical properties. We investigate the $^{12}$CO excitation of twelve dusty, luminous star-forming galaxies at $z \sim 2-4$ by combining observations of the $^{12}$CO from $J_{\rm up} = 1$ to $J_{\rm up} = 8$. The spectral line energy distribution (SLED) has a similar shape to NGC 253, M82, and local ULIRGs, with much stronger excitation than the Milky Way inner disc. By combining with resolved dust continuum sizes from high-resolution $870$-$\mu$m ALMA observations and dust mass measurements determined from multi-wavelength SED fitting, we measure the relationship between the $^{12}$CO SLED and probable physical drivers of excitation: star-formation efficiency, the average intensity of the radiation field $\langle U\rangle$, and the star-formation rate surface density. The primary driver of high-$J_{\rm up}$ $^{12}$CO excitation in star-forming galaxies is star-formation rate surface density. We use the ratio of the CO($3-2$) and CO($6-5$) line fluxes to infer the CO excitation in each source and find that the average ratios for our sample are elevated compared to observations of low-redshift, less actively star-forming galaxies and agree well with predictions from numerical models that relate the ISM excitation to the star-formation rate surface density. The significant scatter in the line ratios of a factor $\approx 3$ within our sample likely reflects intrinsic variations in the ISM properties which may be caused by other effects on the excitation of the molecular gas, such as cosmic ray ionization rates and mechanical heating through turbulence dissipation., Comment: Accepted for publication in MNRAS; 17 pages, 7 figures

An established method measuring the hydrogen ionisation fraction in shock excited gas is the BE99 method, which utilises six bright forbidden emission lines of [SII]6716, 6731, [NII]6548, 6583, and [OI]6300, 6363. We aim to extent the BE99 method by including more emission lines in the blue and near-infrared part of the spectrum ($\lambda$ = 3500-11000A), and considering higher hydrogen ionisation fractions ($x_e > 0.3$). In addition, we investigate how a non-equilibrium state of the gas and the presence of extinction influence the BE99 technique. We find that plenty additional emission line ratios can in principle be exploited as extended curves (or stripes) in the ($x_e, T_e$)-diagram. If the BE99 equilibrium is reached and extinction is corrected for, all stripes overlap in one location in the ($x_e, T_e$)-diagram indicating the existing gas parameters. The application to the Par Lup 3-4 outflow shows that the classical BE99 lines together with the [NI]5198+5200 lines do not meet in one locationin the ($x_e, T_e$)-diagram. This indicates that the gas parameters derived from the classical BE99 method are not fully consistent with other observed line ratios. A multi-line approach is necessary to determine the gas parameters. From our analysis we derive $n_e \sim$ 45 000 cm^-3 - 53000 cm^-3 , $T_e$ = 7600K - 8000K, and $x_e \sim$ 0.027 - 0.036 for the Par Lup 3-4 outflow. For the 244-440 Proplyd we were able to use the line ratios of [SII]6716+6731, [OI]6300+6363, and [OII]7320, 7330 in the BE99 diagram to estimate the ionisation fraction at knot E3 ($x_e = 0.58 \pm 0.05$). In conclusion, exploiting new line ratios reveals more insights on the state of the gas. Our analysis indicates, however, that a multi-line approach is more robust in deriving gas parameters, especially for high density gas., Comment: 12 pages, 10 figures, accepted to A&A

We study the structure, interstellar absorption, color-magnitude diagrams, kinematics, and dynamical state of embedded star clusters in the star-forming region associated with the giant molecular cloud G174+2.5. Our investigation is based on photometric data from the UKIDSS Galactic Plane Survey catalog and astrometric data from the Gaia DR3 catalogs. First, we recover all the known embedded clusters and candidate clusters in the region using surface density maps. Then, for the detected clusters, we determine their general parameters: the center positions, radii, number of stars, and reddening. To evaluate the reddening, we use both the NICEST algorithm and the Q-method. Both methods produce consistent extinction maps in the regions of the four studied clusters. However, the Q-method yields a much smaller color scatter in the CMD. For four clusters in particular (S235~North-West, S235~A-B-C, S235~Central, and S235~East1+East2), we were able to compute individual membership probabilities, the cluster distances, the cluster masses, and their average proper motions. By building on these results, we have studied the clusters' kinematics and dynamics. Moreover, we estimate the mass of the gas component and the star formation efficiency (SFE) in the regions of these four clusters. Finally, we provide an estimate of the total energy of the stellar and gas components in the area of these four clusters to determine whether the clusters are bound (here we consider a `cluster' as the system `stars + gas'). The gravitational bound strongly depends on the region for which we estimate the gas mass. If we consider the mass of the entire cloud, all these four clusters turn out to be bound., Comment: 30 pages, 16 figures, accepted for publication in the Astronomical Journal

Context. Quasar outflows are key players in the feedback processes that influence the evolution of galaxies and the intergalactic medium. The chemical abundance of these outflows provides crucial insights into their origin and impact. Aims. To determine the absolute abundances of nitrogen and sulfur and the physical conditions of the outflow seen in quasar 3C298. Methods. We analyze archival spectral data from the Hubble Space Telescope (HST) for 3C298. We measure Ionic column densities from the absorption troughs and compare the results to photoionization predictions made by the Cloudy code for three different spectral energy distributions (SED), including MF87, UVsoft, and HE0238 SEDs. We also calculate the ionic column densities of excited and ground states of N iii to estimate the electron number density and location of the outflow using the Chianti atomic database. Results. The MF87, UVsoft, and HE0238 SEDs yield nitrogen and sulfur abundances at super-solar, solar, and sub-solar values, respectively, with a spread of 0.4 to 3 times solar. Additionally, we determined an electron number density of log(ne) greater than 3.3 cm-3, with the outflow possibly extending up to a maximum distance of 2.8 kpc. Conclusions. Our results indicate solar metallicity within a 60 percent uncertainty range, driven by variations in the chosen SED and photoionization models. This study underscores the importance of SEDs impact on determining chemical abundances in quasars outflows. These findings highlight the necessity of considering a wider range of possible abundances, spanning from sub solar to super solar values., Comment: Accepted for publication in A&A- 11 pages, 10 figures

We investigate how stellar feedback from the first stars (Population III) distributes metals through the interstellar and intergalactic medium using the star-by-star cosmological hydrodynamics simulation, Aeos. We find that energy injected from the supernovae of the first stars is enough to expel a majority of gas and injected metals beyond the virial radius of halos with mass $M_* \lesssim 10^7$ M$_\odot$, regardless of the number of supernovae. This prevents self-enrichment and results in a non-monotonic increase in metallicity at early times. Most minihalos ($M \gtrsim 10^5 \, \rm M_\odot$) do not retain significant fractions of the yields produced within their virial radii until they have grown to halo masses of $M \gtrsim 10^7 \, \rm M_\odot$. The loss of metals to regions well beyond the virial radius delays the onset of enriched star formation and extends the period that Population III star formation can persist. We also explore the contributions of different nucleosynthetic channels to 10 individual elements. On the timescale of the simulation (lowest redshift $z=14.3$), enrichment is dominated by core-collapse supernovae for all elements, but with a significant contribution from asymptotic giant branch winds to the s-process elements, which are normally thought to only be important at late times. In this work, we establish important mechanisms for early chemical enrichment which allows us to apply Aeos in later epochs to trace the evolution of enrichment during the complete transition from Population III to Population II stars., Comment: 16 pages, 12 figures, 1 animated figure (see source files), submitted to ApJ

Massive black holes (MBHs) exist in the Galactic center (GC) and other nearby galactic nuclei. As natural outcome of galaxy mergers, some MBHs may have a black hole (BH) companion. In this paper, assuming that the MBHs in the GC and some nearby galaxies are in binaries with orbital periods ranging from months to years (gravitational-wave frequency $\sim1-100$\,nHz), we investigate the detectability of gravitational-waves from these binary MBHs (BBHs) and constraints on the parameter space for the existence of BBHs in the GC, LMC, M31, M32, and M87, that may be obtained by current/future pulsar timing array (PTA) observations. We find that a BBH in the GC, if any, can be revealed by the Square Kilometer Array PTA (SKA-PTA) if its mass ratio $q\gtrsim10^{-4}-10^{-3}$ and semimajor axis $a\sim20-10^3$\,AU. The existence of a BH companion of the MBH can be revealed by SKA-PTA with $\sim20$-year observations in M31 if $q\gtrsim10^{-4}$ and $a\sim10^2-10^4$\,AU or in M87 if $q\gtrsim10^{-5}$ and $a\sim10^3-2\times10^4$\,AU, but not in LMC and M32 if $q\ll1$. If a number of milli-second stable pulsars with distances $\lesssim0.1-1$\,pc away from the central MBH in the GC, LMC, M32, or M31, can be detected in future and applied to PTAs, the BH companion with mass even down to $\sim100M_\odot$, close to stellar masses, can be revealed by such PTAs. Future PTAs are expected to provide an independent way to reveal BBHs and low-mass MBH companions in the GC and nearby galaxies, improving our understandings of the formation and evolution of MBHs and galaxies., Comment: 20 pages, 13 figures, accepted for publication in ApJ

Massive stars lead the ionization and mechanical feedback within young star-forming regions. The Large Magellanic Cloud (LMC) is an ideal galaxy for studying individual massive stars and quantifying their feedback contribution to the environment. We analyze eight exemplary targets in LMC N11 B from the Hubble UV Legacy Library of Young Stars as Essential Standards (ULLYSES) program, using novel spectra from HST (COS and STIS) in the UV, and from VLT (X-shooter) in the optical. We model the spectra of early to late O-type stars by using state-of-the-art PoWR atmosphere models. We determine the stellar and wind parameters (e.g., $T_\star$, $\log g$, $L_{\star}$, $\dot{M}$, $v_\infty$) of the analyzed objects, chemical abundances (C, N, O), ionizing and mechanical feedback ($Q_\mathrm{H}$, $Q_\mathrm{He{\small{I}}}$, $Q_\mathrm{He{\small{II}}}$, $L_\mathrm{mec}$) and X-rays. We report ages of $2-4.5$ Myr and masses of $30-60$ $M_\odot$ for the analyzed stars in N11 B, consistent with a scenario of sequential star formation. We note that the observed wind-momentum luminosity relation is consistent with theoretical predictions. We detect nitrogen enrichment in most of the stars, up to a factor of seven. However, we do not find a correlation between nitrogen enrichment and projected rotational velocity. Finally, based on their spectral type, we estimate the total ionizing photons injected from the O-type stars in N11 B into its environment. We report $\log$ ($\sum$ $Q_\mathrm{H}$)$=50.5$ ph s$^{-1}$, $\log$ ($\sum$ $Q_\mathrm{He{\small{I}}}$)$=49.6$ ph s$^{-1}$ and $\log$ ($\sum$ $Q_\mathrm{He{\small{II}}}$)$=44.4$ ph s$^{-1}$, consistent with the total ionizing budget in N11., Comment: Accepted in Astronomy and Astrophysics

In this paper, we evaluate the potential of multiband gravitational wave observations to constrain the properties of static dark matter spikes around intermediate-mass ratio inspirals. The influence of dark matter on the orbital evolution of the compact binary is incorporated as a correction to the inspiral Newtonian gravitational waveform. We show that the observations from the proposed space-based detector GWSat, sensitive within the deci-Hz frequency band, when combined with that of the third-generation ground-based detectors like the Einstein Telescope and Cosmic Explorer, will produce significantly improved error estimates for all parameters. In particular, our results demonstrate that the joint multiband approach substantially refines the bounds on the dark matter spike parameters-namely, the power-law index and spike density-by factors of approximately $10^6$ and $10^3$, respectively, compared to observations employing only third-generation gravitational wave detectors., Comment: 13 pages, 5 figures

Accurate calculations of starspot spectra are essential for multiple applications in astronomy. The current standard is to represent starspot spectra by spectra of stars that are cooler than the quiet star regions. This implies approximating a starspot as a non-magnetic 1D structure in radiative-convective equilibrium, parametrizing convective energy transport by mixing length theory. It is the inhibition of convection by the starspot magnetic field that is emulated by using a lower spot temperature relative to the quiet stellar regions. Here, we take a different approach avoiding the approximate treatment of convection and instead self-consistently accounting for the interaction between matter, radiation, and the magnetic field. We simulate spots on G2V, K0V, M0V stars with the 3D radiative magnetohydrodynamics code MURaM and calculate spectra ($R \approx 500$ from 250~nm to 6000~nm) using ray-by-ray radiative transfer with the MPS-ATLAS code. We find that the 1D models fail to return accurate umbral and penumbral spectra on K0V and M0V stars where convective and radiative transfer of energy is simultaneously important over a broad range of atmospheric heights rendering mixing length theory inaccurate. However, 1D models work well for G2V stars, where both radiation and convection significantly contribute to energy transfer only in a narrow region near the stellar surface. Quantitatively, the 1D approximation leads to errors longward of 500 nm of about 50\% for both umbral and penumbral flux contrast relative to quiet star regions on M0V stars, and less than 2\% (for umbrae) and 10\% (for penumbrae) for G2V stars., Comment: Accepted for publication in the Astrophysical Journal Letters

Studying the distribution and properties of ionised gas in outflows driven by AGN is crucial for understanding the feedback mechanisms at play in extragalactic environments. In this study, we explore the connection between ionised outflows traced by rest-frame UV absorption and optical emission lines in GS133, a Compton thick AGN at z = 3.47. We combine observations from the JWST NIRSpec Integral Field Spectrograph (IFS) with archival VLT VIMOS long-slit spectroscopic data, as part of the GA-NIFS project. We perform a multi-component kinematic decomposition of the UV and optical line profiles to derive the physical properties of the absorbing and emitting gas in GS133. Our kinematic decomposition reveals two distinct components in the optical lines. The first component likely traces a rotating disk with a dynamical mass of 2e10 Msun. The second component corresponds to a galaxy-wide, bi-conical outflow, with a velocity of 1000 km/s and an extension of 3 kpc. The UV absorption lines show two outflow components, with bulk velocities v_out = -900 km/s and -1900 km/s, respectively. This characterises GS133 as a mini-BAL system. Balmer absorption lines with similar velocities are tentatively detected in the NIRSpec spectrum. Both photoionisation models and outflow energetics suggest that the ejected absorbing gas is located at 1-10 kpc from the AGN. We use 3D gas kinematic modelling to infer the orientation of the [O III] bi-conical outflow, and find that a portion of the emitting gas resides along our line of sight, suggesting that [O III] and absorbing gas clouds are partially mixed in the outflow. The derived mass-loading factor (i.e. the mass outflow rate divided by the SFR) of 1-10, and the kinetic coupling efficiency (i.e. the kinetic power divided by LAGN) of 0.1-1% per cent suggest that the outflow in GS133 provides significant feedback on galactic scales., Comment: 18 pages, 13 figures, submitted to Astronomy & Astrophysics. Comments are welcome

Numerical QUMOND-simulations of star clusters orbiting in a Galactic disk potential show that the leading tidal arm of open star clusters contains tendentially more members than the trailing arm. However, these type of simulations are performed by solving the field-equations of QUMOND and already become non-practical for star cluster masses at around 5000 Msun. Nearby star clusters have masses of 1000 Msun or ~1000 particles and less/fewer and can currently not be simulated reliably in field-theoretical formulations of MOND. In order to handle particle numbers below the QUMOND-limit the star cluster is simulated in Milgrom-law dynamics (MLD): Milgrom's law is postulated to be valid for discrete systems in vectorial form. In order to suppress the Newtonisation of compact subsystems in the star cluster the gravitational force is softened below particle distances of 0.001 pc ~206 AU. Thus, MLD can only be considered as an approximation of a full MOND-theoretical description of discrete systems which are internally in the MOND regime. The MLD equations of motion are integrated by the standard Hermite scheme generally applied to Newtonian N-body systems, which is extended to solve for the accelerations and jerks associated with Milgrom's law. It is found that the tidal tails of a low-mass star cluster are populated asymmetrically in the MLD-treatment, very similar to the QUMOND simulations of the higher-mass star clusters. In the MLD-simulations the leading tail hosts up to twice as many members than the trailing arm and the low-mass open star cluster dissolves approximately 25% faster than in the respective Newtonian case., Comment: accepted for publication in A&A

We present radial profiles of luminosity-weighted age, $age_L$, and $\Delta \Sigma_{SFR}$ for various populations of high- and low- mass central and satellite galaxies in the TNG100 cosmological simulation. Using these profiles, we investigate the impact of intrinsic and environmental factors on the radial distribution of star formation. For both central galaxies and satellites, we investigate the effects of black hole mass, cumulative AGN feedback energy, morphology, halo mass, and local galaxy overdensity on the profiles. In addition, we investigate the dependence of radial profiles of the satellite galaxies as a function of the redshifts at which they joined their hosts, as well as the net change in star-forming gas mass since the satellites joined their host. We find that high-mass ($M_*>10^{10.5} M_{\odot}$) central and satellite galaxies show evidence of inside-out quenching driven by AGN feedback. Effects from environmental processes only become apparent in averaged profiles at extreme halo masses and local overdensities. We find that the dominant quenching process for low-mass galaxies ($M_*<10^{10} M_{\odot}$) is environmental, generally occurring at low halo mass and high local galaxy overdensity for low-mass central galaxies and at high host halo masses for low-mass satellite galaxies. Overall, we find that environmental processes generally drive quenching from the outside-in., Comment: submitted to ApJ, revised in response to reviewer's comments 38 pages, 16 figures in main text; 3 appendices with 20 pages, 17 figures

{We present the radial metallicity gradients within the Galactic thin disc population through main-sequence stars selected on the chemical plane using GALAH DR3 accompanied with Gaia DR3 astrometric data. The [Fe/H], [$\alpha$/Fe] and [Mg/H] radial gradients are estimated for guiding radius as $-0.074\pm 0.006$, $+0.004\pm0.002$, $-0.074\pm0.006$ dex kpc$^{-1}$ and for the traceback early orbital radius as $-0.040\pm0.002$, $+0.003\pm 0.001$, $-0.039\pm 0.002$ dex kpc$^{-1}$ for 66,545 thin-disc stars, respectively. Alteration of the chemical structure within the Galactic disc caused by the radial orbital variations complicates results for the radial metallicity gradient. The effect of radial orbital variations on the metallicity gradients as a function on time indicates the following results: (i) The presence of a gradient along the disc throughout the time for which the model provides similar prediction, (ii) the radial orbital variations becomes more pronounced with the age of the stellar population and (iii) the effect of radial orbital variations on the metallicity gradients is minimal. The effect of radial orbital variations is found to be at most 6\% which does not statistically affect the radial gradient results. These findings contribute to a better understanding of the chemical evolution within the Galactic disc and provide an important basis for further research., Comment: 19 pages, including 15 figures & 1 table, accepted for publication in Astronomische Nachrichten. We would like to especially thank the editorial board for their contribution to the publication of the study

CHANCES, the CHileAN Cluster galaxy Evolution Survey, will study the evolution of galaxies in and around ${\sim}$150 massive galaxy clusters, from the local universe out to z=0.45. CHANCES will use the new 4MOST Spectroscopic Survey Facility on the VISTA 4m telescope to obtain spectra for ${\sim}$500,000 galaxies with magnitudes $r_\mathrm{AB} < 20.5$, providing comprehensive spectroscopic coverage of each cluster out to $5r_{200}$. Its wide and deep scope will trace massive and dwarf galaxies from the surrounding filaments and groups to the cores of galaxy clusters, enabling the study of galaxy pre-processing and the role of the evolving environment on galaxy evolution. In this paper we present and characterize the sample of clusters and superclusters to be targeted by CHANCES. We used literature catalogues based on X-ray emission and Sunyaev-Zel'dovich effect to define the cluster sample in a homogeneous way, with attention to cluster mass and redshift, as well as the availability of ancillary data. We calibrated literature mass estimates from various surveys against each other and provide an initial mass estimate for each cluster, which we used to define the radial extent of the 4MOST coverage. We also present an initial assessment of the structure surrounding these clusters based on the redMaPPer red-sequence algorithm as a preview of some of the science CHANCES will enable., Comment: 11 pages, 9 figures, plus references and appendix containing catalog tables, submitted to A&A

We report the discovery of five galaxy candidates at redshifts between $15.9

We present a detailed study of the large-scale shock front in Stephan's Quintet, a byproduct of past and ongoing interactions. Using integral-field spectroscopy from the new William Herschel Telescope Enhanced Area Velocity Explorer (WEAVE), recent 144 MHz observations from the LOFAR Two-metre Sky Survey (LoTSS), and archival data from the Very Large Array and James Webb Space Telescope (JWST), we obtain new measurements of key shock properties and determine its impact on the system. Harnessing the WEAVE large integral field unit's (LIFU) field of view (90 $\times$ 78 arcsec$^{2}$), spectral resolution ($R\sim2500$) and continuous wavelength coverage across the optical band, we perform robust emission line modeling and dynamically locate the shock within the multi-phase intergalactic medium (IGM) with higher precision than previously possible. The shocking of the cold gas phase is hypersonic, and comparisons with shock models show that it can readily account for the observed emission line ratios. In contrast, we demonstrate that the shock is relatively weak in the hot plasma visible in X-rays (with Mach number of $\mathcal{M} \sim 2 - 4$), making it inefficient at producing the relativistic particles needed to explain the observed synchrotron emission. Instead, we propose that it has led to an adiabatic compression of the medium, which has increased the radio luminosity ten-fold. Comparison of the Balmer line-derived extinction map with the molecular gas and hot dust observed with JWST suggests that pre-existing dust may have survived the collision, allowing the condensation of H$_{2}$ - a key channel for dissipating the shock energy., Comment: 23 pages, 15 figures, accepted for publication in MNRAS

The intracluster medium (ICM) in the far outskirts (r $>$ 2-3 R$_{200}$) of galaxy clusters interfaces with the intergalactic medium (IGM) and is theorized to comprise diffuse, multiphase gas. This medium may hold vital clues to clusters' thermodynamic evolution and far-reaching impacts on infalling, future cluster galaxies. The diffuse outskirts of clusters are well-suited for quasar absorption line observations, capable of detecting gas to extremely low column densities. We analyze 18 QSO spectra observed with the Cosmic Origins Spectrograph aboard the Hubble Space Telescope whose lines of sight trace the gaseous environments of 26 galaxy clusters from within R$_{200}$ to 6 R$_{200}$ in projection. We measure the dN/dz and covering fraction of H I and O VI associated with the foreground clusters as a function of normalized impact parameter. We find the dN/dz for H I is consistent with the IGM field value for all impact parameter bins, with an intriguing slight elevation between 2 and 3 R$_{200}$. The dN/dz for O VI is also consistent with the field value (within 3$\sigma$) for all impact parameter bins, with potential elevations in dN/dz both within 1-2 R$_{200}$ and beyond 4 R$_{200}$ at $>2\sigma$. We propose physical scenarios that may give rise to these tentative excesses, such as a buildup of neutral gas at the outer accretion shock front and a signature of the warm-hot IGM. We do not find a systematic excess of potentially associated galaxies near the sightlines where O VI is detected; thus, the detected O VI does not have a clear circumgalactic origin., Comment: Submitted to ApJ, 33 pages, 11 primary figures, 2 tables, 25 supplementary figures

Green Peas (GPs) and blueberry galaxies (BBs) are thought to be local analogs ($z<$0.1) of high redshift Ly$\alpha$ emitters. H I study of these can help us understand the star formation in the primordial Universe. In this Letter, we present the results of H I 21 cm study of 28 high specific star formation rate (sSFR $\gtrsim$10$^{-8}$ yr$^{-1}$) BBs at $z\lesssim$0.05 with the Five-hundred-meter Aperture Spherical radio Telescope. We report significant H I detection towards two BBs namely J1026+0426 and J1132+0809, and discuss possible H I contribution from neighboring galaxies. The median 3$\sigma$ upper limit of $\sim$2.0$\times$10$^{8}$ M$_{\odot}$ was obtained on H I mass for galaxies with nondetections. We find BBs tend to have lower H I-to-stellar mass ratio or gas fraction ($f_{\rm HI}$) than expected from $f_{\rm HI}$-sSFR and $f_{\rm HI}$-$M_{\ast}$ relations for main-sequence galaxies. The BBs also have a median 3$\sigma$ upper limit on H I gas depletion time scale ($\tau_{\rm HI}$) $\sim$0.5 Gyr, about 1 order of magnitude lower than $\tau_{\rm HI}$ for local main-sequence galaxies. We find a significantly low H I detection rate of 2/28 (7.1$^{+9.4}_{-4.6}$ \%) towards these galaxies, which is similar to previous H I studies of low redshift GPs of high ionization parameter indicator, O32 $\equiv$O[{\sc iii}]$\lambda$5007/O[{\sc ii}]$\lambda$3727 ratios $\gtrsim$10., Comment: 10 pages, 4 figures, 1 table, accepted for publication in the Astrophysical Journal Letters (ApJL), this is proof corrected version

The measurement of magnetic fields in cosmic web filaments can be used to reveal the magnetogenesis of the Universe. In previous work, we produced first estimates of the field strength and its redshift evolution using the Faraday Rotation Measure (RM) catalogue of extragalactic background sources at low frequency obtained with LOFAR observations. Here we refine our analysis by selecting sources with low Galactic RM, which reduces its residual contamination. We also conduct a comprehensive analysis of the different contributions to the extragalactic RMs along the line of sight, and confirm that they are dominated by the cosmic filaments component, with only 21 percent originating in galaxy clusters and the circumgalactic medium (CGM) of galaxies. We find a possible hint of a shock at the virial radius of massive galaxies. We also find that the fractional polarization of background sources might be a valuable CGM tracer. The newly selected RMs have a steeper evolution with redshift than previously found. The field strength in filaments ($B_f$) and its evolution are estimated assuming $B_f$ evolves as a power-law $B_f=B_{f,0}\,(1+z)^\alpha$. Our analysis finds an average strength at $z=0$ of $B_{f,0} =11$--15~nG, with an error of 4 nG, and a slope $\alpha=2.3$--$2.6 \pm 0.5$, which is steeper than what we previously found. The comoving field has a slope of $\beta=$ [0.3, 0.6$]\pm 0.5$ that is consistent with being invariant with redshift. Primordial magnetogenesis scenarios are favoured by our data, together with a sub-dominant astrophysical-origin RM component increasing with redshift., Comment: Accepted for publication on A&A

We investigate observations of circumbinary disks (CBD), to find evidence for an equilibrium eccentricity predicted by current binary accretion theory. Although stellar binary demographics in the Milky Way show no evidence for a preferred eccentricity for binary systems, we show that actively accreting systems lie on a predicted equilibrium eccentricity curve. We constrain our sample to only systems that have well defined orbital parameters (e.g,. eccentricity, mass-ratio, inclination angle). We find observations are consistent with theory for stellar binaries that are aligned with the disk and that are separated enough that tidal circularization is negligible. This suggests that eccentricity in these systems evolves after the dissipation of the CBD, given the flat eccentricity distribution of binary systems in the Milky Way., Comment: 5 pages, 2 figures

Galaxies come in different sizes and morphologies, and these differences are thought to correlate with properties of their underlying dark matter halos. However, identifying the specific halo property that controls the galaxy size is a challenging task, especially because most halo properties are dependent on one another. In this work, we demonstrate this challenge by studying how the galaxy-halo size relations impact the galaxy clustering signals. We investigate the reason that a simple linear relation model, which prescribes that the galaxy size is linearly proportional to the dark matter halo's virial radius, can still produce clustering signals that match the observational data reasonably well. We find that this simple linear relation model for galaxy sizes, when combined with the subhalo abundance matching technique, introduces an implicit dependence on the halo formation history. As a result, the effect of halo assembly bias enters the resulting galaxy clustering, especially at lower stellar masses, producing a clustering signal that resembles the observed one. At higher stellar masses, the effect of halo assembly bias weakens and is partially canceled out by the effect of halo bias, and the clustering of large and small galaxies becomes more similar. Our study confirms that the information of halo formation history must play a role in determining galaxy sizes to match the observed clustering signals, but also highlights the challenge of identifying a particular halo property that controls galaxy sizes through constraints from galaxy clustering alone., Comment: 9 pages, 7 figures. To be submitted to the Open Journal of Astrophysics

Pulsar timing array observations have found evidence for an isotropic gravitational wave background with the Hellings-Downs angular correlations, expected from general relativity. This interpretation hinges on the measured shape of the angular correlations, which is predominately quadrupolar under general relativity. Here we explore a more flexible parameterization: we expand the angular correlations into a sum of Legendre polynomials and use a Bayesian analysis to constrain their coefficients with the 15-year pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). When including Legendre polynomials with multipoles $\ell \geq 2$, we only find a significant signal in the quadrupole with an amplitude consistent with general relativity and non-zero at the $\sim 95\%$ confidence level and a Bayes factor of 200. When we include multipoles $\ell \leq 1$, the Bayes factor evidence for quadrupole correlations decreases by more than an order of magnitude due to evidence for a monopolar signal at approximately 4 nHz which has also been noted in previous analyses of the NANOGrav 15-year data. Further work needs to be done in order to better characterize the properties of this monopolar signal and its effect on the evidence for quadrupolar angular correlations., Comment: 15 pages, 6 figures

X-ray quasi-periodic eruptions (QPEs) are intense soft X-ray bursts from the nuclei of nearby low-mass galaxies typically lasting about one hour and repeating every few. Their physical origin remains debated, although so-called impacts models in which a secondary orbiting body pierces through the accretion disc around the primary supermassive black hole (SMBH) in an extreme mass-ratio inspiral (EMRI) system are considered promising. In this work, we study the QPE timing properties of GSN 069, the first galactic nucleus in which QPEs were identified, primarily focusing on Observed minus Calculated (O-C) diagrams. The O-C data in GSN 069 are consistent with a super-orbital modulation on tens of days whose properties do not comply with the impacts model. We suggest that rigid precession of a misaligned accretion disc or, alternatively, the presence of a second SMBH forming a sub-milliparsec binary with the inner EMRI is needed to reconcile the model with the data. In both cases, the quiescent accretion disc emission should also be modulated on similar timescales. Current X-ray monitoring indicates that this might be the case, although a longer baseline of higher-cadence observations is needed to confirm the tentative X-ray flux periodicity on firm statistical grounds. Future dedicated monitoring campaigns will be crucial to test the overall impacts plus modulation model in GSN 069, and to distinguish between the two proposed modulating scenarios. If our interpretation is correct, QPEs in GSN 069 represent the first electromagnetic detection of a short-period EMRI system in an external galaxy, opening the way to future multi-messenger astronomical observations. [abridged], Comment: Accepted for publication in A&A

The Surface Photometry and Accurate Rotation Curves (SPARC) database has provided the community with mass models for 175 nearby galaxies, allowing different research teams to test different dark matter models, galaxy evolution models, and modified gravity theories. Extensive tests, however, are hampered by the somewhat heterogeneous nature of the HI rotation curves and the limited sample size of SPARC. To overcome these limitations, we are working on BIG-SPARC, a new database that consists of about 4000 galaxies with HI datacubes from public telescope archives (APERTIF, ASKAP, ATCA, GMRT, MeerKAT, VLA, and WSRT) and near infrared photometry from WISE. For these galaxies, we will provide homogeneously derived HI rotation curves, surface brightness profiles, and mass models. BIG-SPARC is expected to increase the size of its predecessor by a factor of more than 20. This is a necessary step to prepare for the additional order of magnitude increase in sample size expected from ongoing and future HI surveys with the Square Kilometre Array (SKA) and its pathfinders, Comment: To be published in the proceedings of IAU Symposium 392 "Neutral Hydrogen in and around Galaxies in the SKA Era", 6 pages, 3 figures, 1 table

We introduce TPCNet, a neural network predictor that combines Convolutional and Transformer architectures with Positional encodings, for neutral atomic hydrogen (HI) spectral analysis. Trained on synthetic datasets, our models predict cold neutral gas fraction ($f_\text{CNM}$) and HI opacity correction factor ($R_\text{HI}$) from emission spectra based on the learned relationships between the desired output parameters and observables (optically-thin column density and peak brightness). As a follow-up to Murray et al. (2020)'s shallow Convolutional Neural Network (CNN), we construct deep CNN models and compare them to TPCNet models. TPCNet outperforms deep CNNs, achieving a 10% average increase in testing accuracy, algorithmic (training) stability, and convergence speed. Our findings highlight the robustness of the proposed model with sinusoidal positional encoding applied directly to the spectral input, addressing perturbations in training dataset shuffling and convolutional network weight initializations. Higher spectral resolutions with increased spectral channels offer advantages, albeit with increased training time. Diverse synthetic datasets enhance model performance and generalization, as demonstrated by producing $f_\text{CNM}$ and $R_\text{HI}$ values consistent with evaluation ground truths. Applications of TPCNet to observed emission data reveal strong agreement between the predictions and Gaussian decomposition-based estimates (from emission and absorption surveys), emphasizing its potential in HI spectral analysis., Comment: This paper has 27 pages, 27 figures and two tables. The work has been accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal

While WISE is the largest, best quality infrared all-sky survey to date, a smaller coverage mission, Spitzer, was designed to have better sensitivity and spatial resolution at similar wavelengths. Confusion and contamination in WISE data result in discrepancies between them. We present a novel approach to work with WISE measurements with the goal of maintaining both its high coverage and vast amount of data while taking full advantage of the higher sensitivity and spatial resolution of Spitzer. We have applied machine learning (ML) techniques to a complete WISE data sample of open cluster members, using a training set of paired data from high-quality Spitzer Enhanced Imaging Products (SEIP), MIPS and IRAC, and allWISE catalogs, W1 (3.4 {\mu}m) to W4 (22 {\mu}m) bands. We have tested several ML regression models with the aim of predicting mid-infrared fluxes at MIPS1 (24 {\mu}m) and IRAC4 (8 {\mu}m) bands from WISE fluxes and quality flags. In addition, to improve the prediction quality, we have implemented feature selection techniques to remove irrelevant WISE variables. We have notably enhanced WISE detection capabilities, mostly at lowest magnitudes, which previously showed the largest discrepancies with Spitzer. In our particular case, extremely randomized trees was found to be the best algorithm to predict mid-infrared fluxes from WISE variables. We have tested our results in the SED of members of IC 348. We show discrepancies in the measurements of Spitzer and WISE and demonstrate the good concordance of our predicted fluxes with the real ones. ML is a fast and powerful tool that can be used to find hidden relationships between datasets, as the ones that exist between WISE and Spitzer fluxes. We believe this approach could be employed for other samples from the allWISE catalog with SEIP positional counterparts, and in other astrophysical studies with analogous discrepancies., Comment: 13 pages, 10 figures

We present a new iterative deblending method to separate the host galaxy (HG) and their Active Galactic Nuclei (AGN) emission with the use of Integral Field spectroscopic (IFS) data. The method decomposes the resolved HG emission from the unresolved AGN emission by modelling the two-dimensional surface brightness (SB) profile of the point-spread function (PSF) and the two-dimensional SB HG continuum simultaneously per each monochromatic slide. Our method does not require any prior information about the observed SB profile or a detailed fitting of the PSF, making it ideal for the automatic analysis of large galaxy samples. In this work, we test the quality of our method, its advantages, and its disadvantages. We test our method by using a set of IFS mock data cubes to quantify the reliability of our deblending process and further compare our method with the {\sc QDeblend3D} analysis tool. Furthermore, we applied our method to three data cubes selected from the MaNGA survey according to the dominance of either its HG or its AGN. We show that our deblending method is capable of disengaging the bright, nonresolved AGN emission from the HG continuum and its narrow emission lines. However, the decoupling depends on how well the IFS spatially resolves the PSF, and on the relative flux intensity of the HG-AGN. Therefore, the method is ideal for disentangling the bright-flux contribution from AGN-dominated spectra., Comment: 25 pages, 18 figures, accepted for publication on the MNRAS

Understanding the interstellar and potentially circumstellar extinction in the sight lines of classical T Tauri stars is an important ingredient for constructing reliable spectral energy distributions, which catalyze protoplanetary disk chemistry, for example. Therefore, some attempts of measuring $A_{V}$ toward individual stars have been made using partly different wavelength regimes and different underlying assumptions. We used strong lines of Ly{\alpha} fluorescent H2 and derived the extinction based on the assumption of optically thin transitions. We investigated a sample of 72 classical T Tauri stars observed with the Hubble Space Telescope in the framework of the ULLYSES program. We computed $A_{V}$ and $R_{V}$ values for the 34 objects with sufficient data quality and an additionally $A_{V}$ value for the canonical $R_{V}$ = 3.1 value. Our results agree largely with values obtained from optical data. Moreover, we confirm the degeneracy between $A_{V}$ and $R_{V}$ and present possibilities to break this. Finally, we discuss whether the assumption of optical thin lines is valid., Comment: 12 pages, 6 figures, accepted to A&A

This study aims to investigate whether the environment and the nuclear activity of a particular galaxy influence the ageing and quenching at the transition stage of the galaxy evolution using the volume-limited sample constructed from the twelve release of the Sloan Digital Sky Survey. To this end, the galaxies were classified into isolated and non-isolated environments and then each subsample was further classified according to their nuclear activity using the WHAN diagnostic diagram, and ageing diagram to obtain ageing and quenching galaxies. The ageing and quenching galaxies at the transition stage were selected for the rest of the analysis. Using the star formation rate and the $u-r$ colour-stellar mass diagrams, the study revealed a significant change of $0.03$ dex in slope and $0.30$ dex in intercept for ageing galaxies and an insignificant change of $0.02$ dex in slope and $0.12$ dex in intercept of the star formation main sequence between isolated and non-isolated quenching galaxies. Further, a more significant change in the number of ageing galaxies above, within and below the main sequence and the green valley was observed. On the other hand, an insignificant change in the number of quenching galaxies above, within and below the main sequence and the green valley was observed. The study concludes that ageing depends on the environment and the dependence is influenced by the nuclear activity of a particular galaxy while quenching does not depend on the environment and this independence is not influenced by the nuclear activity., Comment: 14 pages and 6 figures. arXiv admin note: text overlap with arXiv:2405.00481