Search

This paper addresses, for the first time, a key aspect of the phenomenology of Compact Symmetric Objects (CSOs) -- the characteristics of their radio spectra. We present a radio-spectrum description of a complete sample of high luminosity CSOs (CSO-2s), which shows that they exhibit the \textit{complete} range of spectral types, including flat-spectrum sources ($\alpha \ge -0.5$), steep-spectrum sources ($\alpha < -0.5$), and peaked-spectrum sources. We show that there is no clear correlation between spectral type and size, but there is a correlation between the high-frequency spectral index and both object type and size. We also show that, to avoid biasing the data and to understand the various classes of jetted-AGN involved, the complete range of spectral types should be included in studying the general phenomenology of compact jetted-AGN, and that complete samples must be used, selected over a wide range of frequencies. We discuss examples that demonstrate these points. We find that the high-frequency spectral indices of CSO-2s span $-1.3 <\alpha_{\rm hi} < -0.3$, and hence that radio spectral signatures cannot be used to discriminate definitively between CSO-2s, binary galactic nuclei, and millilensed objects, unless they have $\alpha_{\rm hi} >-0.3$.

We investigate the ability of the Euclid telescope to detect galaxy-scale gravitational lenses. To do so, we perform a systematic visual inspection of the $0.7\,\rm{deg}^2$ Euclid ERO data towards the Perseus cluster using both the high-resolution VIS $I_{\scriptscriptstyle\rm E}$ band, and the lower resolution NISP bands. We inspect every extended source brighter than magnitude $23$ in $I_{\scriptscriptstyle\rm E}$ with $41$ expert human classifiers. This amounts to $12\,086$ stamps of $10^{\prime\prime}\,\times\,10^{\prime\prime}$. We find $3$ grade A and $13$ grade B candidates. We assess the validity of these $16$ candidates by modelling them and checking that they are consistent with a single source lensed by a plausible mass distribution. Five of the candidates pass this check, five others are rejected by the modelling and six are inconclusive. Extrapolating from the five successfully modelled candidates, we infer that the full $14\,000\,{\rm deg}^2$ of the Euclid Wide Survey should contain $100\,000^{+70\,000}_{-30\,000}$ galaxy-galaxy lenses that are both discoverable through visual inspection and have valid lens models. This is consistent with theoretical forecasts of $170\,000$ discoverable galaxy-galaxy lenses in Euclid. Our five modelled lenses have Einstein radii in the range $0.\!\!^{\prime\prime}68\,<\,\theta_\mathrm{E}\,<1.\!\!^{\prime\prime}24$, but their Einstein radius distribution is on the higher side when compared to theoretical forecasts. This suggests that our methodology is likely missing small Einstein radius systems. Whilst it is implausible to visually inspect the full Euclid data set, our results corroborate the promise that Euclid will ultimately deliver a sample of around $10^5$ galaxy-scale lenses., Comment: 21 pages, 20 figures, submitted to A&A

Owens Valley Radio Observatory (OVRO) observations of supermassive black hole binary (SMBHB) candidate PKS~2131$-$021 revealed, for the first time, six likely characteristics of the phenomenology exhibited by SMBHB in blazars, of which the most unexpected and critical is sinusoidal flux density variations. We have now identified a second blazar, PKS~J0805$-$0111, showing significant sinusoidal variations, with an observed period that translates to $1.422 \pm 0.005$ yr in the rest frame of the $z = 1.388$ object. We generate $10^6$ simulated light curves to reproduce the radio variability characteristics of PKS~J0805$-$0111, and show that the global probability, considering the \textit{look-elsewhere effect}, indicates that the observed periodicity can be attributed to the red noise tail of the power spectral density, with a $p_0$ value of $7.8 \times 10^{-5}$ (i.e. 3.78$\sigma$). PKS J0805$-$0111 displays all six characteristics observed in PKS 2131$-$021. Taking into account the well-defined OVRO sample size, the false positive probability $\sim 0.22$, but the rare behavior makes this a strong SMBHB candidate. The discovery of a second SMBHB candidate exhibiting these rare characteristics reveals that PKS~2131$-$021 is not a unique, isolated case. With these two strong cases we are clearly seeing only the tip of the iceberg. We estimate that the number of SMBHB candidates amongst blazars $\sim$ 1 in 100.

The Carbon monOxide Mapping Array Project (COMAP) Pathfinder survey continues to demonstrate the feasibility of line-intensity mapping using high-redshift carbon monoxide (CO) line emission traced at cosmological scales. The latest COMAP Pathfinder power spectrum analysis is based on observations through the end of Season 2, covering the first three years of Pathfinder operations. We use our latest constraints on the CO(1-0) line-intensity power spectrum at $z\sim3$ to update corresponding constraints on the cosmological clustering of CO line emission and thus the cosmic molecular gas content at a key epoch of galaxy assembly. We first mirror the COMAP Early Science interpretation, considering how Season 2 results translate to limits on the shot noise power of CO fluctuations and the bias of CO emission as a tracer of the underlying dark matter distribution. The COMAP Season 2 results place the most stringent limits on the CO tracer bias to date, at $\langle{Tb}\rangle<4.8$ $\mu$K. These limits narrow the model space significantly compared to previous CO line-intensity mapping results while maintaining consistency with small-volume interferometric surveys of resolved line candidates. The results also express a weak preference for CO emission models used to guide fiducial forecasts from COMAP Early Science, including our data-driven priors. We also consider directly constraining a model of the halo-CO connection, and show qualitative hints of capturing the total contribution of faint CO emitters through the improved sensitivity of COMAP data. With continued observations and matching improvements in analysis, the COMAP Pathfinder remains on track for a detection of cosmological clustering of CO emission., Comment: 9 pages + bibliography and appendices (13 pages total); 9 figures, 1 table; v2 reflects minor changes made for version submitted to A&A, with no changes to top-line results

Single-dish far-infrared (far-IR) and sub-millimetre (sub-mm) point source catalogues and their connections with catalogues at other wavelengths are of paramount importance. However, due to the large mismatch in spatial resolution, cross-matching galaxies at different wavelengths is challenging. This work aims to develop the next-generation deblended far-IR and sub-mm catalogues and present the first application in the COSMOS field. Our progressive deblending used the Bayesian probabilistic framework known as XID+. The deblending started from the Spitzer/MIPS 24 micron data, using an initial prior list composed of sources selected from the COSMOS2020 catalogue and radio catalogues from the VLA and the MeerKAT surveys, based on spectral energy distribution modelling which predicts fluxes of the known sources at the deblending wavelength. To speed up flux prediction, we made use of a neural network-based emulator. After deblending the 24 micron data, we proceeded to the Herschel PACS (100 & 160 micron) and SPIRE wavebands (250, 350 & 500 micron). Each time we constructed a tailor-made prior list based on the predicted fluxes of the known sources. Using simulated far-IR and sub-mm sky, we detailed the performance of our deblending pipeline. After validation with simulations, we then deblended the real observations from 24 to 500 micron and compared with blindly extracted catalogues and previous versions of deblended catalogues. As an additional test, we deblended the SCUBA-2 850 micron map and compared our deblended fluxes with ALMA measurements, which demonstrates a higher level of flux accuracy compared to previous results.We publicly release our XID+ deblended point source catalogues. These deblended long-wavelength data are crucial for studies such as deriving the fraction of dust-obscured star formation and better separation of quiescent galaxies from dusty star-forming galaxies., Comment: 23 pages, 30 figures, accepted for publication in A&A. Catalogues can be downloaded from https://hedam.lam.fr/HELP/dataproducts/dmu26/dmu26_XID+COSMOS2024/

Thin stellar streams originating from globular clusters are among the most sensitive tracers of low-mass dark-matter subhalos. Joint analysis of the entire population of stellar streams will place the most robust constraints on the dark-matter subhalo mass function, and therefore the nature of dark matter. Here we use a hierarchical model of globular cluster formation to forecast the total number, masses and radial distribution of dissolved globular cluster in Milky Way-like galaxies. Furthermore, we generate mock stellar streams from these progenitors' orbital histories taking into account the clusters' formation and accretion time, mass, and metallicity. Out of $\sim$10,000 clusters more massive than $10^4$ M$_{\odot}$, $\sim$9000 dissolved in the central bulge and are fully phase-mixed at the present, while the remaining $\sim$1000 survive as coherent stellar streams. This suggests that the current census of $\sim$80 globular cluster streams in the Milky Way is severely incomplete. Beyond 15 kpc from the Galactic center we are missing hundreds of streams, of which the vast majority are from accreted GCs. Deep Rubin photometry $(g\lesssim27.5)$ would be able to detect these streams, even the most distant ones beyond $> 75$ kpc. We also find that M31 will have an abundance of streams at galactocentric radii of 30-100 kpc. We conclude that future surveys will find a multitude of stellar streams from globular clusters which can be used for dark matter subhalo searches., Comment: 21 pages, 13 figures, submitted to AAS journals, comments welcome

Aims. This work aims to reproduce the time before or after a merger event of merging galaxies from the IllustrisTNG cosmological simulation using machine learning. Methods. Images of merging galaxies were created in the u, g, r, and i bands from IllustrisTNG. The merger times were determined using the time difference between the last simulation snapshot where the merging galaxies were tracked as two galaxies and the first snapshot where the merging galaxies were tracked as a single galaxy. This time was then further refined using simple gravity simulations. These data were then used to train a residual network (ResNet50), a Swin Transformer (Swin), a convolutional neural network (CNN), and an autoencoder (using a single latent neuron) to reproduce the merger time. The full latent space of the autoencoder was also studied to see if it reproduces the merger time better than the other methods. This was done by reducing the latent space dimensions using Isomap, linear discriminant analysis (LDA), neighbourhood components analysis, sparse random projection, truncated singular value decomposition and uniform manifold approximation and projection. Results. The CNN is the best of all the neural networks. The performance of the autoencoder was close to the CNN, with Swin close behind the autoencoder. ResNet50 performed the worst. The LDA dimensionality reduction performed the best of the six methods used. The exploration of the full latent space produced worse results than the single latent neuron of the autoencoder. For the test data set, we found a median error of 190 Myr, comparable to the time separation between snapshots in IllustrisTNG. Galaxies more than $\approx$ 625 Myr before a merger have poorly recovered merger times, as well as galaxies more than $\approx$ 125 Myr after a merger event., Comment: 16 pages, 20 figures, 1 appendix, Accepted for publication in Astronomy & Astrophysics

We present the first degree-scale tomography map of the dusty magnetized interstellar medium (ISM) from stellar polarimetry and distance measurements. We used the RoboPol polarimeter at Skinakas Observatory to conduct a survey of starlight polarization in a region of the sky of 4 square degrees. We propose a Bayesian method to decompose the stellar-polarization source field along the distance to invert the 3D volume occupied by the observed stars. We used it to obtain the first 3D map of the dusty magnetized ISM. Specifically, we produced a tomography map of the orientation of the plane-of-sky (POS) component of the magnetic field threading the diffuse, dusty regions responsible for the stellar polarization. For the targeted region centered on Galactic coordinates $(l,b) \approx (103.3^\circ, 22.3^\circ)$, we identified several ISM clouds. Most of the lines of sight intersect more than one cloud. A very nearby component was detected in the foreground of a dominant component from which most of the polarization signal comes. Farther clouds, with a distance of up to 2~kpc, were similarly detected. Some of them likely correspond to intermediate-velocity clouds seen in HI spectra in this region of the sky. We found that the orientation of the POS component of the magnetic field changes along distance for most of the lines of sight. Our study demonstrates that starlight polarization data coupled to distance measures have the power to reveal the great complexity of the dusty magnetized ISM in 3D and, in particular, to provide local measurements of the POS component of the magnetic field. This demonstrates that the inversion of large data volumes, as expected from the PASIPHAE survey, will provide the necessary means to move forward in the modeling of the Galactic magnetic field and of the dusty magnetized ISM as a contaminant in observations of the cosmic microwave background polarization., Comment: Preprint of the accepted version at Astronomy & Astrophysics. The 3D map obtained in this paper can be visualized online at https://pasiphae.science/visualization and a video featuring it is accessible at https://www.youtube.com/watch?v=dB_6J1zhmPI

We present the first systematic investigation of supervised scaling laws outside of an ImageNet-like context - on images of galaxies. We use 840k galaxy images and over 100M annotations by Galaxy Zoo volunteers, comparable in scale to Imagenet-1K. We find that adding annotated galaxy images provides a power law improvement in performance across all architectures and all tasks, while adding trainable parameters is effective only for some (typically more subjectively challenging) tasks. We then compare the downstream performance of finetuned models pretrained on either ImageNet-12k alone vs. additionally pretrained on our galaxy images. We achieve an average relative error rate reduction of 31% across 5 downstream tasks of scientific interest. Our finetuned models are more label-efficient and, unlike their ImageNet-12k-pretrained equivalents, often achieve linear transfer performance equal to that of end-to-end finetuning. We find relatively modest additional downstream benefits from scaling model size, implying that scaling alone is not sufficient to address our domain gap, and suggest that practitioners with qualitatively different images might benefit more from in-domain adaption followed by targeted downstream labelling., Comment: 10+6 pages, 12 figures. Appendix C2 based on arxiv:2206.11927. Code, demos, documentation at https://github.com/mwalmsley/zoobot

Various galaxy merger detection methods have been applied to diverse datasets. However, it is difficult to understand how they compare. We aim to benchmark the relative performance of machine learning (ML) merger detection methods. We explore six leading ML methods using three main datasets. The first one (the training data) consists of mock observations from the IllustrisTNG simulations and allows us to quantify the performance metrics of the detection methods. The second one consists of mock observations from the Horizon-AGN simulations, introduced to evaluate the performance of classifiers trained on different, but comparable data. The third one consists of real observations from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) survey. For the binary classification task (mergers vs. non-mergers), all methods perform reasonably well in the domain of the training data. At $0.1

Aims. In this work, we intend to examine how environment influences the merger fraction, from the low density field environment to higher density groups and clusters. We also aim to study how the properties of a group or cluster, as well as the position of a galaxy in the group or cluster, influences the merger fraction. Methods. We identified galaxy groups and clusters in the North Ecliptic Pole using a friends-of-friends algorithm and the local density. Once identified, we determined the central galaxies, group radii, velocity dispersions, and group masses of these groups and clusters. Merging systems were identified with a neural network as well as visually. With these, we examined how the merger fraction changes as the local density changes for all galaxies as well as how the merger fraction changes as the properties of the groups or clusters change. Results. We find that the merger fraction increases as local density increases and decreases as the velocity dispersion increases, as is often found in literature. A decrease in merger fraction as the group mass increases is also found. We also find groups with larger radii have higher merger fractions. The number of galaxies in a group does not influence the merger fraction. Conclusions. The decrease in merger fraction as group mass increases is a result of the link between group mass and velocity dispersion. Hence, this decrease of merger fraction with increasing mass is a result of the decrease of merger fraction with velocity dispersion. The increasing relation between group radii and merger fraction may be a result of larger groups having smaller velocity dispersion at a larger distance from the centre or larger groups hosting smaller, infalling groups with more mergers. However, we do not find evidence of smaller groups having higher merger fractions., Comment: 15 pages, 10 figures, 8 tables, 2 appendices, accepted for publication in Astronomy & Astrophysics

Galaxy mergers play a crucial role in galaxy evolution. However, the correlation between mergers and the local environment of galaxies is not fully understood. We aim to address the question of whether galaxy mergers prefer denser or less dense environments by quantifying the spatial clustering of mergers and non-mergers. We use two different indicators to classify mergers and non-mergers - classification based on a deep learning technique ($f$) and non-parametric measures of galaxy morphology, Gini-$M_{20}$ ($g$). We used a set of galaxy samples in the redshift range $0.1 < z < 0.15$ from the Galaxy and Mass Assembly (GAMA) survey with a stellar mass cut of $\log (M_{\star}/M_{\odot} ) > 9.5$. We measured and compared the two-point correlation function (2pCF) of mergers and non-mergers classified using the two merger indicators $f$ and $g$. We measured the marked correlation function (MCF), in which the galaxies are weighted by $f$ to probe the environmental dependence of galaxy mergers. We do not observe a statistically significant difference between the clustering strengths of mergers and non-mergers obtained using 2pCF. However, using the MCF measurements with $f$ as a mark, we observe an anti-correlation between the likelihood of a galaxy being a merger and its environment. Our results emphasise the advantage of MCF over 2pCF in probing the environmental correlations. Based on the MCF measurements, we conclude that the galaxy mergers prefer to occur in the under-dense environments on scales $> 50 \, h^{-1} \mathrm{kpc}$ of the large-scale structure (LSS). We attribute this observation to the high relative velocities of galaxies in the densest environments that prevent them from merging., Comment: 13 pages, 9 figures, matches the version published in A&A

Large future sky surveys, such as the LSST, will provide optical photometry for billions of objects. This paper aims to construct a proxy for the far ultraviolet attenuation (AFUVp) from the optical data alone, enabling the rapid estimation of the star formation rate (SFR) for galaxies that lack UV or IR data. To mimic LSST observations, we use the deep panchromatic optical coverage of the SDSS Photometric Catalogue DR~12, complemented by the estimated physical properties for the SDSS galaxies from the GALEX-SDSS-WISE Legacy Catalog (GSWLC) and inclination information obtained from the SDSS DR7. We restricted our sample to the 0.025-0.1 z-spec range and investigated relations among surface brightness, colours, and dust attenuation in the far UV range for star-forming galaxies obtained from the spectral energy distribution (SED). {Dust attenuation is best correlated with (u-r) colour and the surface brightness in the u band ($\rm \mu_{u}$). We provide a dust attenuation proxy for galaxies on the star-forming main sequence, which can be used for the LSST or any other type of broadband optical survey. The mean ratio between the catalogue values of SFR and those estimated using optical-only SDSS data with the AFUVp prior calculated as $\Delta$SFR=log(SFR$_{\tiny{\mbox{this work}}}$/SFR$_{\tiny{}\texttt{GSWLC}}$) is found to be less than 0.1~dex, while runs without priors result in an SFR overestimation larger than 0.3~dex. The presence or absence of theAFUVp has a negligible influence on the stellar mass estimation (with $\Delta$M$_{star}$ in the range from 0 to $-0.15$ dex). Forthcoming deep optical observations of the LSST Deep Drilling Fields, which also have multi-wavelength data, will enable one to calibrate the obtained relation for higher redshift galaxies and, possibly, extend the study towards other types of galaxies, such as early-type galaxies off the main sequence., Comment: 18 pages, accepted for publication in Astronomy and Astrophysics

We used 5487 star-forming galaxies at redshift z~0.63 extracted from the VIPERS and 143774 comparison galaxies in the local Universe from the GSWLC catalog. We employed two families of methods: parametric and non-parametric. In the former approaches, we compared the FMR projections plagued by observational biases on differently constructed control samples at various redshifts. Then, the metallicity difference between different redshifts in M*-SFR bins. In the latter approach, we related the metallicity and the normalized sSFR. The methodologies implemented to construct fair, complete samples for studying the MZR and the FMR produced consistent results showing a small, but still statistically significant evolution of both relations up to z~0.63. In particular, we observed a systematic trend where the median metallicity of the sample at z=0.63 is lower than that of the local sample at the same M* and SFR. The average difference in the metallicity of the low and intermediate redshifts is approximately 1.8 times the metallicity standard deviation of the median, of the intermediate redshift sample, in M*-SFR bins. We confirmed this result using the Kolmogorov-Smirnov test. When we applied the M*-completeness criterion to catalogs, the metallicity difference in redshifts decreased to approximately 0.96 times the metallicity standard deviation of the median, thus not statistically significant. This result may be dominated by the limited parameter space, being the lower M* galaxies where the difference is larger out from the analysis. A careful reading of the results, and their underlying selection criteria, are crucial in studies of the mass-metallicity and FMRs., Comment: Accepted for publication in A&A, 22 pages, 20 figures

We present a near-IR survey of the Trapezium Cluster and inner Orion Nebula using the NASA/ESA/CSA James Webb Space Telescope. The survey with the NIRCam instrument covers 10.9 x 7.5 arcminutes (~1.25 x 0.85 pc) in twelve wide-, medium-, and narrow-band filters from 1-5 microns and is diffraction-limited at all wavelengths, providing a maximum spatial resolution of 0.063 arcsec at 2 microns, corresponding to ~25 au at Orion. The suite of filters chosen was designed to address a number of scientific questions including the form of the extreme low-mass end of the IMF into the planetary-mass range to 1 Jupiter mass and below; the nature of ionised and non-ionised circumstellar disks and associated proplyds in the near-IR with a similar resolution to prior HST studies; to examine the large fragmented outflow from the embedded BN-KL region at very high resolution and fidelity; and to search for new jets and outflows from young stars in the Trapezium Cluster and the Orion Molecular Cloud 1 behind. In this paper, we present a description of the design of the observational programme, explaining the rationale for the filter set chosen and the telescope and detector modes used to make the survey; the reduction of the data using the JWST pipeline and other tools; the creation of large colour mosaics covering the region; and an overview of the discoveries made in the colour images and in the individual filter mosaics. Highlights include the discovery of large numbers of free-floating planetary-mass candidates as low as 0.6 Jupiter masses, a significant fraction of which are in wide binaries; new emission phenomena associated with the explosive outflow from the BN-KL region; and a mysterious "dark absorber" associated with a number of disparate features in the region, but which is seen exclusively in the F115W filter. Further papers will examine those discoveries and others in more detail., Comment: 39 pages, 18 figures, 1 table

The physical mechanisms responsible for bar formation and destruction in galaxies remain a subject of debate. While we have gained valuable insight into how bars form and evolve from isolated idealized simulations, in the cosmological domain, galactic bars evolve in complex environments with mergers, gas accretion events, in presence of turbulent Inter Stellar Medium (ISM) with multiple star formation episodes, in addition to coupling to their host galaxies' dark matter halos. We investigate bar formation in 13 Milky Way-mass galaxies from the FIRE-2 (Feedback in Realistic Environments) cosmological zoom-in simulations. 8 of the 13 simulated galaxies form bars at some point during their history: three from tidal interactions and five from internal evolution of the disk. The bars in FIRE-2 are generally shorter than the corotation radius (mean bar radius $\sim 1.53$ kpc), have a wide range of pattern speeds (36--97 km s$^{-1}$kpc$^{-1}$), and live for a wide range of dynamical times (2--160 bar rotations). We find that bar formation in FIRE-2 galaxies is influenced by satellite interactions and the stellar-to-dark matter mass ratio in the inner galaxy, but neither is a sufficient condition for bar formation. Bar formation is more likely to occur, and the bars formed are stronger and longer-lived, if the disks are kinematically cold; galaxies with high central gas fractions and/or vigorous star formation, on the other hand, tend to form weaker bars. In the case of the FIRE-2 galaxies these properties combine to produce ellipsoidal bars with strengths $A_2/A_0 \sim$ 0.1--0.2., Comment: 36 pages (including Appendix), 22 figures, to be submitted to the Astrophysical Journal

Low mass galaxy pair fractions are understudied, and it is unclear whether low mass pair fractions evolve in the same way as more massive systems over cosmic time. In the era of JWST, Roman, and Rubin, selecting galaxy pairs in a self-consistent way will be critical to connect observed pair fractions to cosmological merger rates across all mass scales and redshifts. Utilizing the Illustris TNG100 simulation, we create a sample of physically associated low mass ($\rm 10^81$. At $z=0$, the low mass major (1:4 mass ratio) pair fraction is 4$\times$ lower than high mass pairs, consistent with findings for cosmological merger rates. We show that separation limits that vary with the mass and redshift of the system, such as scaling by the virial radius of the host halo ($r_{\mathrm{sep}}< 1 R_{\rm vir}$), are critical for recovering pair fraction differences between low mass and high mass systems. Alternatively, static physical separation limits applied equivalently to all galaxy pairs do not recover the differences between low and high mass pair fractions, even up to separations of $300$ kpc. Finally, we place isolated mass-analogs of Local Group galaxy pairs, i.e., Milky Way (MW)--M31, MW--LMC, LMC--SMC, in a cosmological context, showing that isolated analogs of LMC--SMC-mass pairs and low-separation ($<50$ kpc) MW--LMC-mass pairs are $2-3\times$ more common at $z\gtrsim2-3$., Comment: 18 pages, 5 figures, accepted for publication in ApJ, updated Feb 2024 to reflect accepted version

Bright galaxies at sub-millimetre wavelengths from Herschel are now well known to be predominantly strongly gravitationally lensed. The same models that successfully predicted this strongly lensed population also predict about one percent of faint $450{\mu}$m-selected galaxies from deep James Clerk Maxwell Telescope (JCMT) surveys will also be strongly lensed. Follow-up ALMA campaigns have so far found one potential lens candidate, but without clear compelling evidence e.g. from lensing arcs. Here we report the discovery of a compelling gravitational lens system confirming the lensing population predictions, with a $z_{s} = 3.4 {\pm} 0.4$ submm source lensed by a $z_{spec} = 0.360$ foreground galaxy within the COSMOS field, identified through public JWST imaging of a $450{\mu}$m source in the SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES) catalogue. These systems will typically be well within the detectable range of future wide-field surveys such as Euclid and Roman, and since sub-millimetre galaxies are predominantly very red at optical/near-infrared wavelengths, they will tend to appear in near-infrared channels only. Extrapolating to the Euclid-Wide survey, we predict tens of thousands of strongly lensed near-infrared galaxies. This will be transformative for the study of dusty star-forming galaxies at cosmic noon, but will be a contaminant population in searches for strongly lensed ultra-high-redshift galaxies in Euclid and Roman., Comment: 9 pages, 5 figures, 4 tables; added figure, updated Discussion section, and minor revisions; accepted for publication in MNRAS

This work aims to determine how the galaxy main sequence (MS) changes using seven different commonly used methods to select the star-forming galaxies within VIPERS data over $0.5 \leq z < 1.2$. The form and redshift evolution of the MS will then be compared between selection methods. The star-forming galaxies were selected using widely known methods: a specific star-formation rate (sSFR), Baldwin, Phillips and Terlevich (BPT) diagram, 4000\AA\ spectral break (D4000) cut and four colour-colour cuts: NUVrJ, NUVrK, u-r, and UVJ. The main sequences were then fitted for each of the seven selection methods using a Markov chain Monte Carlo forward modelling routine, fitting both a linear main sequence and a MS with a high-mass turn-over to the star-forming galaxies. This was done in four redshift bins of $0.50 \leq z < 0.62$, $0.62 \leq z < 0.72$, $0.72 \leq z < 0.85$, and $0.85 \leq z < 1.20$. The slopes of all star-forming samples were found to either remain constant or increase with redshift, and the scatters were approximately constant. There is no clear redshift dependency of the presence of a high-mass turn-over for the majority of samples, with the NUVrJ and NUVrK being the only samples with turn-overs only at low redshift. No samples have turn-overs at all redshifts. Star-forming galaxies selected with sSFR and u-r are the only samples to have no high-mass turn-over in all redshift bins. The normalisation of the MS increases with redshift, as expected. The scatter around the MS is lower than the $\approx$0.3~dex typically seen in MS studies for all seven samples. The lack, or presence, of a high-mass turn-over is at least partially a result of the method used to select star-forming galaxies. However, whether a turn-over should be present or not is unclear., Comment: 20 pages, 3 appendices, 14 figures, 5 tables, accepted for publication in Astronomy & Astrophysics

As our nearest large neighbor, the Andromeda Galaxy provides a unique laboratory for investigating galaxy formation and the distribution and substructure properties of dark matter in a Milky Way-like galaxy. Here, we propose an initial 2-epoch ($\Delta t\approx 5$yr), 2-band Roman survey of the entire halo of Andromeda, covering 500 square degrees, which will detect nearly every red giant star in the halo (10$\sigma$ detection in F146, F062 of 26.5, 26.1AB mag respectively) and yield proper motions to $\sim$25 microarcsec/year (i.e., $\sim$90 km/s) for all stars brighter than F146 $\approx 23.6$ AB mag (i.e., reaching the red clump stars in the Andromeda halo). This survey will yield (through averaging) high-fidelity proper motions for all satellites and compact substructures in the Andromeda halo and will enable statistical searches for clusters in chemo-dynamical space. Adding a third epoch during the extended mission will improve these proper motions by $\sim t^{-1.5}$, to $\approx 11$ km/s, but this requires obtaining the first epoch in Year 1 of Roman operations. In combination with ongoing and imminent spectroscopic campaigns with ground-based telescopes, this Roman survey has the potential to yield full 3-d space motions of $>$100,000 stars in the Andromeda halo, including (by combining individual measurements) robust space motions of its entire globular cluster and most of its dwarf galaxy satellite populations. It will also identify high-velocity stars in Andromeda, providing unique information on the processes that create this population. These data offer a unique opportunity to study the immigration history, halo formation, and underlying dark matter scaffolding of a galaxy other than our own., Comment: Submitted in response to the call for Roman Space Telescope Core Community Survey white papers

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

We analyze pre-explosion near- and mid-infrared (IR) imaging of the site of SN 2023ixf in the nearby spiral galaxy M101 and characterize the candidate progenitor star. The star displays compelling evidence of variability with a possible period of $\approx$1000 days and an amplitude of $\Delta m \approx 0.6$ mag in extensive monitoring with the Spitzer Space Telescope since 2004, likely indicative of radial pulsations. Variability consistent with this period is also seen in the near-IR $J$ and $K_{s}$ bands between 2010 and 2023, up to just 10 days before the explosion. Beyond the periodic variability, we do not find evidence for any IR-bright pre-supernova outbursts in this time period. The IR brightness ($M_{K_s} = -10.7$ mag) and color ($J-K_{s} = 1.6$ mag) of the star suggest a luminous and dusty red supergiant. Modeling of the phase-averaged spectral energy distribution (SED) yields constraints on the stellar temperature ($T_{\mathrm{eff}} = 3500_{-1400}^{+800}$ K) and luminosity ($\log L/L_{\odot} = 5.1\pm0.2$). This places the candidate among the most luminous Type II supernova progenitors with direct imaging constraints, with the caveat that many of these rely only on optical measurements. Comparison with stellar evolution models gives an initial mass of $M_{\mathrm{init}} = 17\pm4 M_{\odot}$. We estimate the pre-supernova mass-loss rate of the star between 3 and 19 yr before explosion from the SED modeling at $\dot M \approx 3\times10^{-5}$ to $3\times10^{-4} M_{\odot}$ yr$^{-1}$ for an assumed wind velocity of $v_w = 10$ km s$^{-1}$, perhaps pointing to enhanced mass loss in a pulsation-driven wind., Comment: 13 pages, 5 figures, published in ApJL, replacement with revisions to match published version

Although low surface brightness galaxies (LSBs) contribute a large fraction to the number density of galaxies, their properties are still poorly known. LSBs are often considered dust poor, based only on a few studies. We use, for the first time, a large sample of LSBs and high surface brightness galaxies (HSBs) with deep observational data to study their dust properties as a function of surface brightness. Our sample consists of 1631 optically selected galaxies at $z < 0.1$ from the North Ecliptic Pole (NEP) wide field. We use the large set of data available in this field, from UV to FIR. We measured the optical size and the surface brightness of the targets, and analyzed their spectral energy distribution using the CIGALE fitting code. We found that the specific star formation rate and specific infrared luminosity (total infrared luminosity per stellar mass) remain mostly flat as a function of surface brightness for both LSBs and HSBs that are star-forming but decline steeply for the quiescent galaxies. The majority of LSBs in our sample have negligible dust attenuation (A$_{V} < 0.1$ mag), except for about 4% of them that show significant attenuation with a mean A$_{V}$ of 0.8 mag. We found that these LSBs also have a high $\textit{r}$-band mass-to-light ratio ($M/L_r>3$ M$_{\odot}$/L$_{\odot}$), and show similarity to the extreme giant LSBs from the literature, indicating a possibly higher dust attenuation in giant LSBs as well. This work provides a large catalog of LSBs and HSBs with detailed measurements of their several optical and infrared physical properties. Our results suggest that the dust content of LSBs is more varied than previously thought, with some of them having significant attenuation making them fainter than their intrinsic value. This will have serious implications for the observation and analysis of LSBs with current/upcoming surveys like JWST and LSST., Comment: 16 pages, 12 figures, accepted for publication in A&A

Stellar streams form through the tidal disruption of satellite galaxies or globular clusters orbiting a host galaxy. Globular cluster streams are exciting since they are thin (dynamically cold) and, therefore sensitive to perturbations from low-mass subhalos. Since the subhalo mass function differs depending on the dark matter composition, these gaps can provide unique constraints on dark matter models. However, current samples are limited to the Milky Way. With its large field of view, deep imaging sensitivity, and high angular resolution, the upcoming Nancy Grace Roman Space Telescope (Roman) presents a unique opportunity to increase the number of observed streams and gaps significantly. This paper presents a first exploration of the prospects for detecting gaps in streams in M31 and other nearby galaxies with resolved stars. We simulate the formation of gaps in a Palomar-5-like stream and generate mock observations of these gaps with background stars in M31 and the foreground Milky Way stellar fields. We assess Roman's ability to detect gaps out to 10 Mpc through visual inspection and with the gap-finding tool ${\texttt{FindTheGap}}$. We conclude that gaps of $\approx 1.5$ kpc in streams that are created from subhalos of masses $\geq5 \times 10^6$M$_{\odot}$ are detectable within a 2-3 Mpc volume in exposures of 1000s to 1 hour. This volume contains $\approx 150$ galaxies, including $\approx 8$ galaxies with luminosities $>10^{9}~$L$_{\odot}$. Large samples of stream gaps in external galaxies will open up a new era of statistical analyses of gap characteristics in stellar streams and help constrain dark matter models., Comment: ApJ version

The Andromeda galaxy (M31) is our closest neighbouring spiral galaxy, making it an ideal target for studying the physics of the interstellar medium in a galaxy very similar to our own. Using new observations of M31 at 4.76GHz by the C-Band All-Sky Survey (C-BASS), and all available radio data at $1^\circ$ resolution, we produce the integrated spectrum and put new constraints on the synchrotron spectral index and anomalous microwave emission (AME) from M31. We use aperture photometry and spectral modelling to fit for the integrated spectrum of M31, and subtract a comprehensive model of nearby background radio sources. The AME in M31 is detected at $3\sigma$ significance with a peak near 30GHz and flux density $0.27\pm0.09$Jy. The synchrotron spectral index of M31 is flatter than our own Galaxy at $\alpha = -0.66 \pm 0.03$ with no strong evidence of spectral curvature. The emissivity of AME, averaged over the total emission from M31 is lower than typical AME sources in our Galaxy, implying that AME is not uniformly distributed throughout M31 and instead is likely confined to sub-regions -- this will need to be confirmed using future higher resolution observations around 20--30GHz., Comment: 16 pages, 6 figures. Accepted by MNRAS

Compact Symmetric Objects (CSOs) form a distinct class of jetted active galactic nuclei (jetted-AGN). We examine a carefully selected sample of 54 CSOs, and confirm that there are two unrelated classes: an edge-dimmed, low-luminosity class (CSO 1), and an edge-brightened, high-luminosity class (CSO 2). Using statistically significant blind tests, we show that CSO 2s themselves consist of two morphologically distinct classes: CSO 2.0, having prominent hot-spots at the leading edges of narrow jets and/or narrow lobes; and CSO 2.2, without prominent hot-spots, and with broad jets and/or lobes. An intermediate class, CSO 2.1, exhibits mixed properties. The four classes occupy different, overlapping, portions of the luminosity-size plane, with the sizes of largest CSOs being $\sim 500$pc. We advance the hypothesis that CSO 2.0s are young and evolve through CSO 2.1s into CSO 2.2s, which are old (up to $\sim 5000$ yr). Thus CSOs do not evolve into larger types of jetted-AGN, but spend their whole life cycle as CSOs. The radio emission region energies in the CSO 2s we have studied range from $\sim 10^{-4}\, M_\odot {c}^2$ to $\sim 7 \, M_\odot {c}^2$. We show that the transient nature of CSO 2s, and their birthrate, can be explained through ignition in the tidal disruption events of giant stars. We also consider possibilities for tapping the spin energy of the supermassive black hole, and tapping the energy of the accretion disk, in a manner similar to, but not the same as, that which occurs in dwarf novae. Our results demonstrate conclusively that CSOs constitute a large family of AGN in which we have thus far studied only the brightest. More comprehensive radio studies, with higher sensitivity, resolution, and dynamic range, will revolutionize our understanding of AGN and the central engines that power them., Comment: 41 pages, 13 figures, 10 tables

Compact Symmetric Objects (CSOs) are jetted Active Galactic Nuclei (AGN) with overall projected size < 1 kpc. The classification was introduced to distinguish these objects from the majority of known compact jetted-AGN, where the observed emission is relativistically boosted towards the observer. The original classification criteria for CSOs were: (i) evidence of emission on both sides of the center of activity, and (ii) overall size < 1 kpc. However some relativistically boosted objects with jet axes close to the line of sight appear symmetric and have been mis-classified as CSOs in the literature, thereby undermining the CSO classification. We introduce two new CSO classification criteria based on (iii) flux density variability, and (iv) the apparent velocity of components moving along the jets. As a first step towards creating a comprehensive catalog of "bona fide" CSOs, in this paper we identify 79 bona fide CSOs that meet our expanded CSO selection criteria. This sample of bona fide CSOs can be used for astrophysical studies of CSOs without fear of contamination by objects incorrectly identified as CSOs. We define three complete sub-samples of the 79 CSOs, which are suitable for statistical tests, and show that the fraction of CSOs in flux density limited samples with $\rm S_{5~GHz}$ > 700 mJy is between 6.4% and 8.2%., Comment: 21 pages, 7 figures, 3 tables

Galaxy mergers are hugely important in our current dark matter cosmology. These powerful events cause the disruption of the merging galaxies, pushing the gas, stars and dust of the galaxies resulting in changes to morphologies. This disruption can also cause more extreme events inside the galaxies: periods of extreme star formation rates and the rapid increase in active galactic nuclei activity. Hence, to better understand what goes on in these rare events, we need to be able to identify statistically large samples. In the last few years, the growth of artificial intelligence techniques has seen application to identifying galaxy mergers. These techniques have shown to be highly accurate and their application has grown beyond academic studies of ``can we?'' to deeper scientific use. However, these classifications are not without their problems. In this proceedings, we will explore how galaxy merger classification can be improved by adding pre-extracted galaxy morphologies alongside the traditional imaging data. This demonstrates that current neural networks are not extracting all the information from the images they are given. It will also explore how the resulting samples of rare objects could be highly contaminated. This has a knock on impact on the upcoming large scale surveys like Euclid and Rubin-LSST., Comment: 4 pages, 3 figures, proceedings for EAS 2022 S11, to be published in Memorie della SAIt

Observations have found black holes spanning ten orders of magnitude in mass across most of cosmic history. The Kerr black hole solution is however provisional as its behavior at infinity is incompatible with an expanding universe. Black hole models with realistic behavior at infinity predict that the gravitating mass of a black hole can increase with the expansion of the universe independently of accretion or mergers, in a manner that depends on the black hole's interior solution. We test this prediction by considering the growth of supermassive black holes in elliptical galaxies over $0

We report on analysis using the James Webb Space Telescope (JWST) to identify a candidate progenitor star of the Type II-plateau supernova SN 2022acko in the nearby, barred spiral galaxy NGC 1300. To our knowledge, our discovery represents the first time JWST has been used to localize a progenitor system in pre-explosion archival Hubble Space Telescope (HST) images. We astrometrically registered a JWST NIRCam image from 2023 January, in which the SN was serendipitously captured, to pre-SN HST F160W and F814W images from 2017 and 2004, respectively. An object corresponding precisely to the SN position has been isolated with reasonable confidence. That object has a spectral energy distribution (SED) and overall luminosity consistent with a single-star model having an initial mass possibly somewhat less than the canonical 8 Msun theoretical threshold for core collapse (although masses as high as 9 Msun for the star are also possible); however, the star's SED and luminosity are inconsistent with that of a super-asymptotic giant branch star which might be a forerunner of an electron-capture SN. The properties of the progenitor alone imply that SN 2022acko is a relatively normal SN II-P, albeit most likely a low-luminosity one. The progenitor candidate should be confirmed with follow-up HST imaging at late times, when the SN has sufficiently faded. This potential use of JWST opens a new era of identifying SN progenitor candidates at high spatial resolution., Comment: 8 pages, substantial changes from v1, to appear in MNRAS

Recent observations have revealed a trove of unexpected morphological features in many of the Milky Way's stellar streams. Explanations for such features include time-dependent deformations of the Galactic gravitational potential, local disruptions induced by dark matter substructure, and special configurations of the streams' progenitors. In this paper, we study how these morphologies can also arise in certain static, non-spherical gravitational potentials that host a subset of resonantly-trapped orbit families. The transitions, or separatrices, between these orbit families mark abrupt discontinuities in the orbital structure of the potential. We develop a novel numerical approach for measuring the libration frequencies of resonant and near-resonant orbits, and apply it to study the evolution of stellar streams on these orbits. We reveal two distinct morphological features that arise in streams on near-resonant orbits: fans, that come about due to a large spread in the libration frequencies near a separatrix; and bifurcations, that arise when a separatrix splits the orbital distribution of the stellar stream between two (or more) distinct orbit families. We demonstrate that these effects can arise in some Milky Way streams for certain choices of the dark matter halo potential, and discuss how this might be used to probe and constrain the global shape of the Milky Way's gravitational potential., Comment: 24 pages, 8 figures. Submitted to ApJ; comments welcome!

The assembly of stellar and supermassive black hole (SMBH) mass in elliptical galaxies since $z\sim1$ can help to diagnose the origins of locally-observed correlations between SMBH mass and stellar mass. We therefore construct three samples of elliptical galaxies, one at $z\sim0$ and two at $0.7\lesssim z \lesssim2.5$, and quantify their relative positions in the $M_{BH}-M_*$ plane. Using a Bayesian analysis framework, we find evidence for translational offsets in both stellar mass and SMBH mass between the local sample and both higher redshift samples. The offsets in stellar mass are small, and consistent with measurement bias, but the offsets in SMBH mass are much larger, reaching a factor of seven between $z\sim1$ and $z\sim0$. The magnitude of the SMBH offset may also depend on redshift, reaching a factor of $\sim20$ at $z\sim 2$. The result is robust against variation in the high and low redshift samples and changes in the analysis approach. The magnitude and redshift evolution of the offset are challenging to explain in terms of selection and measurement biases. We conclude that either there is a physical mechanism that preferentially grows SMBHs in elliptical galaxies at $z\lesssim 2$, or that selection and measurement biases are both underestimated, and depend on redshift., Comment: ApJ, accepted

We investigate the variation in the mid-infrared spectral energy distributions of 373 low-redshift ($z<0.4$) star-forming galaxies, which reflects a variety of polycyclic aromatic hydrocarbon (PAH) emission features. The relative strength of PAH emission is parameterized as $q_\mathrm{PAH}$, which is defined as the mass fraction of PAH particles in the total dust mass. With the aid of continuous mid-infrared photometric data points covering 7-24$\mu$m and far-infrared flux densities, $q_\mathrm{PAH}$ values are derived through spectral energy distribution fitting. The correlation between $q_\mathrm{PAH}$ and other physical properties of galaxies, i.e., gas-phase metallicity ($12+\mathrm{log(O/H)}$), stellar mass, and specific star-formation rate (sSFR) are explored. As in previous studies, $q_\mathrm{PAH}$ values of galaxies with high metallicity are found to be higher than those with low metallicity. The strength of PAH emission is also positively correlated with the stellar mass and negatively correlated with the sSFR. The correlation between $q_\mathrm{PAH}$ and each parameter still exists even after the other two parameters are fixed. In addition to the PAH strength, the application of metallicity-dependent gas-to-dust mass ratio appears to work well to estimate gas mass that matches the observed relationship between molecular gas and physical parameters. The result obtained will be used to calibrate the observed PAH luminosity-total infrared luminosity relation, based on the variation of MIR-FIR SED, which is used in the estimation of hidden star formation., Comment: 19 pages, 9 figures, AJ, in press

Aims. We present the application of a fully connected neural network (NN) for galaxy merger identification using exclusively photometric information. Our purpose is not only to test the method's efficiency, but also to understand what merger properties the NN can learn and what their physical interpretation is. Methods. We created a class-balanced training dataset of 5\,860 galaxies split into mergers and non-mergers. The galaxy observations came from SDSS DR6 and were visually identified in Galaxy Zoo. The 2$\,$930 mergers were selected from known SDSS mergers and the respective non-mergers were the closest match in both redshift and $r$ magnitude. The NN architecture was built by testing a different number of layers with different sizes and variations of the dropout rate. We compared input spaces constructed using: the five SDSS filters: $u$, $g$, $r$, $i$, and $z$; combinations of bands, colours, and their errors; six magnitude types; and variations of input normalization. Results. We find that the fibre magnitude errors contribute the most to the training accuracy. Studying the parameters from which they are calculated, we show that the input space built from the sky error background in the five SDSS bands alone leads to 92.64 $\pm$ 0.15 \% training accuracy. We also find that the input normalization, that is to say, how the data are presented to the NN, has a significant effect on the training performance. Conclusions. We conclude that, from all the SDSS photometric information, the sky error background is the most sensitive to merging processes. This finding is supported by an analysis of its five-band feature space by means of data visualization. Moreover, studying the plane of the $g$ and $r$ sky error bands shows that a decision boundary line is enough to achieve an accuracy of 91.59\%., Comment: 16 pages, 20 figures, accepted in A&A on 05/11/2022, final version submitted on 09/12/2022 (also current one in arXiv)

We aim to generate a catalogue of merging galaxies within the 5.4 sq. deg. North Ecliptic Pole over the redshift range $0.0 < z < 0.3$. To do this, imaging data from the Hyper Suprime-Cam are used along with morphological parameters derived from these same data. The catalogue was generated using a hybrid approach. Two neural networks were trained to perform binary merger non-merger classifications: one for galaxies with $z < 0.15$ and another for $0.15 \leq z < 0.30$. Each network used the image and morphological parameters of a galaxy as input. The galaxies that were identified as merger candidates by the network were then visually checked by experts. The resulting mergers will be used to calculate the merger fraction as a function of redshift and compared with literature results. We found that 86.3% of galaxy mergers at $z < 0.15$ and 79.0% of mergers at $0.15 \leq z < 0.30$ are expected to be correctly identified by the networks. Of the 34 264 galaxies classified by the neural networks, 10 195 were found to be merger candidates. Of these, 2109 were visually identified to be merging galaxies. We find that the merger fraction increases with redshift, consistent with literature results from observations and simulations, and that there is a mild star-formation rate enhancement in the merger population of a factor of $1.102 \pm 0.084$., Comment: Accepted to A&A, 26 pages, 20 figures, 8 tables, 3 appendixes, full tables 1 and 4 will be available on CDS

We present new radial velocity measurements from the Magellan and the Anglo-Australian Telescopes for 174 previously known and 122 newly confirmed globular clusters (GCs) around NGC 5128, the nearest accessible massive early-type galaxy at D=3.8 Mpc. Remarkably, 28 of these newly confirmed GCs are at projected radii >50' ($\gtrsim 54$ kpc), extending to $\sim 130$ kpc, in the outer halo where few GCs had been confirmed in previous work. We identify several subsets of GCs that spatially trace halo substructures that are visible in red giant branch star maps of the galaxy. In some cases, these subsets of GCs are kinematically cold, and may be directly associated with and originate from these specific stellar substructures. From a combined kinematic sample of 645 GCs, we see evidence for coherent rotation at all radii, with a higher rotation amplitude for the metal-rich GC subpopulation. Using the tracer mass estimator, we measure a total enclosed mass of $2.5\pm0.3 \times 10^{12} M_{\odot}$ within $\sim 120$ kpc, an estimate that will be sharpened with forthcoming dynamical modeling. The combined power of stellar mapping and GC kinematics makes NGC 5128 an ongoing keystone for understanding galaxy assembly at mass scales inaccessible in the Local Group., Comment: 24 pages, 13 figures, submitted to AAS Journals

We present the first Bayesian method for tomographic decomposition of the plane-of-sky orientation of the magnetic field with the use of stellar polarimetry and distance. This standalone tomographic inversion method presents an important step forward in reconstructing the magnetized interstellar medium (ISM) in 3D within dusty regions. We develop a model in which the polarization signal from the magnetized and dusty ISM is described by thin layers at various distances. Our modeling makes it possible to infer the mean polarization (amplitude and orientation) induced by individual dusty clouds and to account for the turbulence-induced scatter in a generic way. We present a likelihood function that explicitly accounts for uncertainties in polarization and parallax. We develop a framework for reconstructing the magnetized ISM through the maximization of the log-likelihood using a nested sampling method. We test our Bayesian inversion method on mock data taking into account realistic uncertainties from Gaia and as expected for the optical polarization survey PASIPHAE according to the currently planned observing strategy. We demonstrate that our method is effective at recovering the cloud properties as soon as the polarization induced by a cloud to its background stars is higher than $\sim 0.1\%$ for the adopted survey exposure time and level of systematic uncertainty. Our method makes it possible to recover not only the mean polarization properties but also to characterize the intrinsic scatter, thus creating new ways to characterize ISM turbulence and the magnetic field strength. Finally, we apply our method to an existing data set of starlight polarization with known line-of-sight decomposition, demonstrating agreement with previous results and an improved quantification of uncertainties in cloud properties., Comment: Closely matches the published version: A&A 670, A164, 2023 The accompanying code (BISP-1) is available at https://github.com/vpelgrims/Bisp_1

Studies of resolved stellar populations in the Milky Way and nearby galaxies reveal an amazingly detailed and clear picture of galaxy evolution. Within the Local Group, the ability to probe the stellar populations of small and large galaxies opens up the possibility of exploring key questions such as the nature of dark matter, the detailed formation history of different galaxy components, and the role of accretion in galactic formation. Upcoming wide-field surveys promise to extend this ability to all galaxies within 10~Mpc, drastically increasing our capability to decipher galaxy evolution and enabling statistical studies of galaxies' stellar populations. To facilitate the optimum use of these upcoming capabilities we develop a simple formalism to predict the density of resolved stars for an observation of a stellar population at fixed surface brightness and population parameters. We provide an interface to calculate all quantities of interest to this formalism via a public release of the code: \texttt{walter}. This code enables calculation of (i) the expected number density of detected stars, (ii) the exposure time needed to reach certain population features, such as the horizontal branch, and (iii) an estimate of the crowding limit, among other features. These calculations will be very useful for planning surveys with NASA's upcoming Nancy Grace Roman Space Telescope (Roman, formerly WFIRST), which we use for example calculations throughout this work., Comment: 16 pages, 6 figures, submitted to AJ. Comments welcome

We present Hubble Space Telescope (HST) imaging of the site of SN 2015bh in the nearby spiral galaxy NGC 2770 taken between 2017 and 2019, nearly four years after the peak of the explosion. In 2017-2018, the transient fades steadily in optical filters before declining more slowly to $F814W = -7.1$ mag in 2019, $\approx$4 mag below the level of its eruptive luminous blue variable (LBV) progenitor observed with HST in 2008-2009. The source fades at a constant color of $F555W - F814W = 0.4$ mag until 2018, similar to SN 2009ip and consistent with a spectrum dominated by interaction of the ejecta with circumstellar material (CSM). A deep optical spectrum obtained in 2021 lacks signatures of ongoing interaction ($L_{\mathrm{H}\alpha} \lesssim 10^{38}$ erg s$^{-1}$ for broadened emission $\lesssim$2000 km s$^{-1}$), but indicates the presence of a nearby H II region ($\lesssim$300 pc). The color evolution of the fading source makes it unlikely that emission from a scattered-light echo or binary OB companion of the progenitor contributes significantly to the flattening of the late-time light curve. The remaining emission in 2019 may plausibly be attributed an evolved/inflated companion or an unresolved ($\lesssim$3 pc), young stellar cluster. Importantly, the color evolution of SN 2015bh rules out scenarios in which the surviving progenitor is obscured by nascent dust and does not clearly indicate a transition to a hotter, optically faint state. The simplest explanation is that the massive progenitor did not survive. SN 2015bh likely represents a remarkable example of the terminal explosion of a massive star preceded by decades of end-stage eruptive variability., Comment: 15 pages, 5 figures, published in ApJL. Replacement with revisions to match published article

We present a novel method for constraining the length of the Galactic bar using 6D phase space information to directly integrate orbits. We define a pseudo-length for the Galactic bar, named $R_{Freq}$, based on the maximal extent of trapped bar orbits. We find the $R_{Freq}$ measured from orbits is consistent with the $R_{Freq}$ of the assumed potential only when the length of the bar and pattern speed of said potential is similar to the model from which the initial phase-space coordinates of the orbits are derived. Therefore, one can measure the model's or the Milky Way's bar length from 6D phase-space coordinates by determining which assumed potential leads to a self-consistent measured $R_{Freq}$. When we apply this method to $\approx$210,000 stars in APOGEE DR17 and $Gaia$ eDR3 data, we find a consistent result only for potential models with a dynamical bar length of $\approx$3.5 kpc. We find the Milky Way's trapped bar orbits extend out to only $\approx$3.5 kpc, but there is also an overdensity of stars at the end of the bar out to 4.8 kpc which could be related to an attached spiral arm. We also find that the measured orbital structure of the bar is strongly dependent on the properties of the assumed potential., Comment: 15 pages, 8 figures, 2 tables, accepted to MNRAS, comments welcome

In the coming decade, thousands of stellar streams will be observed in the halos of external galaxies. What fundamental discoveries will we make about dark matter from these streams? As a first attempt to look at these questions, we model Magellan/Megacam imaging of the Centaurus A's (Cen A) disrupting dwarf companion Dwarf 3 (Dw3) and its associated stellar stream, to find out what can be learned about the Cen A dark-matter halo. We develop a novel external galaxy stream-fitting technique and generate model stellar streams that reproduce the stream morphology visible in the imaging. We find that there are many viable stream models that fit the data well, with reasonable parameters, provided that Cen A has a halo mass larger than M$_{200}$ $>4.70\times 10^{12}$ M$_{\odot}$. There is a second stream in Cen A's halo that is also reproduced within the context of this same dynamical model. However, stream morphology in the imaging alone does not uniquely determine the mass or mass distribution for the Cen A halo. In particular, the stream models with high likelihood show covariances between the inferred Cen A mass distribution, the inferred Dw3 progenitor mass, the Dw3 velocity, and the Dw3 line-of-sight position. We show that these degeneracies can be broken with radial-velocity measurements along the stream, and that a single radial velocity measurement puts a substantial lower limit on the halo mass. These results suggest that targeted radial-velocity measurements will be critical if we want to learn about dark matter from extragalactic stellar streams., Comment: 11 figures, 23 pages, Accepted to ApJ

We study stellar and black hole mass assembly in a sample of 42 infrared-luminous galaxy mergers at z<0.3 by combining results from radiative transfer modelling with archival measures of molecular gas and black hole mass. The ratios of stellar mass, molecular gas mass, and black hole mass to each other are consistent with those of massive gas-rich galaxies at z<0.3. The advanced mergers may show increased black hole mass to stellar mass ratios, consistent with the transition from AGN to ellipticals and implying substantial black hole mass growth over the course of the merger. Star formation rates are enhanced relative to the local main sequence, by factors of ~100 in the starburst and ~1.8 in the host. The starburst star formation rates appear distinct to star formation in the main sequence at all redshifts up to at least z~5. Starbursts may prefer late-stage mergers, but are observed at any merger stage. We do not find evidence that the starbursts in these low-redshift systems substantially increase the total stellar mass, with a soft upper limit on the stellar mass increase from starburst activity of about a factor of two. In contrast, 12 objects show evidence for super-Eddington accretion, associated with late-stage mergers, suggesting that many AGN in infrared-luminous mergers go through a super-Eddington phase. The super-Eddington phase may increase black hole mass by up to an order of magnitude at an accretion efficiency of 42+/-33% over a period of 44+/-22Myr. Our results imply that super-Eddington accretion is an important black hole growth channel in infrared-luminous galaxies at all redshifts., Comment: MNRAS Accepted

We present the ultraviolet to submillimetre spectral energy distributions (SEDs) of the HERschel Ultra Luminous Infrared Galaxy Survey (HERUS) sample of 42 local ultraluminous infrared galaxies (ULIRGs) and fit them with a Markov chain Monte Carlo (MCMC) code using the CYprus models for Galaxies and their NUclear Spectra (CYGNUS) radiative transfer models for starbursts, active galactic nucleus (AGN) tori and host galaxy. The Spitzer IRS spectroscopy data are included in the fitting. Our bayesian SED fitting method takes comparable time to popular energy balance methods but it is more physically motivated and versatile. All HERUS galaxies harbor high rates of star formation but we also find bolometrically significant AGN in all of the galaxies of the sample. We estimate the correction of the luminosities of the AGN in the ULIRGs due to the anisotropic emission of the torus and find that it could be up to a factor of $\sim10$ for nearly edge-on tori. We present a comparison of our results with the smooth torus model of Fritz et al. and the two-phase models of Siebenmorgen et al. and SKIRTOR. We find that the CYGNUS AGN torus models fit significantly better the SEDs of our sample compared to all other models. We find no evidence that strong AGN appear either at the beginning or end of a starburst episode or that starbursts and AGN affect each other. IRAS 01003-2238 and Mrk 1014 show evidence for dual AGN in their SED fits suggesting a minimum dual AGN fraction in the sample of 5%., Comment: Published in Monthly Notices of the Royal Astronomical Society. The CYGNUS models used in this paper are available at https://arc.euc.ac.cy/cygnus/

We present the multi-wavelength counterparts of 850-$\mu$m selected submillimetre sources over a 2-deg$^2$ field centred on the North Ecliptic Pole. In order to overcome the large beam size (15 arcsec) of the 850-$\mu$m images, deep optical to near-infrared (NIR) photometric data and arcsecond-resolution 20-cm images are used to identify counterparts of submillimetre sources. Among 647 sources, we identify 514 reliable counterparts for 449 sources (69 per cent in number), based either on probabilities of chance associations calculated from positional offsets or offsets combined with the optical-to-NIR colours. In the radio imaging, the fraction of 850-$\mu$m sources having multiple counterparts is 7 per cent. The photometric redshift, infrared luminosity, stellar mass, star-formation rate (SFR), and the AGN contribution to the total infrared luminosity of the identified counterparts are investigated through spectral energy distribution fitting. The SMGs are infrared-luminous galaxies at an average $\langle z\rangle=2.5$ with $\mathrm{log}_{10} (L_\mathrm{IR}/\mathrm{L}_\odot)=11.5-13.5$, with a mean stellar mass of $\mathrm{log}_{10} (M_\mathrm{star}/\mathrm{M}_\odot)=10.90$ and SFR of $\mathrm{log}_{10} (\mathrm{SFR/M_\odot\,yr^{-1}})=2.34$. The SMGs show twice as large SFR as galaxies on the star-forming main sequence, and about 40 per cent of the SMGs are classified as objects with bursty star formation. At $z\ge4$, the contribution of AGN luminosity to total luminosity for most SMGs is larger than 30 per cent. The FIR-to-radio correlation coefficient of SMGs is consistent with that of main-sequence galaxies at $z\simeq2$., Comment: 24 pages, 16 figures. Accepted for publication in MNRAS

Galaxy clusters provide an excellent probe in various research fields in astrophysics and cosmology. However, the number of galaxy clusters detected so far in the $AKARI$ North Ecliptic Pole (NEP) field is limited. In this work, we provide galaxy cluster candidates in the $AKARI$ NEP field with the minimum requisites based only on coordinates and photometric redshift (photo-$z$) of galaxies. We used galaxies detected in 5 optical bands ($g$, $r$, $i$, $z$, and $Y$) by the Subaru Hyper Suprime-Cam (HSC), assisted with $u$-band from Canada-France-Hawaii Telescope (CFHT) MegaPrime/MegaCam, and IRAC1 and IRAC2 bands from the $Spitzer$ space telescope for photo-$z$ estimation. We calculated the local density around every galaxy using the 10$^{th}$-nearest neighbourhood. Cluster candidates were determined by applying the friends-of-friends algorithm to over-densities. 88 cluster candidates containing 4390 member galaxies below redshift 1.1 in 5.4 deg$^2$ have been detected. The reliability of our method was examined through false detection tests, redshift uncertainty tests, and applications on the COSMOS data, giving false detection rates of 0.01 to 0.05 and recovery rate of 0.9 at high richness. 3 X-ray clusters previously observed by $ROSAT$ and $Chandra$ were recovered. The cluster galaxies show higher stellar mass and lower star formation rate (SFR) compared to the field galaxies in two-sample Z-tests. These cluster candidates are useful for environmental studies of galaxy evolution and future astronomical surveys in the NEP, where $AKARI$ has performed unique 9-band mid-infrared photometry for tens of thousands of galaxies., Comment: 19 pages, 16 figures, has been accepted for publication in MNRAS

The C-Band All-Sky Survey (C-BASS) has observed the Galaxy at 4.76GHz with an angular resolution of $0.73^\circ$ full-width half-maximum, and detected Galactic synchrotron emission with high signal-to-noise ratio over the entire northern sky ($\delta > -15^{\circ}$). We present the results of a spatial correlation analysis of Galactic foregrounds at mid-to-high ($b > 10^\circ$) Galactic latitudes using a preliminary version of the C-BASS intensity map. We jointly fit for synchrotron, dust, and free-free components between $20$ and $1000$GHz and look for differences in the Galactic synchrotron spectrum, and the emissivity of anomalous microwave emission (AME) when using either the C-BASS map or the 408MHz all-sky map to trace synchrotron emission. We find marginal evidence for a steepening ($\left<\Delta\beta\right> = -0.06\pm0.02$) of the Galactic synchrotron spectrum at high frequencies resulting in a mean spectral index of $\left<\beta\right> = -3.10\pm0.02$ over $4.76-22.8$GHz. Further, we find that the synchrotron emission can be well modelled by a single power-law up to a few tens of GHz. Due to this, we find that the AME emissivity is not sensitive to changing the synchrotron tracer from the 408MHz map to the 4.76GHz map. We interpret this as strong evidence for the origin of AME being spinning dust emission., Comment: 20 pages, 13 figures, published with MNRAS

Symmetric Achromatic Variability (SAV) is a rare form of radio variability in blazars that has been attributed to gravitational millilensing by a ~$10^2 - 10^5$ $M_\odot$ mass condensate. Four SAVs have been identified between 1980 and 2020 in the long-term radio monitoring data of the blazar PKS 1413+135. We show that all four can be fitted with the same, unchanging, gravitational lens model. If SAV is due to gravitational millilensing, PKS 1413+135 provides a unique system for studying active galactic nuclei with unprecedented microarcsecond resolution, as well as for studying the nature of the millilens itself. We discuss two possible candidates for the putative millilens: a giant molecular cloud hosted in the intervening edge-on spiral galaxy, and an undetected dwarf galaxy with a massive black hole. We find a significant dependence of SAV crossing time on frequency, which could indicate a fast shock moving in a slower underlying flow. We also find tentative evidence for a 989-day periodicity in the SAVs, which, if real, makes possible the prediction of future SAVs: the next three windows for possible SAVs begin in August 2022, May 2025, and February 2028., Comment: Accepted to ApJ

Context. The North Ecliptic Pole (NEP) field provides a unique set of panchromatic data, well suited for active galactic nuclei (AGN) studies. Selection of AGN candidates is often based on mid-infrared (MIR) measurements. Such method, despite its effectiveness, strongly reduces a catalog volume due to the MIR detection condition. Modern machine learning techniques can solve this problem by finding similar selection criteria using only optical and near-infrared (NIR) data. Aims. Aims of this work were to create a reliable AGN candidates catalog from the NEP field using a combination of optical SUBARU/HSC and NIR AKARI/IRC data and, consequently, to develop an efficient alternative for the MIR-based AKARI/IRC selection technique. Methods. A set of supervised machine learning algorithms was tested in order to perform an efficient AGN selection. Best of the models were formed into a majority voting scheme, which used the most popular classification result to produce the final AGN catalog. Additional analysis of catalog properties was performed in form of the spectral energy distribution (SED) fitting via the CIGALE software. Results. The obtained catalog of 465 AGN candidates (out of 33 119 objects) is characterized by 73% purity and 64% completeness. This new classification shows consistency with the MIR-based selection. Moreover, 76% of the obtained catalog can be found only with the new method due to the lack of MIR detection for most of the new AGN candidates. Training data, codes and final catalog are available via the github repository. Final AGN candidates catalog will be also available via the CDS service after publication., Comment: Accepted for publication in Astronomy and Astrophysics

Mergers and tidal interactions between massive galaxies and their dwarf satellites are a fundamental prediction of the Lambda-Cold Dark Matter cosmology. These events are thought to provide important observational diagnostics of nonlinear structure formation. Stellar streams in the Milky Way and Andromeda are spectacular evidence for ongoing satellite disruption. However, constructing a statistically meaningful sample of tidal streams beyond the Local Group has proven a daunting observational challenge, and the full potential for deepening our understanding of galaxy assembly using stellar streams has yet to be realised. Here we introduce the Stellar Stream Legacy Survey, a systematic imaging survey of tidal features associated with dwarf galaxy accretion around a sample of ~3100 nearby galaxies within z~0.02, including about 940 Milky Way analogues. Our survey exploits public deep imaging data from the DESI Legacy Imaging Surveys, which reach surface brightness as faint as ~29 mag/arcsec^2 in the r band. As a proof of concept of our survey, we report the detection and broad-band photometry of 24 new stellar streams in the local Universe. We discuss how these observations can yield new constraints on galaxy formation theory through comparison to mock observations from cosmological galaxy simulations. These tests will probe the present-day mass assembly rate of galaxies, the stellar populations and orbits of satellites, the growth of stellar halos and the resilience of stellar disks to satellite bombardment., Comment: Accepted in Astronomy & Astrophysics, 25 pages, 17 figures, 3 tables

In this work we study the striking case of a narrow blue stream around the NGC 7241 galaxy and its foreground dwarf companion. We want to figure out if the stream was generated by tidal interaction with NGC 7241 or it first interacted with the foreground dwarf companion and later both fell together towards NGC 7241. We use four sets of observations, including a follow-up spectroscopic study with the MEGARA instrument at the 10.4-m Gran Telescopio Canarias. Our data suggest that the compact object we detected in the stream is a foreground Milky Way halo star. Near this compact object we detect emission lines overlapping a bluer and fainter blob of the stream that is clearly visible in both ultra-violet and optical deep images. From its heliocentric systemic radial velocity (Vsyst= 1548.58+/-1.80 km s^-1) and new UV and optical broad-band photometry, we conclude that this over-density could be the actual core of the stream, with an absolute magnitude of M_g ~ -10 and a (g-r) = 0.08 +/- 0.11, consistent with a remnant of a low-mass dwarf satellite undergoing a current episode of star formation. From the width of the stream and assuming a circular orbit, we calculate that the progenitor mass can be the typical of a dwarf galaxy, but it could also be substantially lower if the stream is on a very radial orbit or it was created by tidal interaction with the companion dwarf instead of with NGC 7241. Finally, we find that blue stellar streams containing star formation regions are commonly predicted by high-resolution cosmological simulations of galaxies lighter than the Milky Way. This scenario is consistent with the processes explaining the bursty star formation history of some dwarf satellites, which are followed by a gas depletion and a fast quenching once they enter within the virial radius of their host galaxies for the first time., Comment: 12 pages, 5 figures, accepted for publication in Astronomy & Astrophysics