Showing 1–50 of 238 results for author: Knapen, J H

Deriving the star formation histories of galaxies from spectra with simulation-based inference

Authors: Patricia Iglesias-Navarro, Marc Huertas-Company, Ignacio Martín-Navarro, Johan H. Knapen, Emilie Pernet

Abstract: High-resolution galaxy spectra encode information about the stellar populations within galaxies. The properties of the stars, such as their ages, masses, and metallicities, provide insights into the underlying physical processes that drive the growth and transformation of galaxies over cosmic time. We explore a simulation-based inference (SBI) workflow to infer from optical absorption spectra th… ▽ More High-resolution galaxy spectra encode information about the stellar populations within galaxies. The properties of the stars, such as their ages, masses, and metallicities, provide insights into the underlying physical processes that drive the growth and transformation of galaxies over cosmic time. We explore a simulation-based inference (SBI) workflow to infer from optical absorption spectra the posterior distributions of metallicities and the star formation histories (SFHs) of galaxies (i.e. the star formation rate as a function of time). We generated a dataset of synthetic spectra to train and test our model using the spectroscopic predictions of the MILES stellar population library and non-parametric SFHs. We reliably estimate the mass assembly of an integrated stellar population with well-calibrated uncertainties. Specifically, we reach a score of $0.97\,R^2$ for the time at which a given galaxy from the test set formed $50\%$ of its stellar mass, obtaining samples of the posteriors in only $10^{-4}$\,s. We then applied the pipeline to real observations of massive elliptical galaxies, recovering the well-known relationship between the age and the velocity dispersion, and show that the most massive galaxies ($σ\sim300$ km/s) built up to 90\% of their total stellar masses within $1$\,Gyr of the Big Bang. The inferred properties also agree with the state-of-the-art inversion codes, but the inference is performed up to five orders of magnitude faster. This SBI approach coupled with machine learning and applied to full spectral fitting makes it possible to address large numbers of galaxies while performing a thick sampling of the posteriors. It will allow both the deterministic trends and the inherent uncertainties of the highly degenerated inversion problem to be estimated for large and complex upcoming spectroscopic surveys, such as DESI, WEAVE, or 4MOST. △ Less

Submitted 28 June, 2024; v1 submitted 26 June, 2024; originally announced June 2024.

Comments: 12 pages, 9 figures, 1 table. Accepted for publication on A&A

Retrieval of the physical parameters of galaxies from WEAVE-StePS-like data using machine learning

Authors: J. Angthopo, B. R. Granett, F. La Barbera, M. Longhetti, A. Iovino, M. Fossati, F. R. Ditrani, L. Costantin, S. Zibetti, A. Gallazzi, P. Sánchez-Blázquez, C. Tortora, C. Spiniello, B. Poggianti, A. Vazdekis, M. Balcells, S. Bardelli, C. R. Benn, M. Bianconi, M. Bolzonella, G. Busarello, L. P. Cassarà, E. M. Corsini, O. Cucciati, G. Dalton , et al. (24 additional authors not shown)

Abstract: The WHT Enhanced Area Velocity Explorer (WEAVE) is a new, massively multiplexing spectrograph. This new instrument will be exploited to obtain high S/N spectra of $\sim$25000 galaxies at intermediate redshifts for the WEAVE Stellar Population Survey (WEAVE-StePS). We test machine learning methods for retrieving the key physical parameters of galaxies from WEAVE-StePS-like spectra using both photom… ▽ More The WHT Enhanced Area Velocity Explorer (WEAVE) is a new, massively multiplexing spectrograph. This new instrument will be exploited to obtain high S/N spectra of $\sim$25000 galaxies at intermediate redshifts for the WEAVE Stellar Population Survey (WEAVE-StePS). We test machine learning methods for retrieving the key physical parameters of galaxies from WEAVE-StePS-like spectra using both photometric and spectroscopic information at various S/Ns and redshifts. We simulated $\sim$105000 galaxy spectra assuming SFH with an exponentially declining star formation rate, covering a wide range of ages, stellar metallicities, sSFRs, and dust extinctions. We then evaluated the ability of the random forest and KNN algorithms to correctly predict such parameters assuming no measurement errors. We checked how much the predictive ability deteriorates for different S/Ns and redshifts, finding that both algorithms still accurately estimate the ages and metallicities with low bias. The dispersion varies from 0.08-0.16 dex for ages and 0.11-0.25 dex for metallicity, depending on the redshift and S/N. For dust attenuation, we find a similarly low bias and dispersion. For the sSFR, we find a very good constraining power for star-forming galaxies, log sSFR$\gtrsim$ -11, where the bias is $\sim$ 0.01 dex and the dispersion is $\sim$ 0.10 dex. For more quiescent galaxies, with log sSFR$\lesssim$ -11, we find a higher bias, 0.61-0.86 dex, and a higher dispersion, $\sim$ 0.4 dex, for different S/Ns and redshifts. Generally, we find that the RF outperforms the KNN. Finally, the retrieved sSFR was used to successfully classify galaxies as part of the blue cloud, green valley, or red sequence. We demonstrate that machine learning algorithms can accurately estimate the physical parameters of simulated galaxies even at relatively low S/N=10 per angstrom spectra with available ancillary photometric information. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: 19 pages, 10 + 2 figures, 4 tables, accepted in A&A

Euclid preparation. Detecting globular clusters in the Euclid survey

Authors: Euclid Collaboration, K. Voggel, A. Lançon, T. Saifollahi, S. S. Larsen, M. Cantiello, M. Rejkuba, J. -C. Cuillandre, P. Hudelot, A. A. Nucita, M. Urbano, E. Romelli, M. A. Raj, M. Schirmer, C. Tortora, Abdurro'uf, F. Annibali, M. Baes, P. Boldrini, R. Cabanac, D. Carollo, C. J. Conselice, P. -A. Duc, A. M. N. Ferguson, L. K. Hunt , et al. (247 additional authors not shown)

Abstract: Extragalactic globular clusters (EGCs) are an abundant and powerful tracer of galaxy dynamics and formation, and their own formation and evolution is also a matter of extensive debate. The compact nature of globular clusters means that they are hard to spatially resolve and thus study outside the Local Group. In this work we have examined how well EGCs will be detectable in images from the Euclid… ▽ More Extragalactic globular clusters (EGCs) are an abundant and powerful tracer of galaxy dynamics and formation, and their own formation and evolution is also a matter of extensive debate. The compact nature of globular clusters means that they are hard to spatially resolve and thus study outside the Local Group. In this work we have examined how well EGCs will be detectable in images from the Euclid telescope, using both simulated pre-launch images and the first early-release observations of the Fornax galaxy cluster. The Euclid Wide Survey will provide high-spatial resolution VIS imaging in the broad IE band as well as near-infrared photometry (YE, JE, and HE). We estimate that the galaxies within 100 Mpc in the footprint of the Euclid survey host around 830 000 EGCs of which about 350 000 are within the survey's detection limits. For about half of these EGCs, three infrared colours will be available as well. For any galaxy within 50Mpc the brighter half of its GC luminosity function will be detectable by the Euclid Wide Survey. The detectability of EGCs is mainly driven by the residual surface brightness of their host galaxy. We find that an automated machine-learning EGC-classification method based on real Euclid data of the Fornax galaxy cluster provides an efficient method to generate high purity and high completeness GC candidate catalogues. We confirm that EGCs are spatially resolved compared to pure point sources in VIS images of Fornax. Our analysis of both simulated and first on-sky data show that Euclid will increase the number of GCs accessible with high-resolution imaging substantially compared to previous surveys, and will permit the study of GCs in the outskirts of their hosts. Euclid is unique in enabling systematic studies of EGCs in a spatially unbiased and homogeneous manner and is primed to improve our understanding of many understudied aspects of GC astrophysics. △ Less

Submitted 29 May, 2024; v1 submitted 22 May, 2024; originally announced May 2024.

Comments: Submitted to A&A

Euclid: Early Release Observations -- Globular clusters in the Fornax galaxy cluster, from dwarf galaxies to the intracluster field

Authors: T. Saifollahi, K. Voggel, A. Lançon, Michele Cantiello, M. A. Raj, J. -C. Cuillandre, S. S. Larsen, F. R. Marleau, A. Venhola, M. Schirmer, D. Carollo, P. -A. Duc, A. M. N. Ferguson, L. K. Hunt, M. Kümmel, R. Laureijs, O. Marchal, A. A. Nucita, R. F. Peletier, M. Poulain, M. Rejkuba, R. Sánchez-Janssen, M. Urbano, Abdurro'uf, B. Altieri , et al. (174 additional authors not shown)

Abstract: We present an analysis of Euclid observations of a 0.5 deg$^2$ field in the central region of the Fornax galaxy cluster that were acquired during the performance verification phase. With these data, we investigate the potential of Euclid for identifying GCs at 20 Mpc, and validate the search methods using artificial GCs and known GCs within the field from the literature. Our analysis of artificial… ▽ More We present an analysis of Euclid observations of a 0.5 deg$^2$ field in the central region of the Fornax galaxy cluster that were acquired during the performance verification phase. With these data, we investigate the potential of Euclid for identifying GCs at 20 Mpc, and validate the search methods using artificial GCs and known GCs within the field from the literature. Our analysis of artificial GCs injected into the data shows that Euclid's data in $I_{\rm E}$ band is 80% complete at about $I_{\rm E} \sim 26.0$ mag ($M_{V\rm } \sim -5.0$ mag), and resolves GCs as small as $r_{\rm h} = 2.5$ pc. In the $I_{\rm E}$ band, we detect more than 95% of the known GCs from previous spectroscopic surveys and GC candidates of the ACS Fornax Cluster Survey, of which more than 80% are resolved. We identify more than 5000 new GC candidates within the field of view down to $I_{\rm E}$ mag, about 1.5 mag fainter than the typical GC luminosity function turn-over magnitude, and investigate their spatial distribution within the intracluster field. We then focus on the GC candidates around dwarf galaxies and investigate their numbers, stacked luminosity distribution and stacked radial distribution. While the overall GC properties are consistent with those in the literature, an interesting over-representation of relatively bright candidates is found within a small number of relatively GC-rich dwarf galaxies. Our work confirms the capabilities of Euclid data in detecting GCs and separating them from foreground and background contaminants at a distance of 20 Mpc, particularly for low-GC count systems such as dwarf galaxies. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue `Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

Euclid: Early Release Observations -- Deep anatomy of nearby galaxies

Authors: L. K. Hunt, F. Annibali, J. -C. Cuillandre, A. M. N. Ferguson, P. Jablonka, S. S. Larsen, F. R. Marleau, E. Schinnerer, M. Schirmer, C. Stone, C. Tortora, T. Saifollahi, A. Lançon, M. Bolzonella, S. Gwyn, M. Kluge, R. Laureijs, D. Carollo, M. L. M. Collins, P. Dimauro, P. -A. Duc, D. Erkal, J. M. Howell, C. Nally, E. Saremi , et al. (174 additional authors not shown)

Abstract: Euclid is poised to make significant advances in the study of nearby galaxies in the local Universe. Here we present a first look at 6 galaxies observed for the Nearby Galaxy Showcase as part of the Euclid Early Release Observations acquired between August and November, 2023. These targets, 3 dwarf galaxies (HolmbergII, IC10, NGC6822) and 3 spirals (IC342, NGC2403, NGC6744), range in distance from… ▽ More Euclid is poised to make significant advances in the study of nearby galaxies in the local Universe. Here we present a first look at 6 galaxies observed for the Nearby Galaxy Showcase as part of the Euclid Early Release Observations acquired between August and November, 2023. These targets, 3 dwarf galaxies (HolmbergII, IC10, NGC6822) and 3 spirals (IC342, NGC2403, NGC6744), range in distance from about 0.5 Mpc to 8.8 Mpc. Our assessment of the surface brightness depths in the stacked Euclid images confirms previous estimates in 100 arcsec^2 regions of 1sigma=30.5 mag/arcsec^2 for VIS, but slightly deeper than previous estimates for NISP with 1sigma=29.2-29.4 mag/arcsec^2. By combining Euclid HE, YE, and IE into RGB images, we illustrate the large field-of-view covered by a single Reference Observing Sequence, together with exquisite detail on parsec scales in these nearby galaxies. Radial surface brightness and color profiles demonstrate galaxy colors in agreement with stellar population synthesis models. Standard stellar photometry selection techniques find approximately 1.3 million stars across the 6 galaxy fields. Euclid's resolved stellar photometry allows us to constrain the star-formation histories of these galaxies, by disentangling the distributions of young stars, as well as asymptotic giant branch and red giant branch stellar populations. We finally examine 2 galaxies individually for surrounding satellite systems. Our analysis of the ensemble of dwarf satellites around NGC6744 reveals a new galaxy, EDwC1, a nucleated dwarf spheroidal at the end of a spiral arm. Our new census of the globular clusters around NGC2403 yields 9 new star-cluster candidates, 8 of which with colors indicative of evolved stellar populations. In summary, our investigation of the 6 Showcase galaxies demonstrates that Euclid is a powerful probe of the anatomy of nearby galaxies [abridged]. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: 36 pages; 20 figures in main text; 4 Appendices. Submitted to A&A, as part of the A&A special issue `Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

Euclid: Early Release Observations -- Programme overview and pipeline for compact- and diffuse-emission photometry

Authors: J. -C. Cuillandre, E. Bertin, M. Bolzonella, H. Bouy, S. Gwyn, S. Isani, M. Kluge, O. Lai, A. Lançon, D. A. Lang, R. Laureijs, T. Saifollahi, M. Schirmer, C. Stone, Abdurro'uf, N. Aghanim, B. Altieri, F. Annibali, H. Atek, P. Awad, M. Baes, E. Bañados, D. Barrado, S. Belladitta, V. Belokurov , et al. (240 additional authors not shown)

Abstract: The Euclid ERO showcase Euclid's capabilities in advance of its main mission, targeting 17 astronomical objects, from galaxy clusters, nearby galaxies, globular clusters, to star-forming regions. A total of 24 hours observing time was allocated in the early months of operation, engaging the scientific community through an early public data release. We describe the development of the ERO pipeline t… ▽ More The Euclid ERO showcase Euclid's capabilities in advance of its main mission, targeting 17 astronomical objects, from galaxy clusters, nearby galaxies, globular clusters, to star-forming regions. A total of 24 hours observing time was allocated in the early months of operation, engaging the scientific community through an early public data release. We describe the development of the ERO pipeline to create visually compelling images while simultaneously meeting the scientific demands within months of launch, leveraging a pragmatic, data-driven development strategy. The pipeline's key requirements are to preserve the image quality and to provide flux calibration and photometry for compact and extended sources. The pipeline's five pillars are: removal of instrumental signatures; astrometric calibration; photometric calibration; image stacking; and the production of science-ready catalogues for both the VIS and NISP instruments. We report a PSF with a full width at half maximum of 0.16" in the optical and 0.49" in the three NIR bands. Our VIS mean absolute flux calibration is accurate to about 1%, and 10% for NISP due to a limited calibration set; both instruments have considerable colour terms. The median depth is 25.3 and 23.2 AB mag with a SNR of 10 for galaxies, and 27.1 and 24.5 AB mag at an SNR of 5 for point sources for VIS and NISP, respectively. Euclid's ability to observe diffuse emission is exceptional due to its extended PSF nearly matching a pure diffraction halo, the best ever achieved by a wide-field, high-resolution imaging telescope. Euclid offers unparalleled capabilities for exploring the LSB Universe across all scales, also opening a new observational window in the NIR. Median surface-brightness levels of 29.9 and 28.3 AB mag per square arcsec are achieved for VIS and NISP, respectively, for detecting a 10 arcsec x 10 arcsec extended feature at the 1 sigma level. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Submitted to A&A, 44 pages, 36 figures - Part of the A&A special issue `Euclid on Sky', which contains Euclid key reference papers and first results from the Euclid Early Release Observations

Euclid. I. Overview of the Euclid mission

Authors: Euclid Collaboration, Y. Mellier, Abdurro'uf, J. A. Acevedo Barroso, A. Achúcarro, J. Adamek, R. Adam, G. E. Addison, N. Aghanim, M. Aguena, V. Ajani, Y. Akrami, A. Al-Bahlawan, A. Alavi, I. S. Albuquerque, G. Alestas, G. Alguero, A. Allaoui, S. W. Allen, V. Allevato, A. V. Alonso-Tetilla, B. Altieri, A. Alvarez-Candal, A. Amara, L. Amendola , et al. (1086 additional authors not shown)

Abstract: The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14… ▽ More The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue`Euclid on Sky'

Detection of an intranight optical hard-lag with colour variability in blazar PKS 0735+178

Authors: Callum McCall, Helen E. Jermak, Iain A. Steele, Shiho Kobayashi, Johan H. Knapen, Pablo M. Sánchez-Alarcón

Abstract: Blazars are a highly variable subclass of active galactic nuclei that have been observed to vary significantly during a single night. This intranight variability remains a debated phenomenon, with various mechanisms proposed to explain the behaviour including jet energy density evolution or system geometric changes. We present the results of an intranight optical monitoring campaign of four blazar… ▽ More Blazars are a highly variable subclass of active galactic nuclei that have been observed to vary significantly during a single night. This intranight variability remains a debated phenomenon, with various mechanisms proposed to explain the behaviour including jet energy density evolution or system geometric changes. We present the results of an intranight optical monitoring campaign of four blazars: TXS 0506+056, OJ287, PKS 0735+178, and OJ248 using the Carlos Sánchez Telescope. We detect significant but colourless behaviour in OJ287 and both bluer- and redder-when-brighter colour trends in PKS 0735+178. Additionally, the g band shows a lag of ~10 min with respect to the r,i,z_s bands for PKS 0735+178 on 2023 January 17. This unexpected hard-lag in PKS 0735+178 is not in accordance with the standard synchrotron shock cooling model (which would predict a soft lag) and instead suggests the variability may be a result of changes in the jet's electron energy density distribution, with energy injection from Fermi acceleration processes into a post-shocked medium. △ Less

Submitted 1 February, 2024; originally announced February 2024.

Comments: 19 pages, 19 figures, 5 tables

Strategies for optimal sky subtraction in the low surface brightness regime

Authors: A. E. Watkins, S. Kaviraj, C. C. Collins, J. H. Knapen, L. S. Kelvin, P. -A. Duc, J. Román, J. C. Mihos

Abstract: The low surface brightness (LSB) regime ($μ_{g} \gtrsim 26$ mag arcsec$^{-2}$) comprises a vast, mostly unexplored discovery space, from dwarf galaxies to the diffuse interstellar medium. Accessing this regime requires precisely removing instrumental signatures and light contamination, including, most critically, night sky emission. This is not trivial, as faint astrophysical and instrumental cont… ▽ More The low surface brightness (LSB) regime ($μ_{g} \gtrsim 26$ mag arcsec$^{-2}$) comprises a vast, mostly unexplored discovery space, from dwarf galaxies to the diffuse interstellar medium. Accessing this regime requires precisely removing instrumental signatures and light contamination, including, most critically, night sky emission. This is not trivial, as faint astrophysical and instrumental contamination can bias sky models at the precision needed to characterize LSB structures. Using idealized synthetic images, we assess how this bias impacts two common LSB-oriented sky-estimation algorithms: 1.) masking and parametric modelling, and 2.) stacking and smoothing dithered exposures. Undetected flux limits both methods by imposing a pedestal offset to all derived sky models. Careful, deep masking of fixed sources can mitigate this, but source density always imposes a fundamental limit. Stellar scattered light can contribute $\sim28$--$29$ mag arcsec$^{-2}$ of background flux even in low-density fields; its removal is critical prior to sky estimation. For complex skies, image combining is an effective non-parametric approach, although it strongly depends on observing strategy and adds noise to images on the smoothing kernel scale. Preemptive subtraction of fixed sources may be the only practical approach for robust sky estimation. We thus tested a third algorithm, subtracting a preliminary sky-subtracted coadd from exposures to isolate sky emission. Unfortunately, initial errors in sky estimation propagate through all subsequent sky models, making the method impractical. For large-scale surveys like LSST, where key science goals constrain observing strategy, masking and modelling remains the optimal sky estimation approach, assuming stellar scattered light is removed first. △ Less

Submitted 22 January, 2024; originally announced January 2024.

Comments: 18 pages, 11 figures, accepted for publication in MNRAS

Bar Properties as a Function of Wavelength: A Local Baseline with S4G for High-Redshift Studies

Authors: Karín Menéndez-Delmestre, Thiago S. Gonçalves, Kartik Sheth, Tomás Düringer Jacques de Lima, Taehyun Kim, Dimitri A. Gadotti, Eva Schinnerer, E. Athanassoula, Albert Bosma, Debra Meloy Elmegreen, Johan H. Knapen, Rubens E. G. Machado, Heikki Salo

Abstract: The redshift evolution of bars is an important signpost of the dynamic maturity of disk galaxies. To characterize the intrinsic evolution safe from band-shifting effects, it is necessary to gauge how bar properties vary locally as a function of wavelength. We investigate bar properties in 16 nearby galaxies from the Spitzer Survey of Stellar Structure in Galaxies (S4G) at ultraviolet, optical and… ▽ More The redshift evolution of bars is an important signpost of the dynamic maturity of disk galaxies. To characterize the intrinsic evolution safe from band-shifting effects, it is necessary to gauge how bar properties vary locally as a function of wavelength. We investigate bar properties in 16 nearby galaxies from the Spitzer Survey of Stellar Structure in Galaxies (S4G) at ultraviolet, optical and mid-infrared wavebands. Based on the ellipticity and position angle profiles from fitting elliptical isophotes to the two-dimensional light distribution, we find that both bar length and ellipticity - the latter often used as a proxy for bar strength - increase at bluer wavebands. Bars are 9% longer in the B-band than at 3.6 um. Their ellipticity increases typically by 8% in the B-band, with a significant fraction (>40%) displaying an increase up to 35%. We attribute the increase in bar length to the presence of star forming knots at the end of bars: these regions are brighter in bluer bands, stretching the bar signature further out. The increase in bar ellipticity could be driven by the apparent bulge size: the bulge is less prominent at bluer bands, allowing for thinner ellipses within the bar region. Alternatively, it could be due to younger stellar populations associated to the bar. The resulting effect is that bars appear longer and thinner at bluer wavebands. This indicates that band-shifting effects are significant and need to be corrected for high-redshift studies to reliably gauge any intrinsic evolution of the bar properties with redshift. △ Less

Submitted 7 December, 2023; originally announced December 2023.

Comments: Accepted to MNRAS

Gnuastro: Estimating the Zero Point Magnitude in Astronomical Imaging

Authors: Sepideh Eskandarlou, Mohammad Akhlaghi, Raúl Infante-Sainz, Elham Saremi, Samane Raji, Zahra Sharbaf, Giulia Golini, Zohreh Ghaffari, Johan H. Knapen

Abstract: Calibration of pixel values is a fundamental step for accurate measurements in astronomical imaging. In astronomical jargon this is known as estimating zero point magnitude. Here, we introduce a newly added script in GNU Astronomy Utilities (Gnuastro) version 0.20 for the zero point magnitude estimation, named: astscript-zeropoint. The script offers numerous features, such as the flexibility to us… ▽ More Calibration of pixel values is a fundamental step for accurate measurements in astronomical imaging. In astronomical jargon this is known as estimating zero point magnitude. Here, we introduce a newly added script in GNU Astronomy Utilities (Gnuastro) version 0.20 for the zero point magnitude estimation, named: astscript-zeropoint. The script offers numerous features, such as the flexibility to use either image(s) or a catalog as the reference dataset. Additionally, steps are parallelized to enhance efficiency for big data. Thanks to Gnuastro's minimal dependencies, the script is both flexible and portable. The figures of this research note are reproducible with Maneage, on the Git commit c89275e. △ Less

Submitted 7 December, 2023; originally announced December 2023.

Comments: Supplementary data on Zenodo (https://doi.org/10.5281/zenodo.10256845), project source on Codeberg (https://codeberg.org/gnuastro/papers/src/branch/zeropoint) and archived on Software Heritage (swh:1:dir:8b2d1f63be96de3de03aa3e2bb68fa7fa52df56f; see text for usage)

Preparing for low surface brightness science with the Vera C. Rubin Observatory: A Comparison of Observable and Simulated Intracluster Light Fractions

Authors: Sarah Brough, Syeda Lammim Ahad, Yannick M. Bahe, Amaël Ellien, Anthony H. Gonzalez, Yolanda Jiménez-Teja, Lucas C. Kimmig, Garreth Martin, Cristina Martínez-Lombilla, Mireia Montes, Annalisa Pillepich, Rossella Ragusa, Rhea-Silvia Remus, Chris A. Collins, Johan H. Knapen, J. Chris Mihos

Abstract: Intracluster Light (ICL) provides an important record of the interactions galaxy clusters have undergone. However, we are limited in our understanding by our measurement methods. To address this we measure the fraction of cluster light that is held in the Brightest Cluster Galaxy and ICL (BCG+ICL fraction) and the ICL alone (ICL fraction) using observational methods (Surface Brightness Threshold-S… ▽ More Intracluster Light (ICL) provides an important record of the interactions galaxy clusters have undergone. However, we are limited in our understanding by our measurement methods. To address this we measure the fraction of cluster light that is held in the Brightest Cluster Galaxy and ICL (BCG+ICL fraction) and the ICL alone (ICL fraction) using observational methods (Surface Brightness Threshold-SB, Non-Parametric Measure-NP, Composite Models-CM, Multi-Galaxy Fitting-MGF) and new approaches under development (Wavelet Decomposition-WD) applied to mock images of 61 galaxy clusters (14<log10 M_200c/M_solar <14.5) from four cosmological hydrodynamical simulations. We compare the BCG+ICL and ICL fractions from observational measures with those using simulated measures (aperture and kinematic separations). The ICL fractions measured by kinematic separation are significantly larger than observed fractions. We find the measurements are related and provide equations to estimate kinematic ICL fractions from observed fractions. The different observational techniques give consistent BCG+ICL and ICL fractions but are biased to underestimating the BCG+ICL and ICL fractions when compared with aperture simulation measures. Comparing the different methods and algorithms we find that the MGF algorithm is most consistent with the simulations, and CM and SB methods show the smallest projection effects for the BCG+ICL and ICL fractions respectively. The Ahad (CM), MGF and WD algorithms are best set up to process larger samples, however, the WD algorithm in its current form is susceptible to projection effects. We recommend that new algorithms using these methods are explored to analyse the massive samples that Rubin Observatory's Legacy Survey of Space and Time will provide. △ Less

Submitted 29 November, 2023; originally announced November 2023.

Comments: Submitted for publication in MNRAS, posted to arXiv after responding to two positive rounds of referee comments. Key results in Figs 3, 5, 6 and 11

Constraining the top-light initial mass function in the extended ultraviolet disk of M83

Authors: R. P. V. Rautio, A. E. Watkins, H. Salo, A. Venhola, J. H. Knapen, S. Comerón

Abstract: The universality or non-universality of the initial mass function (IMF) has significant implications for determining star formation rates and star formation histories from photometric properties of stellar populations. We reexamine whether the IMF is deficient in high-mass stars (top-light) in the low-density environment of the outer disk of M83 and constrain the shape of the IMF therein. Using ar… ▽ More The universality or non-universality of the initial mass function (IMF) has significant implications for determining star formation rates and star formation histories from photometric properties of stellar populations. We reexamine whether the IMF is deficient in high-mass stars (top-light) in the low-density environment of the outer disk of M83 and constrain the shape of the IMF therein. Using archival Galaxy Evolution Explorer (GALEX) far ultraviolet (FUV) and near ultraviolet (NUV) data and new deep OmegaCAM narrowband H$α$ imaging, we constructed a catalog of FUV-selected objects in the outer disk of M83. We counted H$α$-bright clusters and clusters that are blue in FUV$-$NUV in the catalog, measured the maximum flux ratio $F_{\mathrm{H}α}/f_{λ\mathrm{FUV}}$ among the clusters, and measured the total flux ratio $ΣF_{\mathrm{H}α}/Σf_{λ\mathrm{FUV}}$ over the catalog. We then compared these measurements to predictions from stellar population synthesis models made with a standard Salpeter IMF, truncated IMFs, and steep IMFs. We also investigated the effect of varying the assumed internal extinction on our results. We are not able to reproduce our observations with models using the standard Salpeter IMF or the truncated IMFs. It is only when assuming an average internal extinction of $0.10 < A_{\mathrm{V}} < 0.15$ in the outer disk stellar clusters that models with steep IMFs ($α> 3.1$) simultaneously reproduce the observed cluster counts, the maximum observed $F_{\mathrm{H}α}/f_{λ\mathrm{FUV}}$, and the observed $ΣF_{\mathrm{H}α}/Σf_{λ\mathrm{FUV}}$. Our results support a non-universal IMF that is deficient in high-mass stars in low-density environments. △ Less

Submitted 6 November, 2023; originally announced November 2023.

Comments: 18 pages, 15 figures, accepted to Astronomy & Astrophysics

Journal ref: A&A 681, A76 (2024)

Stellar metallicity from optical and UV spectral indices: Test case for WEAVE-StePS

Authors: F. R. Ditrani, M. Longhetti, F. La Barbera, A. Iovino, L. Costantin, S. Zibetti, A. Gallazzi, M. Fossati, J. Angthopo, Y. Ascasibar, B. Poggianti, P. Sánchez-Blázquez, M. Balcells, M. Bianconi, M. Bolzonella, L. P. Cassarà, O. Cucciati, G. Dalton, A. Ferré-Mateu, R. García-Benito, B. Granett, M. Gullieuszik, A. Ikhsanova, S. Jin, J. H. Knapen , et al. (13 additional authors not shown)

Abstract: The upcoming generation of optical spectrographs on four meter-class telescopes, with their huge multiplexing capabilities, excellent spectral resolution, and unprecedented wavelength coverage, will provide high-quality spectra for thousands of galaxies. These data will allow us to examine of the stellar population properties at intermediate redshift, an epoch that remains unexplored by large and… ▽ More The upcoming generation of optical spectrographs on four meter-class telescopes, with their huge multiplexing capabilities, excellent spectral resolution, and unprecedented wavelength coverage, will provide high-quality spectra for thousands of galaxies. These data will allow us to examine of the stellar population properties at intermediate redshift, an epoch that remains unexplored by large and deep surveys. We assess our capability to retrieve the mean stellar metallicity in galaxies at different redshifts and S/N, while simultaneously exploiting the UV and optical rest-frame wavelength coverage. The work is based on a comprehensive library of spectral templates of stellar populations, covering a wide range of age and metallicity values and built assuming various SFHs. We simulated realistic observations of a large sample of galaxies carried out with WEAVE at the WHT at different redshifts and S/N values. We measured all the reliable indices on the simulated spectra and on the comparison library. We then adopted a Bayesian approach to obtain the probability distribution of stellar metallicity. The analysis of the spectral indices has shown how some mid-UV indices can provide reliable constraints on stellar metallicity, along with optical indicators. The analysis of the mock observations has shown that even at S/N=10, the metallicity can be derived within 0.3 dex, in particular, for stellar populations older than 2 Gyr. Our results are in good agreement with other theoretical and observational works in the literature and show how the UV indicators can be advantageous in constraining metallicities. This is very promising for the upcoming surveys carried out with new, highly multiplexed, large-field spectrographs, such as StePS at the WEAVE and 4MOST, which will provide spectra of thousands of galaxies covering large spectral ranges at relatively high S/N. △ Less

Submitted 4 August, 2023; originally announced August 2023.

Comments: 16 pages, 14 figures. Submitted 31/03/2023, Accepted 20/07/2023

Journal ref: A&A 677, A93 (2023)

The AMIGA sample of isolated galaxies. XIV. Disc breaks and interactions through ultra-deep optical imaging

Authors: P. M. Sánchez-Alarcón, J. Román, J. H. Knapen, L. Verdes-Montenegro, S. Comerón, R. M. Rich, J. E. Beckman, M. Argudo-Fernández, P. Ramírez-Moreta, J. Blasco, E. Unda-Sanzana, J. Garrido, S. Sánchez-Exposito

Abstract: In the standard cosmological model of galaxy evolution, mergers and interactions play a fundamental role in shaping galaxies. Galaxies that are currently isolated are thus interesting, allowing us to distinguish between internal and external processes affecting the galactic structure. However, current observational limits may obscure crucial information in the low-mass or low-brightness regime. We… ▽ More In the standard cosmological model of galaxy evolution, mergers and interactions play a fundamental role in shaping galaxies. Galaxies that are currently isolated are thus interesting, allowing us to distinguish between internal and external processes affecting the galactic structure. However, current observational limits may obscure crucial information in the low-mass or low-brightness regime. We use optical imaging of a subsample of the AMIGA catalogue of isolated galaxies to explore the impact of different factors on the structure of these galaxies. We study the type of disc break as a function of the degree of isolation and the presence of interaction indicators like tidal streams or plumes only detectable in the low surface brightness regime. We present deep optical imaging of a sample of 25 isolated galaxies. Through careful data processing and analysis techniques, the surface brightness limits achieved are comparable to those to be obtained on the 10-year LSST coadds. The extreme depth of our imaging allows us to study the interaction signatures of 20 galaxies, given that the presence of Galactic cirrus is a strong limiting factor in the characterisation of interactions for the remaining 5 of them. We detect previously unreported interaction features in 8 (40%) galaxies in our sample. We identify 9 galaxies (36%) showing an exponential disc (Type I), 14 galaxies (56%) with down-bending (Type II) profile and only 2 galaxies (8%) with up-bending (Type III) profiles. Isolated galaxies have considerably more purely exponential discs and fewer up-bending surface brightness profiles than field or cluster galaxies. We suggest that major mergers produce up-bending profiles while a threshold in star formation probably forms down-bending profiles. Unperturbed galaxies, evolving slowly with a low star formation rate could cause the high rate of Type I discs in isolated galaxies observed. △ Less

Submitted 5 July, 2023; originally announced July 2023.

Comments: 20 pages, 7 figures, and 2 tables (incl. 2 appendix). Accepted for publication in A&A. Images and profiles will be made publicly available upon publication

Journal ref: A&A 677, A117 (2023)

The truncation of the disk of NGC 4565: Detected up to z=4 kpc, with star formation, and affected by the warp

Authors: Cristina Martinez-Lombilla, Raul Infante-Sainz, Felipe Jimenez-Ibarra, Johan H. Knapen, Ignacio Trujillo, Sebastien Comeron, Alejandro S. Borlaff, Javier Roman

Abstract: Context: The hierarchical model of galaxy formation suggests that galaxies are continuously growing. However, our position inside the Milky Way prevents us from studying the disk edge. Truncations are low surface brightness features located in the disk outskirts of external galaxies. They indicate where the disk brightness abruptly drops and their location is thought to change dynamically. In prev… ▽ More Context: The hierarchical model of galaxy formation suggests that galaxies are continuously growing. However, our position inside the Milky Way prevents us from studying the disk edge. Truncations are low surface brightness features located in the disk outskirts of external galaxies. They indicate where the disk brightness abruptly drops and their location is thought to change dynamically. In previous analyses of Milky Way-like galaxies, truncations were detected up to 3 kpc above the mid-plane but whether they remain present beyond that height remains unclear. Aims: Our goal is to determine whether truncations can be detected above 3 kpc height in the Milky Way-like galaxy NGC 4565, thus establishing the actual disk thickness. We also aim to study how the truncation relates to disk properties such as star formation activity or the warp. Methods: We perform a vertical study of the disk of NGC 4565 edge in unprecedented detail. We explore the truncation radius at different heights above/below the disk mid-plane (0<z<8 kpc) and at different wavelengths. We use new ultra-deep optical data ($μ_{g,\rm{lim}}=30.5$ mag arcsec$^{-2}$; $3 σ$ within $10 \times 10$ arcsec$^{2}$ boxes) in the $g$, $r$ and $i$ broad bands, along with near- and far-ultraviolet, H$α$, and \ion{H}{i} observations. Results: We detect the truncation up to 4 kpc in the $g$, $r$ and $i$ ultra-deep bands which is 1 kpc higher than in any previous study for any galaxy. The radial position of the truncation remains constant up to 3 kpc while higher up it is located at a smaller radius. This result is independent of the wavelength but is affected by the presence of the warp. Conclusions: We propose an inside-out growth scenario for the formation of the disk of NGC 4565. Our results point towards the truncation feature being linked to a star-forming threshold and to the onset of the disk warp. △ Less

Submitted 13 July, 2023; v1 submitted 3 July, 2023; originally announced July 2023.

Comments: 27 pages, 18 figures (incl. 2 appendix); accepted for publication in A&A; Fixed labels in Fig. 6

Journal ref: A&A 678, A62 (2023)

A giant thin stellar stream in the Coma Galaxy Cluster

Authors: Javier Román, R. Michael Rich, Niusha Ahvazi, Laura Sales, Chester Li, Giulia Golini, Ignacio Trujillo, Johan H. Knapen, Reynier F. Peletier, Pablo M. Sánchez-Alarcón

Abstract: The study of dynamically cold stellar streams reveals information about the gravitational potential where they reside and provides important constraints on dark matter properties. However, their intrinsic faintness makes detection beyond Local environments highly challenging. Here we report the detection of an extremely faint stellar stream ($μ_{g,max}=$ 29.5 mag arcsec$^{-2}$) with an extraordina… ▽ More The study of dynamically cold stellar streams reveals information about the gravitational potential where they reside and provides important constraints on dark matter properties. However, their intrinsic faintness makes detection beyond Local environments highly challenging. Here we report the detection of an extremely faint stellar stream ($μ_{g,max}=$ 29.5 mag arcsec$^{-2}$) with an extraordinarily coherent and thin morphology in the Coma Galaxy Cluster. This Giant Coma Stream spans ~510 kpc in length and appears as a free-floating structure located at a projected distance of 0.8 Mpc from the centre of Coma. We do not identify any potential galaxy remnant or core, and the stream structure appears featureless in our data. We interpret the Giant Coma Stream as being a recently accreted, tidally disrupting passive dwarf. Using the Illustris-TNG50 simulation, we identify a case with similar characteristics, showing that, although rare, these types of streams are predicted to exist in $Λ$-CDM. Our work unveils the presence of free-floating, extremely faint and thin stellar streams in galaxy clusters, widening the environmental context for their promising future applications in the study of dark matter properties. △ Less

Submitted 15 October, 2023; v1 submitted 4 May, 2023; originally announced May 2023.

Comments: Accepted for publication in A&A

Journal ref: A&A 679, A157 (2023)

LeMMINGs. VI. Connecting nuclear activity to bulge properties of active and inactive galaxies: radio scaling relations and galaxy environment

Authors: B. T. Dullo, J. H. Knapen, R. J. Beswick, R. D. Baldi, D. R. A. Williams, I. M. McHardy, D. A. Green, A. Gil de Paz, S. Aalto, A. Alberdi, M. K. Argo, H. -R. Klöckner, I. M. Mutie, D. J. Saikia, P. Saikia, I. R. Stevens

Abstract: Multiwavelength studies indicate that nuclear activity and bulge properties are closely related, but the details remain unclear. To study this further, we combine $Hubble~Space~Telescope$ bulge structural and photometric properties with 1.5 GHz, $e$-MERLIN nuclear radio continuum data from the LeMMINGs survey for a large sample of 173 `active' galaxies (LINERs and Seyferts) and `inactive' galaxies… ▽ More Multiwavelength studies indicate that nuclear activity and bulge properties are closely related, but the details remain unclear. To study this further, we combine $Hubble~Space~Telescope$ bulge structural and photometric properties with 1.5 GHz, $e$-MERLIN nuclear radio continuum data from the LeMMINGs survey for a large sample of 173 `active' galaxies (LINERs and Seyferts) and `inactive' galaxies (H IIs and absorption line galaxies, ALGs). Dividing our sample into active and inactive, they define distinct (radio core luminosity)$-$(bulge mass), L_R,core-M_*,bulge, relations, with a mass turnover at M_*, bulge ~ 10^(9.8 +- 0.3) M_sun (supermassive black hole mass M_BH ~ 10^(6.8 +- 0.3) M_sun), which marks the transition from AGN-dominated nuclear radio emission in more massive bulges to that mainly driven by stellar processes in low-mass bulges. None of our 10/173 bulgeless galaxies host an AGN. The AGN fraction increases with increasing M_*, bulge such that f_optical_AGN $\propto$ M_*,bulge^(0.24 +- 0.06) and f_radio_AGN $\propto$ M_*,bulge^(0.24 +- 0.05). Between M_*,bulge ~ 10^8.5 and 10^11.3 M_sun, f_optical_AGN steadily rises from 15 +- 4 to 80 +- 5 per cent. We find that at fixed bulge mass, the radio loudness, nuclear radio activity and the (optical and radio) AGN fraction exhibit no dependence on environment. Radio-loud hosts preferentially possess an early-type morphology than radio-quiet hosts, the two types are however indistinguishable in terms of bulge Sérsic index and ellipticity, while results on the bulge inner logarithmic profile slope are inconclusive. We finally discuss the importance of bulge mass in determining the AGN triggering processes, including potential implications for the nuclear radio emission in nearby galaxies. △ Less

Submitted 13 April, 2023; originally announced April 2023.

Comments: 27 pages, 15 figures, 4 tables, accepted for publication in MNRAS

LeMMINGs. V. Nuclear activity and bulge properties: a detailed multi-component decomposition of $e$-MERLIN Palomar galaxies with $HST$

Authors: B. T. Dullo, J. H. Knapen, R. J. Beswick, R. D. Baldi, D. R. A. Williams, I. M. McHardy, J. S. Gallagher, S. Aalto, M. K. Argo, A. Gil de Paz, H. -R. Klöckner, J. M. Marcaide, C. G. Mundell, I. M. Mutie, P. Saikia

Abstract: [Abridged] We use high-resolution $HST$ imaging and $e$-MERLIN 1.5-GHz observations of galaxy cores from the LeMMINGs survey to investigate the relation between optical structural properties and nuclear radio emission for a large sample of galaxies. We perform accurate, multi-component decompositions of new surface brightness profiles extracted from $HST$ images for 163 LeMMINGs galaxies and fit u… ▽ More [Abridged] We use high-resolution $HST$ imaging and $e$-MERLIN 1.5-GHz observations of galaxy cores from the LeMMINGs survey to investigate the relation between optical structural properties and nuclear radio emission for a large sample of galaxies. We perform accurate, multi-component decompositions of new surface brightness profiles extracted from $HST$ images for 163 LeMMINGs galaxies and fit up to six galaxy components (e.g., bulges, discs, AGN, bars, rings, spiral arms, and nuclear star clusters) simultaneously with Sérsic and/or core-Sérsic models. By adding such decomposition data for 10 LeMMINGs galaxies from our past work, the final sample of 173 nearby galaxies (102 Ss, 42 S0s, 23 Es plus 6 Irr) with bulge stellar mass (typically) M_*, bulge ~ 10^6-10^12.5 M_sun, encompasses all optical spectral classes (LINER, Seyfert, ALG and H II). We show that the bulge mass can be significantly overestimated in many galaxies when components such as bars, rings and spirals are not included in the fits. We additionally implement a Monte Carlo method to determine errors on bulge, disc and other fitted structural parameters. Moving (in the opposite direction) across the Hubble sequence, i.e., from the irregular to elliptical galaxies, we confirm that bulges become larger, more prominent and round. Such bulge dominance is associated with a brighter radio core luminosity. We also find that the radio detection fraction increases with bulge mass. At M_*,bulge > 10^11 M_sun, the radio detection fraction is 77%, declining to 24% for M_bulge < 10^10 M_sun. Furthermore, we observe core-Sérsic bulges tend to be systematically round and to possess high radio core luminosities and boxy-distorted or pure elliptical isophotes. △ Less

Submitted 20 March, 2023; originally announced March 2023.

Comments: 102 pages, 26 figures, 10 tables, accepted for publication in A&A

Journal ref: A&A 675, A105 (2023)

MIGHTEE-\HI: Possible interactions with the galaxy NGC~895

Authors: Brenda Namumba, Javier Román, Jesus Falcon Barroso, Johan H. Knapen, Ianjamasimanana Roger, Elizabeth Naluminsa, Gyula I. G. Jozsa, Marie Korsaga, Natasha Maddox, Brad Frank, Sinenhlanhla Sikhosana, Samuel Legodi, Claude Carignan, Anastasia A. Ponomareva, Tom Jarrett, Danielle Lucero, Oleg M. Smirnov, Thijs van der Hulst, D. J. Pisano, kasia Malek, Lucia Marchetti, Mattia Vaccari, Matt Jarvis, Maarten Baes, Martin Meyer , et al. (7 additional authors not shown)

Abstract: The transformation and evolution of a galaxy is strongly influenced by interactions with its environment. Neutral hydrogen (\HI) is an excellent way to trace these interactions. Here, we present \HI\ observations of the spiral galaxy NGC~895, which was previously thought to be isolated. High-sensitivity \HI\ observations from the MeerKAT large survey project MIGHTEE reveal possible interaction fea… ▽ More The transformation and evolution of a galaxy is strongly influenced by interactions with its environment. Neutral hydrogen (\HI) is an excellent way to trace these interactions. Here, we present \HI\ observations of the spiral galaxy NGC~895, which was previously thought to be isolated. High-sensitivity \HI\ observations from the MeerKAT large survey project MIGHTEE reveal possible interaction features, such as extended spiral arms, and the two newly discovered \HI\ companions, that drive us to change the narrative that it is an isolated galaxy. We combine these observations with deep optical images from the Hyper Suprime Camera to show an absence of tidal debris between NGC 895 and its companions. We do find an excess of light in the outer parts of the companion galaxy MGTH$\_$J022138.1-052631 which could be an indication of external perturbation and thus possible sign of interactions. Our analysis shows that NGC~895 is an actively star-forming galaxy with a SFR of $\mathrm{1.75 \pm 0.09 [M_{\odot}/yr]}$, a value typical for high stellar mass galaxies on the star forming main sequence. It is reasonable to state that different mechanisms may have contributed to the observed features in NGC~895 and this emphasizes the need to revisit the target with more detailed observations. Our work shows the high potential and synergy of using state-of-the-art data in both \HI\ and optical to reveal a more complete picture of galaxy environments. △ Less

Submitted 16 March, 2023; originally announced March 2023.

Comments: 14 pages, 10 figures. Accepted for publication in MNRAS

Identification of tidal features in deep optical galaxy images with Convolutional Neural Networks

Authors: H. Domínguez Sánchez, G. Martin, I. Damjanov, F. Buitrago, M. Huertas-Company, C. Bottrell, M. Bernardi, J. H. Knapen, J. Vega-Ferrero, R. Hausen, E. Kado-Fong, D. Población-Criado, H. Souchereau, O. K. Leste, B. Robertson, B. Sahelices, K. V. Johnston

Abstract: Interactions between galaxies leave distinguishable imprints in the form of tidal features which hold important clues about their mass assembly. Unfortunately, these structures are difficult to detect because they are low surface brightness features so deep observations are needed. Upcoming surveys promise several orders of magnitude increase in depth and sky coverage, for which automated methods… ▽ More Interactions between galaxies leave distinguishable imprints in the form of tidal features which hold important clues about their mass assembly. Unfortunately, these structures are difficult to detect because they are low surface brightness features so deep observations are needed. Upcoming surveys promise several orders of magnitude increase in depth and sky coverage, for which automated methods for tidal feature detection will become mandatory. We test the ability of a convolutional neural network to reproduce human visual classifications for tidal detections. We use as training $\sim$6000 simulated images classified by professional astronomers. The mock Hyper Suprime Cam Subaru (HSC) images include variations with redshift, projection angle and surface brightness ($μ_{lim}$ =26-35 mag arcsec$^{-2}$). We obtain satisfactory results with accuracy, precision and recall values of Acc=0.84, P=0.72 and R=0.85, respectively, for the test sample. While the accuracy and precision values are roughly constant for all surface brightness, the recall (completeness) is significantly affected by image depth. The recovery rate shows strong dependence on the type of tidal features: we recover all the images showing shell features and 87% of the tidal streams; these fractions are below 75% for mergers, tidal tails and bridges. When applied to real HSC images, the performance of the model worsens significantly. We speculate that this is due to the lack of realism of the simulations and take it as a warning on applying deep learning models to different data domains without prior testing on the actual data. △ Less

Submitted 6 March, 2023; originally announced March 2023.

Comments: 13 pages, 10 figures, accepted for publication in MNRAS

A possible signature of the influence of tidal perturbations in dwarf galaxy scaling relations

Authors: A. E. Watkins, H. Salo, S. Kaviraj, C. A. Collins, J. H. Knapen, A. Venhola, J. Román

Abstract: Dwarf galaxies are excellent cosmological probes, because their shallow potential wells make them very sensitive to the key processes that drive galaxy evolution, including baryonic feedback, tidal interactions, and ram pressure stripping. However, some of the key parameters of dwarf galaxies, which help trace the effects of these processes, are still debated, including the relationship between th… ▽ More Dwarf galaxies are excellent cosmological probes, because their shallow potential wells make them very sensitive to the key processes that drive galaxy evolution, including baryonic feedback, tidal interactions, and ram pressure stripping. However, some of the key parameters of dwarf galaxies, which help trace the effects of these processes, are still debated, including the relationship between their sizes and masses. We re-examine the Fornax Cluster dwarf population from the point of view of isomass-radius--stellar mass relations (IRSMRs) using the Fornax Deep Survey Dwarf galaxy Catalogue, with the centrally located (among dwarfs) $3.63 \mathcal{M}_{\odot}$~pc$^{-2}$ isodensity radius defining our fiducial relation. This relation is a powerful diagnostic tool for identifying dwarfs with unusual structure, as dwarf galaxies' remarkable monotonicity in light profile shapes, as a function of stellar mass, reduces the relation's scatter tremendously. By examining how different dwarf properties (colour, tenth-nearest-neighbour distance, etc.) correlate with distance from our fiducial relation, we find a significant population of structural outliers with comparatively lower central mass surface density and larger half-light-radii, residing in locally denser regions in the cluster, albeit with similar red colours. We propose that these faint, extended outliers likely formed through tidal disturbances, which make the dwarfs more diffuse, but with little mass loss. Comparing these outliers with ultra-diffuse galaxies (UDGs), we find that the term UDG lacks discriminatory power; UDGs in the Fornax Cluster lie both on and off of IRSMRs defined at small radii, while IRSMR outliers with masses below $\sim 10^{7.5} \mathcal{M}_{\odot}$ are excluded from the UDG classification due to their small effective radii. △ Less

Submitted 27 February, 2023; originally announced February 2023.

Comments: 16 pages (+2 appendix), 10 figures, accepted for publication in MNRAS

Evidence for globular cluster collapse after a dwarf-dwarf merger: A potential nuclear star cluster in formation

Authors: J. Román, P. M. Sánchez-Alarcón, J. H. Knapen, R. Peletier

Abstract: Direct observational evidence for the creation of nuclear star clusters (NSCs) is needed to support the proposed scenarios for their formation. We analysed the dwarf galaxy UGC 7346, located in the peripheral regions of the Virgo Cluster, to highlight a series of properties that indicate the formation of a NSC caught in its earlier stages. First, we report on remnants of a past interaction in the… ▽ More Direct observational evidence for the creation of nuclear star clusters (NSCs) is needed to support the proposed scenarios for their formation. We analysed the dwarf galaxy UGC 7346, located in the peripheral regions of the Virgo Cluster, to highlight a series of properties that indicate the formation of a NSC caught in its earlier stages. First, we report on remnants of a past interaction in the form of diffuse streams or shells, suggesting a recent merging of two dwarf galaxies with a 1:5 stellar mass ratio. Second, we identify a number of globular cluster (GC) candidates that are broadly compatible in colour with the main component that is both more extended and more massive. Strikingly, we find these GCs candidates to be highly concentrated towards the centre of the galaxy (R$_{GC}$ = 0.41 R$_{e}$). We suggest that the central concentration of the GCs is likely produced by the dynamical friction of this merger. This would make UGC 7346 a unique case of a galaxy caught in the earlier stages of NSC formation. The formation of NSCs due to collapse of GCs by dynamical friction in dwarf mergers would provide a natural explanation of the environmental correlations found for the nucleation fraction for early-type dwarf galaxies, whereby denser environments host galaxies with a higher nucleation fraction. △ Less

Submitted 8 March, 2023; v1 submitted 16 February, 2023; originally announced February 2023.

Comments: Matching published version in A&A Letters

Journal ref: A&A 671, L7 (2023)

WEAVE-StePS. A stellar population survey using WEAVE at WHT

Authors: A. Iovino, B. M. Poggianti, A. Mercurio, M. Longhetti, M. Bolzonella, G. Busarello, M. Gullieuszik, F. LaBarbera, P. Merluzzi, L. Morelli, C. Tortora, D. Vergani, S. Zibetti, C. P. Haines, L. Costantin, F. R. Ditrani, L. Pozzetti, J. Angthopo, M. Balcells, S. Bardelli, C. R. Benn, M. Bianconi, L. P. Cassarà, E. M. Corsini, O. Cucciati , et al. (22 additional authors not shown)

Abstract: The upcoming new generation of optical spectrographs on four-meter-class telescopes will provide valuable opportunities for forthcoming galaxy surveys through their huge multiplexing capabilities, excellent spectral resolution, and unprecedented wavelength coverage. WEAVE is a new wide-field spectroscopic facility mounted on the 4.2 m William Herschel Telescope in La Palma. WEAVE-StePS is one of t… ▽ More The upcoming new generation of optical spectrographs on four-meter-class telescopes will provide valuable opportunities for forthcoming galaxy surveys through their huge multiplexing capabilities, excellent spectral resolution, and unprecedented wavelength coverage. WEAVE is a new wide-field spectroscopic facility mounted on the 4.2 m William Herschel Telescope in La Palma. WEAVE-StePS is one of the five extragalactic surveys that will use WEAVE during its first five years of operations. It will observe galaxies using WEAVE MOS (~950 fibres across a field of view of ~3 deg2 on the sky) in low-resolution mode (R~5000, spanning the wavelength range 3660-9590 AA). WEAVE-StePS will obtain high-quality spectra (S/N ~ 10 per AA at R~5000) for a magnitude-limited (I_AB = 20.5) sample of ~25,000 galaxies, the majority selected at z>=0.3. The survey goal is to provide precise spectral measurements in the crucial interval that bridges the gap between LEGA-C and SDSS data. The wide area coverage of ~25 deg2 will enable us to observe galaxies in a variety of environments. The ancillary data available in each observed field (including X-ray coverage, multi-narrow-band photometry and spectroscopic redshift information) will provide an environmental characterisation for each observed galaxy. This paper presents the science case of WEAVE-StePS, the fields to be observed, the parent catalogues used to define the target sample, and the observing strategy chosen after a forecast of the expected performance of the instrument for our typical targets. WEAVE-StePS will go back further in cosmic time than SDSS, extending its reach to encompass more than ~6 Gyr, nearly half of the age of the Universe. The spectral and redshift range covered by WEAVE-StePS will open a new observational window by continuously tracing the evolutionary path of galaxies in the largely unexplored intermediate-redshift range. △ Less

Submitted 14 February, 2023; originally announced February 2023.

Comments: 15 pages, 9 figures, A&A in press

Journal ref: A&A 672, A87 (2023)

On the nature of disks at high redshift seen by JWST/CEERS with contrastive learning and cosmological simulations

Authors: J. Vega-Ferrero, M. Huertas-Company, L. Costantin, P. G. Pérez-González, R. Sarmiento, J. S. Kartaltepe, A. Pillepich, M. B. Bagley, S. L. Finkelstein, E. J. McGrath, J. H. Knapen, P. Arrabal Haro, E. F. Bell, F. Buitrago, A. Calabrò, A. Dekel, M. Dickinson, H. Domínguez Sánchez, D. Elbaz, H. C. Ferguson, M. Giavalisco, B. W. Holwerda, D. D. Kocesvski, A. M. Koekemoer, V. Pandya , et al. (4 additional authors not shown)

Abstract: Visual inspections of the first optical rest-frame images from JWST have indicated a surprisingly high fraction of disk galaxies at high redshifts. Here, we alternatively apply self-supervised machine learning to explore the morphological diversity at $z \geq 3$. Our proposed data-driven representation scheme of galaxy morphologies, calibrated on mock images from the TNG50 simulation, is shown to… ▽ More Visual inspections of the first optical rest-frame images from JWST have indicated a surprisingly high fraction of disk galaxies at high redshifts. Here, we alternatively apply self-supervised machine learning to explore the morphological diversity at $z \geq 3$. Our proposed data-driven representation scheme of galaxy morphologies, calibrated on mock images from the TNG50 simulation, is shown to be robust to noise and to correlate well with the physical properties of the simulated galaxies, including their 3D structure. We apply the method simultaneously to F200W and F356W galaxy images of a mass-complete sample ($M_*/M_\odot>10^9$) at $ 3 \leq z \leq 6$ from the first JWST/NIRCam CEERS data release. We find that the simulated and observed galaxies do not exactly populate the same manifold in the representation space from contrastive learning. We also find that half the galaxies classified as disks -- either CNN-based or visually -- populate a similar region of the representation space as TNG50 galaxies with low stellar specific angular momentum and non-oblate structure. Although our data-driven study does not allow us to firmly conclude on the true nature of these galaxies, it suggests that the disk fraction at $z \geq 3$ remains uncertain and possibly overestimated by traditional supervised classifications. Deeper imaging and spectroscopic follow-ups as well as comparisons with other simulations will help to unambiguously determine the true nature of these galaxies, and establish more robust constraints on the emergence of disks at very high redshift. △ Less

Submitted 25 October, 2023; v1 submitted 14 February, 2023; originally announced February 2023.

Comments: 28 pages, 22 figures. Accepted for publication in ApJ

First detection of the outer edge of an AGN accretion disc: Very fast multiband optical variability of NGC 4395 with GTC/HiPERCAM and LT/IO:O

Authors: I. M. McHardy, M. Beard, E. Breedt, J. H. Knapen, F. M. Vincentelli, M. Veresvarska, V. S. Dhillon, T. R. Marsh, S. P. Littlefair, K. Horne, R. Glew, M. R. Goad, E. Kammoun, D. Emmanoulopoulos

Abstract: We present fast (~200s sampling) ugriz photometry of the low mass AGN NGC 4395 with the Liverpool Telescope, followed by very fast (3s sampling) us, gs, rs, is and zs simultaneous monitoring with HiPERCAM on the 10.4m GTC. These observations provide the fastest ever AGN multiband photometry and very precise lag measurements. Unlike in all other AGN, gs lags us by a large amount, consistent with di… ▽ More We present fast (~200s sampling) ugriz photometry of the low mass AGN NGC 4395 with the Liverpool Telescope, followed by very fast (3s sampling) us, gs, rs, is and zs simultaneous monitoring with HiPERCAM on the 10.4m GTC. These observations provide the fastest ever AGN multiband photometry and very precise lag measurements. Unlike in all other AGN, gs lags us by a large amount, consistent with disc reprocessing but not with reprocessing in the Broad Line Region (BLR). There is very little increase in lag with wavelength at long wavelengths, indicating an outer edge (Rout) to the reprocessor. We have compared truncated disc reprocessing models to the combined HiPERCAM and previous X-ray/UV lags. For the normally accepted mass of 3.6E5 solar, we obtain reasonable agreement with zero spin, Rout ~1700 Rg, and the DONE physically-motivated temperature-dependent disc colour correction factor (fcol). A smaller mass of 4E4 solar can only be accomodated if fcol=2.4, which is probably unrealistically high. Disc self gravity is probably unimportant in this low mass AGN but an obscuring wind may provide an edge. For the small mass the dust sublimation radius is similar to Rout, so the wind could be dusty. However for the more likely large mass the sublimation radius is further out so the optically-thick base of a line-driven gaseous wind is more likely. The inner edge of the BLR is close to Rout in both cases. These observations provide the first good evidence for a truncated AGN disc and caution that truncation should be included in reverberation lag modelling. △ Less

Submitted 15 December, 2022; originally announced December 2022.

Comments: 19 pages, 27 figures

The wide-field, multiplexed, spectroscopic facility WEAVE: Survey design, overview, and simulated implementation

Authors: Shoko Jin, Scott C. Trager, Gavin B. Dalton, J. Alfonso L. Aguerri, J. E. Drew, Jesús Falcón-Barroso, Boris T. Gänsicke, Vanessa Hill, Angela Iovino, Matthew M. Pieri, Bianca M. Poggianti, D. J. B. Smith, Antonella Vallenari, Don Carlos Abrams, David S. Aguado, Teresa Antoja, Alfonso Aragón-Salamanca, Yago Ascasibar, Carine Babusiaux, Marc Balcells, R. Barrena, Giuseppina Battaglia, Vasily Belokurov, Thomas Bensby, Piercarlo Bonifacio , et al. (190 additional authors not shown)

Abstract: WEAVE, the new wide-field, massively multiplexed spectroscopic survey facility for the William Herschel Telescope, will see first light in late 2022. WEAVE comprises a new 2-degree field-of-view prime-focus corrector system, a nearly 1000-multiplex fibre positioner, 20 individually deployable 'mini' integral field units (IFUs), and a single large IFU. These fibre systems feed a dual-beam spectrogr… ▽ More WEAVE, the new wide-field, massively multiplexed spectroscopic survey facility for the William Herschel Telescope, will see first light in late 2022. WEAVE comprises a new 2-degree field-of-view prime-focus corrector system, a nearly 1000-multiplex fibre positioner, 20 individually deployable 'mini' integral field units (IFUs), and a single large IFU. These fibre systems feed a dual-beam spectrograph covering the wavelength range 366$-$959\,nm at $R\sim5000$, or two shorter ranges at $R\sim20\,000$. After summarising the design and implementation of WEAVE and its data systems, we present the organisation, science drivers and design of a five- to seven-year programme of eight individual surveys to: (i) study our Galaxy's origins by completing Gaia's phase-space information, providing metallicities to its limiting magnitude for $\sim$3 million stars and detailed abundances for $\sim1.5$ million brighter field and open-cluster stars; (ii) survey $\sim0.4$ million Galactic-plane OBA stars, young stellar objects and nearby gas to understand the evolution of young stars and their environments; (iii) perform an extensive spectral survey of white dwarfs; (iv) survey $\sim400$ neutral-hydrogen-selected galaxies with the IFUs; (v) study properties and kinematics of stellar populations and ionised gas in $z<0.5$ cluster galaxies; (vi) survey stellar populations and kinematics in $\sim25\,000$ field galaxies at $0.3\lesssim z \lesssim 0.7$; (vii) study the cosmic evolution of accretion and star formation using $>1$ million spectra of LOFAR-selected radio sources; (viii) trace structures using intergalactic/circumgalactic gas at $z>2$. Finally, we describe the WEAVE Operational Rehearsals using the WEAVE Simulator. △ Less

Submitted 31 October, 2023; v1 submitted 7 December, 2022; originally announced December 2022.

Comments: 41 pages, 27 figures, accepted for publication by MNRAS; updated version including information on individual grants in a revised Acknowledgements section, corrections to the affiliation list, and an updated references list

MaNGIA: 10,000 mock galaxies for stellar population analysis

Authors: Regina Sarmiento, Marc Huertas-Company, Johan H. Knapen, Héctor Ibarra-Medel, Annalisa Pillepich, Sebastián F. Sánchez, Alina Boecker

Abstract: Modern astronomical observations give unprecedented access to the physical properties of nearby galaxies, including spatially resolved stellar populations. However, observations can only give a present-day view of the Universe, whereas cosmological simulations give access to the past record of the processes that galaxies have experienced in their evolution. To connect the events that happened in t… ▽ More Modern astronomical observations give unprecedented access to the physical properties of nearby galaxies, including spatially resolved stellar populations. However, observations can only give a present-day view of the Universe, whereas cosmological simulations give access to the past record of the processes that galaxies have experienced in their evolution. To connect the events that happened in the past with galactic properties as seen today, simulations must be taken to a common ground before being compared to observations. We emulate data from the MaNGA survey, which is the largest integral field spectroscopic galaxy survey to date with its 10,000 nearby galaxies of all types. For this, we use the cosmological simulations TNG50 to generate MaNGIA (Mapping Nearby Galaxies with IllustrisTNG Astrophysics), a mock MaNGA sample of similar size that emulates observations of galaxies for stellar population analysis. We choose TNG galaxies to match the MaNGA sample selection to limit the impact of selection effects. We produce MaNGA-like datacubes from all simulated galaxies, and process these with the pyPipe3D analysis code. This allows us to extract spatially resolved stellar maps. This first paper presents the approach to generate the mock sample and provides an initial exploration of its properties. We show that the stellar populations and kinematics of the simulated MaNGIA galaxies are overall in good agreement with observations. Specific discrepancies, especially in the age and metallicity gradients in low- to intermediate-mass regimes and in massive galaxies' kinematics, require further investigation. We compare our results to other attempts to mock similar observations, all of smaller data sets. Our final dataset will be released with the publication, consisting of >10,000 post-processed data-cubes analysed with pyPipe3D, along with the codes developed to create it. △ Less

Submitted 21 November, 2022; originally announced November 2022.

Comments: Submitted to A&A, 21 pages, 17 figures, 1 table

The edges of galaxies: Tracing the limits of star formation

Authors: N. Chamba, I. Trujillo, J. H. Knapen

Abstract: The outskirts of galaxies have been studied from multiple perspectives for the past few decades. However, it is still unknown if all galaxies have clear-cut edges like everyday objects. We address this question by developing physically motivated criteria to define the edges of galaxies. Based on the gas density threshold required for star formation, we define the edge of a galaxy as the outermost… ▽ More The outskirts of galaxies have been studied from multiple perspectives for the past few decades. However, it is still unknown if all galaxies have clear-cut edges like everyday objects. We address this question by developing physically motivated criteria to define the edges of galaxies. Based on the gas density threshold required for star formation, we define the edge of a galaxy as the outermost radial location associated with a significant drop in either past or ongoing in-situ star formation. We explore $\sim$1000 low-inclination galaxies with a wide range in morphology (dwarfs to ellipticals) and stellar mass ($10^7 M_{\odot} < M_{\star} < 10^{12}M_{\odot}$). The location of the edges of these galaxies ($R_{\rm edge}$) are visually identified as the outermost cut-off or truncation in their radial profiles using deep multi-band optical imaging from the IAC Stripe82 Legacy Project. We find this characteristic feature at the following mean stellar mass density which varies with galaxy morphology: $2.9\pm0.10\,M_{\odot}$/pc$^2$ for ellipticals, $1.1\pm0.04\,M_{\odot}/$pc$^2$ for spirals and $0.6\pm0.03\,M_{\odot}/$pc$^2$ for present-day star forming dwarfs. Additionally, we find that $R_{\rm edge}$ depends on its age (colour) where bluer galaxies have larger $R_{\rm edge}$ at a fixed stellar mass. The resulting stellar mass--size plane using $R_{\rm edge}$ as a physically motivated galaxy size measure has a very narrow intrinsic scatter ($\lesssim 0.06$ dex). These results highlight the importance of new deep imaging surveys to explore the growth of galaxies and trace the limits of star formation in their outskirts. △ Less

Submitted 12 September, 2022; originally announced September 2022.

Comments: 22 pages (including appendix), 11 Figures, accepted for publication in A&A. Definition and concept explained in Section 2. Criteria to identify edges for each morphological type is detailed in Section 5. Key results in Figs. 5-8

Journal ref: A&A 667, A87 (2022)

Cloud-scale Radio Surveys of Star Formation and Feedback in Triangulum Galaxy M 33: VLA Observations

Authors: F. S. Tabatabaei, W. Cotton, E. Schinnerer, R. Beck, A. Brunthaler, K. M. Menten, J. Braine, E. Corbelli, C. Kramer, J. E. Beckman, J. H. Knapen, R. Paladino, E. Koch, A. Camps Farina

Abstract: Studying the interplay between massive star formation and the interstellar medium (ISM) is paramount to understand the evolution of galaxies. Radio continuum (RC) emission serves as an extinction-free tracer of both massive star formation and the energetic components of the interstellar medium. We present a multi-band radio continuum survey of the local group galaxy M 33 down to ~30 pc linear reso… ▽ More Studying the interplay between massive star formation and the interstellar medium (ISM) is paramount to understand the evolution of galaxies. Radio continuum (RC) emission serves as an extinction-free tracer of both massive star formation and the energetic components of the interstellar medium. We present a multi-band radio continuum survey of the local group galaxy M 33 down to ~30 pc linear resolution observed with the Karl G. Jansky Very Large Array (VLA). We calibrate the star-formation rate surface density and investigate the impact of diffuse emission on this calibration using a structural decomposition. Separating the thermal and nonthermal emission components, the correlation between different phases of the interstellar medium and the impact of massive star formation are also investigated. Radio sources with sizes <~ 200 pc constitute about 36% (46%) of the total RC emission at 1.5 GHz (6.3 GHz) in the inner 18' x 18' (or 4kpc x 4kpc) disk of M 33. The nonthermal spectral index becomes flatter with increasing star-formation rate surface density, indicating the escape of cosmic ray electrons {from their birth places}. The magnetic field strength also increases with star-formation rate following a bi-modal relation, indicating that the small-scale turbulent dynamo acts more efficiently at higher luminosities and star-formation rates. Although the correlations are tighter in star-forming regions, the nonthermal emission is correlated also with the more quiescent molecular gas in the ISM. An almost linear molecular star-formation law exists in M 33 when excluding diffuse structures. Massive star formation amplifies the magnetic field and increases the number of high-energy cosmic ray electrons, which can help the onset of winds and outflows. △ Less

Submitted 3 September, 2022; originally announced September 2022.

Comments: 20 pages, 12 figures, accepted for publication in MNRAS journal

Linking star formation thresholds and truncations in the thin and thick disks of the low-mass galaxy UGC 7321

Authors: Simón Díaz-García, Sébastien Comerón, Stéphane Courteau, Aaron E. Watkins, Johan H. Knapen, Javier Román

Abstract: Thin and thick disks are found in most spiral galaxies, yet their formation scenarios remain uncertain. Whether thick disks form through slow or fast, internal or environmental, processes is unclear. The physical origin of outer truncations in thin and thick disks, observed as a drop in optical and near-infrared (NIR) surface brightness profiles, is also a much debated topic. These truncations hav… ▽ More Thin and thick disks are found in most spiral galaxies, yet their formation scenarios remain uncertain. Whether thick disks form through slow or fast, internal or environmental, processes is unclear. The physical origin of outer truncations in thin and thick disks, observed as a drop in optical and near-infrared (NIR) surface brightness profiles, is also a much debated topic. These truncations have been linked to star formation (SF) thresholds in Milky-Way type galaxies, but no such connection has been made for their low-mass counterparts or in thick disks. Our photometric analysis of the edge-on galaxy UGC 7321 offers a possible breakthrough. This well-studied diffuse, isolated, bulgeless, ultra-thin galaxy is thought to be under-evolved both dynamically and in SF. It is an ideal target to disentangle internal effects in the formation of thick disks and truncations. Our axial light profiles from deep far- and near-ultraviolet (UV; GALEX) images, tracing recent SF, and optical (DESI grz) and NIR (Spitzer 3.6 microns) images, tracing old stellar populations, enable a detailed identification of an outer truncation in all probed wavelengths in both the thin and thick disks. After deprojecting to a face-on view, a sharp truncation signature is found at a stellar density of roughly 1.5 solar masses per square parsec, in agreement with theoretical expectations of gas density SF thresholds. The redder colours beyond the truncation radius are indicative of stellar migration towards the outer regions. We thus show that thick disks and truncations can form via internal mechanisms alone, given the pristine nature of UGC 7321. We report the discovery of a truncation at and above the mid-plane of a diffuse galaxy that is linked to a SF threshold; this poses a constraint on physically-motivated disk size measurements among low-mass galaxies. △ Less

Submitted 4 November, 2022; v1 submitted 29 August, 2022; originally announced August 2022.

Comments: Accepted for publication in A&A (August 29, 2022). 7 pages, 3 figures

Journal ref: A&A 667, A109 (2022)

How to plan your astronomy research paper in ten steps

Authors: Nushkia Chamba, Johan H. Knapen, Diane Black

Abstract: Scientific writing is an important skill for a career as a professional astrophysicist. Very few researchers, however, receive any formal training in how to write scientific research papers of high quality in an efficient manner. This paper (Paper I) is the first of a two-part self-help guide in scientific writing to address this skills gap. Paper I focuses on planning your academic research paper… ▽ More Scientific writing is an important skill for a career as a professional astrophysicist. Very few researchers, however, receive any formal training in how to write scientific research papers of high quality in an efficient manner. This paper (Paper I) is the first of a two-part self-help guide in scientific writing to address this skills gap. Paper I focuses on planning your academic research paper in astronomy. We discuss how to crystallise the ideas that underlie a research project, analyse how the paper can be constructed considering the audience and the chosen journal, and give an overview of the publishing process. Paper II is a detailed description of the different sections that make up a research paper in astronomy and shares the best practice in how to write in English. Whether you are a student writing your first paper or an experienced author, you may find the ideas presented here useful. △ Less

Submitted 12 October, 2022; v1 submitted 22 July, 2022; originally announced July 2022.

Comments: 15 pages; Paper I, submitted as a two-part series of Perspectives to Nature Astronomy. Kindly cite and/or acknowledge this work and Paper II (arXiv:2110.05503: 'How to write and develop your astronomy research paper' by Knapen, Chamba & Black) if you use them to plan and write your paper. Comments are very welcome! The Version of Record is available online at https://rdcu.be/cVdtK

Preparing for low surface brightness science with the Vera C. Rubin Observatory: characterisation of tidal features from mock images

Authors: G. Martin, A. E. Bazkiaei, M. Spavone, E. Iodice, J. C. Mihos, M. Montes, J. A. Benavides, S. Brough, J. L. Carlin, C. A. Collins, P. A. Duc, F. A. Gómez, G. Galaz, H. M. Hernández-Toledo, R. A. Jackson, S. Kaviraj, J. H. Knapen, C. Martínez-Lombilla, S. McGee, D. O'Ryan, D. J. Prole, R. M. Rich, J. Román, E. A. Shah, T. K. Starkenburg , et al. (28 additional authors not shown)

Abstract: Tidal features in the outskirts of galaxies yield unique information about their past interactions and are a key prediction of the hierarchical structure formation paradigm. The Vera C. Rubin Observatory is poised to deliver deep observations for potentially of millions of objects with visible tidal features, but the inference of galaxy interaction histories from such features is not straightforwa… ▽ More Tidal features in the outskirts of galaxies yield unique information about their past interactions and are a key prediction of the hierarchical structure formation paradigm. The Vera C. Rubin Observatory is poised to deliver deep observations for potentially of millions of objects with visible tidal features, but the inference of galaxy interaction histories from such features is not straightforward. Utilising automated techniques and human visual classification in conjunction with realistic mock images produced using the NEWHORIZON cosmological simulation, we investigate the nature, frequency and visibility of tidal features and debris across a range of environments and stellar masses. In our simulated sample, around 80 per cent of the flux in the tidal features around Milky Way or greater mass galaxies is detected at the 10-year depth of the Legacy Survey of Space and Time (30-31 mag / sq. arcsec), falling to 60 per cent assuming a shallower final depth of 29.5 mag / sq. arcsec. The fraction of total flux found in tidal features increases towards higher masses, rising to 10 per cent for the most massive objects in our sample (M*~10^{11.5} Msun). When observed at sufficient depth, such objects frequently exhibit many distinct tidal features with complex shapes. The interpretation and characterisation of such features varies significantly with image depth and object orientation, introducing significant biases in their classification. Assuming the data reduction pipeline is properly optimised, we expect the Rubin Observatory to be capable of recovering much of the flux found in the outskirts of Milky Way mass galaxies, even at intermediate redshifts (z<0.2). △ Less

Submitted 7 May, 2022; v1 submitted 15 March, 2022; originally announced March 2022.

Comments: 29 pages, 25 figures; accepted for publication in MNRAS following minor corrections

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 513, Issue 1, June 2022, Pages 1459-1487,

Implications for Galaxy Formation Models from Observations of Globular Clusters around Ultra-Diffuse Galaxies

Authors: Teymoor Saifollahi, Dennis Zaritsky, Ignacio Trujillo, Reynier F. Peletier, Johan H. Knapen, Nicola Amorisco, Michael A. Beasley, Richard Donnerstein

Abstract: We present an analysis of Hubble Space Telescope observations of globular clusters (GCs) in six ultra-diffuse galaxies (UDGs) in the Coma cluster, a sample that represents UDGs with large effective radii ($R_{\rm e}$), and use the results to evaluate competing formation models. We eliminate two significant sources of systematic uncertainty in the determination of the number of GCs, $N_{\rm GC}$ by… ▽ More We present an analysis of Hubble Space Telescope observations of globular clusters (GCs) in six ultra-diffuse galaxies (UDGs) in the Coma cluster, a sample that represents UDGs with large effective radii ($R_{\rm e}$), and use the results to evaluate competing formation models. We eliminate two significant sources of systematic uncertainty in the determination of the number of GCs, $N_{\rm GC}$ by using sufficiently deep observations that (i) reach the turnover of the GC luminosity function and (ii) provide a sufficient number of GCs with which to measure the GC number radial distribution. We find that $N_{\rm GC}$ for these galaxies is on average $\sim$20, which implies an average total mass, $M_{\rm total}$, $\sim$ $10^{11}$ $M_{\odot}$ when applying the relation between $N_{\rm GC}$ and $M_{\rm total}$. This value of $N_{\rm GC}$ lies at the upper end of the range observed for dwarf galaxies of the same stellar mass and is roughly a factor of two larger than the mean. The GC luminosity function, radial profile and average colour are more consistent with those observed for dwarf galaxies than with those observed for the more massive ($L^*$) galaxies, while both the radial and azimuthal GC distributions closely follow those of the stars in the host galaxy. Finally, we discuss why our observations, specifically the GC number and GC distribution around these six UDGs, pose challenges for several of the currently favoured UDG formation models. △ Less

Submitted 3 February, 2022; v1 submitted 27 January, 2022; originally announced January 2022.

Comments: Accepted for publication in MNRAS

Stellar masses, sizes, and radial profiles for 465 nearby early-type galaxies: an extension to the Spitzer Survey of Stellar Structure in Galaxies (S$^{4}$G)

Authors: A. E. Watkins, H. Salo, E. Laurikainen, S. Díaz-García, S. Comerón, J. Janz, A. H. Su, R. Buta, E. Athanassoula, A. Bosma, L. C. Ho, B. W. Holwerda, T. Kim, J. H. Knapen, S. Laine, K. Menéndez-Delmestre, R. F. Peletier, K. Sheth, D. Zaritsky

Abstract: The Spitzer Survey of Stellar Structure in Galaxies (S$^{4}$G) is a detailed study of over 2300 nearby galaxies in the near-infrared (NIR), which has been critical to our understanding of the detailed structures of nearby galaxies. Because the sample galaxies were selected only using radio-derived velocities, however, the survey favored late-type disk galaxies over lenticulars and ellipticals. A f… ▽ More The Spitzer Survey of Stellar Structure in Galaxies (S$^{4}$G) is a detailed study of over 2300 nearby galaxies in the near-infrared (NIR), which has been critical to our understanding of the detailed structures of nearby galaxies. Because the sample galaxies were selected only using radio-derived velocities, however, the survey favored late-type disk galaxies over lenticulars and ellipticals. A follow-up Spitzer survey was conducted to rectify this bias, adding 465 early-type galaxies (ETGs) to the original sample, to be analyzed in a manner consistent with the initial survey. We present the data release of this ETG extension, up to the third data processing pipeline (P3): surface photometry. We produce curves of growth and radial surface brightness profiles (with and without inclination corrections) using reduced and masked Spitzer IRAC 3.6$μ$m and 4.5$μ$m images produced through Pipelines 1 and 2, respectively. From these profiles, we derive the following integrated quantities: total magnitudes, stellar masses, concentration parameters, and galaxy size metrics. We showcase NIR scaling relations for ETGs among these quantities. We examine general trends across the whole S$^{4}$G and ETG extension among our derived parameters, highlighting differences between ETGs and late-type galaxies (LTGs). ETGs are, on average, more massive and more concentrated than LTGs, and also show subtle distinctions among ETG morphological sub-types. We also derive the following scaling relations and compare with previous results in visible light: mass--size (both half-light and isophotal), mass--concentration, mass--surface brightness (central, effective, and within 1 kpc), and mass--color. We find good agreement with previous works, though some relations (e.g., mass--central surface brightness) will require more careful multi-component decompositions to be fully understood. △ Less

Submitted 20 January, 2022; originally announced January 2022.

Comments: 25 pages, 17 figures, accepted for publication in A&A

Journal ref: A&A 660, A69 (2022)

The Seventeenth Data Release of the Sloan Digital Sky Surveys: Complete Release of MaNGA, MaStar and APOGEE-2 Data

Authors: Abdurro'uf, Katherine Accetta, Conny Aerts, Victor Silva Aguirre, Romina Ahumada, Nikhil Ajgaonkar, N. Filiz Ak, Shadab Alam, Carlos Allende Prieto, Andres Almeida, Friedrich Anders, Scott F. Anderson, Brett H. Andrews, Borja Anguiano, Erik Aquino-Ortiz, Alfonso Aragon-Salamanca, Maria Argudo-Fernandez, Metin Ata, Marie Aubert, Vladimir Avila-Reese, Carles Badenes, Rodolfo H. Barba, Kat Barger, Jorge K. Barrera-Ballesteros, Rachael L. Beaton , et al. (316 additional authors not shown)

Abstract: This paper documents the seventeenth data release (DR17) from the Sloan Digital Sky Surveys; the fifth and final release from the fourth phase (SDSS-IV). DR17 contains the complete release of the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, which reached its goal of surveying over 10,000 nearby galaxies. The complete release of the MaNGA Stellar Library (MaStar) accompanies… ▽ More This paper documents the seventeenth data release (DR17) from the Sloan Digital Sky Surveys; the fifth and final release from the fourth phase (SDSS-IV). DR17 contains the complete release of the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, which reached its goal of surveying over 10,000 nearby galaxies. The complete release of the MaNGA Stellar Library (MaStar) accompanies this data, providing observations of almost 30,000 stars through the MaNGA instrument during bright time. DR17 also contains the complete release of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) survey which publicly releases infra-red spectra of over 650,000 stars. The main sample from the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), as well as the sub-survey Time Domain Spectroscopic Survey (TDSS) data were fully released in DR16. New single-fiber optical spectroscopy released in DR17 is from the SPectroscipic IDentification of ERosita Survey (SPIDERS) sub-survey and the eBOSS-RM program. Along with the primary data sets, DR17 includes 25 new or updated Value Added Catalogs (VACs). This paper concludes the release of SDSS-IV survey data. SDSS continues into its fifth phase with observations already underway for the Milky Way Mapper (MWM), Local Volume Mapper (LVM) and Black Hole Mapper (BHM) surveys. △ Less

Submitted 13 January, 2022; v1 submitted 3 December, 2021; originally announced December 2021.

Comments: 40 pages, 8 figures, 6 tables. In press at ApJSS (arxiv v2 corrects some minor typos and updates references)

LeMMINGs IV: The X-ray properties of a statistically-complete sample of the nuclei in active and inactive galaxies from the Palomar sample

Authors: D. R. A. Williams, M. Pahari, R. D. Baldi, I. M. McHardy, S. Mathur, R. J. Beswick, A. Beri, P. Boorman, S. Aalto, A. Alberdi, M. K. Argo, B. T. Dullo, D. M. Fenech, D. A. Green, J. H. Knapen, I. Martí-Vidal, J. Moldon, C. G. Mundell, T. W. B. Muxlow, F. Panessa, M. Pérez-Torres, P. Saikia, F. Shankar, I. R. Stevens, P. Uttley

Abstract: All 280 of the statistically-complete Palomar sample of nearby (<120 Mpc) galaxies dec > 20 degrees have been observed at 1.5 GHz as part of the LeMMINGs e-MERLIN legacy survey. Here, we present Chandra X-ray observations of the nuclei of 213 of these galaxies, including a statistically-complete sub-set of 113 galaxies in the declination range 40 degrees to 65 degrees. We observed galaxies of all… ▽ More All 280 of the statistically-complete Palomar sample of nearby (<120 Mpc) galaxies dec > 20 degrees have been observed at 1.5 GHz as part of the LeMMINGs e-MERLIN legacy survey. Here, we present Chandra X-ray observations of the nuclei of 213 of these galaxies, including a statistically-complete sub-set of 113 galaxies in the declination range 40 degrees to 65 degrees. We observed galaxies of all optical spectral types, including 'active' galaxies (e.g., LINERs and Seyferts) and 'inactive' galaxies like HII galaxies and absorption line galaxies (ALG). The X-ray flux limit of our survey is 1.65$\times$10$^{-14}$~erg s$^{-1}$ cm$^{-2}$ (0.3$-$10 keV). We detect X-ray emission coincident within 2-arcsec of the nucleus in 150/213 galaxies, including 13/14 Seyferts, 68/77 LINERs, 13/22 ALGs and 56/100 HII galaxies, but cannot completely rule out contamination from non-AGN processes in sources with nuclear luminosities <10$^{39}$ erg s$^{-1}$. We construct an X-ray Luminosity function (XLF) and find that the local galaxy XLF, when including all AGN types, can be represented as a single power-law of slope $-0.54 \pm 0.06$. The Eddington ratio of the Seyferts is usually 2-4 decades higher than that of the LINERs, ALGs and HII galaxies, which are mostly detected with Eddington ratios <10$^{-3}$. Using [O III] line measurements and BH masses from the literature, we show that LINERs, HII galaxies and ALGs follow similar correlations to low luminosities, suggesting that some 'inactive' galaxies may harbour AGN. △ Less

Submitted 17 November, 2021; originally announced November 2021.

Comments: 20 pages + 42 pages of online supplementary material 15 figures, 7 tables (additional 150 figures and 6 tables in online supplementary material) Accepted for publication in MNRAS, 15th November 2021

How to write and develop your astronomy research paper

Authors: Johan H. Knapen, Nushkia Chamba, Diane Black

Abstract: Writing is a vital component of a modern career in scientific research. But how to write correctly and effectively is often not included in the training that young astronomers receive from their supervisors and departments. We offer a step-by-step guide to tackle this deficiency, published as a set of two papers. In the first, we addressed how to plan and outline your paper and decide where to pub… ▽ More Writing is a vital component of a modern career in scientific research. But how to write correctly and effectively is often not included in the training that young astronomers receive from their supervisors and departments. We offer a step-by-step guide to tackle this deficiency, published as a set of two papers. In the first, we addressed how to plan and outline your paper and decide where to publish. In the current second paper, we describe the various sections that constitute a typical research paper in astronomy, sharing best practice for the most efficient use of each of them. We also discuss a selection of issues that often cause trouble to writers, from sentence to paragraph structure, the `writing mechanics' used to develop a manuscript. Our two-part guide is aimed primarily at master's and PhD level students who are presented with the daunting task of writing their first scientific paper, but more senior researchers or writing instructors may well find the ideas presented here useful. △ Less

Submitted 10 January, 2023; v1 submitted 11 October, 2021; originally announced October 2021.

Comments: Submitted to Nature Astronomy as Paper 2 in a two-part mini-series of 'Perspectives' articles. Paper 1 in the series is 'How to plan your astronomy research paper in ten steps' by Chamba, Knapen & Black. Kindly cite and/or acknowledge this work and Paper 1 (arXiv:2207.12959) if you use them to plan and write your paper. The Version of Record is available online at https://rdcu.be/cVdtJ

LeMMINGs. III. The e-MERLIN Legacy Survey of the Palomar sample. Exploring the origin of nuclear radio emission in active and inactive galaxies through the [O III] -- radio connection

Authors: R. D. Baldi, D. R. A. Williams, R. J. Beswick, I. McHardy, B. T. Dullo, J. H. Knapen, L. Zanisi, M. K. Argo, S. Aalto, A. Alberdi, W. A. Baan, G. J. Bendo, D. M. Fenech, D. A. Green, H. -R. Klöckner, E. Körding, T. J. Maccarone, J. M. Marcaide, I. Mutie, F. Panessa, M. A. Pérez-Torres, C. Romero-Cañizales, D. J. Saikia, P. Saikia, F. Shankar , et al. (9 additional authors not shown)

Abstract: What determines the nuclear radio emission in local galaxies? We combine optical [O III] line emission, robust black hole (BH) mass estimates, and high-resolution e-MERLIN 1.5-GHz data, from the LeMMINGs survey, of a statistically-complete sample of 280 nearby, optically active (LINER and Seyfert) and inactive HII and Absorption line galaxies [ALG]) galaxies. Using [O III] luminosity (… ▽ More What determines the nuclear radio emission in local galaxies? We combine optical [O III] line emission, robust black hole (BH) mass estimates, and high-resolution e-MERLIN 1.5-GHz data, from the LeMMINGs survey, of a statistically-complete sample of 280 nearby, optically active (LINER and Seyfert) and inactive HII and Absorption line galaxies [ALG]) galaxies. Using [O III] luminosity ($L_{\rm [O~III]}$) as a proxy for the accretion power, local galaxies follow distinct sequences in the optical-radio planes of BH activity, which suggest different origins of the nuclear radio emission for the optical classes. The 1.5-GHz radio luminosity of their parsec-scale cores ($L_{\rm core}$) is found to scale with BH mass ($M_{\rm BH}$) and [O~III] luminosity. Below $M_{\rm BH} \sim$10$^{6.5}$ M$_{\odot}$, stellar processes from non-jetted HII galaxies dominate with $L_{\rm core} \propto M_{\rm BH}^{0.61\pm0.33}$ and $L_{\rm core} \propto L_{\rm [O~III]}^{0.79\pm0.30}$. Above $M_{\rm BH} \sim$10$^{6.5}$ M$_{\odot}$, accretion-driven processes dominate with $L_{\rm core} \propto M_{\rm BH}^{1.5-1.65}$ and $L_{\rm core} \propto L_{\rm [O~III]}^{0.99-1.31}$ for active galaxies: radio-quiet/loud LINERs, Seyferts and jetted HII galaxies always display (although low) signatures of radio-emitting BH activity, with $L_{\rm 1.5\, GHz}\gtrsim$10$^{19.8}$ W Hz$^{-1}$ and $M_{\rm BH}\gtrsim10^{7}$ M$_{\odot}$, on a broad range of Eddington-scaled accretion rates ($\dot{m}$). Radio-quiet and radio-loud LINERs are powered by low-$\dot{m}$ discs launching sub-relativistic and relativistic jets, respectively. Low-power slow jets and disc/corona winds from moderately high to high-$\dot{m}$ discs account for the compact and edge-brightened jets of Seyferts, respectively. Jetted HII galaxies may host weakly active BHs. Fuel-starved BHs and recurrent activity account for ALG properties. [abridged] △ Less

Submitted 13 September, 2021; originally announced September 2021.

Comments: Accpeted for publication on MNRAS (27 pages, 9 Figures, 3 Tables, Appendix A)

Euclid preparation: XVI. Exploring the ultra low-surface brightness Universe with Euclid/VIS

Authors: A. S. Borlaff, P. Gómez-Alvarez, B. Altieri, P. M. Marcum, R. Vavrek, R. Laureijs, R. Kohley, F. Buitrago, J. C. Cuillandre, P. A. Duc, L. M. Gaspar Venancio, A. Amara, S. Andreon, N. Auricchio, R. Azzollini, C. Baccigalupi, A. Balaguera-Antolínez, M. Baldi, S. Bardelli, R. Bender, A. Biviano, C. Bodendorf, D. Bonino, E. Bozzo, E. Branchini , et al. (158 additional authors not shown)

Abstract: While Euclid is an ESA mission specifically designed to investigate the nature of Dark Energy and Dark Matter, the planned unprecedented combination of survey area ($\sim15\,000$ deg$^2$), spatial resolution, low sky-background, and depth also make Euclid an excellent space observatory for the study of the low surface brightness Universe. Scientific exploitation of the extended low surface brightn… ▽ More While Euclid is an ESA mission specifically designed to investigate the nature of Dark Energy and Dark Matter, the planned unprecedented combination of survey area ($\sim15\,000$ deg$^2$), spatial resolution, low sky-background, and depth also make Euclid an excellent space observatory for the study of the low surface brightness Universe. Scientific exploitation of the extended low surface brightness structures requires dedicated calibration procedures yet to be tested. We investigate the capabilities of Euclid to detect extended low surface brightness structure by identifying and quantifying sky background sources and stray-light contamination. We test the feasibility of generating sky flat-fields to reduce large-scale residual gradients in order to reveal the extended emission of galaxies observed in the Euclid Survey. We simulate a realistic set of Euclid/VIS observations, taking into account both instrumental and astronomical sources of contamination, including cosmic rays, stray-light, zodiacal light, ISM, and the CIB, while simulating the effects of the presence of background sources in the FOV. We demonstrate that a combination of calibration lamps, sky flats and self-calibration would enable recovery of emission at a limiting surface brightness magnitude of $μ=29.5^{+0.08}_{-0.27} $ mag arcsec$^{-2}$ ($3σ$, $10\times10$ arcsec$^2$) in the Wide Survey, reaching regions 2 magnitudes deeper in the Deep Surveys. Euclid/VIS has the potential to be an excellent low surface brightness observatory. Covering the gap between pixel-to-pixel calibration lamp flats and self-calibration observations for large scales, the application of sky flat-fielding will enhance the sensitivity of the VIS detector at scales of larger than 1 degree, up to the size of the FOV, enabling Euclid to detect extended surface brightness structures below $μ=31$ mag arcsec$^{-2}$ and beyond. △ Less

Submitted 23 August, 2021; originally announced August 2021.

Comments: 23 pages, 13 figures, Euclid Consortium Key Project, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 657, A92 (2022)

Euclid preparation: XIII. Forecasts for galaxy morphology with the Euclid Survey using Deep Generative Models

Authors: Euclid Collaboration, H. Bretonnière, M. Huertas-Company, A. Boucaud, F. Lanusse, E. Jullo, E. Merlin, D. Tuccillo, M. Castellano, J. Brinchmann, C. J. Conselice, H. Dole, R. Cabanac, H. M. Courtois, F. J. Castander, P. A. Duc, P. Fosalba, D. Guinet, S. Kruk, U. Kuchner, S. Serrano, E. Soubrie, A. Tramacere, L. Wang, A. Amara , et al. (171 additional authors not shown)

Abstract: We present a machine learning framework to simulate realistic galaxies for the Euclid Survey. The proposed method combines a control on galaxy shape parameters offered by analytic models with realistic surface brightness distributions learned from real Hubble Space Telescope observations by deep generative models. We simulate a galaxy field of $0.4\,\rm{deg}^2$ as it will be seen by the Euclid vis… ▽ More We present a machine learning framework to simulate realistic galaxies for the Euclid Survey. The proposed method combines a control on galaxy shape parameters offered by analytic models with realistic surface brightness distributions learned from real Hubble Space Telescope observations by deep generative models. We simulate a galaxy field of $0.4\,\rm{deg}^2$ as it will be seen by the Euclid visible imager VIS and show that galaxy structural parameters are recovered with similar accuracy as for pure analytic Sérsic profiles. Based on these simulations, we estimate that the Euclid Wide Survey will be able to resolve the internal morphological structure of galaxies down to a surface brightness of $22.5\,\rm{mag}\,\rm{arcsec}^{-2}$, and $24.9\,\rm{mag}\,\rm{arcsec}^{-2}$ for the Euclid Deep Survey. This corresponds to approximately $250$ million galaxies at the end of the mission and a $50\,\%$ complete sample for stellar masses above $10^{10.6}\,\rm{M}_\odot$ (resp. $10^{9.6}\,\rm{M}_\odot$) at a redshift $z\sim0.5$ for the wide (resp. deep) survey. The approach presented in this work can contribute to improving the preparation of future high-precision cosmological imaging surveys by allowing simulations to incorporate more realistic galaxies. △ Less

Submitted 10 January, 2022; v1 submitted 25 May, 2021; originally announced May 2021.

Comments: 22 pages, 21 figures (6 in appendinx). Accepted for publication to A&A. Change in the title for Euclid submission chronology coherence. Few changes for clarification

The centres of M83 and the Milky Way: opposite extremes of a common star formation cycle

Authors: Daniel Callanan, Steven N. Longmore, J. M. Diederik Kruijssen, Andreas Schruba, Adam Ginsburg, Mark R. Krumholz, Nate Bastian, Joao Alves, Jonathan D. Henshaw, Johan H. Knapen, Melanie Chevance

Abstract: In the centres of the Milky Way and M83, the global environmental properties thought to control star formation are very similar. However, M83's nuclear star formation rate (SFR), as estimated by synchrotron and H-alpha emission, is an order of magnitude higher than the Milky Way's. To understand the origin of this difference we use ALMA observations of HCN (1-0) and HCO+ (1-0) to trace the dense g… ▽ More In the centres of the Milky Way and M83, the global environmental properties thought to control star formation are very similar. However, M83's nuclear star formation rate (SFR), as estimated by synchrotron and H-alpha emission, is an order of magnitude higher than the Milky Way's. To understand the origin of this difference we use ALMA observations of HCN (1-0) and HCO+ (1-0) to trace the dense gas at the size scale of individual molecular clouds (0.54", 12pc) in the inner ~500 pc of M83, and compare this to gas clouds at similar resolution and galactocentric radius in the Milky Way. We find that both the overall gas distribution and the properties of individual clouds are very similar in the two galaxies, and that a common mechanism may be responsible for instigating star formation in both circumnuclear rings. Given the considerable similarity in gas properties, the most likely explanation for the order of magnitude difference in SFR is time variability, with the Central Molecular Zone (CMZ) currently being at a more quiescent phase of its star formation cycle. We show M83's SFR must have been an order of magnitude higher 5-7 Myr ago. M83's `starburst' phase was highly localised, both spatially and temporally, greatly increasing the feedback efficiency and ability to drive galactic-scale outflows. This highly dynamic nature of star formation and feedback cycles in galaxy centres means (i) modeling and interpreting observations must avoid averaging over large spatial areas or timescales, and (ii) understanding the multi-scale processes controlling these cycles requires comparing snapshots of a statistical sample of galaxies in different evolutionary stages. △ Less

Submitted 20 May, 2021; originally announced May 2021.

Comments: 31 pages, 22 figures, accepted by MNRAS

Capturing the physics of MaNGA galaxies with self-supervised Machine Learning

Authors: Regina Sarmiento, Marc Huertas-Company, Johan H. Knapen, Sebastián F. Sánchez, Helena Domínguez Sánchez, Niv Drory, Jesús Falcón-Barroso

Abstract: As available data sets grow in size and complexity, advanced visualization tools enabling their exploration and analysis become more important. In modern astronomy, integral field spectroscopic galaxy surveys are a clear example of increasing dimensionality and complexity of datasets, which challenge the traditional methods used to extract the physical information they contain. We present the use… ▽ More As available data sets grow in size and complexity, advanced visualization tools enabling their exploration and analysis become more important. In modern astronomy, integral field spectroscopic galaxy surveys are a clear example of increasing dimensionality and complexity of datasets, which challenge the traditional methods used to extract the physical information they contain. We present the use of a novel self-supervised Machine Learning method to visualize the multi-dimensional information on stellar population and kinematics in the MaNGA survey in a two dimensional plane. Our framework is insensitive to non-physical properties such as the size of integral field unit (IFU) and is therefore able to order galaxies according to their resolved physical properties. Using the extracted representations, we study how galaxies distribute based on their resolved and global physical properties. We show that even when using exclusively information about the internal structure, galaxies naturally cluster into two well-known categories from a purely data driven perspective: rotating main-sequence disks and massive slow rotators, hence confirming distinct assembly channels. Low-mass rotation-dominated quenched galaxies appear only as a third cluster if information about the integrated physical properties is preserved, suggesting a mixture of assembly processes for these galaxies without any particular signature in their internal kinematics that distinguishes them from the two main groups. The framework for data exploration is publicly released with this publication, ready to be used with the MaNGA or other integral field data sets. △ Less

Submitted 19 August, 2021; v1 submitted 16 April, 2021; originally announced April 2021.

Comments: 17 pages, 9 figures. Accepted for publication in The Astrophysical Journal

Ultramassive black holes in the most massive galaxies: $M_{\rm BH}-σ$ versus $M_{\rm BH}-R_{\rm b}$

Authors: Bililign T. Dullo, Armando Gil de Paz, Johan H. Knapen

Abstract: [Abridged] We investigate the nature of the relations between black hole (BH) mass ($M_{\rm BH}$) and the central velocity dispersion ($σ$) and, for core-Sérsic galaxies, the size of the depleted core ($R_{\rm b}$). Our sample of 144 galaxies with dynamically determined $M_{\rm BH}$ encompasses 24 core-Sérsic galaxies, thought to be products of gas-poor mergers, and reliably identified based on hi… ▽ More [Abridged] We investigate the nature of the relations between black hole (BH) mass ($M_{\rm BH}$) and the central velocity dispersion ($σ$) and, for core-Sérsic galaxies, the size of the depleted core ($R_{\rm b}$). Our sample of 144 galaxies with dynamically determined $M_{\rm BH}$ encompasses 24 core-Sérsic galaxies, thought to be products of gas-poor mergers, and reliably identified based on high-resolution HST imaging. For core-Sérsic galaxies -- i.e., combining normal-core ($R_{\rm b} < 0.5 $ kpc) and large-core galaxies ($R_{\rm b} \gtrsim 0.5$ kpc), we find that $M_{\rm BH}$ correlates remarkably well with $R_{\rm b}$ such that $M_{\rm BH} \propto R_{\rm b}^{1.20 \pm 0.14}$ (rms scatter in log $M_{\rm BH}$ of $Δ_{\rm rms} \sim 0.29$ dex), confirming previous works on the same galaxies except three new ones. Separating the sample into Sérsic, normal-core and large-core galaxies, we find that Sérsic and normal-core galaxies jointly define a single log-linear $M_{\rm BH}-σ$ relation $M_{\rm BH} \propto σ^{ 4.88 \pm 0.29}$ with $Δ_{\rm rms} \sim 0.47$ dex, however, at the high-mass end large-core galaxies (four with measured $M_{\rm BH}$) are offset upward from this relation by ($2.5-4) \times σ_{\rm s}$, explaining the previously reported steepening of the $M_{\rm BH}-σ$ relation for massive galaxies. Large-core spheroids have magnitudes $M_{V} \le -23.50$ mag, half-light radii Re $>$ 10 kpc and are extremely massive $M_{*} \ge 10^{12}M_{\odot}$. Furthermore, these spheroids tend to host ultramassive BHs ($M_{\rm BH} \ge 10^{10}M_{\odot}$) tightly connected with their $R_{\rm b}$ rather than $σ$. The less popular $M_{\rm BH}-R_{\rm b}$ relation exhibits $\sim$ 62% less scatter in log $M_{\rm BH}$ than the $M_{\rm BH}- σ$ relations. △ Less

Submitted 8 December, 2020; originally announced December 2020.

Comments: 11 pages, 4 figures, accepted (2020, November 27) for publication in ApJ

LeMMINGs. II. The e-MERLIN legacy survey of nearby galaxies. The deepest radio view of the Palomar sample on parsec scale

Authors: R. D. Baldi, D. R. A. Williams, I. M. McHardy, R. J. Beswick, E. Brinks, B. T. Dullo, J. H. Knapen, M. K. Argo, S. Aalto, A. Alberdi, W. A. Baan, G. J. Bendo, S. Corbel, D. M. Fenech, J. S. Gallagher, D. A. Green, R. C. Kennicutt, H. -R. Klöckner, E. Körding, T. J. Maccarone, T. W. B. Muxlow, C. G. Mundell, F. Panessa, A. B. Peck, M. A. Pérez-Torres , et al. (8 additional authors not shown)

Abstract: We present the second data release of high-resolution ($\leq0.2$ arcsec) 1.5-GHz radio images of 177 nearby galaxies from the Palomar sample, observed with the e-MERLIN array, as part of the LeMMINGs (Legacy e-MERLIN Multi-band Imaging of Nearby Galaxy Sample) survey. Together with the 103 targets of the first LeMMINGs data release, this represents a complete sample of 280 local active (LINER and… ▽ More We present the second data release of high-resolution ($\leq0.2$ arcsec) 1.5-GHz radio images of 177 nearby galaxies from the Palomar sample, observed with the e-MERLIN array, as part of the LeMMINGs (Legacy e-MERLIN Multi-band Imaging of Nearby Galaxy Sample) survey. Together with the 103 targets of the first LeMMINGs data release, this represents a complete sample of 280 local active (LINER and Seyfert) and inactive galaxies HII galaxies and Absorption Line Galaxies, ALG). This large program is the deepest radio survey of the local Universe, $\gtrsim$10$^{17.6}$ W Hz$^{-1}$, regardless of the host and nuclear type: we detect radio emission $\gtrsim$0.25 mJy beam$^{-1}$ for 125/280 galaxies (44.6 per cent) with sizes of typically $\lesssim$100 pc. Of those 125, 106 targets show a core which coincides within 1.2 arcsec with the optical nucleus. Although we observed mostly cores, around one third of the detected galaxies features jetted morphologies. The detected radio core luminosities of the sample range between $\sim$10$^{34}$ and 10$^{40}$ erg s$^{-1}$. LINERs and Seyferts are the most luminous sources, whereas HII galaxies are the least. LINERs show FRI-like core-brightened radio structures, while Seyferts reveal the highest fraction of symmetric morphologies. The majority of HII galaxies have single radio core or complex extended structures, which probably conceal a nuclear starburst and/or a weak active nucleus (seven of them show clear jets). ALGs, which are typically found in evolved ellipticals, although the least numerous, exhibit on average the most luminous radio structures, similar to LINERs. △ Less

Submitted 5 November, 2020; originally announced November 2020.

Comments: Accepted for publication on MNRAS (48 pages, 7 figures, 10 tables, Appendix)

The complex multi-component outflow of the Seyfert galaxy NGC 7130

Authors: S. Comerón, J. H. Knapen, C. Ramos Almeida, A. E. Watkins

Abstract: AGN are a key ingredient for understanding galactic evolution. AGN-driven outflows are one of the manifestations of feedback. The AO mode for MUSE at the VLT permits to study the innermost tens of parsecs of nearby AGN in the optical. We present a detailed analysis of the ionised gas in the central regions of NGC 7130, an archetypical composite Seyfert and nuclear starburst galaxy. We achieve an a… ▽ More AGN are a key ingredient for understanding galactic evolution. AGN-driven outflows are one of the manifestations of feedback. The AO mode for MUSE at the VLT permits to study the innermost tens of parsecs of nearby AGN in the optical. We present a detailed analysis of the ionised gas in the central regions of NGC 7130, an archetypical composite Seyfert and nuclear starburst galaxy. We achieve an angular resolution of 0.17$^{\prime\prime}$ (50 pc). We performed a multi-component analysis of the main ISM lines and identified nine kinematic components, six of which correspond to the outflow. The outflow is biconic and has velocities of a few $100\,{\rm km\,s^{-1}}$ with respect to the disc. We decompose the approaching side of the outflow into a broad and a narrow component with typical velocity dispersions below and above $\sim200\,{\rm km\,s^{-1}}$, respectively. The blueshifted narrow component has substructure, in particular a collimated plume aligned with the radio jet, indicating that it may be jet-powered. The redshifted lobe is composed of two Narrow Components and a Broad Component. An additional redshifted component is seen outside the main outflow axis. Line ratio diagnostics indicate that the outflow gas in the main axis is AGN-powered whereas the off-axis component has LINER properties. The ionised gas mass outflow rate is $\dot{M}=1.2\pm0.7\,M_{\odot}\,{\rm yr^{-1}}$ and the kinetic power is $\dot{E}_{\rm kin}=(2.7\pm2.0)\times10^{41}\,{\rm erg\,s^{-1}}$, which corresponds to $F_{\rm kin}=0.12\pm0.09\%$ of the bolometric AGN power. The combination of high angular resolution integral field spectroscopy and a careful multi-component decomposition allows a uniquely detailed view of the outflow in NGC 7130, illustrating that AGN kinematics are more complex than traditionally derived from less sophisticated data and analyses. (abridged) △ Less

Submitted 7 April, 2021; v1 submitted 5 November, 2020; originally announced November 2020.

Comments: Published in A&A. This version of the preprint has been corrected with the modifications appearing in the A&A corrigendum

Journal ref: A&A 645, A130 (2021)

The kinematics of young and old stellar populations in nuclear rings of MUSE TIMER galaxies

Authors: D. Rosado-Belza, J. Falcón-Barroso, J. H. Knapen, A. Bittner, D. A. Gadotti, J. Neumann, A. de Lorenzo-Cáceres, J. Méndez-Abreu, M. Querejeta, I. Martín-Navarro, P. Sánchez-Blázquez, P. R. T. Coelho, M. Martig, G. van de Ven, T. Kim

Abstract: Studying the stellar kinematics of galaxies is a key tool in the reconstruction of their evolution. However, the current measurements of the stellar kinematics are complicated by several factors, including dust extinction and the presence of multiple stellar populations. We use integral field spectroscopic data of four galaxies from the TIMER survey to explore and compare the kinematics measured i… ▽ More Studying the stellar kinematics of galaxies is a key tool in the reconstruction of their evolution. However, the current measurements of the stellar kinematics are complicated by several factors, including dust extinction and the presence of multiple stellar populations. We use integral field spectroscopic data of four galaxies from the TIMER survey to explore and compare the kinematics measured in different spectral regions that are sensitive to distinct stellar populations. We derive the line-of-sight velocity and velocity dispersion of both a young (<2 Gyr) and an old stellar population from the spectral regions around the H$β$ line and the Ca II Triplet. In addition we obtain colour excess, mean age, and metallicity. We report a correlation of the colour excess with the difference in the kinematic parameters of the H$β$ line and the Ca II Triplet range, which are dominated by young and old stellar populations, respectively. Young stellar populations, located primarily in nuclear rings, have higher velocity dispersions than old ones. These differences in the rings are typically 10 km/s in velocity dispersion, but up to a mean value of 24 km/s in the most extreme case. Trends with age exist in the nuclear rings but are less significant than those with dust extinction. We report different degrees of correlation of these trends among the galaxies in the sample, which are related to the size of the Voronoi bins in their rings. No clear trends for the difference of line-of-sight velocity are observed. The absence of these trends can be explained as a consequence of the masking process of the H$β$ line during the kinematic extraction, as confirmed by dedicated simulations. Our study demonstrates that kinematic differences caused by different stellar populations can be identified in the central regions of nearby galaxies even from intermediate resolution spectroscopy. △ Less

Submitted 27 October, 2020; v1 submitted 22 October, 2020; originally announced October 2020.

Comments: 22 pages, 24 figures

Journal ref: A&A 644, A116 (2020)

The distribution of star formation in galactic bars as seen with H$α$ and stacked GALEX UV imaging

Authors: Simón Díaz-García, Facundo D. Moyano, Sébastien Comerón, Johan H. Knapen, Heikki Salo, Alexandre Y. K. Bouquin

Abstract: We investigate the spatial distribution of star formation (SF) within bars of nearby disk galaxies (inclination $< 65^{\circ}$) from the S$^4$G survey. We use archival GALEX far- and near-UV imaging for 772 barred galaxies. We also assemble a compilation of continuum-subtracted H$α$ images for 433 barred galaxies, of which 70 are produced by ourselves from ancillary photometry and MUSE/CALIFA IFU… ▽ More We investigate the spatial distribution of star formation (SF) within bars of nearby disk galaxies (inclination $< 65^{\circ}$) from the S$^4$G survey. We use archival GALEX far- and near-UV imaging for 772 barred galaxies. We also assemble a compilation of continuum-subtracted H$α$ images for 433 barred galaxies, of which 70 are produced by ourselves from ancillary photometry and MUSE/CALIFA IFU data cubes. We employ two complementary approaches: i) the analysis of bar/disk stacks built from co-added UV images of hundreds of galaxies; and ii) the classification of the morphology of ionised regions in galaxies into three main SF classes: A) only circumnuclear SF, B) SF at the bar ends, but not along the bar, and C) SF along the bar. Lenticular galaxies typically belong to SF class A: this is probably related to bar-induced SF quenching. The distribution of SF class B peaks for early- and intermediate-type spirals: this most likely results from the interplay of gas flow, shocks, and enhanced shear in centrally concentrated galaxies with large bar amplitudes. Late-type galaxies are mainly assigned to SF class C: we argue that this is a consequence of low shear. In bar stacks of spirals, the UV emission traces the stellar bars and dominates on their leading side, as witnessed in simulations. For early-types, the central UV emission is $\sim$0.5 mag brighter in strongly barred galaxies, relative to their weakly barred counterparts: this is related to the efficiency of strong bars sweeping the disk gas and triggering central starbursts. We also show that the distributions of SF in inner ringed galaxies are broadly the same in barred and non-barred galaxies, including a UV/H$α$ deficit in the middle part of the bar: this hints at the effect of resonance rings trapping gas. Distinct distributions of SF within bars are reported in galaxies of different morphological types (Abridged). △ Less

Submitted 20 November, 2020; v1 submitted 2 September, 2020; originally announced September 2020.

Comments: Accepted for publication in A&A (September 2, 2020). 18 pages, 14 figures, 2 tables, plus appendices (16 pages, 8 figures, 1 table)

Journal ref: A&A 644, A38 (2020)

An extensive spectroscopic time series of three Wolf-Rayet stars -- II. A search for wind asymmetries in the dust-forming WC7 binary WR137

Authors: N. St-Louis, C. Piaulet, N. D. Richardson, T. Shenar, A. F. J. Moffat, T. Eversberg, G. M. Hill, B. Gauza, J. H. Knapen, J. Kubat, B. Kubatova, D. P. Sablowski, S. Simon-Diaz, F. Bolduan, F. M. Dias, P. Dubreuil, D. Fuchs, T. Garrel, G. Grutzeck, T. Hunger, D. Kusters, M. Langenbrink, R. Leadbeater, D. Li, A. Lopez , et al. (17 additional authors not shown)

Abstract: We present the results of a four-month, spectroscopic campaign of the Wolf-Rayet dust-making binary, WR137. We detect only small-amplitude, random variability in the CIII5696 emission line and its integrated quantities (radial velocity, equivalent width, skewness, kurtosis) that can be explained by stochastic clumps in the wind of the WC star. We find no evidence of large-scale, periodic variation… ▽ More We present the results of a four-month, spectroscopic campaign of the Wolf-Rayet dust-making binary, WR137. We detect only small-amplitude, random variability in the CIII5696 emission line and its integrated quantities (radial velocity, equivalent width, skewness, kurtosis) that can be explained by stochastic clumps in the wind of the WC star. We find no evidence of large-scale, periodic variations often associated with Corotating Interaction Regions that could have explained the observed intrinsic continuum polarization of this star. Our moderately high-resolution and high signal-to-noise average Keck spectrum shows narrow double-peak emission profiles in the Halpha, Hbeta, Hgamma, HeII6678 and HeII5876 lines. These peaks have a stable blue-to-red intensity ratio with a mean of 0.997 and a root-mean-square of 0.004, commensurate with the noise level; no variability is found during the entire observing period. We suggest that these profiles arise in a decretion disk around the O9 companion, which is thus an O9e star. The characteristics of the profiles are compatible with those of other Be/Oe stars. The presence of this disk can explain the constant component of the continuum polarization of this system, for which the angle is perpendicular to the plane of the orbit, implying that the rotation axis of the O9e star is aligned with that of the orbit. It remains to be explained why the disk is so stable within the strong ultraviolet radiation field of the O star. We present a binary evolutionary scenario that is compatible with the current stellar and system parameters. △ Less

Submitted 17 July, 2020; originally announced July 2020.

Comments: 12 pages, 5 Figures

The number of globular clusters around the iconic UDG DF44 is as expected for dwarf galaxies

Authors: Teymoor Saifollahi, Ignacio Trujillo, Michael A. Beasley, Reynier F. Peletier, Johan H. Knapen

Abstract: There is a growing consensus that the vast majority of ultra-diffuse galaxies (UDGs) are dwarf galaxies. However, there remain a few UDGs that seem to be special in terms of their globular cluster (GC) systems. In particular, according to some authors, certain UDGs exhibit large GC populations when compared to expectations from their stellar (or total) mass. Among these special UDGs, DF44 in the C… ▽ More There is a growing consensus that the vast majority of ultra-diffuse galaxies (UDGs) are dwarf galaxies. However, there remain a few UDGs that seem to be special in terms of their globular cluster (GC) systems. In particular, according to some authors, certain UDGs exhibit large GC populations when compared to expectations from their stellar (or total) mass. Among these special UDGs, DF44 in the Coma cluster is one of the better-known examples. DF44 has been claimed to have a relatively high number of GCs, $N_{GC}=74^{+18}_{-18}$, for a stellar mass of only $3\times 10^8$ $M_{ \odot }$ which would indicate a much larger dark halo mass than dwarfs of similar stellar mass. In this paper we revisit this number and, contrary to previous results, find $N_{GC}=21^{+7}_{-9}$ assuming that the distribution of the GCs follows the same geometry as the galaxy. If we assume that the GCs around DF44 are distributed in a (projected) circularly symmetric way and, if we use a less strict criterion for the selection of the GCs, we find $N_{GC}=18^{+23}_{-12}$. Making use of the $M_{\rm GC} - M_{\rm halo}$ relation, this number of GCs suggests a dark matter halo mass of $M_{halo}=1.1^{+0.4}_{-0.5} \times 10^{11} M_{\odot}$, a value which is consistent with the expected total mass for DF44 based on its velocity dispersion, $σ=33^{+3}_{-3}$ km s$^{-1}$. We conclude that the number of GCs around DF44 is as expected for regular dwarf galaxies of similar stellar mass and DF44 is not extraordinary in this respect. △ Less

Submitted 29 September, 2020; v1 submitted 25 June, 2020; originally announced June 2020.

Comments: 16 pages, 13 figures, accepted for publication in MNRAS, minor changes on the text to match the accepted version