Showing 1–50 of 413 results for author: Baes, M

MICONIC: JWST/MIRI MRS observations of the nuclear and circumnuclear regions of Mrk231

Authors: A. Alonso-Herrero, L. Hermosa Muñoz, A. Labiano, P. Guillard, V. A. Buiten, D. Dicken, P. van der Werf, J. Álvarez-Márquez, T. Böker, L. Colina, A. Eckart, M. García-Marín, O. C. Jones, L. Pantoni, P. G. Pérez-González, D. Rouan, M. J. Ward, M. Baes, G. Östlin, P. Royer, G. S. Wright, M. Güdel, Th. Henning, P. -O. Lagage, E. F. van Dishoeck

Abstract: We present JWST/MIRI MRS spatially resolved $\sim 5-28\,μ$m observations of the central ~4-8kpc of the ultraluminous infrared galaxy and broad absorption line quasar Mrk231. These are part of the Mid-Infrared Characterization of Nearby Iconic galaxy Centers (MICONIC) program of the MIRI European Consortium guaranteed time observations. No high excitation lines (i.e., [MgV] at 5.61$μ$m or [NeV] at… ▽ More We present JWST/MIRI MRS spatially resolved $\sim 5-28\,μ$m observations of the central ~4-8kpc of the ultraluminous infrared galaxy and broad absorption line quasar Mrk231. These are part of the Mid-Infrared Characterization of Nearby Iconic galaxy Centers (MICONIC) program of the MIRI European Consortium guaranteed time observations. No high excitation lines (i.e., [MgV] at 5.61$μ$m or [NeV] at 14.32$μ$m) typically associated with the presence of an active galactic nucleus (AGN) are detected in the nuclear region of Mrk231. This is likely due to the intrinsically X-ray weak nature of its quasar. Some intermediate ionization potential lines, for instance, [ArIII] at 8.99$μ$m and [SIV] at 10.51$μ$m, are not detected either, even though they are clearly observed in a star-forming region ~920pc south-east of the AGN. Thus, the strong nuclear mid-infrared (mid-IR) continuum is also in part hampering the detection of faint lines in the nuclear region. The nuclear [NeIII]/[NeII]line ratio is consistent with values observed in star-forming galaxies. Moreover, we resolve for the first time the nuclear starburst in the mid-IR low-excitation line emission (size of ~400pc, FWHM). Several pieces of evidence also indicate that it is partly obscured even at these wavelengths. At the AGN position, the ionized and warm molecular gas emission lines have modest widths (W_80~300km/s). There are, however, weak blueshifted wings reaching velocities v_02~-400km/s in [NeII]. The nuclear starburst is at the center of a large (~8kpc), massive rotating disk with widely-spread, low velocity outflows. Given the high star formation rate of Mrk231, we speculate that part of the nuclear outflows and the large-scale non-circular motions observed in the mid-IR are driven by its powerful nuclear starburst. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: Accepted for publication in Astronomy & Astrophysics

The Size-Mass relation at Rest-Frame $1.5μ$m from JWST/NIRCam in the COSMOS-WEB and PRIMER-COSMOS fields

Authors: Marco Martorano, Arjen van der Wel, Maarten Baes, Eric F. Bell, Gabriel Brammer, Marijn Franx, Angelos Nersesian

Abstract: We present the galaxy stellar mass - size relation in the rest-frame near-IR ($1.5~μ{\text{m}}$) and its evolution with redshift up to $z=2.5$. Sérsic profiles are measured for $\sim$ $26\,000$ galaxies with stellar masses $M_\star > 10^9~{\text{M}}_\odot$ from JWST/NIRCam F277W and F444W imaging provided by the COSMOS-WEB and PRIMER surveys, using coordinates, redshifts, colors and stellar mass e… ▽ More We present the galaxy stellar mass - size relation in the rest-frame near-IR ($1.5~μ{\text{m}}$) and its evolution with redshift up to $z=2.5$. Sérsic profiles are measured for $\sim$ $26\,000$ galaxies with stellar masses $M_\star > 10^9~{\text{M}}_\odot$ from JWST/NIRCam F277W and F444W imaging provided by the COSMOS-WEB and PRIMER surveys, using coordinates, redshifts, colors and stellar mass estimates from the COSMOS2020 catalog. The new rest-frame near-IR effective radii are generally smaller than previously measured rest-frame optical sizes, on average by 0.14~dex, with no significant dependence on redshift. For quiescent galaxies this size offset does not depend on stellar mass, but for star-forming galaxies the offset increases from -0.1~dex at $M_\star = 10^{9.5}~{\text{M}}_\odot$ to -0.25~dex at $M_\star > 10^{11}~{\text{M}}_\odot$. That is, we find that the near-IR stellar mass - size relation for star-forming galaxies is flatter in the rest-frame near-IR than in the rest-frame optical at all redshifts $0.5<z<2.5$. The general pace of size evolution is the same in the near-IR as previously demonstrated in the optical, with slower evolution ($R_{\text{e}} \propto (1+z)^{-0.7}$) for $L^*$~star-forming galaxies and faster evolution ($R_{\text{e}} \propto (1+z)^{-1.3}$) for $L^*$~quiescent galaxies. Massive ($M_\star>10^{11}~{\text{M}}_\odot$) star-forming galaxies evolve in size almost as fast as quiescent galaxies. Low-mass ($M_\star<10^{10}~{\text{M}}_\odot$)~quiescent galaxies evolve as slow as star-forming galaxies. Our main conclusion is that the size evolution narrative as it has emerged over the past two decades does not radically change when accessing with JWST the rest-frame near-IR, a better proxy of the underlying stellar mass distribution. △ Less

Submitted 25 June, 2024; originally announced June 2024.

Comments: Accepted for publication

Self-consistent dynamical models with a finite extent -- IV. Wendland models based on compactly supported radial basis functions

Authors: Maarten Baes

Abstract: We present a new step in our systematic effort to develop self-consistent dynamical models with a finite radial extent. The focus is on models with simple analytical density profiles allowing for analytical calculations of many dynamical properties. In this paper, we introduce a family of models, termed Wendland models, based on compactly supported radial basis functions. The family of models is c… ▽ More We present a new step in our systematic effort to develop self-consistent dynamical models with a finite radial extent. The focus is on models with simple analytical density profiles allowing for analytical calculations of many dynamical properties. In this paper, we introduce a family of models, termed Wendland models, based on compactly supported radial basis functions. The family of models is characterised by a parameter $k$ that controls the smoothness of the transition at the truncation radius. In the limit $k\to\infty$, the Wendland model reduces to a non-truncated model with a Gaussian density profile. For each Wendland model, the density, mass and gravitational potential are simple truncated polynomial functions of radius. Via the SpheCow tool we demonstrate that all Wendland models can be supported by isotropic distribution functions. Surprisingly, the isotropic distribution function exhibits varied behaviour across different Wendland models. Additionally, each model can be supported by a continuum of Osipkov--Merritt orbital structures, ranging from radially anisotropic to completely tangential at the truncation radius. To the best of our knowledge, the Wendland models presented here are the first family of models accommodating both radial and tangential Osipkov--Merritt distribution functions. Using linear superposition, these models can easily be combined to generate Wendland models with even more diverse orbital structures. While the Wendland models are not fully representative of real dynamical systems due to their Gaussian-like density profile, this study lays important groundwork for constructing more realistic models with truncated density profiles that can be supported by a range of orbital structures. △ Less

Submitted 15 June, 2024; originally announced June 2024.

Comments: 12 pages, 7 figures, accepted for publication in MNRAS

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 -- Overview of the Perseus cluster and analysis of its luminosity and stellar mass functions

Authors: J. -C. Cuillandre, M. Bolzonella, A. Boselli, F. R. Marleau, M. Mondelin, J. G. Sorce, C. Stone, F. Buitrago, Michele Cantiello, K. George, N. A. Hatch, L. Quilley, F. Mannucci, T. Saifollahi, R. Sánchez-Janssen, F. Tarsitano, C. Tortora, X. Xu, H. Bouy, S. Gwyn, M. Kluge, A. Lançon, R. Laureijs, M. Schirmer, Abdurro'uf , et al. (177 additional authors not shown)

Abstract: The Euclid ERO programme targeted the Perseus cluster of galaxies, gathering deep data in the central region of the cluster over 0.7 square degree, corresponding to approximately 0.25 r_200. The data set reaches a point-source depth of IE=28.0 (YE, JE, HE = 25.3) AB magnitudes at 5 sigma with a 0.16" and 0.48" FWHM, and a surface brightness limit of 30.1 (29.2) mag per square arcsec. The exception… ▽ More The Euclid ERO programme targeted the Perseus cluster of galaxies, gathering deep data in the central region of the cluster over 0.7 square degree, corresponding to approximately 0.25 r_200. The data set reaches a point-source depth of IE=28.0 (YE, JE, HE = 25.3) AB magnitudes at 5 sigma with a 0.16" and 0.48" FWHM, and a surface brightness limit of 30.1 (29.2) mag per square arcsec. The exceptional depth and spatial resolution of this wide-field multi-band data enable the simultaneous detection and characterisation of both bright and low surface brightness galaxies, along with their globular cluster systems, from the optical to the NIR. This study advances beyond previous analyses of the cluster and enables a range of scientific investigations summarised here. We derive the luminosity and stellar mass functions (LF and SMF) of the Perseus cluster in the Euclid IE band, thanks to supplementary u,g,r,i,z and Halpha data from the CFHT. We adopt a catalogue of 1100 dwarf galaxies, detailed in the corresponding ERO paper. We identify all other sources in the Euclid images and obtain accurate photometric measurements using AutoProf or AstroPhot for 138 bright cluster galaxies, and SourceExtractor for half a million compact sources. Cluster membership for the bright sample is determined by calculating photometric redshifts with Phosphoros. Our LF and SMF are the deepest recorded for the Perseus cluster, highlighting the groundbreaking capabilities of the Euclid telescope. Both the LF and SMF fit a Schechter plus Gaussian model. The LF features a dip at M(IE)=-19 and a faint-end slope of alpha_S = -1.2 to -1.3. The SMF displays a low-mass-end slope of alpha_S = -1.2 to -1.35. These observed slopes are flatter than those predicted for dark matter halos in cosmological simulations, offering significant insights for models of galaxy formation and evolution. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Submitted to A&A, 44 pages, 35 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: 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 -- 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'

Towards an observationally motivated AGN dusty torus model. I. Dust chemical composition from the modeling of Spitzer spectra

Authors: Omar Ulises Reyes-Amador, Jacopo Fritz, Omaira González-Martín, Sundar Srinivasan, Maarten Baes, Enrique Lopez-Rodriguez, Natalia Osorio-Clavijo, Cesar Iván Victoria-Ceballos, Marko Stalevski, C. Ramos Almeida

Abstract: Spectral energy distribution (SED) fitting is one of most commonly used techniques to study the dust properties in Active Galactic Nuclei (AGN). Works implementing this technique commonly use radiative transfer models that assume a variety of dust properties. Despite the key role of this aspect, limited effort has been put forward to explore the chemical composition, the role of different optical… ▽ More Spectral energy distribution (SED) fitting is one of most commonly used techniques to study the dust properties in Active Galactic Nuclei (AGN). Works implementing this technique commonly use radiative transfer models that assume a variety of dust properties. Despite the key role of this aspect, limited effort has been put forward to explore the chemical composition, the role of different optical properties, and the grain size distribution of dust, all of which can have a substantial impact on the theoretical radiative transfer calculations. In this work, we explore the role of the dust chemical composition in the AGN dusty torus through SED fitting to \emph{Spitzer}/IRS spectra of a sample of 49 nearby AGN with silicate features in emission. We implement a mineralogy model including the popular astronomical silicates and a set of oxides and amorphous silicates with different grain sizes. We find that best fits use principally porous alumina, periclase, and olivine. In terms of mass fractions, $\sim99\%$ of the dust is composed of dust grains of size $\rm{0.1 μm}$, with a $<1\%$ contribution from $\rm{3 μm}$ grains. Moreover, the astronomical silicates have a very low occurrence in the best fits, suggesting that they are not the most suited dust species to reproduce the silicate features in our sample. △ Less

Submitted 14 May, 2024; originally announced May 2024.

Comments: 16 pages, 13 figures, Accepted by MNRAS

The many colors of the TNG100 simulation

Authors: Andrea Gebek, Ana Trčka, Maarten Baes, Marco Martorano, Annalisa Pillepich, Anand Utsav Kapoor, Angelos Nersesian, Arjen van der Wel

Abstract: We apply the 3D dust radiative transfer code SKIRT to the low-redshift ($z\leq0.1$) galaxy population in the TNG100 cosmological simulation, the fiducial run of the IllustrisTNG project. We compute global fluxes and spectral energy distributions (SEDs) from the far-ultraviolet to the sub-millimeter for $\approx\,$60 000 galaxies. Our post-processing methodology follows the study of Trčka et al. (2… ▽ More We apply the 3D dust radiative transfer code SKIRT to the low-redshift ($z\leq0.1$) galaxy population in the TNG100 cosmological simulation, the fiducial run of the IllustrisTNG project. We compute global fluxes and spectral energy distributions (SEDs) from the far-ultraviolet to the sub-millimeter for $\approx\,$60 000 galaxies. Our post-processing methodology follows the study of Trčka et al. (2022) of the higher-resolution TNG50 simulation. We verify that TNG100 reproduces observational luminosity functions at low redshifts to excellent precision, unlike TNG50. Additionally, we test the realism of our TNG100 plus SKIRT fluxes by comparing various flux and color relations to data from the GAMA survey. TNG100 broadly reproduces the observed distributions, but we predict ultraviolet colors that are too blue by $\approx\,$0.4 mag, possibly related to the extinction in the star-forming regions subgrid model not being selective enough. Furthermore, we find that the simulated galaxies exhibit mid-infrared fluxes elevated by up to $\approx\,$0.5 mag that we attribute to overly effective stochastic heating of the diffuse dust. All synthetic broadband fluxes and SEDs are made publicly available in three orientations and four apertures, and can readily be used to study TNG100 galaxies in a mock observational fashion. △ Less

Submitted 30 May, 2024; v1 submitted 8 May, 2024; originally announced May 2024.

Comments: Main text 14 pages, 5 figures, accepted to MNRAS. Generated data products are publicly available at www.tng-project.org/gebek24

JWST observations of the Horsehead photon-dominated region I. First results from multi-band near- and mid-infrared imaging

Authors: A. Abergel, K. Misselt, K. D. Gordon, A. Noriega-Crespo, P. Guillard, D. Van De Putte, A. N. Witt, N. Ysard, M. Baes, H. Beuther, P. Bouchet, B. R. Brandl, M. Elyajouri, O. Kannavou, S. Kendrew, P. Klassen, B. Trahin

Abstract: The JWST has captured the sharpest IR images ever taken of the Horsehead nebula, a prototypical moderately irradiated PDR that is fully representative of most of the UV-illuminated molecular gas in the Milky Way and star-forming galaxies. We investigate the impact of FUV radiation of a molecular cloud and constrain the structure of the edge of the PDR and its illumination conditions. We used NIRCa… ▽ More The JWST has captured the sharpest IR images ever taken of the Horsehead nebula, a prototypical moderately irradiated PDR that is fully representative of most of the UV-illuminated molecular gas in the Milky Way and star-forming galaxies. We investigate the impact of FUV radiation of a molecular cloud and constrain the structure of the edge of the PDR and its illumination conditions. We used NIRCam and MIRI to obtain 17 broadband and 6 narrowband maps from 0.7 to 28 $μ$m. We mapped the dust emission, scattered light, and several gas phase lines. We also used HST-WFC3 maps at 1.1 and 1. 6 $μ$m, along with HST-STIS spectroscopic observations of the H$α$ line. We probed the structure of the edge of the Horsehead and resolved its spatial complexity. We detected a network of faint striated features extending perpendicularly to the PDR front into the H\,II region in filters sensitive to nano-grain emission and light scattered by larger grains. This may indeed figure as the first detection of the entrainment of dust particles in the evaporative flow. The map of the 1-0 S(1) line of H$_2$ presents sharp sub-structures on scales as small as 1.5 arcsec. The ionization and dissociation fronts appear at distances 1-2 arcsec behind the edge of the PDR and seem to spatially coincide, indicating a thickness of the neutral atomic layer below 100 au. All broadband maps present strong color variations which can be explained by dust attenuation. Deviations of the emissions in the H$α$, Pa$α,$ and Br$α$ lines also indicate dust attenuation. With a very simple model, we derive the main features of the extinction curve. A small excess of extinction at 3 $μ$m may be attributed to icy H$_2$O mantles onto grains. In all lines of sight crossing the inner regions of the Horsehead, it appears that dust attenuation is non-negligible over the entire spectral range. △ Less

Submitted 24 April, 2024; originally announced April 2024.

Comments: 25 pages, 17 figures

MIGHTEE-HI: HI galaxy properties in the large scale structure environment at z~0.37 from a stacking experiment

Authors: Francesco Sinigaglia, Giulia Rodighiero, Ed Elson, Alessandro Bianchetti, Mattia Vaccari, Natasha Maddox, Anastasia A. Ponomareva, Bradley S. Frank, Matt J. Jarvis, Barbara Catinella, Luca Cortese, Sambit Roychowdhury, Maarten Baes, Jordan D. Collier, Olivier Ilbert, Ali A. Khostovan, Sushma Kurapati, Hengxing Pan, Isabella Prandoni, Sambatriniaina H. A. Rajohnson, Mara Salvato, Srikrishna Sekhar, Gauri Sharma

Abstract: We present the first measurement of HI mass of star-forming galaxies in different large scale structure environments from a blind survey at $z\sim 0.37$. In particular, we carry out a spectral line stacking analysis considering $2875$ spectra of colour-selected star-forming galaxies undetected in HI at $0.23 < z < 0.49$ in the COSMOS field, extracted from the MIGHTEE-HI Early Science datacubes, ac… ▽ More We present the first measurement of HI mass of star-forming galaxies in different large scale structure environments from a blind survey at $z\sim 0.37$. In particular, we carry out a spectral line stacking analysis considering $2875$ spectra of colour-selected star-forming galaxies undetected in HI at $0.23 < z < 0.49$ in the COSMOS field, extracted from the MIGHTEE-HI Early Science datacubes, acquired with the MeerKAT radio telescope. We stack galaxies belonging to different subsamples depending on three different definitions of large scale structure environment: local galaxy overdensity, position inside the host dark matter halo (central, satellite, or isolated), and cosmic web type (field, filament, or knot). We first stack the full star-forming galaxy sample and find a robust HI detection yielding an average galaxy HI mass of $M_{\rm HI}=(8.12\pm 0.75)\times 10^9\, {\rm M}_\odot$ at $\sim 11.8σ$. Next, we investigate the different subsamples finding a negligible difference in $M_{\rm HI}$ as a function of the galaxy overdensity. We report an HI excess compared to the full sample in satellite galaxies ($M_{\rm HI}=(11.31\pm1.22)\times 10^9$, at $\sim 10.2 σ$) and in filaments ($M_{\rm HI}=(11.62\pm 0.90)\times 10^9$. Conversely, we report non-detections for the central and knot galaxies subsamples, which appear to be HI-deficient. We find the same qualitative results also when stacking in units of HI fraction ($f_{\rm HI}$). We conclude that the HI amount in star-forming galaxies at the studied redshifts correlates with the large scale structure environment. △ Less

Submitted 1 March, 2024; originally announced March 2024.

Comments: Accepted for publication in MNRAS. 15 figures, 3 tables

Observational signatures of the dust size evolution in isolated galaxy simulations

Authors: Kosei Matsumoto, Hiroyuki Hirashita, Kentaro Nagamine, Stefan van der Giessen, Leonard E. C. Romano, Monica Relaño, Ilse De Looze, Maarten Baes, Angelos Nersesian, Peter Camps, Kuan-chou Hou, Yuri Oku

Abstract: We aim to provide observational signatures of the dust size evolution in the ISM, in particular exploring indicators of polycyclic aromatic hydrocarbon (PAH) mass fraction ($q_{PAH}$) defined as the mass fraction of PAHs relative to total dust grains. Additionally, we validate our dust evolution model by comparing the observational signatures from our simulations to those from observations. We mod… ▽ More We aim to provide observational signatures of the dust size evolution in the ISM, in particular exploring indicators of polycyclic aromatic hydrocarbon (PAH) mass fraction ($q_{PAH}$) defined as the mass fraction of PAHs relative to total dust grains. Additionally, we validate our dust evolution model by comparing the observational signatures from our simulations to those from observations. We model the evolution of grain size distribution of Milky Way-like and NGC 628-like galaxies representing star-forming galaxies with a hydrodynamic simulation code, GADGET4-OSAKA, which considers dust production and interstellar processing. Furthermore, we perform post-processing dust radiative transfer with SKIRT based on the simulations to predict the observational properties. We find that the intensity ratio between 8 um and 24 um correlates with $q_{PAH}$ and can be used as an indicator of PAH mass fraction. However, this ratio is influenced by the radiation field. As another indicator, we suggest the 8 um-to-total infrared intensity ratio ($νI_ν(8 μm)/I$(TIR)), which tightly correlates with $q_{PAH}$. Furthermore, we explore the spatial evolution of $q_{PAH}$ in the simulated Milky Way-like galaxy using $νI_ν(8 μm)/I$(TIR). We find that the spatially resolved $q_{PAH}$ increases with metallicity at lower metallicity (Z<0.2 Zsun) due to the interplay between accretion and shattering while it decreases with metallicity at higher metallicity (Z>0.2 Zsun) due to coagulation. Finally, we compare the above indicators in the NGC 628-like simulation with those observed in NGC 628. Consequently, our simulation underestimates the PAH mass fraction throughout the entire galaxy. This is probably because PAH is too efficiently lost by coagulation in the interstellar medium in our model, which suggests that the inhibition of coagulation of the PAHs is key to enhancing PAH formation. △ Less

Submitted 4 February, 2024; originally announced February 2024.

Comments: 16 pages, 12 figures, submit to A&A

The TNG50-SKIRT Atlas: wavelength dependence of the effective radius

Authors: Maarten Baes, Aleksandr Mosenkov, Raymond Kelly, Abdurro'uf, Nick Andreadis, Sena Bokona Tulu, Peter Camps, Abdissa Tassama Emana, Jacopo Fritz, Andrea Gebek, Inja Kovacic, Antonio La Marca, Marco Martorano, Angelos Nersesian, Vicente Rodriguez-Gomez, Crescenzo Tortora, Ana Trcka, Bert Vander Meulen, Arjen van der Wel, Lingyu Wang

Abstract: Galaxy sizes correlate with many other important properties of galaxies, and the cosmic evolution of galaxy sizes is an important observational diagnostic for constraining galaxy evolution models. The effective radius is probably the most widely used indicator of galaxy size. We used the TNG50-SKIRT Atlas to investigate the wavelength dependence of the effective radius of galaxies at optical and n… ▽ More Galaxy sizes correlate with many other important properties of galaxies, and the cosmic evolution of galaxy sizes is an important observational diagnostic for constraining galaxy evolution models. The effective radius is probably the most widely used indicator of galaxy size. We used the TNG50-SKIRT Atlas to investigate the wavelength dependence of the effective radius of galaxies at optical and near-infrared (NIR) wavelengths. We find that, on average, the effective radius in every band exceeds the stellar mass effective radius, and that this excess systematically decreases with increasing wavelength. The optical g-band (NIR Ks-band) effective radius is on average 58% (13%) larger than the stellar mass effective radius. Effective radii measured from dust-obscured images are systematically larger than those measured from dust-free images, although the effect is limited (8.7% in the g-band, 2.1% in the Ks-band). We find that stellar population gradients are the dominant factor (about 80%) in driving the wavelength dependence of the effective radius, and that differential dust attenuation is a secondary factor (20%). Comparing our results to recent observational data, we find offsets in the absolute values of the median effective radii, up to 50% for the population of blue galaxies. We find better agreement in the slope of the wavelength dependence of the effective radius, with red galaxies having a slightly steeper slope than green-blue galaxies. Comparing our effective radii with those of galaxies from the Siena Galaxy Atlas in separate bins in z-band absolute magnitude and g-z colour, we find excellent agreement for the reddest galaxies, but again significant offsets for the blue populations: up to 70% for galaxies around Mz=-21.5. This difference in median effective radius for the bluer galaxies is most probably due to (abridged...). △ Less

Submitted 8 January, 2024; originally announced January 2024.

Comments: 11 pages, 10 figures, accepted for publication in A&A

The TNG50-SKIRT Atlas: post-processing methodology and first data release

Authors: Maarten Baes, Andrea Gebek, Ana Trcka, Peter Camps, Arjen van der Wel, Abdurro'uf, Nick Andreadis, Sena Bokona Tulu, Abdissa Tassama Emana, Jacopo Fritz, Raymond Kelly, Inja Kovacic, Antonio La Marca, Marco Martorano, Aleksandr Mosenkov, Angelos Nersesian, Vicente Rodriguez-Gomez, Crescenzo Tortora, Bert Vander Meulen, Lingyu Wang

Abstract: Galaxy morphology is a powerful diagnostic to assess the realism of cosmological hydrodynamical simulations. Determining the morphology of simulated galaxies requires the generation of synthetic images through 3D radiative transfer post-processing that properly accounts for different stellar populations and interstellar dust attenuation. We use the SKIRT code to generate the TNG50-SKIRT Atlas, a s… ▽ More Galaxy morphology is a powerful diagnostic to assess the realism of cosmological hydrodynamical simulations. Determining the morphology of simulated galaxies requires the generation of synthetic images through 3D radiative transfer post-processing that properly accounts for different stellar populations and interstellar dust attenuation. We use the SKIRT code to generate the TNG50-SKIRT Atlas, a synthetic UV to near-infrared broadband image atlas for a complete stellar-mass selected sample of 1154 galaxies extracted from the TNG50 cosmological simulation at $z=0$. The images have a high spatial resolution (100 pc) and a wide field of view (160 kpc). In addition to the dust-obscured images, we also release dust-free images and physical parameter property maps with matching characteristics. As a sanity check and preview application we discuss the UVJ diagram of the galaxy sample. We investigate the effect of dust attenuation on the UVJ diagram and find that it affects both the star-forming and the quiescent galaxy populations. The quiescent galaxy region is polluted by younger and star-forming highly inclined galaxies, while dust attenuation induces a separation in inclination of the star-forming galaxy population, with low-inclination galaxies remaining at the blue side of the diagram and high-inclination galaxies systematically moving towards the red side. This image atlas can be used for a variety of other applications, including galaxy morphology studies and the investigation of local scaling relations. We publicly release the images and parameter maps, and we invite the community to use them. △ Less

Submitted 8 January, 2024; originally announced January 2024.

Comments: 13 pages, 7 figures, accepted for publication in A&A

The discovery of a z=0.7092 OH megamaser with the MIGHTEE survey

Authors: Matt J. Jarvis, Ian Heywood, Sophie M. Jewell, Roger P. Deane, H. -R. Klöckner, Anastasia A. Ponomareva, Natasha Maddox, Andrew J. Baker, Alessandro Bianchetti, Kelley M. Hess, Hayley Roberts, Giulia Rodighiero, Ilaria Ruffa, Francesco Sinigaglia, R. G. Varadaraj, I. H. Whittam, Elizabeth A. K. Adams, Maarten Baes, Eric J. Murphy, Hengxing Pan, Mattia Vaccari

Abstract: We present the discovery of the most distant OH megamaser to be observed in the main lines, using data from the MeerKAT International Giga-Hertz Tiered Extragalactic Exploration (MIGHTEE) survey. At a newly measured redshift of $z = 0.7092$, the system has strong emission in both the 1665MHz ($L \approx 2500$ L$_{\odot}$) and 1667 MHz ($L \approx 4.5\times10^4$ L$_{\odot}$) transitions, with both… ▽ More We present the discovery of the most distant OH megamaser to be observed in the main lines, using data from the MeerKAT International Giga-Hertz Tiered Extragalactic Exploration (MIGHTEE) survey. At a newly measured redshift of $z = 0.7092$, the system has strong emission in both the 1665MHz ($L \approx 2500$ L$_{\odot}$) and 1667 MHz ($L \approx 4.5\times10^4$ L$_{\odot}$) transitions, with both narrow and broad components. We interpret the broad line as a high-velocity-dispersion component of the 1667 MHz transition, with velocity $v \sim 330$km s$^{-1}$ with respect to the systemic velocity. The host galaxy has a stellar mass of $M_{\star} = 2.95 \times 10^{10}$ M$_{\odot}$ and a star-formation rate of SFR = 371 M$_{\odot}$yr$^{-1}$, placing it $\sim 1.5$dex above the main sequence for star-forming galaxies at this redshift, and can be classified as an ultra-luminous infrared galaxy. Alongside the optical imaging data, which exhibits evidence for a tidal tail, this suggests that the OH megamaser arises from a system that is currently undergoing a merger, which is stimulating star formation and providing the necessary conditions for pumping the OH molecule to saturation. The OHM is likely to be lensed, with a magnification factor of $\sim 2.5$, and perhaps more if the maser emitting region is compact and suitably offset relative to the centroid of its host galaxy's optical light. This discovery demonstrates that spectral line mapping with the new generation of radio interferometers may provide important information on the cosmic merger history of galaxies. △ Less

Submitted 7 December, 2023; originally announced December 2023.

Comments: 11 pages, accepted for publication in MNRAS

Voice Conversion for Stuttered Speech, Instruments, Unseen Languages and Textually Described Voices

Authors: Matthew Baas, Herman Kamper

Abstract: Voice conversion aims to convert source speech into a target voice using recordings of the target speaker as a reference. Newer models are producing increasingly realistic output. But what happens when models are fed with non-standard data, such as speech from a user with a speech impairment? We investigate how a recent voice conversion model performs on non-standard downstream voice conversion ta… ▽ More Voice conversion aims to convert source speech into a target voice using recordings of the target speaker as a reference. Newer models are producing increasingly realistic output. But what happens when models are fed with non-standard data, such as speech from a user with a speech impairment? We investigate how a recent voice conversion model performs on non-standard downstream voice conversion tasks. We use a simple but robust approach called k-nearest neighbors voice conversion (kNN-VC). We look at four non-standard applications: stuttered voice conversion, cross-lingual voice conversion, musical instrument conversion, and text-to-voice conversion. The latter involves converting to a target voice specified through a text description, e.g. "a young man with a high-pitched voice". Compared to an established baseline, we find that kNN-VC retains high performance in stuttered and cross-lingual voice conversion. Results are more mixed for the musical instrument and text-to-voice conversion tasks. E.g., kNN-VC works well on some instruments like drums but not on others. Nevertheless, this shows that voice conversion models - and kNN-VC in particular - are increasingly applicable in a range of non-standard downstream tasks. But there are still limitations when samples are very far from the training distribution. Code, samples, trained models: https://rf5.github.io/sacair2023-knnvc-demo/. △ Less

Submitted 12 October, 2023; originally announced October 2023.

Comments: 11 pages, 1 figure, 5 tables. Accepted at SACAIR 2023

IAS/CEA Evolution of Dust in Nearby Galaxies (ICED): the spatially-resolved dust properties of NGC4254

Authors: L. Pantoni, R. Adam, P. Ade, H. Ajeddig, P. André, E. Artis, H. Aussel, M. Baes, A. Beelen, A. Benoît, S. Berta, L. Bing, O. Bourrion, M. Calvo, A. Catalano, M. De Petris, F. -X. Désert, S. Doyle, E. F. C. Driessen, G. Ejlali, F. Galliano, A. Gomez, J. Goupy, A. P. Jones, C. Hanser , et al. (35 additional authors not shown)

Abstract: We present the first preliminary results of the project \textit{ICED}, focusing on the face-on galaxy NGC4254. We use the millimetre maps observed with NIKA2 at IRAM-30m, as part of the IMEGIN Guaranteed Time Large Program, and of a wide collection of ancillary data (multi-wavelength photometry and gas phase spectral lines) that are publicly available. We derive the global and local properties of… ▽ More We present the first preliminary results of the project \textit{ICED}, focusing on the face-on galaxy NGC4254. We use the millimetre maps observed with NIKA2 at IRAM-30m, as part of the IMEGIN Guaranteed Time Large Program, and of a wide collection of ancillary data (multi-wavelength photometry and gas phase spectral lines) that are publicly available. We derive the global and local properties of interstellar dust grains through infrared-to-radio spectral energy distribution fitting, using the hierarchical Bayesian code HerBIE, which includes the grain properties of the state-of-the-art dust model, THEMIS. Our method allows us to get the following dust parameters: dust mass, average interstellar radiation field, and fraction of small grains. Also, it is effective in retrieving the intrinsic correlations between dust parameters and interstellar medium properties. We find an evident anti-correlation between the interstellar radiation field and the fraction of small grains in the centre of NGC4254, meaning that, at strong radiation field intensities, very small amorphous carbon grains are efficiently destroyed by the ultra-violet photons coming from newly formed stars, through photo-desorption and sublimation. We observe a flattening of the anti-correlation at larger radial distances, which may be driven by the steep metallicity gradient measured in NGC4254. △ Less

Submitted 10 October, 2023; originally announced October 2023.

Comments: to appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

Journal ref: EPJ Web of Conferences 293 (2024) 00038

Exploring the interstellar medium of NGC 891 at millimeter wavelengths using the NIKA2 camera

Authors: S. Katsioli, R. Adam, P. Ade, H. Ajeddig, P. André, E. Artis, H. Aussel, M. Baes, A. Beelen, A. Benoît, S. Berta, L. Bing, O. Bourrion, M. Calvo, A. Catalano, C. J. R. Clark, I. De Looze, M. De Petris, F. -X. Désert, S. Doyle, E. F. C. Driessen, G. Ejlali, M. Galametz, F. Galliano, A. Gomez , et al. (39 additional authors not shown)

Abstract: In the framework of the IMEGIN Large Program, we used the NIKA2 camera on the IRAM 30-m telescope to observe the edge-on galaxy NGC 891 at 1.15 mm and 2 mm and at a FWHM of 11.1" and 17.6", respectively. Multiwavelength data enriched with the new NIKA2 observations fitted by the HerBIE SED code (coupled with the THEMIS dust model) were used to constrain the physical properties of the ISM. Emission… ▽ More In the framework of the IMEGIN Large Program, we used the NIKA2 camera on the IRAM 30-m telescope to observe the edge-on galaxy NGC 891 at 1.15 mm and 2 mm and at a FWHM of 11.1" and 17.6", respectively. Multiwavelength data enriched with the new NIKA2 observations fitted by the HerBIE SED code (coupled with the THEMIS dust model) were used to constrain the physical properties of the ISM. Emission originating from the diffuse dust disk is detected at all wavelengths from mid-IR to mm, while mid-IR observations reveal warm dust emission from compact HII regions. Indications of mm excess emission have also been found in the outer parts of the galactic disk. Furthermore, our SED fitting analysis constrained the mass fraction of the small (< 15 Angstrom) dust grains. We found that small grains constitute 9.5% of the total dust mass in the galactic plane, but this fraction increases up to ~ 20% at large distances (|z| > 3 kpc) from the galactic plane. △ Less

Submitted 6 October, 2023; originally announced October 2023.

Comments: To appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

Journal ref: EPJ Web of conferences 293 (2024) 00026

Constraining Millimeter Dust Emission in Nearby Galaxies with NIKA2: the case of NGC2146 and NGC2976

Authors: G. Ejlali, R. Adam, P. Ade, H. Ajeddig, P. André, E. Artis, H. Aussel, M. Baes, A. Beelen, Benoît, S. Berta, L. Bing, O. Bourrion, M. Calvo, A. Catalano, M. De Petris, F. -X. Désert, S. Doyle, E. F. C. Driessen, F. Galliano, A. Gomez, J. Goupy, A. P. Jones, C. Hanser, A. Hughes , et al. (35 additional authors not shown)

Abstract: This study presents the first millimeter continuum mapping observations of two nearby galaxies, the starburst spiral galaxy NGC2146 and the dwarf galaxy NGC2976, at 1.15 mm and 2 mm using the NIKA2 camera on the IRAM 30m telescope, as part of the Guaranteed Time Large Project IMEGIN. These observations provide robust resolved information about the physical properties of dust in nearby galaxies by… ▽ More This study presents the first millimeter continuum mapping observations of two nearby galaxies, the starburst spiral galaxy NGC2146 and the dwarf galaxy NGC2976, at 1.15 mm and 2 mm using the NIKA2 camera on the IRAM 30m telescope, as part of the Guaranteed Time Large Project IMEGIN. These observations provide robust resolved information about the physical properties of dust in nearby galaxies by constraining their FIR-radio SED in the millimeter domain. After subtracting the contribution from the CO line emission, the SEDs are modeled spatially using a Bayesian approach. Maps of dust mass surface density, temperature, emissivity index, and thermal radio component of the galaxies are presented, allowing for a study of the relations between the dust properties and star formation activity (using observations at 24$μ$m as a tracer). We report that dust temperature is correlated with star formation rate in both galaxies. The effect of star formation activity on dust temperature is stronger in NGC2976, an indication of the thinner interstellar medium of dwarf galaxies. Moreover, an anti-correlation trend is reported between the dust emissivity index and temperature in both galaxies. △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: To appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

Journal ref: EPJ Web of conferences 293 (2024) 00016

TODDLERS: A new UV-mm emission library for star-forming regions. 1. Integration with SKIRT and public release

Authors: Anand Utsav Kapoor, Maarten Baes, Arjen van der Wel, Andrea Gebek, Peter Camps, Angelos Nersesian, Sharon E. Meidt, Aaron Smith, Sebastien Vicens, Francesco D'Eugenio, Marco Martorano, Daniela Barrientos, Nina Sanches Sartorio

Abstract: We present and publicly release a new star-forming regions emission library TODDLERS (Time evolution of Observables including Dust Diagnostics and Line Emission from Regions containing young Stars) for the publicly available radiative transfer code SKIRT. The library generation involves the spherical evolution of a homogeneous gas cloud around a young stellar cluster that accounts for stellar feed… ▽ More We present and publicly release a new star-forming regions emission library TODDLERS (Time evolution of Observables including Dust Diagnostics and Line Emission from Regions containing young Stars) for the publicly available radiative transfer code SKIRT. The library generation involves the spherical evolution of a homogeneous gas cloud around a young stellar cluster that accounts for stellar feedback processes including stellar winds, supernovae, and radiation pressure, as well as the gravitational forces on the gas. The semi-analytical evolution model is coupled with the photoionization code Cloudy to calculate time-dependent UV-mm spectral energy distributions (SEDs) from star-forming regions of varying metallicity, star-formation efficiency, birth-cloud density, and mass. The calculated SEDs include the stellar, nebular, and dust continuum emission along with a wide range of emission lines originating from H ii, photo-dissociation, and molecular gas regimes tabulated at high resolution. The SEDs incorporated in SKIRT are generated by calculating a stellar-mass normalized luminosity, which assumes that each emission source is composed of a power-law population of star-forming clouds. When compared to the previous treatment of star-forming regions in SKIRT, TODDLERS shows a better agreement with low-redshift observational data in the IR wavelength range while offering a more comprehensive line-emission support. This paves the way for a variety of applications using simulated galaxies at low and high redshift. △ Less

Submitted 20 October, 2023; v1 submitted 30 September, 2023; originally announced October 2023.

Comments: Accepted for publication in MNRAS on 27/09/23

The stratification of ISM properties in the edge-on galaxy NGC 891 revealed by NIKA2

Authors: S. Katsioli, E. M. Xilouris, C. Kramer, R. Adam, P. Ade, H. Ajeddig, P. André, E. Artis, H. Aussel, M. Baes, A. Beelen, A. Benoît, S. Berta, L. Bing, O. Bourrion, M. Calvo, A. Catalano, C. J. R. Clark, I. De Looze, M. De Petris, F. -X. Désert, S. Doyle, E. F. C. Driessen, G. Ejlali, M. Galametz , et al. (38 additional authors not shown)

Abstract: As the millimeter wavelength range remains a largely unexplored spectral region for galaxies, the IMEGIN large program aims to map the millimeter continuum emission of 22 nearby galaxies at 1.15 and 2 mm. Using the high-resolution maps produced by the NIKA2 camera, we explore the existence of very cold dust and take possible contamination by free-free and synchrotron emission into account. We stud… ▽ More As the millimeter wavelength range remains a largely unexplored spectral region for galaxies, the IMEGIN large program aims to map the millimeter continuum emission of 22 nearby galaxies at 1.15 and 2 mm. Using the high-resolution maps produced by the NIKA2 camera, we explore the existence of very cold dust and take possible contamination by free-free and synchrotron emission into account. We study the IR-to-radio emission coming from different regions along the galactic plane and at large vertical distances. New observations of NGC 891, using the NIKA2 camera on the IRAM 30m telescope, along with a suite of observations at other wavelengths were used to perform a multiwavelength study of the spectral energy distribution in the interstellar medium in this galaxy. This analysis was performed globally and locally, using the advanced hierarchical Bayesian fitting code, HerBIE, coupled with the THEMIS dust model. Our dust modeling is able to reproduce the near-IR to millimeter emission of NGC 891, with the exception of an excess at a level of 25% obtained by the NIKA2 observations in the outermost parts of the disk. The radio continuum and thermal dust emission are distributed differently in the disk and galaxy halo. Different dusty environments are also revealed by a multiwavelength investigation of the emission features. Our detailed decomposition at millimeter and centimeter wavelengths shows that emission at 1 mm is purely originated by dust. Radio components become progressively important with increasing wavelengths. Finally, we find that emission arising from small dust grains accounts for ~ 9.5% of the total dust mass, reaching up to 20% at large galactic latitudes. Shock waves in the outflows that shatter the dust grains might explain this higher fraction of small grains in the halo. △ Less

Submitted 15 September, 2023; originally announced September 2023.

Comments: 15 pages, 10 figures, Accepted for publication in A&A

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

The problem of dust attenuation in photometric decomposition of edge-on galaxies and possible solutions

Authors: Sergey Savchenko, Denis Poliakov, Aleksandr Mosenkov, Anton Smirnov, Alexander Marchuk, Vladimir Il'in, George Gontcharov, Jonah Seguine, Maarten Baes

Abstract: The presence of dust in spiral galaxies affects the ability of photometric decompositions to retrieve the parameters of their main structural components. For galaxies in an edge-on orientation, the optical depth integrated over the line-of-sight is significantly higher than for those with intermediate or face-on inclinations, so it is only natural to expect that for edge-on galaxies, dust attenuat… ▽ More The presence of dust in spiral galaxies affects the ability of photometric decompositions to retrieve the parameters of their main structural components. For galaxies in an edge-on orientation, the optical depth integrated over the line-of-sight is significantly higher than for those with intermediate or face-on inclinations, so it is only natural to expect that for edge-on galaxies, dust attenuation should severely influence measured structural parameters. In this paper, we use radiative transfer simulations to generate a set of synthetic images of edge-on galaxies which are then analysed via decomposition. Our results demonstrate that for edge-on galaxies, the observed systematic errors of the fit parameters are significantly higher than for moderately inclined galaxies. Even for models with a relatively low dust content, all structural parameters suffer offsets that are far from negligible. In our search for ways to reduce the impact of dust on retrieved structural parameters, we test several approaches, including various masking methods and an analytical model that incorporates dust absorption. We show that using such techniques greatly improves the reliability of decompositions for edge-on galaxies. △ Less

Submitted 12 September, 2023; originally announced September 2023.

Comments: 17 pages, 10 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 524, Issue 3, pp.4729-4745 (2023)

Self-consistent dust and non-LTE line radiative transfer with SKIRT

Authors: Kosei Matsumoto, Peter Camps, Maarten Baes, Frederik De Ceuster, Keiichi Wada, Takao Nakagawa, Kentaro Nagamine

Abstract: We introduce Monte Carlo-based non-LTE line radiative transfer calculations in the 3D dust radiative transfer code SKIRT, which was originally set up as a dust radiative transfer code. By doing so, we develop a generic and powerful 3D radiative transfer code that can self-consistently generate spectra with molecular and atomic lines against the underlying continuum. We test the accuracy of the non… ▽ More We introduce Monte Carlo-based non-LTE line radiative transfer calculations in the 3D dust radiative transfer code SKIRT, which was originally set up as a dust radiative transfer code. By doing so, we develop a generic and powerful 3D radiative transfer code that can self-consistently generate spectra with molecular and atomic lines against the underlying continuum. We test the accuracy of the non-LTE line radiative transfer module in the extended SKIRT code using standard benchmarks. We find excellent agreement between the SKIRT results, the published benchmark results, and results obtained using the ray-tracing non-LTE line radiative transfer code MAGRITTE, which validates our implementation. We apply the extended SKIRT code on a 3D hydrodynamic simulation of a dusty AGN torus model and generate multi-wavelength images with CO rotational-line spectra against the underlying dust continuum. We find that the low-J CO emission traces the geometrically thick molecular torus, whereas the higher-J CO lines originate from the gas with high kinetic temperature located in the innermost regions of the torus. Comparing the calculations with and without dust radiative transfer, we find that higher-J CO lines are slightly attenuated by the surrounding cold dust when seen edge-on. This shows that atomic and molecular lines can experience attenuation, an effect that is particularly important for transitions at mid- and near-infrared wavelengths. Therefore, our self-consistent dust and non-LTE line radiative transfer calculations can help interpret the observational data from Herschel, ALMA, and JWST. △ Less

Submitted 5 September, 2023; originally announced September 2023.

Comments: 12 pages, 9 figures, Accepted for publication in A&A

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

Rest-Frame Near-Infrared Radial Light Profiles up to z=3 from JWST/NIRCam: Wavelength Dependence of the Sérsic Index

Authors: Marco Martorano, Arjen van der Wel, Eric F. Bell, Marijn Franx, Katherine E. Whitaker, Angelos Nersesian, Sedona H. Price, Maarten Baes, Katherine A. Suess, Erica J. Nelson, Tim B. Miller, Rachel Bezanson, Gabriel Brammer

Abstract: We examine the wavelength dependence of radial light profiles based on Sérsic index $n$ measurements of 1067 galaxies with M$_*\geq$ 10$^{9.5}$M$_\odot$ and in the redshift range $0.5 < z < 3$. The sample and rest-frame optical light profiles are drawn from CANDELS$+$3D-HST; rest-frame near-infrared light profiles are inferred from CEERS JWST/NIRCam imaging. $n$ shows only weak dependence on wavel… ▽ More We examine the wavelength dependence of radial light profiles based on Sérsic index $n$ measurements of 1067 galaxies with M$_*\geq$ 10$^{9.5}$M$_\odot$ and in the redshift range $0.5 < z < 3$. The sample and rest-frame optical light profiles are drawn from CANDELS$+$3D-HST; rest-frame near-infrared light profiles are inferred from CEERS JWST/NIRCam imaging. $n$ shows only weak dependence on wavelength, regardless of redshift, galaxy mass and type: on average, star-forming galaxies have $n = 1-1.5$ and quiescent galaxies have $n = 3-4$ in the rest-frame optical and near-infrared. The strong correlation at all wavelengths between $n$ and star-formation activity implies a physical connection between the radial stellar mass profile and star-formation activity. The main caveat is that the current sample is too small to discern trends for the most massive galaxies (M$_* > 10^{11}M_\odot$). △ Less

Submitted 22 August, 2023; originally announced August 2023.

Comments: 15 pages, 7 figures, Submitted to ApJ. Comments welcome

Journal ref: ApJ 957 46 (2023)

Self-consistent dynamical models with a finite extent -- III. Truncated power-law spheres

Authors: Maarten Baes, Bert Vander Meulen

Abstract: Fully analytical dynamical models usually have an infinite extent, while real star clusters, galaxies, and dark matter haloes have a finite extent. The standard method for generating dynamical models with a finite extent consists of taking a model with an infinite extent and applying a truncation in binding energy. This method, however, cannot be used to generate models with a pre-set analytical m… ▽ More Fully analytical dynamical models usually have an infinite extent, while real star clusters, galaxies, and dark matter haloes have a finite extent. The standard method for generating dynamical models with a finite extent consists of taking a model with an infinite extent and applying a truncation in binding energy. This method, however, cannot be used to generate models with a pre-set analytical mass density profile. We investigate the self-consistency and dynamical properties of a family of power-law spheres with a general tangential Cuddeford (TC) orbital structure. By varying the density power-law slope $γ$ and the central anisotropy $β_0$, these models cover a wide parameter space in density and anisotropy profiles. We explicitly calculate the phase-space distribution function for various parameter combinations, and interpret our results in terms of the energy distribution of bound orbits. We find that truncated power-law spheres can be supported by a TC orbital structure if and only if $γ\geqslant 2β_0$, which means that the central density slope-anisotropy inequality is both a sufficient and a necessary condition for this family. We provide closed expressions for structural and dynamical properties such as the radial and tangential velocity dispersion profiles, which can be compared against more complex numerical modelling results. This work significantly adds to the available suite of self-consistent dynamical models with a finite extent and an analytical description. △ Less

Submitted 1 August, 2023; originally announced August 2023.

Comments: 12 pages, 4 figures, accepted for publication in MNRAS

Disentanglement in a GAN for Unconditional Speech Synthesis

Authors: Matthew Baas, Herman Kamper

Abstract: Can we develop a model that can synthesize realistic speech directly from a latent space, without explicit conditioning? Despite several efforts over the last decade, previous adversarial and diffusion-based approaches still struggle to achieve this, even on small-vocabulary datasets. To address this, we propose AudioStyleGAN (ASGAN) -- a generative adversarial network for unconditional speech syn… ▽ More Can we develop a model that can synthesize realistic speech directly from a latent space, without explicit conditioning? Despite several efforts over the last decade, previous adversarial and diffusion-based approaches still struggle to achieve this, even on small-vocabulary datasets. To address this, we propose AudioStyleGAN (ASGAN) -- a generative adversarial network for unconditional speech synthesis tailored to learn a disentangled latent space. Building upon the StyleGAN family of image synthesis models, ASGAN maps sampled noise to a disentangled latent vector which is then mapped to a sequence of audio features so that signal aliasing is suppressed at every layer. To successfully train ASGAN, we introduce a number of new techniques, including a modification to adaptive discriminator augmentation which probabilistically skips discriminator updates. We apply it on the small-vocabulary Google Speech Commands digits dataset, where it achieves state-of-the-art results in unconditional speech synthesis. It is also substantially faster than existing top-performing diffusion models. We confirm that ASGAN's latent space is disentangled: we demonstrate how simple linear operations in the space can be used to perform several tasks unseen during training. Specifically, we perform evaluations in voice conversion, speech enhancement, speaker verification, and keyword classification. Our work indicates that GANs are still highly competitive in the unconditional speech synthesis landscape, and that disentangled latent spaces can be used to aid generalization to unseen tasks. Code, models, samples: https://github.com/RF5/simple-asgan/ △ Less

Submitted 25 January, 2024; v1 submitted 4 July, 2023; originally announced July 2023.

Comments: 12 pages, 5 tables, 4 figures. Accepted to IEEE TASLP. arXiv admin note: substantial text overlap with arXiv:2210.05271

Spatially resolved mock observations of stellar kinematics: full radiative transfer treatment of simulated galaxies

Authors: Daniela Barrientos Acevedo, Arjen van der Wel, Maarten Baes, Robert J. J Grand, Anand Utsav Kapoor, Peter Camps, Anna de Graaff, Caroline M. S. Straatman, Rachel Bezanson

Abstract: We present a framework to build realistic mock spectroscopic observations for state-of-the-art hydrodynamical simulations, using high spectral resolution stellar population models and full radiative transfer treatment with SKIRT. As a first application we generate stellar continuum mock observations for the Auriga cosmological zoom simulations emulating integral-field observations from the SAMI Ga… ▽ More We present a framework to build realistic mock spectroscopic observations for state-of-the-art hydrodynamical simulations, using high spectral resolution stellar population models and full radiative transfer treatment with SKIRT. As a first application we generate stellar continuum mock observations for the Auriga cosmological zoom simulations emulating integral-field observations from the SAMI Galaxy survey. We perform spectral fitting on our synthetic cubes and compute the resulting rotation velocity ($V_{\rm{rot}}$) and velocity dispersion within 1$R_{\text{e}}$ ($σ_{\text{e}}$) for a sub-set of the Auriga sample. We find that the kinematics produced by Auriga are in good agreement with the observations from the SAMI Galaxy survey after taking into account the effects of dust and the systematics produced by the observation limitations. We also explore the effects of seeing convolution, inclination, and attenuation on the line-of-sight velocity distribution. For highly inclined galaxies, these effects can lead to an artificial decrease in the measured $V/σ$ by nearly a factor two (after inclination correction). We also demonstrate the utility of our method for high-redshift galaxies by emulating spatially resolved continuum spectra from the LEGA-C survey and, looking forward, E-ELT HARMONI. Our framework represents a crucial link between the ground truth for stellar populations and kinematics in simulations and the observed stellar continuum observations at low and high redshift. △ Less

Submitted 13 June, 2023; originally announced June 2023.

Comments: 17 pages, 11 figures. Accepted for publication in MNRAS

Voice Conversion With Just Nearest Neighbors

Authors: Matthew Baas, Benjamin van Niekerk, Herman Kamper

Abstract: Any-to-any voice conversion aims to transform source speech into a target voice with just a few examples of the target speaker as a reference. Recent methods produce convincing conversions, but at the cost of increased complexity -- making results difficult to reproduce and build on. Instead, we keep it simple. We propose k-nearest neighbors voice conversion (kNN-VC): a straightforward yet effecti… ▽ More Any-to-any voice conversion aims to transform source speech into a target voice with just a few examples of the target speaker as a reference. Recent methods produce convincing conversions, but at the cost of increased complexity -- making results difficult to reproduce and build on. Instead, we keep it simple. We propose k-nearest neighbors voice conversion (kNN-VC): a straightforward yet effective method for any-to-any conversion. First, we extract self-supervised representations of the source and reference speech. To convert to the target speaker, we replace each frame of the source representation with its nearest neighbor in the reference. Finally, a pretrained vocoder synthesizes audio from the converted representation. Objective and subjective evaluations show that kNN-VC improves speaker similarity with similar intelligibility scores to existing methods. Code, samples, trained models: https://bshall.github.io/knn-vc △ Less

Submitted 30 May, 2023; originally announced May 2023.

Comments: 5 page, 1 table, 2 figures. Accepted at Interspeech 2023

MIGHTEE-HI: The first MeerKAT HI mass function from an untargeted interferometric survey

Authors: Anastasia A. Ponomareva, Matt J. Jarvis, Hengxing Pan, Natasha Maddox, Michael G. Jones, Bradley S. Frank, Sambatriniaina H. A. Rajohnson, Wanga Mulaudzi, Martin Meyer, Elizabeth A. K. Adams, Maarten Baes, Kelley M. Hess, Sushma Kurapati, Isabella Prandoni, Francesco Sinigaglia, Kristine Spekkens, Madalina Tudorache, Ian Heywood, Jordan D. Collier, Srikrishna Sekhar

Abstract: We present the first measurement of the HI mass function (HIMF) using data from MeerKAT, based on 276 direct detections from the MIGHTEE Survey Early Science data covering a period of approximately a billion years ($0 \leq z \leq 0.084 $). This is the first HIMF measured using interferometric data over non-group or cluster field, i.e. a deep blank field. We constrain the parameters of the Schechte… ▽ More We present the first measurement of the HI mass function (HIMF) using data from MeerKAT, based on 276 direct detections from the MIGHTEE Survey Early Science data covering a period of approximately a billion years ($0 \leq z \leq 0.084 $). This is the first HIMF measured using interferometric data over non-group or cluster field, i.e. a deep blank field. We constrain the parameters of the Schechter function which describes the HIMF with two different methods: $1/\rm V_{\rm max}$ and Modified Maximum Likelihood (MML). We find a low-mass slope $α=-1.29^{+0.37}_{-0.26}$, `knee' mass $\log_{10}(M_{*}/{\rm M_{\odot}}) = 10.07^{+0.24}_{-0.24}$ and normalisation $\log_{10}(φ_{*}/\rm Mpc^{-3})=-2.34^{+0.32}_{-0.36}$ ($H_0 = 67.4$ kms$^{-1}$ Mpc$^{-1}$) for $1/\rm V_{\rm max}$ and $α=-1.44^{+0.13}_{-0.10}$, `knee' mass $\log_{10}(M_{*}/{\rm M_{\odot}}) = 10.22^{+0.10}_{-0.13}$ and normalisation $\log_{10}(φ_{*}/\rm Mpc^{-3})=-2.52^{+0.19}_{-0.14}$ for MML. When using $1/\rm V_{\rm max}$ we find both the low-mass slope and `knee' mass to be consistent within $1σ$ with previous studies based on single-dish surveys. The cosmological mass density of HI is found to be slightly larger than previously reported: $Ω_{\rm HI}=5.46^{+0.94}_{-0.99} \times 10^{-4}h^{-1}_{67.4}$ from $1/\rm V_{\rm max}$ and $Ω_{\rm HI}=6.31^{+0.31}_{-0.31} \times 10^{-4}h^{-1}_{67.4}$ from MML but consistent within the uncertainties. We find no evidence for evolution of the HIMF over the last billion years. △ Less

Submitted 25 April, 2023; originally announced April 2023.

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

X-ray radiative transfer in full 3D with SKIRT

Authors: Bert Vander Meulen, Peter Camps, Marko Stalevski, Maarten Baes

Abstract: Models of active galactic nuclei (AGN) suggest that their circumnuclear media are complex with clumps and filaments, while recent observations hint towards polar extended structures of gas and dust, as opposed to the classical torus paradigm. The X-ray band forms an interesting window to study these circumnuclear media in great detail. In this work, we extended the radiative transfer code SKIRT wi… ▽ More Models of active galactic nuclei (AGN) suggest that their circumnuclear media are complex with clumps and filaments, while recent observations hint towards polar extended structures of gas and dust, as opposed to the classical torus paradigm. The X-ray band forms an interesting window to study these circumnuclear media in great detail. In this work, we extended the radiative transfer code SKIRT with the X-ray processes that govern the broadband X-ray spectra of obscured AGN, to study the structure of AGN circumnuclear media in full 3D, based on their reflected X-ray emission. We extended the SKIRT code with Compton scattering on free electrons, photo-absorption and fluorescence by cold atomic gas, scattering on bound electrons, and extinction by dust. This includes a novel treatment of extreme-forward scattering by dust, and a detailed description of anomalous Rayleigh scattering. To verify our X-ray implementation, we performed the first dedicated benchmark of X-ray torus models, comparing five X-ray radiative transfer codes. Finally, we illustrated the 3D nature of the code by producing synthetic X-ray images and spectra of clumpy torus models. SKIRT forms a powerful new tool to model AGN circumnuclear media in full 3D from X-ray to millimetre wavelengths, and is now publicly available. In the X-ray regime, we find an excellent agreement with the simulation results of the MYTorus and RefleX codes, which validates our X-ray implementation. We find some discrepancies with other codes, which motivates the need for a robust framework that can handle non-linear 3D radiative transfer effects. The new X-ray functionalities of the SKIRT code allow for uncomplicated access to a broad suite of 3D X-ray models for AGN that can easily be tested and modified. This will be particularly useful with the advent of X-ray microcalorimeter observations. △ Less

Submitted 20 April, 2023; originally announced April 2023.

Comments: 20 pages, 20 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 674, A123 (2023)

Simulated observations of star formation regions: infrared evolution of globally collapsing clouds

Authors: Jesús M. Jáquez-Domínguez, Roberto Galván-Madrid, Jacopo Fritz, Manuel Zamora-Avilés, Peter Camps, Gustavo Bruzual, Maarten Baes, Yuxin Lin, Enrique Vázquez-Semadeni

Abstract: The direct comparison between hydrodynamical simulations and observations is needed to improve the physics included in the former and test biases in the latter. Post-processing radiative transfer and synthetic observations are now the standard way to do this. We report on the first application of the \texttt{SKIRT} radiative transfer code to simulations of a star-forming cloud. The synthetic obser… ▽ More The direct comparison between hydrodynamical simulations and observations is needed to improve the physics included in the former and test biases in the latter. Post-processing radiative transfer and synthetic observations are now the standard way to do this. We report on the first application of the \texttt{SKIRT} radiative transfer code to simulations of a star-forming cloud. The synthetic observations are then analyzed following traditional observational workflows. We find that in the early stages of the simulation, stellar radiation is inefficient in heating dust to the temperatures observed in Galactic clouds, thus the addition of an interstellar radiation field is necessary. The spectral energy distribution of the cloud settles rather quickly after $\sim3$ Myr of evolution from the onset of star formation, but its morphology continues to evolve for $\sim8$ Myr due to the expansion of \textsc{Hii} regions and the respective creation of cavities, filaments, and ridges. Modeling synthetic \textit{Herschel} fluxes with 1- or 2-component modified black bodies underestimates total dust masses by a factor of $\sim2$. Spatially-resolved fitting recovers up to about $70\%$ of the intrinsic value. This ``missing mass'' is located in a very cold dust component with temperatures below $10$ K, which does not contribute appreciably to the far-infrared flux. This effect could bias real observations if such dust exists in large amounts. Finally, we tested observational calibrations of the SFR based on infrared fluxes and concluded that they are in agreement when compared to the intrinsic SFR of the simulation averaged over $\sim100$ Myr. △ Less

Submitted 10 April, 2023; originally announced April 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

The atomic-to-molecular hydrogen transition in the TNG50 simulation: Using realistic UV fields to create spatially resolved HI maps

Authors: Andrea Gebek, Maarten Baes, Benedikt Diemer, W. J. G. de Blok, Dylan Nelson, Anand Utsav Kapoor, Peter Camps, Omphile Rabyang, Lerothodi Leeuw

Abstract: Cold gas in galaxies provides a crucial test to evaluate the realism of cosmological hydrodynamical simulations. To extract the atomic and molecular hydrogen properties of the simulated galaxy population, postprocessing methods taking the local UV field into account are required. We improve upon previous studies by calculating realistic UV fields with the dust radiative transfer code SKIRT to mode… ▽ More Cold gas in galaxies provides a crucial test to evaluate the realism of cosmological hydrodynamical simulations. To extract the atomic and molecular hydrogen properties of the simulated galaxy population, postprocessing methods taking the local UV field into account are required. We improve upon previous studies by calculating realistic UV fields with the dust radiative transfer code SKIRT to model the atomic-to-molecular transition in TNG50, the highest-resolution run of the IllustrisTNG suite. Comparing integrated quantities such as the HI mass function, we study to what detail the UV field needs to be modelled in order to calculate realistic cold gas properties. We then evaluate new, spatially resolved comparisons for cold gas in galaxies by exploring synthetic maps of atomic hydrogen at redshift zero and compare them to 21-cm observations of local galaxies from the WHISP survey. In terms of non-parametric morphologies, we find that TNG50 HI maps are less concentrated than their WHISP counterparts (median $ΔC\approx0.3$), due in part to central HI deficits related to the ejective character of supermassive black hole feedback in TNG. In terms of the HI column density distribution function, we find discrepancies between WHISP and IllustrisTNG that depend on the total HI abundance in these datasets as well as the postprocessing method. To fully exploit the synergy between cosmological simulations and upcoming deep HI/H2 data, we advocate the use of accurate methods to estimate the UV radiation field and to generate mock maps. △ Less

Submitted 13 March, 2023; originally announced March 2023.

Comments: 21 pages, 10 figures (main text). Accepted to MNRAS

Nonparametric galaxy morphology from stellar and nebular emission with the CALIFA sample

Authors: Angelos Nersesian, Stefano Zibetti, Francesco D'Eugenio, Maarten Baes

Abstract: We present a nonparametric morphology analysis of the stellar continuum and nebular emission lines for a sample of local galaxies. We explore the dependence of the various morphological parameters on wavelength and morphological type. Our goal is to quantify the difference in morphology between the stellar and nebular components. We derive the nonparametric morphological indicators of 364 galaxies… ▽ More We present a nonparametric morphology analysis of the stellar continuum and nebular emission lines for a sample of local galaxies. We explore the dependence of the various morphological parameters on wavelength and morphological type. Our goal is to quantify the difference in morphology between the stellar and nebular components. We derive the nonparametric morphological indicators of 364 galaxies from the CALIFA Survey. To calculate those indicators, we apply the StatMorph package on the high-quality integral field spectroscopic data cubes, as well as to the most prominent nebular emission-line maps, namely [OIII]$λ$5007, H$α$, and [NII]$λ$6583. We show that the physical size of galaxies, M$_{20}$ index, and concentration have a strong gradient from blue to red optical wavelengths. We find that the light distribution of the nebular emission is less concentrated than the stellar continuum. A comparison between the nonparametric indicators and the galaxy physical properties revealed a very strong correlation of the concentration with the specific star-formation rate and morphological type. Furthermore, we explore how the galaxy inclination affects our results. We find that edge-on galaxies show a more rapid change in physical size and concentration with increasing wavelength due to the increase in optical free path. We conclude that the apparent morphology of galaxies originates from the pure stellar distribution, but the morphology of the ISM presents differences with respect to the morphology of the stellar component. Our analysis also highlights the importance of dust attenuation and galaxy inclination in the measurement of nonparametric morphological indicators, especially in the the wavelength range 4000-5000 Å. △ Less

Submitted 3 March, 2023; originally announced March 2023.

Journal ref: A&A 673, A63 (2023)

Paradigmatic examples for testing models of optical light polarization by spheroidal dust

Authors: C. Peest, R. Siebenmorgen, F. Heymann, T. Vannieuwenhuyse, M. Baes

Abstract: We present a general framework on how the polarization of radiation due to scattering, dichroic extinction, and birefringence of aligned spheroidal dust grains can be implemented and tested in 3D Monte Carlo radiative transfer (MCRT) codes. We derive a methodology for solving the radiative transfer equation governing the changes of the Stokes parameters in dust-enshrouded objects. We utilize the M… ▽ More We present a general framework on how the polarization of radiation due to scattering, dichroic extinction, and birefringence of aligned spheroidal dust grains can be implemented and tested in 3D Monte Carlo radiative transfer (MCRT) codes. We derive a methodology for solving the radiative transfer equation governing the changes of the Stokes parameters in dust-enshrouded objects. We utilize the Müller matrix, and the extinction, scattering, linear, and circular polarization cross sections of spheroidal grains as well as electrons. An established MCRT code is used and its capabilities are extended to include the Stokes formalism. We compute changes in the polarization state of the light by scattering, dichroic extinction, and birefringence on spheroidal grains. The dependency of the optical depth and the albedo on the polarization is treated. The implementation of scattering by spheroidal grains both for random walk steps as well as for directed scattering (peel-off) are described. The observable polarization of radiation of the objects is determined through an angle binning method for photon packages leaving the model space as well as through an inverse ray-tracing routine for the generation of images. We present paradigmatic examples for which we derive analytical solutions of the optical light polarization by spheroidal dust particles. These tests are suited for benchmark verification of MCpol and other such codes, and allow to quantify the numerical precision reached. We demonstrate that MCpol is in excellent agreement to within 0.1% of the Stokes parameters when compared to the analytical solutions. △ Less

Submitted 26 February, 2023; originally announced February 2023.

Comments: A&A accepted by Editor

Journal ref: A&A 673, A112 (2023)

Self-consistent dynamical models with a finite extent -- II. Radially truncated models

Authors: Maarten Baes

Abstract: Galaxies, dark matter haloes, and star clusters have a finite extent, yet most simple dynamical models have an infinite extent. The default method to generate dynamical models with a finite extent is to apply an energy truncation to the distribution function, but this approach is not suited to construct models with a preset density profile and it imposes unphysical constraints on the orbit populat… ▽ More Galaxies, dark matter haloes, and star clusters have a finite extent, yet most simple dynamical models have an infinite extent. The default method to generate dynamical models with a finite extent is to apply an energy truncation to the distribution function, but this approach is not suited to construct models with a preset density profile and it imposes unphysical constraints on the orbit population. We investigate whether it is possible to construct simple dynamical models for spherical systems with a preset density profile with a finite extent, and ideally with a different range of orbital structures. We systematically investigate the consistency of radially truncated dynamical models, and demonstrate that no spherical models with a discontinuous density truncation can be supported by an ergodic orbital structure. On the other hand, we argue that many radially truncated models can be supported by a tangential Osipkov-Merritt orbital structure that becomes completely tangential at the truncation radius. We formulate a consistency hypothesis for radially truncated models with such an orbital structure, and test it using an analytical example and the numerical exploration of a large model parameter space using the SpheCow code. We physically interpret our results in terms of the occupancy of bound orbits, and we discuss possible extensions of the tangential Osipkov-Merritt orbital structure that can support radially truncated models. △ Less

Submitted 10 January, 2023; originally announced January 2023.

Comments: Accepted for publication in MNRAS

Effects of dust sources on dust attenuation properties in IllustrisTNG galaxies at $z\sim 7$

Authors: Yuan-Ming Hsu, Hiroyuki Hirashita, Yen-Hsing Lin, Peter Camps, Maarten Baes

Abstract: Dust emission from high-redshift galaxies gives us a clue to the origin and evolution of dust in the early Universe. Previous studies have shown that different sources of dust (stellar dust production and dust growth in dense clouds) predict different ultraviolet (UV) extinction curves for galaxies at $z\sim 7$ but that the observed attenuation curves depend strongly on the geometry of dust and st… ▽ More Dust emission from high-redshift galaxies gives us a clue to the origin and evolution of dust in the early Universe. Previous studies have shown that different sources of dust (stellar dust production and dust growth in dense clouds) predict different ultraviolet (UV) extinction curves for galaxies at $z\sim 7$ but that the observed attenuation curves depend strongly on the geometry of dust and star distributions. Thus, we perform radiative transfer calculations under the dust-stars geometries computed by a cosmological hydrodynamic simulation (IllustrisTNG). This serves to investigate the dust attenuation curves predicted from `realistic' geometries. We choose objects with stellar mass and star formation rate appropriate for Lyman break galaxies at $z\sim 7$. We find that the attenuation curves are very different from the original extinction curves in most of the galaxies. This makes it difficult to constrain the dominant dust sources from the observed attenuation curves. We further include infrared dust emission in the analysis and plot the infrared excess (IRX)-UV spectral slope ($β$) diagram. We find that different sources of dust cause different IRX-$β$ relations for the simulated galaxies. In particular, if dust growth is the main source of dust, a variation of dust-to-metal ratio causes a more extended sequence with smaller IRX in the IRX-$β$ diagram. Thus, the comprehensive analysis of the abundances of dust and metals, the UV slope, and the dust emission could provide a clue to the dominant dust sources in the Universe. △ Less

Submitted 9 December, 2022; originally announced December 2022.

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

The resolved scaling relations in DustPedia: Zooming in on the local Universe

Authors: Viviana Casasola, Simone Bianchi, Laura Magrini, Aleksandr V. Mosenkov, Francesco Salvestrini, Maarten Baes, Francesco Calura, Letizia P. Cassara', Christopher J. R. Clark, Edvige Corbelli, Jacopo Fritz, Frederic Galliano, Elisabetta Liuzzo, Suzanne Madden, Angelos Nersesian, Francesca Pozzi, Sambit Roychowdhury, Ivano Baronchelli, Matteo Bonato, Carlotta Gruppioni, Lara Pantoni

Abstract: We perform a homogeneous analysis of an unprecedented set of spatially resolved scaling relations (SRs) between ISM components and other properties in the range of scales 0.3-3.4 kpc. We also study some ratios: dust-to-stellar, dust-to-gas, and dust-to-metal. We use a sample of 18 large, spiral, face-on DustPedia galaxies. All the SRs are moderate/strong correlations except the dust-HI SR that doe… ▽ More We perform a homogeneous analysis of an unprecedented set of spatially resolved scaling relations (SRs) between ISM components and other properties in the range of scales 0.3-3.4 kpc. We also study some ratios: dust-to-stellar, dust-to-gas, and dust-to-metal. We use a sample of 18 large, spiral, face-on DustPedia galaxies. All the SRs are moderate/strong correlations except the dust-HI SR that does not exist or is weak for most galaxies. The SRs do not have a universal form but each galaxy is characterized by distinct correlations, affected by local processes and galaxy peculiarities. The SRs hold starting from 0.3 kpc, and if a breaking down scale exists it is < 0.3 kpc. By evaluating all galaxies at 3.4 kpc, differences due to peculiarities of individual galaxies are cancelled out and the corresponding SRs are consistent with those of whole galaxies. By comparing subgalactic and global scales, the most striking result emerges from the SRs involving ISM components: the dust-total gas SR is a good correlation at all scales, while the dust-H2 and dust-HI SRs are good correlations at subkpc/kpc and total scales, respectively. For the other explored SRs, there is a good agreement between small and global scales and this may support the picture where the main physical processes regulating the properties and evolution of galaxies occur locally. Our results are consistent with the hypothesis of self-regulation of the SF process. The analysis of subgalactic ratios shows that they are consistent with those derived for whole galaxies, from low to high z, supporting the idea that also these ratios could be set by local processes. Our results highlight the heterogeneity of galaxy properties and the importance of resolved studies on local galaxies in the context of galaxy evolution. They also provide observational constraints to theoretical models and updated references for high-z studies. △ Less

Submitted 28 October, 2022; originally announced October 2022.

Comments: 42 pages, 11 figures and 5 tables in the main text, 2 figures and 1 table in Appendix. Accepted for publication in A&A

Journal ref: A&A 668, A130 (2022)

TransFusion: Transcribing Speech with Multinomial Diffusion

Authors: Matthew Baas, Kevin Eloff, Herman Kamper

Abstract: Diffusion models have shown exceptional scaling properties in the image synthesis domain, and initial attempts have shown similar benefits for applying diffusion to unconditional text synthesis. Denoising diffusion models attempt to iteratively refine a sampled noise signal until it resembles a coherent signal (such as an image or written sentence). In this work we aim to see whether the benefits… ▽ More Diffusion models have shown exceptional scaling properties in the image synthesis domain, and initial attempts have shown similar benefits for applying diffusion to unconditional text synthesis. Denoising diffusion models attempt to iteratively refine a sampled noise signal until it resembles a coherent signal (such as an image or written sentence). In this work we aim to see whether the benefits of diffusion models can also be realized for speech recognition. To this end, we propose a new way to perform speech recognition using a diffusion model conditioned on pretrained speech features. Specifically, we propose TransFusion: a transcribing diffusion model which iteratively denoises a random character sequence into coherent text corresponding to the transcript of a conditioning utterance. We demonstrate comparable performance to existing high-performing contrastive models on the LibriSpeech speech recognition benchmark. To the best of our knowledge, we are the first to apply denoising diffusion to speech recognition. We also propose new techniques for effectively sampling and decoding multinomial diffusion models. These are required because traditional methods of sampling from acoustic models are not possible with our new discrete diffusion approach. Code and trained models are available: https://github.com/RF5/transfusion-asr △ Less

Submitted 14 October, 2022; originally announced October 2022.

Comments: 12 pages, 4 figures, 1 table. Accepted at SACAIR 2022

GAN You Hear Me? Reclaiming Unconditional Speech Synthesis from Diffusion Models

Authors: Matthew Baas, Herman Kamper

Abstract: We propose AudioStyleGAN (ASGAN), a new generative adversarial network (GAN) for unconditional speech synthesis. As in the StyleGAN family of image synthesis models, ASGAN maps sampled noise to a disentangled latent vector which is then mapped to a sequence of audio features so that signal aliasing is suppressed at every layer. To successfully train ASGAN, we introduce a number of new techniques,… ▽ More We propose AudioStyleGAN (ASGAN), a new generative adversarial network (GAN) for unconditional speech synthesis. As in the StyleGAN family of image synthesis models, ASGAN maps sampled noise to a disentangled latent vector which is then mapped to a sequence of audio features so that signal aliasing is suppressed at every layer. To successfully train ASGAN, we introduce a number of new techniques, including a modification to adaptive discriminator augmentation to probabilistically skip discriminator updates. ASGAN achieves state-of-the-art results in unconditional speech synthesis on the Google Speech Commands dataset. It is also substantially faster than the top-performing diffusion models. Through a design that encourages disentanglement, ASGAN is able to perform voice conversion and speech editing without being explicitly trained to do so. ASGAN demonstrates that GANs are still highly competitive with diffusion models. Code, models, samples: https://github.com/RF5/simple-asgan/. △ Less

Submitted 11 October, 2022; originally announced October 2022.

Comments: 6 pages, 2 figures, 2 tables. Accepted at IEEE SLT 2022

MIGHTEE-HI: The HI mass-stellar mass relation over the last billion years

Authors: Hengxing Pan, Matt J. Jarvis, Mario G. Santos, Natasha Maddox, Bradley S. Frank, Anastasia A. Ponomareva, Isabella Prandoni, Sushma Kurapati, Maarten Baes, Pavel E. Mancera Piña, Giulia Rodighiero, Martin J. Meyer, Romeel Davé, Gauri Sharma, Sambatriniaina H. A. Rajohnson, Nathan J. Adams, Rebecca A. A. Bowler, Francesco Sinigaglia, Thijs van der Hulst, Peter W. Hatfield, Srikrishna Sekhar, Jordan D. Collier

Abstract: We study the $M_{\rm HI}-M_{\star}$ relation over the last billion years using the MIGHTEE-HI sample. We first model the upper envelope of the $M_{\rm HI}-M_{\star}$ relation with a Bayesian technique applied to a total number of 249 HI-selected galaxies, without binning the datasets, while taking account of the intrinsic scatter. We fit the envelope with both linear and non-linear models, and fin… ▽ More We study the $M_{\rm HI}-M_{\star}$ relation over the last billion years using the MIGHTEE-HI sample. We first model the upper envelope of the $M_{\rm HI}-M_{\star}$ relation with a Bayesian technique applied to a total number of 249 HI-selected galaxies, without binning the datasets, while taking account of the intrinsic scatter. We fit the envelope with both linear and non-linear models, and find that the non-linear model is preferred over the linear one with a measured transition stellar mass of $\log_{10}(M_\star$/$M_{\odot})$ = $9.15\pm0.87$, beyond which the slope flattens. This finding supports the view that the lack of HI gas is ultimately responsible for the decreasing star formation rate observed in the massive main-sequence galaxies. For spirals alone, which are biased towards the massive galaxies in our sample, the slope beyond the transition mass is shallower than for the full sample, indicative of distinct gas processes ongoing for the spirals/high-mass galaxies from other types with lower stellar masses. We then create mock catalogues for the MIGHTEE-HI detections and non-detections with two main galaxy populations of late- and early-type galaxies to measure the underlying $M_{\rm HI}-M_{\star}$ relation. We find that the turnover in this relation persists whether considering the two galaxy populations as a whole or separately. We note that an underlying linear relation could mimic this turnover in the observed scaling relation, but a model with a turnover is strongly preferred. Measurements on the logarithmic average of HI masses against the stellar mass are provided as a benchmark for future studies. △ Less

Submitted 27 July, 2023; v1 submitted 10 October, 2022; originally announced October 2022.

Comments: Revised manuscript accepted for publication in MNRAS

The Einasto model for dark matter haloes

Authors: Maarten Baes

Abstract: Context: The Einasto model has become one of the most popular models for describing the density profile of dark matter haloes. There have been relatively few comprehensive studies on the dynamical structure of the Einasto model, mainly because only a limited number of properties can be calculated analytically. Aims: We want to systematically investigate the photometric and dynamical structure of t… ▽ More Context: The Einasto model has become one of the most popular models for describing the density profile of dark matter haloes. There have been relatively few comprehensive studies on the dynamical structure of the Einasto model, mainly because only a limited number of properties can be calculated analytically. Aims: We want to systematically investigate the photometric and dynamical structure of the family of Einasto models over the entire model parameter space. Methods: We used the SpheCow code to explore the properties of the Einasto model. We systematically investigated how the most important properties change as a function of the Einasto index $n$. We considered both isotropic models and radially anisotropic models with an Osipkov-Merritt orbital structure. Results: We find that all Einasto models with $n<\tfrac12$ have a formal isotropic or Osipkov-Merritt distribution function that is negative in parts of phase space, and hence cannot be supported by such orbital structures. On the other hand, all models with larger values of $n$ can be supported by an isotropic orbital structure, or by an Osipkov-Merritt anisotropy, as long as the anisotropy radius is larger than a critical value. This critical anisotropy radius is a decreasing function of $n$, indicating that less centrally concentrated models allow for a larger degree of radial anisotropy. Conclusions: Studies of the structure and dynamics of models for galaxies and dark matter haloes should not be restricted to completely analytical models. Numerical codes such as SpheCow can help open up the range of models that are systematically investigated. This applies to the Einasto model discussed here, but also to other proposed models for dark matter haloes, including different extensions to the Einasto model. △ Less

Submitted 8 September, 2022; originally announced September 2022.

Comments: 11 pages, accepted for publication in A&A

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

UV to submillimetre luminosity functions of TNG50 galaxies

Authors: Ana Trcka, Maarten Baes, Peter Camps, Anand Utsav Kapoor, Dylan Nelson, Annalisa Pillepich, Daniela Barrientos, Lars Hernquist, Federico Marinacci, Mark Vogelsberger

Abstract: We apply the radiative transfer (RT) code SKIRT on a sample of ~14000 low-redshift (z<= 0.1) galaxies extracted from the TNG50 simulation to enable an apples-to-apples comparison with observations. The RT procedure is calibrated via comparison of a subsample of TNG50 galaxies with the DustPedia observational sample: we compare several luminosity and colour scaling relations and spectral energy dis… ▽ More We apply the radiative transfer (RT) code SKIRT on a sample of ~14000 low-redshift (z<= 0.1) galaxies extracted from the TNG50 simulation to enable an apples-to-apples comparison with observations. The RT procedure is calibrated via comparison of a subsample of TNG50 galaxies with the DustPedia observational sample: we compare several luminosity and colour scaling relations and spectral energy distributions in different specific SFR bins. We consistently derive galaxy luminosity functions for the TNG50 simulation in 14 broadband filters from UV to submillimetre wavelengths and investigate the effects of the aperture, orientation, radiative transfer recipe, and numerical resolution. We find that, while our TNG50+RT fiducial model agrees well with the observed luminosity functions at the knee (+/- 0.04 dex typical agreement), the TNG50+RT luminosity functions evaluated within 5R_1/2 are generally higher than observed at both the faint and bright ends, by 0.004 (total IR)-0.27 (UKIDSS H) dex and 0.12 (SPIRE250)-0.8 (GALEX FUV) dex, respectively. A change in the aperture does affect the bright end of the luminosity function, easily by up to 1 dex depending on the choice. However, we also find that the galaxy luminosity functions of a worse-resolution run of TNG50 (TNG50-2, with 8 times worse mass resolution than TNG50, similar to TNG100) are in better quantitative agreement with observational constraints. Finally, we publicly release the photometry for the TNG50 sample in 53 broadbands from FUV to submillimetre, in three orientations and four apertures, as well as galaxy spectral energy distributions. △ Less

Submitted 12 August, 2022; originally announced August 2022.

Comments: The full galaxy SEDs and the broadband fluxes at z <= 0.1 are made available at https://www.tng- project.org/trcka22 The luminosity function plots are presented on http://luminosity-functions.herokuapp.com/

Monte Carlo radiative transfer with explicit absorption to simulate absorption, scattering, and stimulated emission

Authors: Maarten Baes, Peter Camps, Kosei Matsumoto

Abstract: Context: The Monte Carlo method is probably the most widely used approach to solve the radiative transfer problem, especially in a general 3D geometry. The physical processes of emission, absorption, and scattering are easily incorporated in the Monte Carlo framework. Net stimulated emission, or absorption with a negative cross section, does not fit this method, however. Aims: We explore alteratio… ▽ More Context: The Monte Carlo method is probably the most widely used approach to solve the radiative transfer problem, especially in a general 3D geometry. The physical processes of emission, absorption, and scattering are easily incorporated in the Monte Carlo framework. Net stimulated emission, or absorption with a negative cross section, does not fit this method, however. Aims: We explore alterations to the standard photon packet life cycle in Monte Carlo radiative transfer that allow the treatment of net stimulated emission without loss of generality or efficiency. Methods: We present the explicit absorption technique that allows net stimulated emission to be handled efficiently. It uses the scattering rather than the extinction optical depth along a photon packet's path to randomly select the next interaction location, and offers a separate, deterministic treatment of absorption. We implemented the technique in a special-purpose Monte Carlo code for a two-stream 1D radiative transfer problem and in the fully featured 3D code SKIRT, and we studied its overall performance using quantitative statistical tests. Results: Our special-purpose code is capable of recovering the analytical solutions to the two-stream problem in all regimes, including the one of strong net stimulated emission. The implementation in SKIRT is straightforward, as the explicit absorption technique easily combines with the variance reduction and acceleration techniques already incorporated. In general, explicit absorption tends to improve the efficiency of the Monte Carlo routine in the regime of net absorption. Conclusions: Explicit absorption allows the treatment of net stimulated emission in Monte Carlo radiative transfer, it interfaces smoothly with... (abridged) △ Less

Submitted 12 August, 2022; originally announced August 2022.

Comments: 14 pages, accepted for publication in A&A

Journal ref: A&A 666, A101 (2022)

MIGHTEE-HI: Evolution of HI scaling relations of star-forming galaxies at $z<0.5$

Authors: Francesco Sinigaglia, Giulia Rodighiero, Ed Elson, Mattia Vaccari, Natasha Maddox, Bradley S. Frank, Matt J. Jarvis, Tom Oosterloo, Romeel Davé, Mara Salvato, Maarten Baes, Sabine Bellstedt, Laura Bisigello, Jordan D. Collier, Robin H. W. Cook, Luke J. M. Davies, Jacinta Delhaize, Simon P. Driver, Caroline Foster, Sushma Kurapati, Claudia del P. Lagos, Christopher Lidman, Pavel E. Mancera Piña, Martin J. Meyer, K. Moses Mogotsi , et al. (11 additional authors not shown)

Abstract: We present the first measurements of HI galaxy scaling relations from a blind survey at $z>0.15$. We perform spectral stacking of 9023 spectra of star-forming galaxies undetected in HI at $0.23<z<0.49$, extracted from MIGHTEE-HI Early Science datacubes, acquired with the MeerKAT radio telescope. We stack galaxies in bins of galaxy properties ($M_*$, SFR, and sSFR, with… ▽ More We present the first measurements of HI galaxy scaling relations from a blind survey at $z>0.15$. We perform spectral stacking of 9023 spectra of star-forming galaxies undetected in HI at $0.23<z<0.49$, extracted from MIGHTEE-HI Early Science datacubes, acquired with the MeerKAT radio telescope. We stack galaxies in bins of galaxy properties ($M_*$, SFR, and sSFR, with ${\rm sSFR}\equiv M_*/{\rm SFR}$), obtaining $\gtrsim 5σ$ detections in most cases, the strongest HI-stacking detections to date in this redshift range. With these detections, we are able to measure scaling relations in the probed redshift interval, finding evidence for a moderate evolution from the median redshift of our sample $z_{\rm med}\sim 0.37$ to $z\sim 0$. In particular, low-$M_*$ galaxies ($\log_{10}(M_*/{\rm M_\odot})\sim 9$) experience a strong HI depletion ($\sim 0.5$ dex in $\log_{10}(M_{\rm HI}/{\rm M}_\odot)$), while massive galaxies ($\log_{10}(M_*/{\rm M_\odot})\sim 11$) keep their HI mass nearly unchanged. When looking at the star formation activity, highly star-forming galaxies evolve significantly in $M_{\rm HI}$ ($f_{\rm HI}$, where $f_{\rm HI}\equiv M_{\rm}/M_*$) at fixed SFR (sSFR), while at the lowest probed SFR (sSFR) the scaling relations show no evolution. These findings suggest a scenario in which low-$M_*$ galaxies have experienced a strong HI depletion during the last $\sim4$ Gyr, while massive galaxies have undergone a significant HI replenishment through some accretion mechanism, possibly minor mergers. Interestingly, our results are in good agreement with the predictions of the SIMBA simulation. We conclude that this work sets novel important observational constraints on galaxy scaling relations. △ Less

Submitted 1 August, 2022; originally announced August 2022.

Comments: 12 pages, 5 figures. Accepted for publication in ApJ Letters on 30 July 2022

The distribution of dust in edge-on galaxies: I. The global structure

Authors: Aleksandr V. Mosenkov, Pavel A. Usachev, Zacory Shakespear, Jacob Guerrette, Maarten Baes, Simone Bianchi, Emmanuel M. Xilouris, George A. Gontcharov, Vladimir B. Il'in, Alexander A. Marchuk, Sergey S. Savchenko, Anton A. Smirnov

Abstract: In this first paper in a series we present a study of the global dust emission distribution in nearby edge-on spiral galaxies. Our sample consists of 16 angularly large and 13 less spatially resolved galaxies selected from the DustPedia sample. To explore the dust emission distribution, we exploit the Herschel photometry in the range 100-500 $μ$m. We employ Sérsic and three-dimensional disc models… ▽ More In this first paper in a series we present a study of the global dust emission distribution in nearby edge-on spiral galaxies. Our sample consists of 16 angularly large and 13 less spatially resolved galaxies selected from the DustPedia sample. To explore the dust emission distribution, we exploit the Herschel photometry in the range 100-500 $μ$m. We employ Sérsic and three-dimensional disc models to fit the observed two-dimensional profiles of the galaxies. Both approaches give similar results. Our analysis unequivocally states the case for the presence of extraplanar dust in between 6 to 10 large galaxies. The results reveal that both the disc scale length and height increase as a function of wavelength between 100 and 500 $μ$m. The dust disc scale height positively correlates with the dust disc scale length, similar to what is observed for the stellar discs. We also find correlations between the scale lengths and scale heights in the near- and far-infrared which suggest that the stellar discs and their dust counterparts are tightly connected. Furthermore, the intrinsic flattening of the dust disc is inversely proportional to the maximum rotation velocity and the dust mass of the galaxy: more massive spiral galaxies host, on average, relatively thinner dust discs. Also, there is a tendency for the dust-to-stellar scale height ratio to decrease with the dust mass and rotation velocity. We conclude that low-mass spiral galaxies host a diffuse, puffed-up dust disc with a thickness similar to that of the stellar disc. △ Less

Submitted 23 July, 2022; originally announced July 2022.

Comments: 22 pages, 14 figures, accepted for publication in MNRAS

Euclid preparation: XXIII. Derivation of galaxy physical properties with deep machine learning using mock fluxes and H-band images

Authors: Euclid Collaboration, L. Bisigello, C. J. Conselice, M. Baes, M. Bolzonella, M. Brescia, S. Cavuoti, O. Cucciati, A. Humphrey, L. K. Hunt, C. Maraston, L. Pozzetti, C. Tortora, S. E. van Mierlo, N. Aghanim, N. Auricchio, M. Baldi, R. Bender, C. Bodendorf, D. Bonino, E. Branchini, J. Brinchmann, S. Camera, V. Capobianco, C. Carbone , et al. (174 additional authors not shown)

Abstract: Next generation telescopes, like Euclid, Rubin/LSST, and Roman, will open new windows on the Universe, allowing us to infer physical properties for tens of millions of galaxies. Machine learning methods are increasingly becoming the most efficient tools to handle this enormous amount of data, because they are often faster and more accurate than traditional methods. We investigate how well redshift… ▽ More Next generation telescopes, like Euclid, Rubin/LSST, and Roman, will open new windows on the Universe, allowing us to infer physical properties for tens of millions of galaxies. Machine learning methods are increasingly becoming the most efficient tools to handle this enormous amount of data, because they are often faster and more accurate than traditional methods. We investigate how well redshifts, stellar masses, and star-formation rates (SFR) can be measured with deep learning algorithms for observed galaxies within data mimicking the Euclid and Rubin/LSST surveys. We find that Deep Learning Neural Networks and Convolutional Neutral Networks (CNN), which are dependent on the parameter space of the training sample, perform well in measuring the properties of these galaxies and have a better accuracy than methods based on spectral energy distribution fitting. CNNs allow the processing of multi-band magnitudes together with $H_{\scriptscriptstyle\rm E}$-band images. We find that the estimates of stellar masses improve with the use of an image, but those of redshift and SFR do not. Our best results are deriving i) the redshift within a normalised error of less than 0.15 for 99.9$\%$ of the galaxies with S/N>3 in the $H_{\scriptscriptstyle\rm E}$-band; ii) the stellar mass within a factor of two ($\sim0.3 \rm dex$) for 99.5$\%$ of the considered galaxies; iii) the SFR within a factor of two ($\sim0.3 \rm dex$) for $\sim$70$\%$ of the sample. We discuss the implications of our work for application to surveys as well as how measurements of these galaxy parameters can be improved with deep learning. △ Less

Submitted 4 January, 2023; v1 submitted 29 June, 2022; originally announced June 2022.

Comments: accepted for publication in MNRAS, 21 pages, 22 figures, 6 tables

Dust emissivity in resolved spiral galaxies

Authors: S. Bianchi, V. Casasola, E. Corbelli, F. Galliano, L. Magrini, A. Nersesian, F. Salvestrini, M. Baes, L. P. Cassara, C. J. R. Clark, I. De Looze, A. P. Jones, S. C. Madden, A. Mosenkov, N. Ysard

Abstract: Context: The far-infrared (FIR) and sub-millimeter (submm) emissivity of the Milky Way (MW) cirrus is an important benchmark for dust grain models. Dust masses in other galaxies are generally derived from the FIR/submm using the emission properties of these MW-calibrated models. Aims: We seek to derive the FIR/submm emissivity in nine nearby spiral galaxies to check its compatibility with MW cirru… ▽ More Context: The far-infrared (FIR) and sub-millimeter (submm) emissivity of the Milky Way (MW) cirrus is an important benchmark for dust grain models. Dust masses in other galaxies are generally derived from the FIR/submm using the emission properties of these MW-calibrated models. Aims: We seek to derive the FIR/submm emissivity in nine nearby spiral galaxies to check its compatibility with MW cirrus measurements. Methods: We obtained values of the emissivity at 70 to 500 um, using maps of dust emission from the Herschel satellite and of gas surface density from the THINGS and HERACLES surveys on a scale generally corresponding to 440 pc. We studied the variation of the emissivity with the surface brightness ratio I(250um)/I(500um), a proxy for the intensity of the interstellar radiation field heating the dust. Results: We find that the average value of the emissivity agrees with MW estimates for pixels sharing the same color as the cirrus, namely, for I(250um)/I(500um) = 4.5. For I(250um)/I(500um) > 5, the measured emissivity is instead up to a factor ~2 lower than predicted from MW dust models heated by stronger radiation fields. Regions with higher I(250um)/I(500um) are preferentially closer to the galactic center and have a higher overall (stellar+gas) surface density and molecular fraction. The results do not depend strongly on the adopted CO-to-molecular conversion factor and do not appear to be affected by the mixing of heating conditions. Conclusions: Our results confirm the validity of MW dust models at low density, but are at odds with predictions for grain evolution in higher density environments. If the lower-than-expected emissivity at high I(250um)/I(500um) is the result of intrinsic variations in the dust properties, it would imply an underestimation of the dust mass surface density of up to a factor ~2 when using current dust models. △ Less

Submitted 27 June, 2022; originally announced June 2022.

Comments: 16 pages, 11 figures, A&A accepted

Journal ref: A&A 664, A187 (2022)

Putting GPT-3's Creativity to the (Alternative Uses) Test

Authors: Claire Stevenson, Iris Smal, Matthijs Baas, Raoul Grasman, Han van der Maas

Abstract: AI large language models have (co-)produced amazing written works from newspaper articles to novels and poetry. These works meet the standards of the standard definition of creativity: being original and useful, and sometimes even the additional element of surprise. But can a large language model designed to predict the next text fragment provide creative, out-of-the-box, responses that still solv… ▽ More AI large language models have (co-)produced amazing written works from newspaper articles to novels and poetry. These works meet the standards of the standard definition of creativity: being original and useful, and sometimes even the additional element of surprise. But can a large language model designed to predict the next text fragment provide creative, out-of-the-box, responses that still solve the problem at hand? We put Open AI's generative natural language model, GPT-3, to the test. Can it provide creative solutions to one of the most commonly used tests in creativity research? We assessed GPT-3's creativity on Guilford's Alternative Uses Test and compared its performance to previously collected human responses on expert ratings of originality, usefulness and surprise of responses, flexibility of each set of ideas as well as an automated method to measure creativity based on the semantic distance between a response and the AUT object in question. Our results show that -- on the whole -- humans currently outperform GPT-3 when it comes to creative output. But, we believe it is only a matter of time before GPT-3 catches up on this particular task. We discuss what this work reveals about human and AI creativity, creativity testing and our definition of creativity. △ Less

Submitted 10 June, 2022; originally announced June 2022.

Comments: 5 pages, 6 figures, accepted at the International Conference on Computational Creativity (ICCC) 2022 as a Short Paper. See https://osf.io/vmk3c/ for data, analyses and code