Showing 1–50 of 132 results for author: Malek, k

Methodology to Deploy CNN-Based Computer Vision Models on Immersive Wearable Devices

Authors: Kaveh Malek, Fernando Moreu

Abstract: Convolutional Neural Network (CNN) models often lack the ability to incorporate human input, which can be addressed by Augmented Reality (AR) headsets. However, current AR headsets face limitations in processing power, which has prevented researchers from performing real-time, complex image recognition tasks using CNNs in AR headsets. This paper presents a method to deploy CNN models on AR headset… ▽ More Convolutional Neural Network (CNN) models often lack the ability to incorporate human input, which can be addressed by Augmented Reality (AR) headsets. However, current AR headsets face limitations in processing power, which has prevented researchers from performing real-time, complex image recognition tasks using CNNs in AR headsets. This paper presents a method to deploy CNN models on AR headsets by training them on computers and transferring the optimized weight matrices to the headset. The approach transforms the image data and CNN layers into a one-dimensional format suitable for the AR platform. We demonstrate this method by training the LeNet-5 CNN model on the MNIST dataset using PyTorch and deploying it on a HoloLens AR headset. The results show that the model maintains an accuracy of approximately 98%, similar to its performance on a computer. This integration of CNN and AR enables real-time image processing on AR headsets, allowing for the incorporation of human input into AI models. △ Less

Submitted 28 June, 2024; originally announced July 2024.

Comments: 10 pages 8 figures 4300 words

Immersive Robot Programming Interface for Human-Guided Automation and Randomized Path Planning

Authors: Kaveh Malek, Claus Danielson, Fernando Moreu

Abstract: Researchers are exploring Augmented Reality (AR) interfaces for online robot programming to streamline automation and user interaction in variable manufacturing environments. This study introduces an AR interface for online programming and data visualization that integrates the human in the randomized robot path planning, reducing the inherent randomness of the methods with human intervention. The… ▽ More Researchers are exploring Augmented Reality (AR) interfaces for online robot programming to streamline automation and user interaction in variable manufacturing environments. This study introduces an AR interface for online programming and data visualization that integrates the human in the randomized robot path planning, reducing the inherent randomness of the methods with human intervention. The interface uses holographic items which correspond to physical elements to interact with a redundant manipulator. Utilizing Rapidly Random Tree Star (RRT*) and Spherical Linear Interpolation (SLERP) algorithms, the interface achieves end-effector s progression through collision-free path with smooth rotation. Next, Sequential Quadratic Programming (SQP) achieve robot s configurations for this progression. The platform executes the RRT* algorithm in a loop, with each iteration independently exploring the shortest path through random sampling, leading to variations in the optimized paths produced. These paths are then demonstrated to AR users, who select the most appropriate path based on the environmental context and their intuition. The accuracy and effectiveness of the interface are validated through its implementation and testing with a seven Degree-OF-Freedom (DOF) manipulator, indicating its potential to advance current practices in robot programming. The validation of this paper include two implementations demonstrating the value of human-in-the-loop and context awareness in robotics. △ Less

Submitted 4 June, 2024; originally announced June 2024.

Comments: 10 pages, 13 figures

MSC Class: 90; 37; 70 ACM Class: I.2.9; H.5.2; I.2.10

Novel Data Models for Inter-operable LCA Frameworks

Authors: Kourosh Malek, Max Dreger, Zirui Tang, Qingshi Tu

Abstract: Life cycle assessment (LCA) plays a critical role in assessing the environmental impacts of a product, technology, or service throughout its entire life cycle. Nonetheless, many existing LCA tools and methods lack adequate metadata management, which can hinder their further development and wide adoption. In the example of LCA for clean energy technologies, metadata helps monitor data and the envir… ▽ More Life cycle assessment (LCA) plays a critical role in assessing the environmental impacts of a product, technology, or service throughout its entire life cycle. Nonetheless, many existing LCA tools and methods lack adequate metadata management, which can hinder their further development and wide adoption. In the example of LCA for clean energy technologies, metadata helps monitor data and the environment that holds the integrity of the energy assets and sustainability of the materials sources across their entire value chains. Ontologizing metadata, i.e. a common vocabulary and language to connect multiple data sources, as well as implementing AI-aware data management, can have long-lasting, positive, and accelerating effects along with collecting and utilizing quality data from different sources and across the entire data lifecycle. The integration of ontologies in life cycle assessments has garnered significant attention in recent years. We synthesized the existing literature on ontologies for LCAs, providing insights into this interdisciplinary field's evolution, current state, and future directions. We also proposed the framework for a suitable data model and the workflow thereof to warrant the alignment with existing ontologies, practical frameworks, and industry standards. △ Less

Submitted 16 May, 2024; originally announced May 2024.

The LOFAR Two-metre Sky Survey: The nature of the faint source population and SFR-radio luminosity relation using Prospector

Authors: Soumyadeep Das, Daniel J. B. Smith, Paul Haskell, Martin J. Hardcastle, Philip N. Best, Kenneth J. Duncan, Marina I. Arnaudova, Shravya Shenoy, Rohit Kondapally, Rachel K. Cochrane, Alyssa B. Drake, Gülay Gürkan, Katarzyna Małek, Leah K. Morabito, Isabella Prandoni

Abstract: Spectral energy distribution (SED) fitting has been extensively used to determine the nature of the faint radio source population. Recent efforts have combined fits from multiple SED-fitting codes to account for the host galaxy and any active nucleus that may be present. We show that it is possible to produce similar-quality classifications using a single energy-balance SED fitting code, Prospecto… ▽ More Spectral energy distribution (SED) fitting has been extensively used to determine the nature of the faint radio source population. Recent efforts have combined fits from multiple SED-fitting codes to account for the host galaxy and any active nucleus that may be present. We show that it is possible to produce similar-quality classifications using a single energy-balance SED fitting code, Prospector, to model up to 26 bands of UV$-$far-infrared aperture-matched photometry for $\sim$31,000 sources in the ELAIS-N1 field from the LOFAR Two-Metre Sky Survey (LoTSS) Deep fields first data release. One of a new generation of SED-fitting codes, Prospector accounts for potential contributions from radiative active galactic nuclei (AGN) when estimating galaxy properties, including star formation rates (SFRs) derived using non-parametric star formation histories. Combining this information with radio luminosities, we classify 92 per cent of the radio sources as a star-forming galaxy, high/low-excitation radio galaxy, or radio-quiet AGN and study the population demographics as a function of 150 MHz flux density, luminosity, SFR, stellar mass, redshift and apparent $r$-band magnitude. Finally, we use Prospector SED fits to investigate the SFR$-$150 MHz luminosity relation for a sample of $\sim$$133,000~3.6~μ$m-selected $z<1$ sources, finding that the stellar mass dependence is significantly weaker than previously reported, and may disappear altogether at $\log_{10} (\mathrm{SFR}/M_\odot~\mathrm{yr}^{-1}) > 0.5$. This approach makes it significantly easier to classify radio sources from LoTSS and elsewhere, and may have important implications for future studies of star-forming galaxies at radio wavelengths. △ Less

Submitted 2 May, 2024; originally announced May 2024.

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

Tracing the evolutionary pathways of dust and cold gas in high-z quiescent galaxies with SIMBA

Authors: G. Lorenzon, D. Donevski, K. Lisiecki, C. Lovell, M. Romano, D. Narayanan, R. Davé, A. Man, K. E. Whitaker, A. Nanni, A. Long, M. M. Lee, Junais, K. Małek, G. Rodighiero, Q. Li

Abstract: Recent discoveries of copious amounts of dust in quiescent galaxies (QGs) at high redshifts ($z\gtrsim 1-2$) challenge the conventional view that these objects have poor interstellar medium (ISM) in proportion to their stellar mass. We use the SIMBA cosmological simulation to explore the evolution of dust and cold gas content in QGs in relation to the quenching processes affecting them. We track t… ▽ More Recent discoveries of copious amounts of dust in quiescent galaxies (QGs) at high redshifts ($z\gtrsim 1-2$) challenge the conventional view that these objects have poor interstellar medium (ISM) in proportion to their stellar mass. We use the SIMBA cosmological simulation to explore the evolution of dust and cold gas content in QGs in relation to the quenching processes affecting them. We track the changes in the ISM dust abundance across the evolutionary history of QGs identified at $0 \lesssim z \lesssim2$ in the field and cluster environments. The QGs quench via diverse pathways, both rapid and slow, and exhibit a wide range of times elapsed between the quenching event and cold gas removal (from $\sim650$ Myr to $\sim8$ Gyr). We find that quenching modes attributed to the feedback from active galactic nuclei (AGN) do not affect dust and cold gas within the same timescales. Remarkably, QGs may replenish their dust content in the quenched phase primarily due to internal processes and marginally by external factors such as minor mergers. The key mechanism for re-formation of dust is prolonged grain growth on gas-phase metals, it is effective within $\sim100$ Myr after the quenching event, and rapidly increases the dust-to-gas mass ratio in QGs above the standard values ($δ_{\rm DGR}\gtrsim1/100$). As a result, despite heavily depleted cold gas reservoirs, roughly half of QGs maintain little evolution in their ISM dust with stellar age within the first 2 Gyr following the quenching. Overall, we predict that relatively dusty QGs ($M_{\rm dust}/M_{\star}\gtrsim10^{-3}-10^{-4}$) arise from both fast and slow quenchers, and are prevalent in systems of intermediate and low stellar masses ($9<\log(M_{\star}/M_{\odot})<10.5$). This prediction poses an immediate quest for observational synergy between e.g., James Webb Space Telescope (JWST) and the Atacama Large Millimeter Array (ALMA). △ Less

Submitted 16 April, 2024; originally announced April 2024.

Comments: 21 pages, 10 figures, 1 table, submitted to A&A

Effects of galaxy environment on merger fraction

Authors: W. J. Pearson, D. J. D. Santos, T. Goto, T. -C. Huang, S. J. Kim, H. Matsuhara, A. Pollo, S. C. -C. Ho, H. S. Hwang, K. Małek, T. Nakagawa, M. Romano, S. Serjeant, L. Suelves, H. Shim, G. J. White

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

Submitted 18 March, 2024; originally announced March 2024.

Comments: 15 pages, 10 figures, 8 tables, 2 appendices, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 686, A94 (2024)

Radio-only and Radio-to-far-ultraviolet Spectral Energy Distribution Modeling of 14 ULIRGs: Insights into the Global Properties of Infrared Bright Galaxies

Authors: Subhrata Dey, Arti Goyal, Katarzyna Małek, Tanio Díaz-Santos

Abstract: We present detailed spectral energy distribution (SED) modeling of 14 local ultraluminous infrared galaxies (ULIRGs) with outstanding photometric data from the literature covering the ultraviolet--infrared (FIR) and radio bands ($\sim$50 MHz to $\sim$30 GHz). We employ the CIGALE SED fitting code to model the ultraviolet--FIR--radio SED. For the radio-only SED modeling, we use the UltraNest packag… ▽ More We present detailed spectral energy distribution (SED) modeling of 14 local ultraluminous infrared galaxies (ULIRGs) with outstanding photometric data from the literature covering the ultraviolet--infrared (FIR) and radio bands ($\sim$50 MHz to $\sim$30 GHz). We employ the CIGALE SED fitting code to model the ultraviolet--FIR--radio SED. For the radio-only SED modeling, we use the UltraNest package, leveraging its nested sampling algorithm. Combining the results from our previous study on 11 luminous infrared galaxies (LIRGs), we discuss the global astrophysical properties of a sample of 25 starburst galaxies (z $<$ 0.5). Their radio spectra are frequently characterized by bends and turnovers, with no indication of ULIRGs exhibiting more complicated SEDs than LIRGs despite showing more signs of interactions. Including radio measurements in the CIGALE modeling constrained the dust luminosity and star formation rate (SFR) estimates by more than 1 order of magnitude better than previously reported for starburst galaxies. We show that total and nonthermal radio luminosity at 1.4 and 4.8 GHz frequencies can be good estimators of recent SFRs for all LIRGs and those ULIRGS with an insignificant influence of active galactic nuclei. A weaker but still significant correlation is observed between radio SFRs at 1.4 GHz and old (averaged over 100 Myr) SFRs based on SED modeling, indicative of multiple episodes of starburst activity during their lifetime. The thermal radio luminosity at 4.8 GHz is a better tracer of recent star formation than the thermal luminosity at 1.4 GHz. Statistically, our modeled nonthermal radio spectral indices do not significantly correlate with redshift, stellar mass, SFR, specific SFR, and dust mass. △ Less

Submitted 7 May, 2024; v1 submitted 16 February, 2024; originally announced February 2024.

Comments: Published in the ApJ

Journal ref: ApJ 966 61 (2024)

Attenuation proxy hidden in surface brightness-colour diagrams. A new strategy for the LSST era

Authors: K. Małek, Junais, A. Pollo, M. Boquien, V. Buat, S. Salim, S. Brough, R. Demarco, A. W. Graham, M. Hamed, J. R. Mullaney, M. Romano, C. Sifón, M. Aravena, J. A. Benavides, I. Busà, D. Donevski, O. Dorey, H. M. Hernandez-Toledo, A. Nanni, W. J. Pearson, F. Pistis, R. Ragusa, G. Riccio, J. Román

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

Submitted 1 February, 2024; v1 submitted 23 January, 2024; originally announced January 2024.

Comments: 18 pages, accepted for publication in Astronomy and Astrophysics

A comparative study of the fundamental metallicity relation: the impact of methodology on its observed evolution

Authors: F. Pistis, A. Pollo, M. Figueira, D. Vergani, M. Hamed, K. Małek, A. Durkalec, D. Donevski, S. Salim, A. Iovino, W. J. Pearson, M. Romano, M. Scodeggio

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

Submitted 1 December, 2023; originally announced December 2023.

Comments: Accepted for publication in A&A, 22 pages, 20 figures

Shedding Light on Low Surface Brightness Galaxies in Dark Energy Survey with Transformers

Authors: H. Thuruthipilly, Junais, A. Pollo, U. Sureshkumar, M. Grespan, P. Sawant, K. Malek, A. Zadrozny

Abstract: Low surface brightness galaxies (LSBGs) which are defined as galaxies that are fainter than the night sky, play a crucial role in understanding galaxy evolution and cosmological models. Upcoming large-scale surveys like Rubin Observatory Legacy Survey of Space and Time (LSST) and Euclid are expected to observe billions of astronomical objects. In this context, using semi-automatic methods to ident… ▽ More Low surface brightness galaxies (LSBGs) which are defined as galaxies that are fainter than the night sky, play a crucial role in understanding galaxy evolution and cosmological models. Upcoming large-scale surveys like Rubin Observatory Legacy Survey of Space and Time (LSST) and Euclid are expected to observe billions of astronomical objects. In this context, using semi-automatic methods to identify LSBGs would be a highly challenging and time-consuming process and demand automated or machine learning-based methods to overcome this challenge. We study the use of transformer models in separating LSBGs from artefacts in the data from the Dark Energy Survey (DES) data release 1. Using the transformer models, we then search for new LSBGs from the DES that the previous searches may have missed. Properties of the newly found LSBGs are investigated, along with an analysis of the properties of the total LSBG sample in DES. We identified 4,083 new LSBGs in DES, adding an additional $\sim17\% $ to the LSBGs already known in DES. This also increased the number density of LSBGs in DES to 5.5 deg$^{-2}$. We performed a clustering analysis of the LSBGs in DES using an angular two-point auto-correlation function and found that LSBGs cluster more strongly than their high surface brightness counterparts. We associated 1310 LSBGs with galaxy clusters and identified 317 among them as ultra-diffuse galaxies (UDGs). We found that these cluster LSBGs are getting bluer and larger in size towards the edge of the clusters when compared with those in the centre. Transformer models have the potential to be on par with convolutional neural networks as state-of-the-art algorithms in analysing astronomical data. △ Less

Submitted 20 October, 2023; originally announced October 2023.

Comments: 23 pages, 24 Figures and 4 Tables. Accepted for publication in A&A

MUSE observations of the giant low surface brightness galaxy Malin 1: Numerous HII regions, star formation rate, metallicity, and dust attenuation

Authors: Junais, P. M. Weilbacher, B. Epinat, S. Boissier, G. Galaz, E. J. Johnston, T. H. Puzia, P. Amram, K. Małek

Abstract: Giant low-surface brightness (GLSB) galaxies are an extreme class of objects with very faint and extended gas-rich disks. Malin 1 is the largest GLSB galaxy known to date, but its formation is still poorly understood. We use VLT/MUSE IFU spectroscopic observations of Malin 1 to reveal, for the first time, the presence of H$α$ emission distributed across numerous regions along its disk, up to radia… ▽ More Giant low-surface brightness (GLSB) galaxies are an extreme class of objects with very faint and extended gas-rich disks. Malin 1 is the largest GLSB galaxy known to date, but its formation is still poorly understood. We use VLT/MUSE IFU spectroscopic observations of Malin 1 to reveal, for the first time, the presence of H$α$ emission distributed across numerous regions along its disk, up to radial distances of $\sim$100 kpc. We made an estimate of the dust attenuation using the Balmer decrement and found that Malin 1 has a mean H$α$ attenuation of 0.36 mag. We observe a steep decline in the star formation rate surface density ($Σ_{\rm SFR}$) within the inner 20 kpc, followed by a shallow decline in the extended disk. Similarly, the gas phase metallicity we estimated shows a steep gradient in the inner 20 kpc, followed by a flattening of the metallicity in the extended disk with a relatively high value of $\sim$0.6 $Z_{\odot}$. We found that the normalized abundance gradient of the inner disk is similar to values found in normal galaxies but with an extreme value in the extended disk. A comparison of the star formation rate surface density and gas surface density shows that, unlike normal disk galaxies or other LSBs, Malin 1 exhibits a very low star formation efficiency. Owing to the detection of emission lines over a large part of the disk of Malin 1, this work sheds light on the star formation processes in this unique galaxy, highlighting its extended star-forming disk, dust attenuation, almost flat metallicity distribution in the outer disk, and exceptionally low star-formation efficiency. Our findings contribute to a more detailed understanding of the formation of the giant disk of Malin 1 and also constrain possible proposed scenarios on the nature of GLSB galaxies in general. △ Less

Submitted 18 October, 2023; originally announced October 2023.

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

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

Decoding the IRX-β dust attenuation relation in star-forming galaxies at intermediate redshift

Authors: M. Hamed, F. Pistis, M. Figueira, K. Małek, A. Nanni, V. Buat, A. Pollo, D. Vergani, M. Bolzonella, Junais, J. Krywult, T. Takeuchi, G. Riccio, T. Moutard

Abstract: We aim to understand what drives the IRX-βdust attenuation relation at intermediate redshift (0.5 < z < 0.8) in star-forming galaxies. We investigate the role of various galaxy properties in shaping this observed relation. We use robust [O ii] λ3727, [O iii] λλ4959, 5007, and Hβline detections of our statistical sample of 1049 galaxies to estimate the gas-phase metallicities. We derive key physica… ▽ More We aim to understand what drives the IRX-βdust attenuation relation at intermediate redshift (0.5 < z < 0.8) in star-forming galaxies. We investigate the role of various galaxy properties in shaping this observed relation. We use robust [O ii] λ3727, [O iii] λλ4959, 5007, and Hβline detections of our statistical sample of 1049 galaxies to estimate the gas-phase metallicities. We derive key physical properties that are necessary to study galaxy evolution, such as the stellar masses and the star formation rates, using the spectral energy distribution fitting tool CIGALE. Equivalently, we study the effect of galaxy morphology (mainly the Sérsic index n and galaxy inclination) on the observed IRX-βscatter. We also investigate the role of the environment in shaping dust attenuation in our sample. We find a strong correlation of the IRX-βrelation on gas-phase metallicity in our sample, and also strong correlation with galaxy compactness characterized by the Sérsic indexes. Correlations are also seen with stellar masses, specific star formation rates and the stellar ages of our sources. Metallicity strongly correlates with the IRX-βscatter, this also results from the older stars and higher masses at higher beta values. Galaxies with higher metallicities show higher IRX and higher beta values. The correlation with specific dust mass strongly shifts the galaxies away from the IRX-βrelation towards lower \b{eta} values. We find that more compact galaxies witness a larger amount of attenuation than less compact galaxies. There is a subtle variation in the dust attenuation scatter between edge-on and face-on galaxies, but the difference is not statistically significant. Galaxy environments do not significantly affect dust attenuation in our sample of star-forming galaxies at intermediate redshift. △ Less

Submitted 4 September, 2023; originally announced September 2023.

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

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

Influence of star-forming galaxy selection on the galaxy main sequence

Authors: W. J. Pearson, F. Pistis, M. Figueira, K. Małek, T. Moutard, D. Vergani, A. Pollo

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

Submitted 31 August, 2023; originally announced August 2023.

Comments: 20 pages, 3 appendices, 14 figures, 5 tables, accepted for publication in Astronomy & Astrophysics

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

Cosmic evolution of radio-AGN feedback: confronting models with data

Authors: R. Kondapally, P. N. Best, M. Raouf, N. L. Thomas, R. Davé, S. S. Shabala, H. J. A. Röttgering, M. J. Hardcastle, M. Bonato, R. K. Cochrane, K. Małek, L. K. Morabito, I. Prandoni, D. J. B. Smith

Abstract: Radio-mode feedback is a key ingredient in galaxy formation and evolution models, required to reproduce the observed properties of massive galaxies in the local Universe. We study the cosmic evolution of radio-AGN feedback out to $z\sim2.5$ using a sample of 9485 radio-excess AGN. We combine the evolving radio luminosity functions with a radio luminosity scaling relationship to estimate AGN jet ki… ▽ More Radio-mode feedback is a key ingredient in galaxy formation and evolution models, required to reproduce the observed properties of massive galaxies in the local Universe. We study the cosmic evolution of radio-AGN feedback out to $z\sim2.5$ using a sample of 9485 radio-excess AGN. We combine the evolving radio luminosity functions with a radio luminosity scaling relationship to estimate AGN jet kinetic powers and derive the cosmic evolution of the kinetic luminosity density, $Ω_{\rm{kin}}$ (i.e. the volume-averaged heating output). Compared to all radio-AGN, low-excitation radio galaxies (LERGs) dominate the feedback activity out to $z\sim2.5$, with both these populations showing a constant heating output of $Ω_{\rm{kin}} \approx 4-5 \times 10^{32}\,\rm{W\,Mpc^{-3}}$ across $0.5 < z < 2.5$. We compare our observations to predictions from semi-analytical and hydrodynamical simulations, which broadly match the observed evolution in $Ω_{\rm{kin}}$, although their absolute normalisation varies. Comparison to the Semi-Analytic Galaxy Evolution (SAGE) model suggests that radio-AGN may provide sufficient heating to offset radiative cooling losses, providing evidence for a self-regulated AGN feedback cycle. We integrate the kinetic luminosity density across cosmic time to obtain the kinetic energy density output from AGN jets throughout cosmic history to be $\sim 10^{50}\,\rm{J\,Mpc^{-3}}$. Compared to AGN winds, the kinetic energy density from AGN jets dominates the energy budget at $z \lesssim 2$; this suggests that AGN jets play an important role in AGN feedback across most of cosmic history. △ Less

Submitted 20 June, 2023; originally announced June 2023.

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

Variation of optical and infrared properties of galaxies with their surface brightness

Authors: Junais, K. Małek, S. Boissier, W. J. Pearson, A. Pollo, A. Boselli, M. Boquien, D. Donevski, T. Goto, M. Hamed, S. J. Kim, J. Koda, H. Matsuhara, G. Riccio, M. Romano

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

Submitted 10 June, 2023; originally announced June 2023.

Comments: 16 pages, 12 figures, accepted for publication in A&A

Journal ref: A&A 676, A41 (2023)

Environments of red nuggets at z~0.7 from the VIPERS survey

Authors: M. Siudek, K. Lisiecki, J. Krywult, D. Donevski, C. P. Haines, A. Karska, K. Małek, T. Moutard, A. Pollo

Abstract: Red ultra-compact massive galaxies, called red nuggets were formed at high redshifts ($\rm{z\sim2-3}$). Survivors of red nuggets, known as relics, observed at lower redshifts ($\rm{z<2}$) are believed to remain almost unchanged since their formation. For the first time, we verify the environmental properties of red nuggets at intermediate redshift ($0.5<\rm{z}<0.9$ ) using 42 red, massive (… ▽ More Red ultra-compact massive galaxies, called red nuggets were formed at high redshifts ($\rm{z\sim2-3}$). Survivors of red nuggets, known as relics, observed at lower redshifts ($\rm{z<2}$) are believed to remain almost unchanged since their formation. For the first time, we verify the environmental properties of red nuggets at intermediate redshift ($0.5<\rm{z}<0.9$ ) using 42 red, massive ($\rm{log(M_{star}/M_{\odot}) \geq 10.9}$) and ultra-compact ($\rm{R_{e}}<1.5$ kpc) from the VIMOS Public Extragalactic Redshift Survey (VIPERS). We found that the increasing fraction of red galaxies, when moving to denser environments, is driven by the red massive normal-size galaxies. Red nuggets, similarly to red intermediate-mass ($\rm{10.4\lesssim log(M_{star}/M_{\odot})<10.9}$) ultra-compact galaxies, are found in various types of environments, with consistent (within $1σ$) fractions across all local densities. Analysis of red nugget stellar ages suggests that relics are preferably found in high-density regions while quiescent red nuggets are overabundant in low-density environments. We speculate that red nuggets have survived to lower redshifts via two channels: i) in low-density environments where the fraction of red nuggets decreases as time passes due to (very) limited merger activity, ii) in high-density environments, where the number of red nuggets drops at higher redshift due to merger activity and is preserved at lower redshift as the high velocities of clusters prevent them from being cannibalised. Even more, the fraction of red nuggets in clusters may increase due to the addition of red massive normal-size galaxies deprived of their envelopes with cosmic time. △ Less

Submitted 12 June, 2023; v1 submitted 5 June, 2023; originally announced June 2023.

Comments: Accepted for publication in MNRAS, 15 pages, 8 figures, 3 appendices

The LOFAR Two-metre Sky Survey Deep Fields Data Release 1: V. Survey description, source classifications and host galaxy properties

Authors: P. N. Best, R. Kondapally, W. L. Williams, R. K. Cochrane, K. J. Duncan, C. L. Hale, P. Haskell, K. Malek, I. McCheyne, D. J. B. Smith, L. Wang, A. Botteon, M. Bonato, M. Bondi, G. Calistro Rivera, F. Gao, G. Gurkan, M. J. Hardcastle, M. J. Jarvis, B. Mingo, H. Miraghaei, L. K. Morabito, D. Nisbet, I. Prandoni, H. J. A. Rottgering , et al. (4 additional authors not shown)

Abstract: Source classifications, stellar masses and star formation rates are presented for 80,000 radio sources from the first data release of the Low Frequency Array Two-metre Sky Survey (LoTSS) Deep Fields, which represents the widest deep radio survey ever undertaken. Using deep multi-wavelength data spanning from the ultraviolet to the far-infrared, spectral energy distribution (SED) fitting is carried… ▽ More Source classifications, stellar masses and star formation rates are presented for 80,000 radio sources from the first data release of the Low Frequency Array Two-metre Sky Survey (LoTSS) Deep Fields, which represents the widest deep radio survey ever undertaken. Using deep multi-wavelength data spanning from the ultraviolet to the far-infrared, spectral energy distribution (SED) fitting is carried out for all of the LoTSS-Deep host galaxies using four different SED codes, two of which include modelling of the contributions from an active galactic nucleus (AGN). Comparing the results of the four codes, galaxies that host a radiative AGN are identified, and an optimised consensus estimate of the stellar mass and star-formation rate for each galaxy is derived. Those galaxies with an excess of radio emission over that expected from star formation are then identified, and the LoTSS-Deep sources are divided into four classes: star-forming galaxies, radio-quiet AGN, and radio-loud high-excitation and low-excitation AGN. Ninety-five per cent of the sources can be reliably classified, of which more than two-thirds are star-forming galaxies, ranging from normal galaxies in the nearby Universe to highly-starbursting systems at z>4. Star-forming galaxies become the dominant population below 150-MHz flux densities of about 1 mJy, accounting for 90 per cent of sources at a 150-MHz flux density of 100 microJy. Radio-quiet AGN comprise around 10 per cent of the overall population. Results are compared against the predictions of the SKADS and T-RECS radio sky simulations, and improvements to the simulations are suggested. △ Less

Submitted 9 May, 2023; originally announced May 2023.

Comments: Accepted for publication in MNRAS. Catalogues available at www.lofar-surveys.org/deepfields.html. 27 pages

The slippery slope of dust attenuation curves: Correlation of dust attenuation laws with star-to-dust compactness up to z = 4

Authors: M. Hamed, K. Małek, V. Buat, Junais, L. Ciesla, D. Donevski, G. Riccio, M. Figueira

Abstract: Aims. We investigate dust attenuation of 122 heavily dust-obscured galaxies detected with the Atacama Large Millimeter Array (ALMA) and Herschel in the COSMOS field. We search for correlations between dust attenuation recipes and the variation of physical parameters, mainly the effective radii of galaxies, their star formation rates (SFR), and stellar masses, and aim to understand which of the com… ▽ More Aims. We investigate dust attenuation of 122 heavily dust-obscured galaxies detected with the Atacama Large Millimeter Array (ALMA) and Herschel in the COSMOS field. We search for correlations between dust attenuation recipes and the variation of physical parameters, mainly the effective radii of galaxies, their star formation rates (SFR), and stellar masses, and aim to understand which of the commonly used laws best describes dust attenuation in dusty star-forming galaxies at high redshift. Methods. We make use of the extensive photometric coverage of the COSMOS data combined with highly-resolved dust continuum maps from ALMA. We use CIGALE to estimate various physical properties of these dusty objects, mainly their SFR, their stellar masses and their attenuation. We infer galaxy effective radii (Re) using GALFIT in the Y band of HSC and ALMA continuum maps. We use these radii to investigate the relative compactness of the dust continuum and the extension of the rest-frame UV/optical Re(y)/Re(ALMA). Results. We find that the physical parameters calculated from our models strongly depend on the assumption of dust attenuation curve. As expected, the most impacted parameter is the stellar mass, which leads to a change in the "starburstiness" of the objects. We find that taking into account the relative compactness of star-to-dust emission prior to SED fitting is crucial, especially when studying dust attenuation of dusty star-forming galaxies. Shallower attenuation curves did not show a clear preference of compactness with attenuation, while the Calzetti attenuation curve preferred comparable spatial extent of unattenuated stellar light and dust emission. The evolution of the Re(UV)/Re(ALMA) ratio with redshift peeks around the cosmic noon in our sample of DSFGs, showing that this compactness is correlated with the cosmic SFR density of these dusty sources. △ Less

Submitted 26 April, 2023; originally announced April 2023.

Comments: 18 pages, 13 figures, accepted for publication in A&A. Abstract abridged for arXiv submission

Journal ref: A&A 674, A99 (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 Art of Measuring Physical Parameters in Galaxies: A Critical Assessment of Spectral Energy Distribution Fitting Techniques

Authors: Camilla Pacifici, Kartheik G. Iyer, Bahram Mobasher, Elisabete da Cunha, Viviana Acquaviva, Denis Burgarella, Gabriela Calistro Rivera, Adam C. Carnall, Yu-Yen Chang, Nima Chartab, Kevin C. Cooke, Ciaran Fairhurst, Jeyhan Kartaltepe, Joel Leja, Katarzyna Malek, Brett Salmon, Marianna Torelli, Alba Vidal-Garcia, Mederic Boquien, Gabriel G. Brammer, Michael J. I. Brown, Peter L. Capak, Jacopo Chevallard, Chiara Circosta, Darren Croton , et al. (30 additional authors not shown)

Abstract: The study of galaxy evolution hinges on our ability to interpret multi-wavelength galaxy observations in terms of their physical properties. To do this, we rely on spectral energy distribution (SED) models which allow us to infer physical parameters from spectrophotometric data. In recent years, thanks to the wide and deep multi-waveband galaxy surveys, the volume of high quality data have signifi… ▽ More The study of galaxy evolution hinges on our ability to interpret multi-wavelength galaxy observations in terms of their physical properties. To do this, we rely on spectral energy distribution (SED) models which allow us to infer physical parameters from spectrophotometric data. In recent years, thanks to the wide and deep multi-waveband galaxy surveys, the volume of high quality data have significantly increased. Alongside the increased data, algorithms performing SED fitting have improved, including better modeling prescriptions, newer templates, and more extensive sampling in wavelength space. We present a comprehensive analysis of different SED fitting codes including their methods and output with the aim of measuring the uncertainties caused by the modeling assumptions. We apply fourteen of the most commonly used SED fitting codes on samples from the CANDELS photometric catalogs at z~1 and z~3. We find agreement on the stellar mass, while we observe some discrepancies in the star formation rate (SFR) and dust attenuation results. To explore the differences and biases among the codes, we explore the impact of the various modeling assumptions as they are set in the codes (e.g., star formation histories, nebular, dust, and AGN models) on the derived stellar masses, SFRs, and A_V values. We then assess the difference among the codes on the SFR-stellar mass relation and we measure the contribution to the uncertainties by the modeling choices (i.e., the modeling uncertainties) in stellar mass (~0.1dex), SFR (~0.3dex), and dust attenuation (~0.3mag). Finally, we present some resources summarizing best practices in SED fitting. △ Less

Submitted 4 December, 2022; originally announced December 2022.

Comments: 25 pages, 11 figures. Accepted for publication in ApJ

Unsupervised classification reveals new evolutionary pathways

Authors: M. Siudek, K. Lisiecki, M. Mezcua, K. Małek, A. Pollo, J. Krywult, A. Karska, Junais

Abstract: While we already seem to have a general scenario of the evolution of different types of galaxies, a complete and satisfactory understanding of the processes that led to the formation of all the variety of today's galaxy types is still beyond our reach. To solve this problem, we need both large datasets reaching high redshifts and novel methodologies for dealing with them. The VIPERS survey statist… ▽ More While we already seem to have a general scenario of the evolution of different types of galaxies, a complete and satisfactory understanding of the processes that led to the formation of all the variety of today's galaxy types is still beyond our reach. To solve this problem, we need both large datasets reaching high redshifts and novel methodologies for dealing with them. The VIPERS survey statistical power, which observed $\sim90,000$ galaxies at $z > 0.5$, and the application of an unsupervised clustering algorithm allowed us to distinguish 12 galaxy classes. Studies of their environmental dependence indicate that this classification may actually reflect different galaxy evolutionary paths. For instance, a class of the most passive red galaxies gathers galaxies $\sim20\%$ smaller than other red galaxies of a similar stellar mass, revealing the first sample of red nuggets at intermediate redshift. On the other end, a class of blue dwarf galaxies is composed mainly of AGN, challenging commonly used mid-infrared AGN selections. △ Less

Submitted 21 November, 2022; originally announced November 2022.

Comments: Accepted for publication in the ML4ASTRO (ICML 2022) proceeding book

SFR estimations from z=0 to z=0.9 -- A comparison of SFR calibrators for star-forming galaxies

Authors: M. Figueira, A. Pollo, K. Małek, V. Buat, M. Boquien, F. Pistis, L. P. Cassarà, D. Vergani, M. Hamed, S. Salim

Abstract: Using VIPERS, we estimated a set of SFR based on photometric and spectroscopic data. We used, as estimators, photometric bands from ultraviolet to mid-infrared, and spectral lines. Assuming a reference SFR obtained from the spectral energy distribution reconstructed with Code Investigating GALaxy Emission, we estimated the reliability of each band as an SFR tracer. We used GSWLC to trace the depen… ▽ More Using VIPERS, we estimated a set of SFR based on photometric and spectroscopic data. We used, as estimators, photometric bands from ultraviolet to mid-infrared, and spectral lines. Assuming a reference SFR obtained from the spectral energy distribution reconstructed with Code Investigating GALaxy Emission, we estimated the reliability of each band as an SFR tracer. We used GSWLC to trace the dependence of these SFR calibrators with redshift. The far and near UV, u-band and 24-$μ$m bands, as well as $L_{TIR}$, are found to be good SFR tracers up to $z\sim0.9$ with a strong dependence on the attenuation prescription used for the bluest bands (scatter of SFR of 0.26, 0.14, 0.15, 0.23, and 0.24dex for VIPERS, and 0.25, 0.24, 0.09, 0.12, and 0.12dex for GSWLC). The 8-$μ$m band provides only a rough estimate of the SFR as it depends on metallicity and polycyclic aromatic hydrocarbon properties (scatter of 0.23dex for VIPERS). We estimated the scatter of rest-frame luminosity estimations from CIGALE to be 0.26, 0.14, 0.12, 0.15, and 0.20dex for FUV, NUV, ugriz, K$_{\mathrm{s}}$, and 8-24$μ$m-$L_{\mathrm{TIR}}$). At intermediate redshift, the H$β$ line is a reliable SFR tracer (scatter of 0.19dex) and the [OII] line gives an equally good estimation when the metallicity from the $R_{23}$ parameter is taken into account (0.17 for VIPERS and 0.20dex for GSWLC). A calibration based on [OIII] retrieves the SFR only when additional information such as the metallicity or the ionization parameter of galaxies are used (0.26 for VIPERS and 0.20dex for GSWLC), diminishing its usability as a direct SFR tracer. Based on rest-frame luminosities estimated with CIGALE, we propose our own set of calibrations from FUV, NUV, u-band, 8, 24$μ$m, $L_{TIR}$, H$β$, [OII], and [OIII]. △ Less

Submitted 9 September, 2022; originally announced September 2022.

Comments: 35 pages, 28 figures, o be published in A&A

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

The first catalogue of spectroscopically confirmed red nuggets at z~0.7 from the VIPERS survey. Linking high-z red nuggets and local relics

Authors: Krzysztof Lisiecki, Katarzyna Małek, Małgorzata Siudek, Agnieszka Pollo, Janusz Krywult, Agata Karska, Junais

Abstract: 'Red nuggets' are a rare population of passive compact massive galaxies thought to be the first massive galaxies that formed in the Universe. First found at $z \sim 3$, they are even less abundant at lower redshifts, and it is believed that with time they mostly transformed through mergers into today's giant ellipticals. Those red nuggets which managed to escape this fate can serve as unique labor… ▽ More 'Red nuggets' are a rare population of passive compact massive galaxies thought to be the first massive galaxies that formed in the Universe. First found at $z \sim 3$, they are even less abundant at lower redshifts, and it is believed that with time they mostly transformed through mergers into today's giant ellipticals. Those red nuggets which managed to escape this fate can serve as unique laboratories to study the early evolution of massive galaxies. In this paper, we aim to make use of the VIMOS Public Extragalactic Redshift Survey to build the largest up-to-date catalogue of spectroscopically confirmed red nuggets at the intermediate redshift $0.5<z<1.0$. Starting from a catalogue of nearly 90 000 VIPERS galaxies we select sources with stellar masses $M_{star} > 8\times10^{10}$ $\rm{M}_{\odot}$ and effective radii $R_\mathrm{e}<1.5$ kpc. Among them, we select red, passive galaxies with old stellar population based on colour--colour NUVrK diagram, star formation rate values, and verification of their optical spectra. Verifying the influence of the limit of the source compactness on the selection, we found that the sample size can vary even up to two orders of magnitude, depending on the chosen criterion. Using one of the most restrictive criteria with additional checks on their spectra and passiveness, we spectroscopically identified only 77 previously unknown red nuggets. The resultant catalogue of 77 red nuggets is the largest such catalogue built based on the uniform set of selection criteria above the local Universe. Number density calculated on the final sample of 77 VIPERS passive red nuggets per comoving Mpc$^3$ increases from 4.7$\times10^{-6}$ at $z \sim 0.61$ to $9.8 \times 10^{-6}$ at $z \sim 0.95$, which is higher than values estimated in the local Universe, and lower than the ones found at $z>2$. It fills the gap at intermediate redshift. △ Less

Submitted 10 August, 2022; v1 submitted 9 August, 2022; originally announced August 2022.

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

Journal ref: A&A 669, A95 (2023)

A Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE).XIII. The role of ram-pressure stripping in transforming the diffuse and ultra-diffuse galaxies in the Virgo cluster

Authors: Junais, S. Boissier, A. Boselli, L. Ferrarese, P. Côté, S. Gwyn, J. Roediger, S. Lim, E. W. Peng, J. -C. Cuillandre, A. Longobardi, M. Fossati, G. Hensler, J. Koda, J. Bautista, M. Boquien, K. Małek, P. Amram, Y. Roehlly

Abstract: Low-surface-brightness galaxies (LSBs) contribute to a significant fraction of all the galaxies in the Universe. Ultra-diffuse galaxies (UDGs) form a subclass of LSBs that has attracted a lot of attention in recent years (although its definition may vary between studies). Although UDGs are found in large numbers in galaxy clusters, groups, and in the field, their formation and evolution are still… ▽ More Low-surface-brightness galaxies (LSBs) contribute to a significant fraction of all the galaxies in the Universe. Ultra-diffuse galaxies (UDGs) form a subclass of LSBs that has attracted a lot of attention in recent years (although its definition may vary between studies). Although UDGs are found in large numbers in galaxy clusters, groups, and in the field, their formation and evolution are still very much debated. Using a comprehensive set of multiwavelength data from the NGVS (optical), VESTIGE (H$α$ narrowband), and GUViCS (UV) surveys, we studied a sample of 64 diffuse galaxies and UDGs in the Virgo cluster to investigate their formation history. We analyzed the photometric colors and surface-brightness profiles of these galaxies and then compared them to models of galaxy evolution, including ram-pressure stripping (RPS) events to infer any possible strong interactions with the hot cluster gas in the past. While our sample consists mainly of red LSBs, which is typical in cluster environments, we found evidence of a color variation with the cluster-centric distance. Blue, HI-bearing, star-forming diffuse galaxies are found at larger distances from the cluster center than the rest of the sample. The comparison of our models with multifrequency observations suggests that most of the galaxies of the sample might have undergone a strong RPS event in their lifetime, on average 1.6 Gyr ago (with a large dispersion, and RPS still ongoing for some of them). This process resulted in the transformation of initially gas-rich diffuse blue galaxies into gas-poor and red ones that form the dominant population now, the more extreme UDGs having undergone the process in a more distant past on average. The RPS in dense environments could be one of the major mechanisms for the formation of the large number of quiescent UDGs we observe in galaxy clusters. △ Less

Submitted 4 August, 2022; originally announced August 2022.

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

Low frequency radio continuum imaging and SED modeling of 11 LIRGs: radio-only and FUV to radio bands

Authors: Subhrata Dey, Arti Goyal, Katarzyna Małek, Timothy J. Galvin, Nicholas Seymour, Tanio Díaz Santos, Julia Piotrowska, Vassilis Charmandaris

Abstract: We present the detailed analysis of 11 local luminous infrared galaxies (LIRGs) from ultraviolet through far-infrared to radio ($\sim$70 MHz to $\sim$15 GHz) bands. We derive the astrophysical properties through spectral energy distribution (SED) modeling using the Code Investigating GALaxy Emission (CIGALE) and UltraNest codes. The radio SEDs include our new observations at 325 and 610 MHz from t… ▽ More We present the detailed analysis of 11 local luminous infrared galaxies (LIRGs) from ultraviolet through far-infrared to radio ($\sim$70 MHz to $\sim$15 GHz) bands. We derive the astrophysical properties through spectral energy distribution (SED) modeling using the Code Investigating GALaxy Emission (CIGALE) and UltraNest codes. The radio SEDs include our new observations at 325 and 610 MHz from the GMRT and the measurements from public archives. Our main results are (1) radio SEDs show turnovers and bends, (2) the synchrotron spectral index of the fitted radio spectra ranges between $-$0.5 and $-$1.7, and (3) the infrared luminosity, dust mass, dust temperature, stellar mass, star-formation rates (SFRs) and AGN fraction obtained from CIGALE falls in the range exhibited by galaxies of the same class. The ratio of 60$μ$m infrared and 1.4 GHz radio luminosity, the 1.4 GHz thermal fraction, and emission measure range between 2.1 and 2.9, 0.1% and 10%, 0.02 and 269.5$\times$10$^{6}$ cm$^{-6}$ pc, respectively. We conclude that the turnovers seen in the radio SEDs are due to free-free absorption; this is supported by the low AGN fraction derived from the CIGALE analysis. The decomposed 1.4 GHz thermal and nonthermal radio luminosities allowed us to compute the star formation rate (SFR) using scaling relations. A positive correlation is observed between the SFR$_{IR}$ obtained 10 Myr ago (compared to 100 Myr ago) and 1.4 GHz radio (total and nonthermal) because similar synchrotron lifetimes are expected for typical magnetic field strengths observed in these galaxies ($\approx$50$μ$G). △ Less

Submitted 18 July, 2022; originally announced July 2022.

Comments: ApJ accepted. Comments are welcome

COSMOS2020: Manifold Learning to Estimate Physical Parameters in Large Galaxy Surveys

Authors: I. Davidzon, K. Jegatheesan, O. Ilbert, S. de la Torre, S. K. Leslie, C. Laigle, S. Hemmati, D. C. Masters, D. Blanquez-Sese, O. B. Kauffmann, G. E. Magdis, K. Małek, H. J. McCracken, B. Mobasher, A. Moneti, D. B. Sanders, M. Shuntov, S. Toft, J. R. Weaver

Abstract: We present a novel method to estimate galaxy physical properties from spectral energy distributions (SEDs), alternate to template fitting techniques and based on self-organizing maps (SOM) to learn the high-dimensional manifold of a photometric galaxy catalog. The method has been previously tested with hydrodynamical simulations in Davidzon et al. (2019) while here is applied to real data for the… ▽ More We present a novel method to estimate galaxy physical properties from spectral energy distributions (SEDs), alternate to template fitting techniques and based on self-organizing maps (SOM) to learn the high-dimensional manifold of a photometric galaxy catalog. The method has been previously tested with hydrodynamical simulations in Davidzon et al. (2019) while here is applied to real data for the first time. It is crucial for its implementation to build the SOM with a high quality, panchromatic data set, which we elect to be the "COSMOS2020" galaxy catalog. After the training and calibration steps with COSMOS2020, other galaxies can be processed through SOM to obtain an estimate of their stellar mass and star formation rate (SFR). Both quantities result to be in good agreement with independent measurements derived from more extended photometric baseline, and also their combination (i.e., the SFR vs. stellar mass diagram) shows a main sequence of star forming galaxies consistent with previous studies. We discuss the advantages of this method compared to traditional SED fitting, highlighting the impact of having, instead of the usual synthetic templates, a collection of empirical SEDs built by the SOM in a "data-driven" way. Such an approach also allows, even for extremely large data sets, an efficient visual inspection to identify photometric errors or peculiar galaxy types. Considering in addition the computational speed of this new estimator, we argue that it will play a valuable role in the analysis of oncoming large-area surveys like Euclid or the Legacy Survey of Space and Time at the Vera Cooper Rubin Telescope. △ Less

Submitted 13 June, 2022; originally announced June 2022.

Comments: to appear on A&A

Journal ref: A&A 665, A34 (2022)

The fundamental metallicity relation from SDSS (z~0) to VIPERS (z~0.7) Data selection or evolution

Authors: F. Pistis, A. Pollo, M. Scodeggio, M. Figueira, A. Durkalec, K. Małek, A. Iovino, D. Vergani, S. Salim

Abstract: Aims. We focus on how to homogenize the comparison between galaxy samples having different characteristics. We check the projections of the fundamental metallicity relation (FMR) and the evolution of these projections between a sample selected at $z\sim0$ (SDSS) and $z\sim0.7$ (VIPERS). We check, in particular, whether and to what extent selection criteria can affect the results. Methods. We check… ▽ More Aims. We focus on how to homogenize the comparison between galaxy samples having different characteristics. We check the projections of the fundamental metallicity relation (FMR) and the evolution of these projections between a sample selected at $z\sim0$ (SDSS) and $z\sim0.7$ (VIPERS). We check, in particular, whether and to what extent selection criteria can affect the results. Methods. We checked the influence of different biases introduced either by physical constraints (evolution of the luminosity function and differences in the fraction of blue galaxies) or data selection (the signal-to-noise ratio and quality of the spectra) on the FMR and its projections. To separate the differences occurring due to the physical evolution of galaxies with redshift from the false evolution mimed by these biases, we first analyzed the effects of these biases individually on the SDSS sample, and next, starting from the SDSS data, we built a VIPERS-equivalent $z \sim 0$ sample, replicating the main characteristics of VIPERS sample at $z\sim0.7$ for a fair comparison. Results. We found that the FMR projections are all sensitive to biases introduced by the selection on S/N and the quality flags of the emission line measurements in the spectra, especially the $\left[\text{O{\,\sc{iii}}}\right]λ4959$ line. The exception is the metallicity versus the sSFR plane which is insensitive to these biases. The evolution of the luminosity function introduces a bias only in the plane metallicity versus the star formation rate (SFR) while the fraction of blue galaxies has no impact on results. △ Less

Submitted 6 June, 2022; originally announced June 2022.

Comments: Accepted for publication in Astronomy & Astrophysics, 18 pages, 19 figures, 3 Tables

Journal ref: A&A 663, A162 (2022)

Star Formation and AGN Feedback in the Local Universe: Combining LOFAR and MaNGA

Authors: C. R. Mulcahey, S. K. Leslie, T. M. Jackson, J. E. Young, I. Prandoni, M. J. Hardcastle, N. Roy, K. Małek, M. Magliocchetti, M. Bonato, H. J. A. Röttgering, A. Drabent

Abstract: The effect of Active Galactic Nuclei (AGN) on their host galaxies -- in particular their levels of star formation -- remains one of the key outstanding questions of galaxy evolution. Successful cosmological models of galaxy evolution require a fraction of energy released by an AGN to be redistributed into the interstellar medium to reproduce the observed stellar mass and luminosity function and to… ▽ More The effect of Active Galactic Nuclei (AGN) on their host galaxies -- in particular their levels of star formation -- remains one of the key outstanding questions of galaxy evolution. Successful cosmological models of galaxy evolution require a fraction of energy released by an AGN to be redistributed into the interstellar medium to reproduce the observed stellar mass and luminosity function and to prevent the formation of over-massive galaxies. Observations have confirmed that the radio-AGN population is energetically capable of heating and redistributing gas at all phases, however, direct evidence of AGN enhancing or quenching star formation remains rare. With modern, deep radio surveys and large integral field spectroscopy (IFS) surveys, we can detect fainter synchrotron emission from AGN jets and accurately probe the star-forming properties of galaxies, respectively. In this paper, we combine data from the LOw Frequency ARray Two-meter Sky Survey with data from one of the largest optical IFS surveys, Mapping Nearby Galaxies at Apache Point Observatory to probe the star-forming properties of 307 local (z $<$ 0.15) galaxies that host radio-detected AGN (RDAGN). We compare our results to a robust control sample of non-active galaxies that each match the stellar mass, redshift, visual morphology, and inclination of a RDAGN host. We find that RDAGN and control galaxies have broad SFR distributions, typically lie below the star-forming main-sequence, and have negative stellar light-weighted age gradients. These results indicate that AGN selected based on their current activity are not responsible for suppressing their host galaxies' star formation. Rather, our results support the maintenance mode role that radio AGN are expected to have in the local Universe. △ Less

Submitted 2 June, 2022; originally announced June 2022.

Comments: Accepted for publication in A&A. 24 pages, 14 figures

Journal ref: A&A 665, A144 (2022)

Shaping physical properties of galaxy subtypes in the VIPERS survey: environment matters

Authors: M. Siudek, K. Malek, A. Pollo, A. Iovino, C. P. Haines, M. Bolzonella, O. Cucciati, A. Gargiulo, B. Granett, J. Krywult, T. Moutard, M. Scodeggio

Abstract: Aims. This study aims to explore the relationship between the physical properties of different galaxy subclasses and their environment based on the analysis of 31 631 VIMOS Public Extragalactic Redshift Survey (VIPERS) galaxies observed at 0.5 < z < 0.9. Methods. We use the results of an unsupervised clustering algorithm to distinguish 11 subclasses of VIPERS galaxies based on the multi-dimensio… ▽ More Aims. This study aims to explore the relationship between the physical properties of different galaxy subclasses and their environment based on the analysis of 31 631 VIMOS Public Extragalactic Redshift Survey (VIPERS) galaxies observed at 0.5 < z < 0.9. Methods. We use the results of an unsupervised clustering algorithm to distinguish 11 subclasses of VIPERS galaxies based on the multi-dimensional feature space defined by rest-frame UV to NIR colours presented in Siudek et al (2018a). We investigate the relationship between the properties of these subclasses of galaxies and their local environment, defined as the galaxy density contrast derived from the 5th nearest neighbour technique. Results. We confirm that the galaxy population-density relation is already in place at z ~ 0.9, with the blue galaxy fraction decreasing with density, compensated by an increase of the red fraction. On average red galaxies in the high-density environment are larger by 28% than the ones in low-density environments. In particular, we find one group of galaxies, subclass C3, whose increase of size with time can be explained mainly as the result of mergers; for other red subclasses, mergers would not seem to play a major role (subclass C2) or play a negligible role (subclass C1). The properties of the green galaxies (subclasses C4-6) depend on whether their stellar mass is above or below a transition mass. Low-mass green galaxies appear to have grown through secular processes, while in high-mass green galaxies mass assembly appears to be dominated by mergers. When it comes to blue galaxies, the trend of decreasing fraction with denser environments seen for the group as a whole (subclasses C7-11) is found to be driven mostly by one group of galaxies, subclass C10. These are compact low-mass galaxies with high sSFRs, that are preferentially found in low-density environments. △ Less

Submitted 29 May, 2022; originally announced May 2022.

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

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

Comparison of the star formation of X-ray selected active galactic nuclei in eFEDS with star-forming galaxies

Authors: G. Mountrichas, V. Buat, G. Yang, M. Boquien, D. Burgarella, L. Ciesla, K. Malek, R. Shirley

Abstract: We use $\sim 1800$ X-ray Active Galactic Nuclei (AGN) in the eROSITA Final Equatorial-Depth Survey (eFEDS), that span over two orders of magnitude in X-ray luminosity, $\rm L_{X,2-10keV} \approx 10^{43-45}\,ergs^{-1}$, and compare their star-formation rate (SFR) relative to that of non-AGN star-forming systems, at $\rm 0.5<z<1.5$. For that purpose, we compile a reference galaxy catalogue with… ▽ More We use $\sim 1800$ X-ray Active Galactic Nuclei (AGN) in the eROSITA Final Equatorial-Depth Survey (eFEDS), that span over two orders of magnitude in X-ray luminosity, $\rm L_{X,2-10keV} \approx 10^{43-45}\,ergs^{-1}$, and compare their star-formation rate (SFR) relative to that of non-AGN star-forming systems, at $\rm 0.5<z<1.5$. For that purpose, we compile a reference galaxy catalogue with $\sim 17000$ sources. Both samples have the same photometric coverage, from optical to far-infrared. We construct the spectral energy distributions (SEDs) of all sources and fit them using the CIGALE code, utilizing the same templates and parametric grid for both samples. We account for the mass incompleteness and exclude quiescent systems from both datasets. These allow us to compare the SFR of the two populations in a uniform manner, minimising systematic effects. Based on our analysis, AGN at low and moderate L$_X$ ($\rm L_{X,2-10keV} < 10^{44}\,ergs^{-1}$), have SFR that is lower, or at most, equal to that of star-forming galaxies, in agreement with previous studies. The large number of luminous X-ray AGN available in the eFEDS dataset, enable us to expand our investigations at higher L$_X$ to test previous, tentative results. At $\rm L_{X,2-10keV} > 10^{44.2}\,ergs^{-1}$, the SFR of AGN appears enhanced, by $\sim 30\%$, compared to that of star-forming sources, for systems with stellar mass, $\rm 10.5 < log\,[M_*(M_\odot)] < 11.5$, confirming indications found in previous studies. The most massive sources $\rm log\,[M_*(M_\odot)] > 11.5$, present a flat SFR$_{norm}$-L$_X$ relation up to $\rm L_{X,2-10keV} \sim 10^{44.5}\,ergs^{-1}$, with SFR similar to that of star-forming galaxies. However, at higher L$_X$ ($\rm L_{X,2-10keV} \sim 10^{45}\,ergs^{-1}$), we find indications that the SFR of these massive AGN hosts may be enhanced compared to that of non-AGN systems. △ Less

Submitted 5 May, 2022; v1 submitted 3 May, 2022; originally announced May 2022.

Comments: Accepted for publication in A&A. 14 pages, 7 figures, 2 tables

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

Cosmic evolution of low-excitation radio galaxies in the LOFAR Two-meter Sky Survey Deep Fields

Authors: R. Kondapally, P. N. Best, R. K. Cochrane, J. Sabater, K. J. Duncan, M. J. Hardcastle, P. Haskell, B. Mingo, H. J. A. Röttgering, D. J. B. Smith, W. L. Williams, M. Bonato, G. Calistro Rivera, F. Gao, C. L. Hale, K. Małek, G. K. Miley, I. Prandoni, L. Wang

Abstract: Feedback from low-excitation radio galaxies (LERGs) plays a key role in the lifecycle of massive galaxies in the local Universe; their evolution, and the impact of these active galactic nuclei on early galaxy evolution, however, remain poorly understood. We use a sample of 10481 LERGs from the first data release of the LOFAR Two-meter Sky Survey Deep Fields, covering $\sim$ 25 deg$^2$, to present… ▽ More Feedback from low-excitation radio galaxies (LERGs) plays a key role in the lifecycle of massive galaxies in the local Universe; their evolution, and the impact of these active galactic nuclei on early galaxy evolution, however, remain poorly understood. We use a sample of 10481 LERGs from the first data release of the LOFAR Two-meter Sky Survey Deep Fields, covering $\sim$ 25 deg$^2$, to present the first measurement of the evolution of the radio luminosity function (LF) of LERGs out to $z\sim2.5$; this shows relatively mild evolution. We split the LERGs into those hosted by quiescent and star-forming galaxies, finding a new dominant population of LERGs hosted by star-forming galaxies at high redshifts. The incidence of LERGs in quiescent galaxies shows a steep dependence on stellar-mass out to $z \sim1.5$, consistent with local Universe measurements of accretion occurring from cooling of hot gas haloes. The quiescent-LERGs dominate the LFs at $z<1$, showing a strong decline in space density with redshift, tracing that of the available host galaxies, while there is an increase in the characteristic luminosity. The star-forming LERG LF increases with redshift, such that this population dominates the space densities at most radio-luminosities by $z \sim 1$. The incidence of LERGs in star-forming galaxies shows a much weaker stellar-mass dependence, and increases with redshift, suggesting a different fuelling mechanism compared to their quiescent counterparts, potentially associated with the cold gas supply present in the star-forming galaxies. △ Less

Submitted 22 April, 2022; v1 submitted 15 April, 2022; originally announced April 2022.

Comments: 26 pages, 16 figures, accepted for publication in MNRAS. LF data available at: https://github.com/rohitk-10/AGN_LF_Kondapally22 . Author list error corrected in updated version

The LOFAR Two-metre Sky Survey -- V. Second data release

Authors: T. W. Shimwell, M. J. Hardcastle, C. Tasse, P. N. Best, H. J. A. Röttgering, W. L. Williams, A. Botteon, A. Drabent, A. Mechev, A. Shulevski, R. J. van Weeren, L. Bester, M. Brüggen, G. Brunetti, J. R. Callingham, K. T. Chyży, J. E. Conway, T. J. Dijkema, K. Duncan, F. de Gasperin, C. L. Hale, M. Haverkorn, B. Hugo, N. Jackson, M. Mevius , et al. (81 additional authors not shown)

Abstract: In this data release from the LOFAR Two-metre Sky Survey (LoTSS) we present 120-168MHz images covering 27% of the northern sky. Our coverage is split into two regions centred at approximately 12h45m +44$^\circ$30' and 1h00m +28$^\circ$00' and spanning 4178 and 1457 square degrees respectively. The images were derived from 3,451hrs (7.6PB) of LOFAR High Band Antenna data which were corrected for th… ▽ More In this data release from the LOFAR Two-metre Sky Survey (LoTSS) we present 120-168MHz images covering 27% of the northern sky. Our coverage is split into two regions centred at approximately 12h45m +44$^\circ$30' and 1h00m +28$^\circ$00' and spanning 4178 and 1457 square degrees respectively. The images were derived from 3,451hrs (7.6PB) of LOFAR High Band Antenna data which were corrected for the direction-independent instrumental properties as well as direction-dependent ionospheric distortions during extensive, but fully automated, data processing. A catalogue of 4,396,228 radio sources is derived from our total intensity (Stokes I) maps, where the majority of these have never been detected at radio wavelengths before. At 6" resolution, our full bandwidth Stokes I continuum maps with a central frequency of 144MHz have: a median rms sensitivity of 83$μ$Jy/beam; a flux density scale accuracy of approximately 10%; an astrometric accuracy of 0.2"; and we estimate the point-source completeness to be 90% at a peak brightness of 0.8mJy/beam. By creating three 16MHz bandwidth images across the band we are able to measure the in-band spectral index of many sources, albeit with an error on the derived spectral index of +/-0.2 which is a consequence of our flux-density scale accuracy and small fractional bandwidth. Our circular polarisation (Stokes V) 20" resolution 120-168MHz continuum images have a median rms sensitivity of 95$μ$Jy/beam, and we estimate a Stokes I to Stokes V leakage of 0.056%. Our linear polarisation (Stokes Q and Stokes U) image cubes consist of 480 x 97.6 kHz wide planes and have a median rms sensitivity per plane of 10.8mJy/beam at 4' and 2.2mJy/beam at 20"; we estimate the Stokes I to Stokes Q/U leakage to be approximately 0.2%. Here we characterise and publicly release our Stokes I, Q, U and V images in addition to the calibrated uv-data. △ Less

Submitted 23 February, 2022; originally announced February 2022.

Comments: 23 figures, 1 table and 29 pages. The catalogues, images and uv-data associated with this data release are publicly available via https://lofar-surveys.org/

North Ecliptic Pole merging galaxy catalogue

Authors: W. J. Pearson, L. E. Suelves, S. C. -C. Ho, N. Oi, S. Brough, B. W. Holwerda, A. M. Hopkins, T. -C. Huang, H. S. Hwang, L. S. Kelvin, S. J. Kim, Á. R. López-Sánchez, K. Małek, C. Pearson, A. Poliszczuk, A. Pollo, V. Rodriguez-Gomez, H. Shim, Y. Toba, L. Wang

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

Submitted 22 February, 2022; originally announced February 2022.

Comments: Accepted to A&A, 26 pages, 20 figures, 8 tables, 3 appendixes, full tables 1 and 4 will be available on CDS

Radio AGN in nearby dwarf galaxies: the important role of AGN in dwarf-galaxy evolution

Authors: F. Davis, S. Kaviraj, M. J. Hardcastle, G. Martin, R. A. Jackson, K. Kraljic, K. Malek, S. Peirani, D. J. B. Smith, M. Volonteri, L. Wang

Abstract: We combine deep optical and radio data, from the Hyper Suprime-Cam and the Low-Frequency Array (LOFAR) respectively, to study 78 radio AGN in nearby (z<0.5) dwarf galaxies. Comparison to a control sample, matched in stellar mass and redshift, indicates that the AGN and controls reside in similar environments, show similar star-formation rates (which trace gas availability) and exhibit a comparable… ▽ More We combine deep optical and radio data, from the Hyper Suprime-Cam and the Low-Frequency Array (LOFAR) respectively, to study 78 radio AGN in nearby (z<0.5) dwarf galaxies. Comparison to a control sample, matched in stellar mass and redshift, indicates that the AGN and controls reside in similar environments, show similar star-formation rates (which trace gas availability) and exhibit a comparable incidence of tidal features (which indicate recent interactions). We explore the AGN properties by combining the predicted gas conditions in dwarfs from a cosmological hydrodynamical simulation with a Monte-Carlo suite of simulated radio sources, based on a semi-analytical model for radio-galaxy evolution. In the subset of LOFAR-detectable simulated sources, which have a similar distribution of radio luminosities as our observed AGN, the median jet powers, ages and accretion rates are $\sim 10^{35}$ W, $\sim 5$ Myr and $\sim 10^{-3.4}$ M$_{\odot}$ yr$^{-1}$ respectively. The median mechanical energy output of these sources is $\sim100$ times larger than the median binding energy expected in dwarf gas reservoirs, making AGN feedback plausible. Since special circumstances (in terms of environment, gas availability and interactions) are not necessary for the presence of AGN, and the central gas masses are predicted to be an order of magnitude larger than that required to fuel the AGN, AGN triggering in dwarfs is likely to be stochastic and a common phenomenon. Together with the plausibility of energetic feedback, this suggests that AGN could be important drivers of dwarf-galaxy evolution, as is the case in massive galaxies. △ Less

Submitted 24 January, 2022; originally announced January 2022.

Comments: Accepted for publication in MNRAS

AGN in the ULIRG HE 0435-5304

Authors: Krzysztof Hryniewicz, Małgorzata Bankowicz, Katarzyna Małek, Aleksander Herzig, Agnieszka Pollo

Abstract: HE 0435-5304 from Hamburg European Southern Observatory survey is a quasar that appears in the literature with two conflicting redshift values: $\sim 1.2$ and $\sim 0.4$. It was used in the studies of the intergalactic medium through fitting of the narrow absorption lines in its ultraviolet (UV) spectrum. This source is also known historically as a luminous infrared galaxy. We present optical spec… ▽ More HE 0435-5304 from Hamburg European Southern Observatory survey is a quasar that appears in the literature with two conflicting redshift values: $\sim 1.2$ and $\sim 0.4$. It was used in the studies of the intergalactic medium through fitting of the narrow absorption lines in its ultraviolet (UV) spectrum. This source is also known historically as a luminous infrared galaxy. We present optical spectra of HE 0435-5304, aiming to precisely measure its redshift and to study its physical properties. In particular, properties of its active nucleus, which is studied in the context of the source being identified here as an ultra-luminous infrared galaxy, allow us to place this quasar in the context of the general population. Fitting the spectra, we focused on modeling H$β$ and [O III] lines. Based on these, we derived the virial black hole mass, bolometric luminosity, and Eddington ratio of the active galactic nucleus (AGN). Additionally, we performed broad band photometry fitting which allows us to quantify host galaxy parameters. The improved redshift value of HE 0435-5304 is estimated to $0.42788 \pm 0.00027$ based on the [O II] line, which is mostly consistent with the narrowest components of the other emission lines. The source was found to be a relatively massive and luminous AGN whose host galaxy is actively forming stars. Although its stellar population seems to be heavily obscured, we did not find evidence for significant obscuration of the nucleus. We conclude that the AGN HE 0435-5304 is a rather prominent iron emitter from the extreme type-A population very close to the narrow-line Seyfert 1 group. The fact that the width of the H$β$ line appears to be systematically growing in its broadest component with time may suggest that this AGN is changing its broad line region. Due to the influence of atmospheric effects, this finding is uncertain. △ Less

Submitted 12 January, 2022; originally announced January 2022.

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

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

Fitting AGN/galaxy X-ray-to-radio SEDs with CIGALE and improvement of the code

Authors: Guang Yang, Médéric Boquien, William N. Brandt, Véronique Buat, Denis Burgarella, Laure Ciesla, Bret D. Lehmer, Katarzyna E. Małek, George Mountrichas, Casey Papovich, Estelle Pons, Marko Stalevski, Patrice Theulé, Shifu Zhu

Abstract: Modern and future surveys effectively provide a panchromatic view for large numbers of extragalactic objects. Consistently modeling these multiwavelength survey data is a critical but challenging task for extragalactic studies. The Code Investigating GALaxy Emission (CIGALE) is an efficient PYTHON code for spectral energy distribution (SED) fitting of galaxies and active galactic nuclei (AGNs). Re… ▽ More Modern and future surveys effectively provide a panchromatic view for large numbers of extragalactic objects. Consistently modeling these multiwavelength survey data is a critical but challenging task for extragalactic studies. The Code Investigating GALaxy Emission (CIGALE) is an efficient PYTHON code for spectral energy distribution (SED) fitting of galaxies and active galactic nuclei (AGNs). Recently, a major extension of CIGALE (named X-CIGALE) has been developed to account for AGN/galaxy X-ray emission and improve AGN modeling at UV-to-IR wavelengths. Here, we apply X-CIGALE to different samples, including COSMOS spectroscopic type 2 AGNs, CDF-S X-ray detected normal galaxies, SDSS quasars, and COSMOS radio objects. From these tests, we identify several weaknesses of X-CIGALE and improve the code accordingly. These improvements are mainly related to AGN intrinsic X-ray anisotropy, X-ray binary emission, AGN accretion-disk SED shape, and AGN radio emission. These updates improve the fit quality and allow new interpretation of the results, based on which we discuss physical implications. For example, we find that AGN intrinsic X-ray anisotropy is moderate, and can be modeled as $L_X(θ) \propto 1+\cos θ$, where $θ$ is the viewing angle measured from the AGN axis. We merge the new code into the major branch of CIGALE, and publicly release this new version as CIGALE v2022.0 on https://cigale.lam.fr △ Less

Submitted 10 January, 2022; originally announced January 2022.

Comments: ApJ accepted. The code is available at https://cigale.lam.fr/news/

The Type II AGN-host galaxy connection: insights from the VVDS and VIPERS surveys

Authors: G. Vietri, B. Garilli, M. Polletta, S. Bisogni, L. P. Cassarà, P. Franzetti, M. Fumana, A. Gargiulo, D. Maccagni, C. Mancini, M. Scodeggio, A. Fritz, K. Malek, G. Manzoni, A. Pollo, M. Siudek, D. Vergani, G. Zamorani, A. Zanichelli

Abstract: We present a study of optically-selected Type II AGN at 0.5 < z < 0.9 from the VIPERS and VVDS surveys, to investigate the connection between AGN activity and physical properties of their host galaxies. The host stellar mass is estimated through spectral energy distribution fitting with the CIGALE code, and star formation rates are derived from the [OII]$λ$3727 $Å$ line luminosity. We find that 49… ▽ More We present a study of optically-selected Type II AGN at 0.5 < z < 0.9 from the VIPERS and VVDS surveys, to investigate the connection between AGN activity and physical properties of their host galaxies. The host stellar mass is estimated through spectral energy distribution fitting with the CIGALE code, and star formation rates are derived from the [OII]$λ$3727 $Å$ line luminosity. We find that 49% of the AGN host galaxies are on or above the main sequence (MS), 40% lie in the sub-MS locus, and 11% in the quiescent locus. Using the [OIII]$λ$5007 $Å$ line luminosity as a proxy of the AGN power, we find that at fixed AGN power Type II AGN host galaxies show a bimodal behaviour: systems with host galaxy stellar mass <10$^{10}$ M$_{\odot}$, reside along the MS or in the starbursts locus (high-SF Type II AGN), while systems residing in massive host-galaxies (>10$^{10}$ M$_{\odot}$) show a lower level of star formation (low-SF Type II AGN). At all stellar masses, the offset from the MS is positively correlated with the AGN power. We interpret this correlation as evidence of co-evolution between the AGN and the host, possibly due to the availability of cold gas. In the most powerful AGN with host galaxies below the MS we find a hint, though weak, of asymmetry in the [OIII] line profile, likely due to outflowing gas, consistent with a scenario in which AGN feedback removes the available gas and halts the star formation in the most massive hosts. △ Less

Submitted 16 November, 2021; originally announced November 2021.

Comments: Accepted for publication in A&A

Journal ref: A&A 659, A129 (2022)

The LOFAR Two-metre Sky Survey Deep fields: A new analysis of low-frequency radio luminosity as a star-formation tracer in the Lockman Hole region

Authors: M. Bonato, I. Prandoni, G. De Zotti, P. N. Best, M. Bondi, G. Calistro Rivera, R. K. Cochrane, G. Gürkan, P. Haskell, R. Kondapally, M. Magliocchetti, S. K. Leslie, K. Malek, H. J. A. Röttgering, D. J. B. Smith, C. Tasse, L. Wang

Abstract: We have exploited LOFAR deep observations of the Lockman Hole field at 150 MHz to investigate the relation between the radio luminosity of star-forming galaxies (SFGs) and their star formation rates (SFRs), as well as its dependence on stellar mass and redshift. The adopted source classification, SFRs and stellar masses are consensus estimates based on a combination of four different SED fitting m… ▽ More We have exploited LOFAR deep observations of the Lockman Hole field at 150 MHz to investigate the relation between the radio luminosity of star-forming galaxies (SFGs) and their star formation rates (SFRs), as well as its dependence on stellar mass and redshift. The adopted source classification, SFRs and stellar masses are consensus estimates based on a combination of four different SED fitting methods. We note a flattening of radio spectra of a substantial minority of sources below $\sim 1.4 $ GHz. Such sources have thus a lower "radio-loudness" level at 150 MHz than expected from extrapolations from 1.4 GHz using the average spectral index. We found a weak trend towards a lower $\hbox{SFR}/L_{150 \rm MHz}$ ratio for higher stellar mass, $M_\star$. We argue that such a trend may account for most of the apparent redshift evolution of the $L_{150 \rm MHz}/\hbox{SFR}$ ratio, in line with previous work. Our data indicate a weaker evolution than found by some previous analyses. We also find a weaker evolution with redshift of the specific star formation rate than found by several (but not all) previous studies. Our radio selection provides a view of the distribution of galaxies in the $\hbox{SFR}$-$M_\star$ plane complementary to that of optical/near-IR selection. It suggests a higher uniformity of the star formation history of galaxies than implied by some analyses of optical/near-IR data. We have derived luminosity functions at 150 MHz of both SFGs and radio-quiet (RQ) AGN at various redshifts. Our results are in very good agreement with the T-RECS simulations and with literature estimates. We also present explicit estimates of SFR functions of SFGs and RQ AGN at several redshifts derived from our radio survey data. △ Less

Submitted 14 September, 2021; originally announced September 2021.

Comments: 18 pages, 9 figures, 2 tables, A&A accepted

Journal ref: A&A 656, A48 (2021)

The viewing angle in AGN SED models, a data-driven analysis

Authors: Andrés Felipe Ramos Padilla, Lingyu Wang, Katarzyna Małek, Andreas Efstathiou, Guang Yang

Abstract: The validity of the unified active galactic nuclei (AGN) model has been challenged in the last decade, especially when different types of AGNs are considered to only differ in the viewing angle to the torus. We aim to assess the importance of the viewing angle in classifying different types of Seyfert galaxies in spectral energy distribution (SED) modelling. We retrieve photometric data from publi… ▽ More The validity of the unified active galactic nuclei (AGN) model has been challenged in the last decade, especially when different types of AGNs are considered to only differ in the viewing angle to the torus. We aim to assess the importance of the viewing angle in classifying different types of Seyfert galaxies in spectral energy distribution (SED) modelling. We retrieve photometric data from publicly available astronomical databases: CDS and NED, to model SEDs with X-CIGALE in a sample of 13 173 Seyfert galaxies located at redshift range from $z=0$ to $z=3.5$, with a median redshift of $z\approx0.2$. We assess whether the estimated viewing angle from the SED models reflects different Seyfert classifications. Two AGN models with either a smooth or clumpy torus structure are adopted in this paper. We find that the viewing angle in Type-1 AGNs is better constrained than in Type-2 AGNs. Limiting the viewing angles representing these two types of AGNs do not affect the physical parameter estimates such as star-formation rate (SFR) or AGN fractional contribution ($f_{\rm{AGN}}$). In addition, the viewing angle is not the most discriminating physical parameter to differentiate Seyfert types. We suggest that the observed and intrinsic AGN disc luminosity can: i) be used in $z<0.5$ studies to distinguish between Type-1 and Type-2 AGNs, and ii) explain the probable evolutionary path between these AGN types. Finally, we propose the use of X-CIGALE for AGN galaxy classification tasks. All data from the 13 173 SED fits are available at https://doi.org/10.5281/zenodo.5221764 △ Less

Submitted 29 November, 2021; v1 submitted 24 August, 2021; originally announced August 2021.

Comments: 19 pages, 10 figures, 6 tables (plus appendix). Accepted for publication in MNRAS

Environmental Effects on AGN activity via Extinction-free Mid-Infrared Census

Authors: Daryl Joe D. Santos, Tomotsugu Goto, Seong Jin Kim, Ting-Wen Wang, Simon C. -C. Ho, Tetsuya Hashimoto, Ting-Chi Huang, Ting-Yi Lu, Alvina Y. L. On, Yi-Hang Valerie Wong, Tiger Yu-Yang Hsiao, Agnieszka Pollo, Matthew A. Malkan, Takamitsu Miyaji, Yoshiki Toba, Ece Kilerci-Eser, Katarzyna Małek, Ho Seong Hwang, Woong-Seob Jeong, Hyunjin Shim, Chris Pearson, Artem Poliszczuk, Bo Han Chen

Abstract: How does the environment affect active galactic nucleus (AGN) activity? We investigated this question in an extinction-free way, by selecting 1120 infrared galaxies in the $AKARI$ North Ecliptic Pole Wide field at redshift $z$ $\leq$ 1.2. A unique feature of the $AKARI$ satellite is its continuous 9-band infrared (IR) filter coverage, providing us with an unprecedentedly large sample of IR spectra… ▽ More How does the environment affect active galactic nucleus (AGN) activity? We investigated this question in an extinction-free way, by selecting 1120 infrared galaxies in the $AKARI$ North Ecliptic Pole Wide field at redshift $z$ $\leq$ 1.2. A unique feature of the $AKARI$ satellite is its continuous 9-band infrared (IR) filter coverage, providing us with an unprecedentedly large sample of IR spectral energy distributions (SEDs) of galaxies. By taking advantage of this, for the first time, we explored the AGN activity derived from SED modelling as a function of redshift, luminosity, and environment. We quantified AGN activity in two ways: AGN contribution fraction (ratio of AGN luminosity to the total IR luminosity), and AGN number fraction (ratio of number of AGNs to the total galaxy sample). We found that galaxy environment (normalised local density) does not greatly affect either definitions of AGN activity of our IRG/LIRG samples (log ${\rm L}_{\rm TIR}$ $\leq$ 12). However, we found a different behavior for ULIRGs (log ${\rm L}_{\rm TIR}$ $>$ 12). At our highest redshift bin (0.7 $\lesssim$ z $\lesssim$ 1.2), AGN activity increases with denser environments, but at the intermediate redshift bin (0.3 $\lesssim$ z $\lesssim$ 0.7), the opposite is observed. These results may hint at a different physical mechanism for ULIRGs. The trends are not statistically significant (p $\geq$ 0.060 at the intermediate redshift bin, and p $\geq$ 0.139 at the highest redshift bin). Possible different behavior of ULIRGs is a key direction to explore further with future space missions (e.g., $JWST$, $Euclid$, $SPHEREx$). △ Less

Submitted 16 August, 2021; originally announced August 2021.

Comments: Accepted for publication in MNRAS. A summary video is available at https://www.youtube.com/watch?v=7Y_a0kJkLI4&ab_channel=NthuCosmology

The Nature of Hyperluminous Infrared Galaxies

Authors: F. Gao, L. Wang, A. Efstathiou, K. Małek, P. N. Best, M. Bonato, D. Farrah, R. Kondapally, I. McCheyne, H. J. A. Röttgering

Abstract: We make use of multi-wavelength data of a large hyperluminous infrared (HLIRG) sample to derive their main physical properties, e.g., stellar mass, star-formation rate (SFR), volume density, contribution to the cosmic stellar mass density and to the cosmic SFR density. We also study the black hole (BH) growth rate and its relationship with the SFR of the host galaxy. We select 526 HLIRGs in three… ▽ More We make use of multi-wavelength data of a large hyperluminous infrared (HLIRG) sample to derive their main physical properties, e.g., stellar mass, star-formation rate (SFR), volume density, contribution to the cosmic stellar mass density and to the cosmic SFR density. We also study the black hole (BH) growth rate and its relationship with the SFR of the host galaxy. We select 526 HLIRGs in three deep fields (Bo$ö$tes, Lockman-Hole, ELAIS-N1) and adopt two spectral energy distribution (SED) fitting codes, CIGALE, which assumes energy balance, and CYGNUS, which is based on radiative transfer models and does not adopt energy balance principle. We use two different active galactic nucleus (AGN) models in CIGALE and three AGN models in CYGNUS to compare the results estimated using different SED fitting codes and different AGN models. The stellar mass, total IR luminosity and AGN luminosity agree well between different models with a typical median offset of 0.1 dex. The SFR estimates show the largest dispersions (up to 0.5 dex). This dispersion has an impact on the subsequent analysis, which may suggest that previous contradictory results could partly be due to different choices of methods. HLIRGs are ultra-massive galaxies with 99% of them having stellar masses larger than $10^{11} M_{\odot}$. Our results reveal a higher space density of ultra-massive galaxies than found in previous surveys or predicted by simulations. We find that HLIRGs contribute more to the cosmic SFR density as redshift increases. In terms of BH growth, the two SED fitting methods provide different results. We can see a clear trend in which SFR decreases as AGN luminosity increases when using CYGNUS estimates, possibly implying quenching by AGN, while this trend is much weaker when using CIGALE estimates. This difference is also influenced by the dispersion between SlFR estimates obtained by the two codes. △ Less

Submitted 2 August, 2021; v1 submitted 19 July, 2021; originally announced July 2021.

Comments: 21 pages, 15 figures. Accepted in A&A

Journal ref: A&A 654, A117 (2021)

Getting ready for the LSST data -- estimating the physical properties of $z<2.5$ main sequence galaxies

Authors: Gabriele Riccio, Katarzyna Małek, Ambra Nanni, Mederic Boquien, Veronique Buat, Denis Burgarella, Darko Donevski, Mahmoud Hamed, Peter Hurley, Raphael Shirley, Agnieszka Pollo

Abstract: In this work we study how to employ the upcoming Legacy Survey of Space and Time (LSST) data to constrain physical properties of normal, star forming galaxies. We use simulated LSST data and existing real observations to test the estimations of the physical properties of galaxies, such as star formation rate (SFR), stellar mass ($M_{star}$), and dust luminosity ($L_{dust}$). We focus on normal sta… ▽ More In this work we study how to employ the upcoming Legacy Survey of Space and Time (LSST) data to constrain physical properties of normal, star forming galaxies. We use simulated LSST data and existing real observations to test the estimations of the physical properties of galaxies, such as star formation rate (SFR), stellar mass ($M_{star}$), and dust luminosity ($L_{dust}$). We focus on normal star-forming galaxies, as they form the majority of the galaxy population in the universe and therefore are more likely to be observed by the LSST. We perform a simulation of LSST observations and uncertainties of 50,385 real galaxies within redshift range $0<z<2.5$. In order to achieve this goal, we used the unique multi-wavelength data from the Herschel Extragalactic Legacy Project (HELP) survey. Our analysis focus on two fields: ELAIS-N1 and COSMOS. To obtain galaxy physical properties we fit their Spectral Energy Distributions (SEDs) using the Code Investigating GALaxy Emission (CIGALE). We compare the main galaxy physical properties obtained from the fit of the observed multi-wavelength photometry of galaxies (from UV to FIR) to the ones obtained from the simulated LSST optical measurements only. The stellar masses estimated based on the LSST measurements are in agreement with the full UV-FIR SED estimations, as they depend mainly on the UV and optical emission, well covered by LSST in the considered redshift range. We obtain a clear overestimation of SFR, $L_{dust}$, $M_{dust}$ estimated with LSST only, highly correlated with redshift. We investigate the cause of this overestimation and we conclude that it is related to an overestimation of the dust attenuation, both UV and NIR. We find that it is necessary to employ auxiliary rest-frame mid-infrared observations, simulated UV observations, or FUV attenuation (AFUV)- Mstar relation, to correct the overestimation. △ Less

Submitted 23 June, 2021; originally announced June 2021.

Comments: Accepted to A&A

Journal ref: A&A 653, A107 (2021)

The role of AGN and its obscuration on the position of the host galaxy relative to Main Sequence

Authors: G. Mountrichas, V. Buat, G. Yang, M. Boquien, D. Burgarella, L. Ciesla, K. Malek, R. Shirley

Abstract: We use X-ray Active Galactic Nuclei (AGN) observed by the Chandra X-ray Observatory within the 9.3 deg$^2$ Bo$\rm \ddot{o}$tes field of the NDWFS to study whether there is a correlation between X-ray luminosity (L$_X$) and star formation rate (SFR) of the host galaxy, at $\rm 0.5<z<2.0$, with respect to the position of the galaxy to the main sequence (SFR$_{norm}$). About half of the sources in th… ▽ More We use X-ray Active Galactic Nuclei (AGN) observed by the Chandra X-ray Observatory within the 9.3 deg$^2$ Bo$\rm \ddot{o}$tes field of the NDWFS to study whether there is a correlation between X-ray luminosity (L$_X$) and star formation rate (SFR) of the host galaxy, at $\rm 0.5<z<2.0$, with respect to the position of the galaxy to the main sequence (SFR$_{norm}$). About half of the sources in the X-ray sample have spectroscopic redshifts. We also construct a reference galaxy catalogue. For both datasets, we use photometric data from optical to the far infrared, compiled by the HELP project and apply spectral energy distribution (SED) fitting, using the X-CIGALE code. We exclude quiescent sources from both the X-ray and the reference samples. We also account for the mass completeness of our dataset, in different redshifts bins. Our analysis highlights the importance of studying the SFR-L$_X$ relation, in a uniform manner, taking into account the systematics and selection effects. Our results suggest that, in less massive galaxies ($\rm log\,[M_*(M_\odot)] \sim 11$), AGN enhances the SFR of the host galaxy by $\sim 50\%$ compared to non AGN systems. A flat relation is observed for the most massive galaxies. SFR$_{norm}$ does not evolve with redshift. The results, although tentative, are consistent with a scenario in which, in less massive systems, both AGN and star formation (SF) are fed by cold gas, supplied by a merger event. In more massive galaxies, the flat relation could be explained by a different SMBH fuelling mechanism that is decoupled from the star formation of the host galaxy (e.g. hot diffuse gas). Finally, we compare the host galaxy properties of X-ray absorbed and unabsorbed sources. Our results show no difference which suggests that X-ray absorption is not linked with the properties of the galaxy. △ Less

Submitted 20 June, 2021; originally announced June 2021.

Comments: Accepted for publication in A&A

Journal ref: A&A 653, A74 (2021)

HELP: The Herschel Extragalactic Legacy Project

Authors: R. Shirley, K. Duncan, M. C. Campos Varillas, P. D. Hurley, K. Malek, Y. Roehlly, M. W. L. Smith, H. Aussel, T. Bakx, V. Buat, D. Burgarella, N. Christopher, S. Duivenvoorden, S. Eales, A. Efstathiou, E. A. Gonzalez Solares, M. Griffin, M. Jarvis, B. Lo Faro, L. Marchetti, I. McCheyne, A. Papadopoulos, K. Penner, E. Pons, M. Prescott , et al. (7 additional authors not shown)

Abstract: We present the Herschel Extragalactic Legacy Project (HELP). This project collates, curates, homogenises, and creates derived data products for most of the premium multi-wavelength extragalactic data sets. The sky boundaries for the first data release cover 1270 deg2 defined by the Herschel SPIRE extragalactic survey fields; notably the Herschel Multi-tiered Extragalactic Survey (HerMES) and the H… ▽ More We present the Herschel Extragalactic Legacy Project (HELP). This project collates, curates, homogenises, and creates derived data products for most of the premium multi-wavelength extragalactic data sets. The sky boundaries for the first data release cover 1270 deg2 defined by the Herschel SPIRE extragalactic survey fields; notably the Herschel Multi-tiered Extragalactic Survey (HerMES) and the Herschel Atlas survey (H-ATLAS). Here, we describe the motivation and principal elements in the design of the project. Guiding principles are transparent or "open" methodologies with care for reproducibility and identification of provenance. A key element of the design focuses around the homogenisation of calibration, meta data and the provision of information required to define the selection of the data for statistical analysis. We apply probabilistic methods that extract information directly from the images at long wavelengths, exploiting the prior information available at shorter wavelengths and providing full posterior distributions rather than maximum likelihood estimates and associated uncertainties as in traditional catalogues. With this project definition paper we provide full access to the first data release of HELP; Data Release 1 (DR1), including a monolithic map of the largest SPIRE extragalactic field at 385 deg2 and 18 million measurements of PACS and SPIRE fluxes. We also provide tools to access and analyse the full HELP database. This new data set includes far-infrared photometry, photometric redshifts, and derived physical properties estimated from modelling the spectral energy distributions. △ Less

Submitted 12 May, 2021; originally announced May 2021.

Active galactic nuclei catalog from the AKARI NEP Wide field

Authors: Artem Poliszczuk, Agnieszka Pollo, Katarzyna Małek, Anna Durkalec, William J. Pearson, Tomotsugu Goto, Seong Jin Kim, Matthew Malkan, Nagisa Oi, Simon C. -C. Ho, Hyunjin Shim, Chris Pearson, Ho Seong Hwang, Yoshiki Toba, Eunbin Kim

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

Submitted 27 April, 2021; originally announced April 2021.

Comments: Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 651, A108 (2021)

A hyperluminous obscured quasar at a redshift of z ~ 4.3

Authors: Andreas Efstathiou, Katarzyna Malek, Denis Burgarella, Peter Hurley, Seb Oliver, Veronique Buat, Raphael Shirley, Steven Duivenvoorden, Vicky Papadopoulou Lesta, Duncan Farrah, Kenneth J. Duncan, Maria del Carmen Campos Varillas

Abstract: In this work we report the discovery of the hyperluminous galaxy HELP_J100156.75+022344.7 at the photometric redshift of z ~ 4.3. The galaxy was discovered in the Cosmological Evolution Survey (COSMOS) field, one of the fields studied by the Herschel Extragalactic Legacy Project (HELP). We present the spectral energy distribution (SED) of the galaxy and fit it with the CYprus models for Galaxies a… ▽ More In this work we report the discovery of the hyperluminous galaxy HELP_J100156.75+022344.7 at the photometric redshift of z ~ 4.3. The galaxy was discovered in the Cosmological Evolution Survey (COSMOS) field, one of the fields studied by the Herschel Extragalactic Legacy Project (HELP). We present the spectral energy distribution (SED) of the galaxy and fit it with the CYprus models for Galaxies and their NUclear Spectra (CYGNUS) multi-component radiative transfer models. We find that its emission is dominated by an obscured quasar with a predicted total 1-1000um luminosity of $3.91^{+1.69}_{-0.55} \times 10^{13} L_\odot$ and an active galactic nucleus (AGN) fraction of ~89%. We also fit HELP_J100156.75+022344.7 with the Code Investigating GALaxy Emission (CIGALE) code and find a similar result. This is only the second z > 4 hyperluminous obscured quasar discovered to date. The discovery of HELP_J100156.75+022344.7 in the ~ 2deg^2 COSMOS field implies that a large number of obscured hyperluminous quasars may lie in the HELP fields which cover ~ 1300deg^2. If this is confirmed, tension between supermassive black hole evolution models and observations will be alleviated. We estimate the space density of objects like HELP_J100156.75+022344.7 at z ~ 4.5 to be $\sim 1.8 \times 10^{-8}$Mpc$^{-3}$. This is slightly higher than the space density of coeval hyperluminous optically selected quasars suggesting that the obscuring torus in z > 4 quasars may have a covering factor $\gtrsim 50\%$. △ Less

Submitted 8 April, 2022; v1 submitted 14 March, 2021; originally announced March 2021.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society Letters. The CYGNUS models used in this work are available at https://arc.euc.ac.cy/cygnus-project-arc/

Journal ref: 2021MNRAS503L..11E

The LOFAR LBA Sky Survey I. survey description and preliminary data release

Authors: F. de Gasperin, W. L. Williams, P. Best, M. Bruggen, G. Brunetti, V. Cuciti, T. J. Dijkema, M. J. Hardcastle, M. J. Norden, A. Offringa, T. Shimwell, R. van Weeren, D. Bomans, A. Bonafede, A. Botteon, J. R. Callingham, R. Cassano, K. T. Chyzy, K. L. Emig, H. Edler, M. Haverkorn, G. Heald, V. Heesen, M. Iacobelli, H. T. Intema , et al. (16 additional authors not shown)

Abstract: LOFAR is the only radio telescope that is presently capable of high-sensitivity, high-resolution (<1 mJy/b and <15") observations at ultra-low frequencies (<100 MHz). To utilise these capabilities, the LOFAR Surveys Key Science Project is undertaking a large survey to cover the entire northern sky with Low Band Antenna (LBA) observations. The LOFAR LBA Sky Survey (LoLSS) aims to cover the entire n… ▽ More LOFAR is the only radio telescope that is presently capable of high-sensitivity, high-resolution (<1 mJy/b and <15") observations at ultra-low frequencies (<100 MHz). To utilise these capabilities, the LOFAR Surveys Key Science Project is undertaking a large survey to cover the entire northern sky with Low Band Antenna (LBA) observations. The LOFAR LBA Sky Survey (LoLSS) aims to cover the entire northern sky with 3170 pointings in the frequency range 42-66 MHz, at a resolution of 15 arcsec and at a sensitivity of 1 mJy/b. Here we outline the survey strategy, the observational status, the current calibration techniques, and briefly describe several scientific motivations. We also describe the preliminary public data release. The preliminary images were produced using a fully automated pipeline that aims to correct all direction-independent effects in the data. Whilst the direction-dependent effects, such as those from the ionosphere, are not yet corrected, the images presented in this work are still 10 times more sensitive than previous surveys available at these low frequencies. The preliminary data release covers 740 sqdeg around the HETDEX spring field region at a resolution of 47" with a median noise level of 5 mJy/b. The images and the catalogue with 25,247 sources are publicly released. We demonstrate that the system is capable of reaching an rms noise of 1 mJy/b and the resolution of 15" once direction-dependent effects are corrected for. LoLSS will provide the ultra-low-frequency information for hundreds of thousands of radio sources, providing critical spectral information and producing a unique dataset that can be used for a wide range of science topics such as: the search for high redshift galaxies and quasars, the study of the magnetosphere of exoplanets, and the detection of the oldest populations of cosmic-rays in galaxies, clusters of galaxies, and from AGN activity. △ Less

Submitted 18 February, 2021; originally announced February 2021.

Comments: 19 pages, 14 figures, Accepted A&A, catalogue and images on www.lofar-surveys.org

Synergies between low- and intermediate-redshift galaxy populations revealed with unsupervised machine learning

Authors: Sebastian Turner, Małgorzata Siudek, Samir Salim, Ivan K. Baldry, Agnieszka Pollo, Steven N. Longmore, Katarzyna Małek, Chris A. Collins, Paulo J. Lisboa, Janusz Krywult, Thibaud Moutard, Daniela Vergani, Alexander Fritz

Abstract: The colour bimodality of galaxies provides an empirical basis for theories of galaxy evolution. However, the balance of processes that begets this bimodality has not yet been constrained. A more detailed view of the galaxy population is needed, which we achieve in this paper by using unsupervised machine learning to combine multi-dimensional data at two different epochs. We aim to understand the c… ▽ More The colour bimodality of galaxies provides an empirical basis for theories of galaxy evolution. However, the balance of processes that begets this bimodality has not yet been constrained. A more detailed view of the galaxy population is needed, which we achieve in this paper by using unsupervised machine learning to combine multi-dimensional data at two different epochs. We aim to understand the cosmic evolution of galaxy subpopulations by uncovering substructures within the colour bimodality. We choose a clustering algorithm that models clusters using only the most discriminative data available, and apply it to two galaxy samples: one from the second edition of the GALEX-SDSS-WISE Legacy Catalogue (GSWLC-2; $z \sim 0.06$), and the other from the VIMOS Public Extragalactic Redshift Survey (VIPERS; $z \sim 0.65$). We cluster within a nine-dimensional feature space defined purely by rest-frame ultraviolet-through-near-infrared colours. Both samples are similarly partitioned into seven clusters, breaking down into four of mostly star-forming galaxies (including the vast majority of green valley galaxies) and three of mostly passive galaxies. The separation between these two families of clusters suggests differences in the evolution of their galaxies, and that these differences are strongly expressed in their colours alone. The samples are closely related, with star-forming/green-valley clusters at both epochs forming morphological sequences, capturing the gradual internally-driven growth of galaxy bulges. At high stellar masses, this growth is linked with quenching. However, it is only in our low-redshift sample that additional, environmental processes appear to be involved in the evolution of low-mass passive galaxies. △ Less

Submitted 6 June, 2021; v1 submitted 9 February, 2021; originally announced February 2021.

Comments: Accepted for publication in MNRAS. 22 pages, 15 figures

Multiwavelength dissection of a massive heavily dust-obscured galaxy and its blue companion at z $\sim $2

Authors: M. Hamed, L. Ciesla, M. Béthermin, K. Małek, E. Daddi, M. T. Sargent, R. Gobat

Abstract: In this work we study a system of two galaxies, Astarte and Adonis, at z $\sim $2 when the Universe was undergoing its peak of star formation activity. Astarte is a dusty star-forming galaxy at the massive-end of the main sequence (MS) and Adonis is a less-massive, bright in ultraviolet (UV), companion galaxy with an optical spectroscopic redshift. We analyse the physical properties of this system… ▽ More In this work we study a system of two galaxies, Astarte and Adonis, at z $\sim $2 when the Universe was undergoing its peak of star formation activity. Astarte is a dusty star-forming galaxy at the massive-end of the main sequence (MS) and Adonis is a less-massive, bright in ultraviolet (UV), companion galaxy with an optical spectroscopic redshift. We analyse the physical properties of this system, and probe the gas mass of Astarte with its ALMA CO emission, to investigate whether this ultra-massive galaxy is quenching or not. We use CIGALE - a spectral energy distribution modeling code - to derive the key physical properties of Astarte and Adonis, mainly their star formation rates (SFRs), stellar masses, and dust luminosities. We inspect the variation of the physical parameters depending on the assumed dust attenuation law. We also estimate the molecular gas mass of Astarte from its CO emission, using different $α_{CO}$ and transition ratios ($r_{31}$) and discuss the implication of the various assumptions on the gas mass derivation. We find that Astarte exhibits a MS-like star formation activity, while Adonis is undergoing a strong starburst (SB) phase. The molecular gas mass of Astarte is far below the gas fraction of typical star-forming galaxies at z=2. This low gas content and high SFR, result in a depletion time of $0.22\pm0.07$ Gyrs, slightly shorter than what is expected for a MS galaxy at this redshift. The CO luminosity versus the total IR luminosity suggests a MS-like activity assuming a galactic conversion factor and a low transition ratio. The SFR of Astarte is of the same order using different attenuation laws, unlike its stellar mass that increases using shallow attenuation laws. We discuss these properties and suggest that Astarte might be experiencing a recent decrease of star formation activity and is quenching through the MS following a SB epoch. △ Less

Submitted 19 January, 2021; originally announced January 2021.

Comments: 10 pages, 7 figures, accepted for publication in A&A. Abstract abridged for arXiv submission

Journal ref: A&A 646, A127 (2021)

The bright end of the infrared luminosity functions and the abundance of hyperluminous infrared galaxies

Authors: L. Wang, F. Gao, P. N. Best, K. Duncan, M. J. Hardcastle, R. Kondapally, K. Malek, I. McCheyne, J. Sabater, T. Shimwell, C. Tasse, M. Bonato, M. Bondi, R. K. Cochrane, D. Farrah, G. Gurkan, P. Haskell, W. J. Pearson, I. Prandoni, H. J. A. Rottgering, D. J. B. Smith, M. Vaccari, W. L. Williams

Abstract: We provide the most accurate estimate yet of the bright end of the infrared (IR) luminosity functions (LFs) and the abundance of hyperluminous IR galaxies (HLIRGs) with IR luminosities > 10^13 L_solar, thanks to the combination of the high sensitivity, angular resolution, and large area of the LOFAR Deep Fields, which probes an unprecedented dynamic range of luminosity and volume. We cross-match H… ▽ More We provide the most accurate estimate yet of the bright end of the infrared (IR) luminosity functions (LFs) and the abundance of hyperluminous IR galaxies (HLIRGs) with IR luminosities > 10^13 L_solar, thanks to the combination of the high sensitivity, angular resolution, and large area of the LOFAR Deep Fields, which probes an unprecedented dynamic range of luminosity and volume. We cross-match Herschel sources and LOFAR sources in Bootes (8.63 deg^2), Lockman Hole (10.28 deg^2), and ELAIS-N1 (6.74 deg^2) with rms sensitivities of around 32, 22, and 20 mJy per beam, respectively. We divide the matched samples into unique and multiple categories. For the multiple matches, we de-blend the Herschel fluxes using the LOFAR positions and the 150-MHz flux densities as priors. We perform spectral energy distribution (SED) fitting, combined with multi-wavelength counterpart identifications and photometric redshift estimates, to derive IR luminosities. The depth of the LOFAR data allows us to identify highly complete (around 92% completeness) samples of bright Herschel sources with a simple selection based on the 250 micron flux density (45, 40, and 35 mJy in Bootes, Lockman Hole, and ELAIS-N1, respectively). Most of the bright Herschel sources fall into the unique category (i.e. a single LOFAR counterpart). For the multiple matches, there is excellent correspondence between the radio emission and the far-IR emission. We find a good agreement in the IR LFs with a previous study out to z around 6 which used de-blended Herschel data. Our sample gives the strongest and cleanest indication to date that the population of HLIRGs has surface densities of around 5 to 18 / deg^2 (with variations due to a combination of the applied flux limit and cosmic variance) and an uncertainty of a factor of 2. In comparison, the GALFORM semi-analytic model significantly under-predicts the abundance of HLIRGs. △ Less

Submitted 17 November, 2020; originally announced November 2020.

Comments: 16 pages, 11 figures, accepted for publication in A&A as part of the LOFAR Deep Fields Paper Splash