-
Probabilistic and progressive deblended far-infrared and sub-millimetre point source catalogues I. Methodology and first application in the COSMOS field
Authors:
Lingyu Wang,
Antonio La Marca,
Fangyou Gao,
William J. Pearson,
Berta Margalef-Bentabol,
Matthieu Béthermin,
Longji Bing,
James Donnellan,
Peter D. Hurley,
Seb J. Oliver,
Catherine L. Hale,
Matt J. Jarvis,
Lucia Marchetti,
Mattia Vaccari,
Imogen H. Whittam
Abstract:
Single-dish far-infrared (far-IR) and sub-millimetre (sub-mm) point source catalogues and their connections with catalogues at other wavelengths are of paramount importance. However, due to the large mismatch in spatial resolution, cross-matching galaxies at different wavelengths is challenging. This work aims to develop the next-generation deblended far-IR and sub-mm catalogues and present the fi…
▽ More
Single-dish far-infrared (far-IR) and sub-millimetre (sub-mm) point source catalogues and their connections with catalogues at other wavelengths are of paramount importance. However, due to the large mismatch in spatial resolution, cross-matching galaxies at different wavelengths is challenging. This work aims to develop the next-generation deblended far-IR and sub-mm catalogues and present the first application in the COSMOS field. Our progressive deblending used the Bayesian probabilistic framework known as XID+. The deblending started from the Spitzer/MIPS 24 micron data, using an initial prior list composed of sources selected from the COSMOS2020 catalogue and radio catalogues from the VLA and the MeerKAT surveys, based on spectral energy distribution modelling which predicts fluxes of the known sources at the deblending wavelength. To speed up flux prediction, we made use of a neural network-based emulator. After deblending the 24 micron data, we proceeded to the Herschel PACS (100 & 160 micron) and SPIRE wavebands (250, 350 & 500 micron). Each time we constructed a tailor-made prior list based on the predicted fluxes of the known sources. Using simulated far-IR and sub-mm sky, we detailed the performance of our deblending pipeline. After validation with simulations, we then deblended the real observations from 24 to 500 micron and compared with blindly extracted catalogues and previous versions of deblended catalogues. As an additional test, we deblended the SCUBA-2 850 micron map and compared our deblended fluxes with ALMA measurements, which demonstrates a higher level of flux accuracy compared to previous results.We publicly release our XID+ deblended point source catalogues. These deblended long-wavelength data are crucial for studies such as deriving the fraction of dust-obscured star formation and better separation of quiescent galaxies from dusty star-forming galaxies.
△ Less
Submitted 28 May, 2024;
originally announced May 2024.
-
Determining the time before or after a galaxy merger event
Authors:
W. J. Pearson,
V. Rodriguez-Gomez,
S. Kruk,
B. Margalef-Bentabol
Abstract:
Aims. This work aims to reproduce the time before or after a merger event of merging galaxies from the IllustrisTNG cosmological simulation using machine learning.
Methods. Images of merging galaxies were created in the u, g, r, and i bands from IllustrisTNG. The merger times were determined using the time difference between the last simulation snapshot where the merging galaxies were tracked as…
▽ More
Aims. This work aims to reproduce the time before or after a merger event of merging galaxies from the IllustrisTNG cosmological simulation using machine learning.
Methods. Images of merging galaxies were created in the u, g, r, and i bands from IllustrisTNG. The merger times were determined using the time difference between the last simulation snapshot where the merging galaxies were tracked as two galaxies and the first snapshot where the merging galaxies were tracked as a single galaxy. This time was then further refined using simple gravity simulations. These data were then used to train a residual network (ResNet50), a Swin Transformer (Swin), a convolutional neural network (CNN), and an autoencoder (using a single latent neuron) to reproduce the merger time. The full latent space of the autoencoder was also studied to see if it reproduces the merger time better than the other methods. This was done by reducing the latent space dimensions using Isomap, linear discriminant analysis (LDA), neighbourhood components analysis, sparse random projection, truncated singular value decomposition and uniform manifold approximation and projection.
Results. The CNN is the best of all the neural networks. The performance of the autoencoder was close to the CNN, with Swin close behind the autoencoder. ResNet50 performed the worst. The LDA dimensionality reduction performed the best of the six methods used. The exploration of the full latent space produced worse results than the single latent neuron of the autoencoder. For the test data set, we found a median error of 190 Myr, comparable to the time separation between snapshots in IllustrisTNG. Galaxies more than $\approx$ 625 Myr before a merger have poorly recovered merger times, as well as galaxies more than $\approx$ 125 Myr after a merger event.
△ Less
Submitted 17 April, 2024;
originally announced April 2024.
-
Galaxy merger challenge: A comparison study between machine learning-based detection methods
Authors:
B. Margalef-Bentabol,
L. Wang,
A. La Marca,
C. Blanco-Prieto,
D. Chudy,
H. Domínguez-Sánchez,
A. D. Goulding,
A. Guzmán-Ortega,
M. Huertas-Company,
G. Martin,
W. J. Pearson,
V. Rodriguez-Gomez,
M. Walmsley,
R. W. Bickley,
C. Bottrell,
C. Conselice,
D. O'Ryan
Abstract:
Various galaxy merger detection methods have been applied to diverse datasets. However, it is difficult to understand how they compare. We aim to benchmark the relative performance of machine learning (ML) merger detection methods. We explore six leading ML methods using three main datasets. The first one (the training data) consists of mock observations from the IllustrisTNG simulations and allow…
▽ More
Various galaxy merger detection methods have been applied to diverse datasets. However, it is difficult to understand how they compare. We aim to benchmark the relative performance of machine learning (ML) merger detection methods. We explore six leading ML methods using three main datasets. The first one (the training data) consists of mock observations from the IllustrisTNG simulations and allows us to quantify the performance metrics of the detection methods. The second one consists of mock observations from the Horizon-AGN simulations, introduced to evaluate the performance of classifiers trained on different, but comparable data. The third one consists of real observations from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) survey. For the binary classification task (mergers vs. non-mergers), all methods perform reasonably well in the domain of the training data. At $0.1<z<0.3$, precision and recall range between $\sim$70\% and 80\%, both of which decrease with increasing $z$ as expected (by $\sim$5\% for precision and $\sim$10\% for recall at $0.76<z<1.0$). When transferred to a different domain, the precision of all classifiers is only slightly reduced, but the recall is significantly worse (by $\sim$20-40\% depending on the method). Zoobot offers the best overall performance in terms of precision and F1 score. When applied to real HSC observations, all methods agree well with visual labels of clear mergers but can differ by more than an order of magnitude in predicting the overall fraction of major mergers. For the multi-class classification task to distinguish pre-, post- and non-mergers, none of the methods offer a good performance, which could be partly due to limitations in resolution and depth of the data. With the advent of better quality data (e.g. JWST and Euclid), it is important to improve our ability to detect mergers and distinguish between merger stages.
△ Less
Submitted 15 April, 2024; v1 submitted 22 March, 2024;
originally announced March 2024.
-
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.
-
Do galaxy mergers prefer under-dense environments?
Authors:
U. Sureshkumar,
A. Durkalec,
A. Pollo,
W. J. Pearson,
D. J. Farrow,
A. Narayanan,
J. Loveday,
E. N. Taylor,
L. E. Suelves
Abstract:
Galaxy mergers play a crucial role in galaxy evolution. However, the correlation between mergers and the local environment of galaxies is not fully understood. We aim to address the question of whether galaxy mergers prefer denser or less dense environments by quantifying the spatial clustering of mergers and non-mergers. We use two different indicators to classify mergers and non-mergers - classi…
▽ More
Galaxy mergers play a crucial role in galaxy evolution. However, the correlation between mergers and the local environment of galaxies is not fully understood. We aim to address the question of whether galaxy mergers prefer denser or less dense environments by quantifying the spatial clustering of mergers and non-mergers. We use two different indicators to classify mergers and non-mergers - classification based on a deep learning technique ($f$) and non-parametric measures of galaxy morphology, Gini-$M_{20}$ ($g$). We used a set of galaxy samples in the redshift range $0.1 < z < 0.15$ from the Galaxy and Mass Assembly (GAMA) survey with a stellar mass cut of $\log (M_{\star}/M_{\odot} ) > 9.5$. We measured and compared the two-point correlation function (2pCF) of mergers and non-mergers classified using the two merger indicators $f$ and $g$. We measured the marked correlation function (MCF), in which the galaxies are weighted by $f$ to probe the environmental dependence of galaxy mergers. We do not observe a statistically significant difference between the clustering strengths of mergers and non-mergers obtained using 2pCF. However, using the MCF measurements with $f$ as a mark, we observe an anti-correlation between the likelihood of a galaxy being a merger and its environment. Our results emphasise the advantage of MCF over 2pCF in probing the environmental correlations. Based on the MCF measurements, we conclude that the galaxy mergers prefer to occur in the under-dense environments on scales $> 50 \, h^{-1} \mathrm{kpc}$ of the large-scale structure (LSS). We attribute this observation to the high relative velocities of galaxies in the densest environments that prevent them from merging.
△ Less
Submitted 30 May, 2024; v1 submitted 28 February, 2024;
originally announced February 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.
-
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.
-
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.
-
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.
-
Pitfalls of AI classification of rare objects: Galaxy Mergers
Authors:
W. J. Pearson,
L. E. Suelves,
S. C. -C. Ho,
N. Oi,
NEP Team,
GAMA Team
Abstract:
Galaxy mergers are hugely important in our current dark matter cosmology. These powerful events cause the disruption of the merging galaxies, pushing the gas, stars and dust of the galaxies resulting in changes to morphologies. This disruption can also cause more extreme events inside the galaxies: periods of extreme star formation rates and the rapid increase in active galactic nuclei activity. H…
▽ More
Galaxy mergers are hugely important in our current dark matter cosmology. These powerful events cause the disruption of the merging galaxies, pushing the gas, stars and dust of the galaxies resulting in changes to morphologies. This disruption can also cause more extreme events inside the galaxies: periods of extreme star formation rates and the rapid increase in active galactic nuclei activity. Hence, to better understand what goes on in these rare events, we need to be able to identify statistically large samples.
In the last few years, the growth of artificial intelligence techniques has seen application to identifying galaxy mergers. These techniques have shown to be highly accurate and their application has grown beyond academic studies of ``can we?'' to deeper scientific use. However, these classifications are not without their problems.
In this proceedings, we will explore how galaxy merger classification can be improved by adding pre-extracted galaxy morphologies alongside the traditional imaging data. This demonstrates that current neural networks are not extracting all the information from the images they are given. It will also explore how the resulting samples of rare objects could be highly contaminated. This has a knock on impact on the upcoming large scale surveys like Euclid and Rubin-LSST.
△ Less
Submitted 22 February, 2023;
originally announced February 2023.
-
Metallicity-PAH Relation of MIR-selected Star-forming Galaxies in AKARI North Ecliptic Pole-wide Survey
Authors:
Hyunjin Shim,
Ho Seong Hwang,
Woong-Seob Jeong,
Yoshiki Toba,
Minjin Kim,
Dohyeong Kim,
Hyunmi Song,
Tetsuya Hashimoto,
Takao Nakagawa,
Ambra Nanni,
William J. Pearson,
Toshinobu Takagi
Abstract:
We investigate the variation in the mid-infrared spectral energy distributions of 373 low-redshift ($z<0.4$) star-forming galaxies, which reflects a variety of polycyclic aromatic hydrocarbon (PAH) emission features. The relative strength of PAH emission is parameterized as $q_\mathrm{PAH}$, which is defined as the mass fraction of PAH particles in the total dust mass. With the aid of continuous m…
▽ More
We investigate the variation in the mid-infrared spectral energy distributions of 373 low-redshift ($z<0.4$) star-forming galaxies, which reflects a variety of polycyclic aromatic hydrocarbon (PAH) emission features. The relative strength of PAH emission is parameterized as $q_\mathrm{PAH}$, which is defined as the mass fraction of PAH particles in the total dust mass. With the aid of continuous mid-infrared photometric data points covering 7-24$μ$m and far-infrared flux densities, $q_\mathrm{PAH}$ values are derived through spectral energy distribution fitting. The correlation between $q_\mathrm{PAH}$ and other physical properties of galaxies, i.e., gas-phase metallicity ($12+\mathrm{log(O/H)}$), stellar mass, and specific star-formation rate (sSFR) are explored. As in previous studies, $q_\mathrm{PAH}$ values of galaxies with high metallicity are found to be higher than those with low metallicity. The strength of PAH emission is also positively correlated with the stellar mass and negatively correlated with the sSFR. The correlation between $q_\mathrm{PAH}$ and each parameter still exists even after the other two parameters are fixed. In addition to the PAH strength, the application of metallicity-dependent gas-to-dust mass ratio appears to work well to estimate gas mass that matches the observed relationship between molecular gas and physical parameters. The result obtained will be used to calibrate the observed PAH luminosity-total infrared luminosity relation, based on the variation of MIR-FIR SED, which is used in the estimation of hidden star formation.
△ Less
Submitted 8 December, 2022;
originally announced December 2022.
-
Merger identification through photometric bands, colours, and their errors
Authors:
L. E. Suelves,
W. J. Pearson,
A. Pollo
Abstract:
Aims. We present the application of a fully connected neural network (NN) for galaxy merger identification using exclusively photometric information. Our purpose is not only to test the method's efficiency, but also to understand what merger properties the NN can learn and what their physical interpretation is. Methods. We created a class-balanced training dataset of 5\,860 galaxies split into mer…
▽ More
Aims. We present the application of a fully connected neural network (NN) for galaxy merger identification using exclusively photometric information. Our purpose is not only to test the method's efficiency, but also to understand what merger properties the NN can learn and what their physical interpretation is. Methods. We created a class-balanced training dataset of 5\,860 galaxies split into mergers and non-mergers. The galaxy observations came from SDSS DR6 and were visually identified in Galaxy Zoo. The 2$\,$930 mergers were selected from known SDSS mergers and the respective non-mergers were the closest match in both redshift and $r$ magnitude. The NN architecture was built by testing a different number of layers with different sizes and variations of the dropout rate. We compared input spaces constructed using: the five SDSS filters: $u$, $g$, $r$, $i$, and $z$; combinations of bands, colours, and their errors; six magnitude types; and variations of input normalization. Results. We find that the fibre magnitude errors contribute the most to the training accuracy. Studying the parameters from which they are calculated, we show that the input space built from the sky error background in the five SDSS bands alone leads to 92.64 $\pm$ 0.15 \% training accuracy. We also find that the input normalization, that is to say, how the data are presented to the NN, has a significant effect on the training performance. Conclusions. We conclude that, from all the SDSS photometric information, the sky error background is the most sensitive to merging processes. This finding is supported by an analysis of its five-band feature space by means of data visualization. Moreover, studying the plane of the $g$ and $r$ sky error bands shows that a decision boundary line is enough to achieve an accuracy of 91.59\%.
△ Less
Submitted 9 December, 2022; v1 submitted 14 November, 2022;
originally announced November 2022.
-
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.
-
Optically-detected galaxy cluster candidates in the $AKARI$ North Ecliptic Pole field based on photometric redshift from Subaru Hyper Suprime-Cam
Authors:
T. -C. Huang,
H. Matsuhara,
T. Goto,
D. J. D. Santos,
S. C. -C. Ho,
S. J. Kim,
T. Hashimoto,
Hiroyuki Ikeda,
Nagisa Oi,
M. A. Malkan,
W. J. Pearson,
A. Pollo,
S. Serjeant,
H. Shim,
T. Miyaji,
H. S. Hwang,
A. Durkalec,
A. Poliszczuk,
T. R. Greve,
C. Pearson,
Y. Toba,
D. Lee,
H. K. Kim,
S. Toft,
W. -S. Jeong
, et al. (1 additional authors not shown)
Abstract:
Galaxy clusters provide an excellent probe in various research fields in astrophysics and cosmology. However, the number of galaxy clusters detected so far in the $AKARI$ North Ecliptic Pole (NEP) field is limited. In this work, we provide galaxy cluster candidates in the $AKARI$ NEP field with the minimum requisites based only on coordinates and photometric redshift (photo-$z$) of galaxies. We us…
▽ More
Galaxy clusters provide an excellent probe in various research fields in astrophysics and cosmology. However, the number of galaxy clusters detected so far in the $AKARI$ North Ecliptic Pole (NEP) field is limited. In this work, we provide galaxy cluster candidates in the $AKARI$ NEP field with the minimum requisites based only on coordinates and photometric redshift (photo-$z$) of galaxies. We used galaxies detected in 5 optical bands ($g$, $r$, $i$, $z$, and $Y$) by the Subaru Hyper Suprime-Cam (HSC), assisted with $u$-band from Canada-France-Hawaii Telescope (CFHT) MegaPrime/MegaCam, and IRAC1 and IRAC2 bands from the $Spitzer$ space telescope for photo-$z$ estimation. We calculated the local density around every galaxy using the 10$^{th}$-nearest neighbourhood. Cluster candidates were determined by applying the friends-of-friends algorithm to over-densities. 88 cluster candidates containing 4390 member galaxies below redshift 1.1 in 5.4 deg$^2$ have been detected. The reliability of our method was examined through false detection tests, redshift uncertainty tests, and applications on the COSMOS data, giving false detection rates of 0.01 to 0.05 and recovery rate of 0.9 at high richness. 3 X-ray clusters previously observed by $ROSAT$ and $Chandra$ were recovered. The cluster galaxies show higher stellar mass and lower star formation rate (SFR) compared to the field galaxies in two-sample Z-tests. These cluster candidates are useful for environmental studies of galaxy evolution and future astronomical surveys in the NEP, where $AKARI$ has performed unique 9-band mid-infrared photometry for tens of thousands of galaxies.
△ Less
Submitted 21 July, 2021;
originally announced July 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.
-
Hidden depths in the local Universe: The Stellar Stream Legacy Survey
Authors:
David Martinez-Delgado,
Andrew P. Cooper,
Javier Roman,
Annalisa Pillepich,
Denis Erkal,
Sarah Pearson,
John Moustakas,
Chervin F. P. Laporte,
Seppo Laine,
Mohammad Akhlaghi,
Dustin Lang,
Dmitry Makarov,
Alejandro S. Borlaff,
Giuseppe Donatiello,
William J. Pearson,
Juan Miro-Carretero,
Jean-Charles Cuillandre,
Helena Dominguez,
Santi Roca-Fabrega,
Carlos S. Frenk,
Judy Schmidt,
Maria A. Gomez-Flechoso,
Rafael Guzman,
Noam I. Libeskind,
Arjun Dey
, et al. (4 additional authors not shown)
Abstract:
Mergers and tidal interactions between massive galaxies and their dwarf satellites are a fundamental prediction of the Lambda-Cold Dark Matter cosmology. These events are thought to provide important observational diagnostics of nonlinear structure formation. Stellar streams in the Milky Way and Andromeda are spectacular evidence for ongoing satellite disruption. However, constructing a statistica…
▽ More
Mergers and tidal interactions between massive galaxies and their dwarf satellites are a fundamental prediction of the Lambda-Cold Dark Matter cosmology. These events are thought to provide important observational diagnostics of nonlinear structure formation. Stellar streams in the Milky Way and Andromeda are spectacular evidence for ongoing satellite disruption. However, constructing a statistically meaningful sample of tidal streams beyond the Local Group has proven a daunting observational challenge, and the full potential for deepening our understanding of galaxy assembly using stellar streams has yet to be realised. Here we introduce the Stellar Stream Legacy Survey, a systematic imaging survey of tidal features associated with dwarf galaxy accretion around a sample of ~3100 nearby galaxies within z~0.02, including about 940 Milky Way analogues. Our survey exploits public deep imaging data from the DESI Legacy Imaging Surveys, which reach surface brightness as faint as ~29 mag/arcsec^2 in the r band. As a proof of concept of our survey, we report the detection and broad-band photometry of 24 new stellar streams in the local Universe. We discuss how these observations can yield new constraints on galaxy formation theory through comparison to mock observations from cosmological galaxy simulations. These tests will probe the present-day mass assembly rate of galaxies, the stellar populations and orbits of satellites, the growth of stellar halos and the resilience of stellar disks to satellite bombardment.
△ Less
Submitted 20 December, 2022; v1 submitted 13 April, 2021;
originally announced April 2021.
-
Galaxy And Mass Assembly: Group and field galaxy morphologies in the star-formation rate - stellar mass plane
Authors:
W. J. Pearson,
L. Wang,
S. Brough,
B. W. Holwerda,
A. M. Hopkins,
J. Loveday
Abstract:
We study the environment in which a galaxy lies (i.e. field or group) and its connection with the morphology of the galaxy. This is done by examining the distribution of parametric and non-parametric statistics across the star-formation rate (SFR) - stellar mass (M$_{\star}$) plane and studying how these distributions change with the environment in the local universe ($z<0.15$).
We determine the…
▽ More
We study the environment in which a galaxy lies (i.e. field or group) and its connection with the morphology of the galaxy. This is done by examining the distribution of parametric and non-parametric statistics across the star-formation rate (SFR) - stellar mass (M$_{\star}$) plane and studying how these distributions change with the environment in the local universe ($z<0.15$).
We determine the concentration (C), Gini, M$_{20}$, asymmetry, Gini-M$_{20}$ bulge statistic (GMB), 50\% light radius ($r_{50}$), total Sérsic index, and bulge Sérsic index ($n_{Bulge}$) for galaxies from the Galaxy and Mass Assembly (GAMA) survey using optical images from the Kilo Degree Survey. We determine the galaxy environment using the GAMA group catalogue and split the galaxies into field or group galaxies. The group galaxies are further divided by the group halo mass (M$_{h}$) - $11\leq\mathrm{log(M}_{h}/\mathrm{M}_\odot)<12$, $12\leq\mathrm{log(M}_{h}/\mathrm{M}_\odot)<13$, and $13\leq\mathrm{log(M}_{h}/\mathrm{M}_\odot)<14$ - and into central and satellite galaxies. The galaxies in each of these samples are then placed onto the SFR-M$_{\star}$ plane, and each parameter is used as a third dimension. We fit the resulting distributions for each parameter in each sample using two two-dimensional Gaussian distributions: one for star-forming galaxies and one for quiescent galaxies. The coefficients of these Gaussian fits are then compared between environments.
Using C and $r_{50}$, we find that galaxies typically become larger as the group mass increases. This change is greater for larger galaxies. There is no indication that galaxies are typically more or less clumpy as the environment changes. Using GMB and $n_{Bulge}$, we see that the star-forming galaxies do not become more bulge or disk dominated as the group mass changes. Asymmetry does not appear to be greatly influenced by environment.
△ Less
Submitted 11 January, 2021;
originally announced January 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.
-
Towards a consistent framework of comparing galaxy mergers in observations and simulations
Authors:
L. Wang,
W. J. Pearson,
V. Rodriguez-Gomez
Abstract:
Aims. We aim to perform consistent comparisons between observations and simulations on the mass dependence of the galaxy major merger fraction at low redshift over an unprecedentedly wide range of stellar masses (10^9 to 10^12 solar masses).
Methods. We first carry out forward modelling of ideal synthetic images of major mergers and non-mergers selected from the Next Generation Illustris Simulat…
▽ More
Aims. We aim to perform consistent comparisons between observations and simulations on the mass dependence of the galaxy major merger fraction at low redshift over an unprecedentedly wide range of stellar masses (10^9 to 10^12 solar masses).
Methods. We first carry out forward modelling of ideal synthetic images of major mergers and non-mergers selected from the Next Generation Illustris Simulations (IllustrisTNG) to include major observational effects. We then train deep convolutional neural networks (CNNs) using realistic mock observations of galaxy samples from the simulations. Subsequently, we apply the trained CNNs to real the Kilo-Degree Survey (KiDS) images of galaxies selected from the Galaxy And Mass Assembly (GAMA) survey. Based on the major merger samples, which are detected in a consistent manner in the observations and simulations, we determine the dependence of major merger fraction on stellar mass at z around 0.15 and make comparisons between the two.
Results. The detected major merger fraction in the GAMA/KiDS observations has a fairly mild decreasing trend with increasing stellar mass over the mass range 10^9 < M_sun < M_star < 10^11.5 M_sun. There is good agreement in the mass dependence of the major merger fraction in the GAMA/KiDS observations and the IllustrisTNG simulations over 10^9.5 M_sun < M_star < 10^10.5 M_sun. However, the observations and the simulations show some differences at M_star > 10^10.5M_sun, possibly due to the supermassive blackhole feedback in its low-accretion state in the simulations which causes a sharp transition in the quenched fractions at this mass scale. The discrepancy could also be due to the relatively small volume of the simulations and/or differences in how stellar masses are measured in simulations and observations.
△ Less
Submitted 7 September, 2020;
originally announced September 2020.
-
Mergers Do Trigger AGNs out to z $\sim$ 0.6
Authors:
F. Gao,
L. Wang,
W. J. Pearson,
Y. A. Gordon,
B. W. Holwerda,
A. M. Hopkins,
M. J. I. Brown,
J. Bland-Hawthorn,
M. S. Owers
Abstract:
Aims. The fueling and feedback of Active Galactic Nuclei (AGNs) are important in understanding the co-evolution between black holes and host galaxies. Mergers are thought to have the capability to bring gas inwards and ignite nuclear activity, especially for more powerful AGNs. However, there is still significant ongoing debate on whether mergers can trigger AGNs and, if they do, whether mergers a…
▽ More
Aims. The fueling and feedback of Active Galactic Nuclei (AGNs) are important in understanding the co-evolution between black holes and host galaxies. Mergers are thought to have the capability to bring gas inwards and ignite nuclear activity, especially for more powerful AGNs. However, there is still significant ongoing debate on whether mergers can trigger AGNs and, if they do, whether mergers are a significant triggering mechanism. Methods. We select a low-redshift ($0.005<z<0.1$) sample from the Sloan Digital Sky Survey (SDSS) and a high-redshift ($0 < z < 0.6$) sample from the Galaxy And Mass Assembly (GAMA) survey. We take advantage of the convolutional neural network technique to identify mergers. We use mid-infrared (MIR) color cut and optical emission line diagnostics to classify AGNs. We also include Low Excitation Radio Galaxies (LERGs) to investigate the connection between mergers and low accretion rate AGNs. Results. We find that AGNs are more likely to be found in mergers than non-mergers, with an AGN excess up to $1.81\pm{0.16}$, suggesting that mergers can trigger AGNs. We also find the fraction of mergers in AGNs is higher than that in non-AGN controls, for both MIR and optically selected AGNs, as well as LERGs, with values between $16.40\pm{0.5}\%$ and $39.23\pm{2.10}\%$, implying a non-negligible to potentially significant role of mergers in triggering AGNs. This merger fraction in AGNs increases as stellar mass increases which supports the idea that mergers are more important for triggering AGNs in more massive galaxies. In terms of merger fraction as a function of AGN power we find a positive trend for MIR selected AGNs and a complex trend for optically selected AGNs, which we interpret under an evolutionary scenario proposed by previous studies. In addition, obscured MIR selected AGNs are more likely to be hosted in mergers than unobscured MIR selected AGNs.
△ Less
Submitted 1 April, 2020;
originally announced April 2020.
-
Effect of galaxy mergers on star formation rates
Authors:
W. J. Pearson,
L. Wang,
M. Alpaslan,
I. Baldry,
M. Bilicki,
M. J. I. Brown,
M. W. Grootes,
B. W. Holwerda,
T. D. Kitching,
S. Kruk,
F. F. S. van der Tak
Abstract:
Galaxy mergers and interactions are an integral part of our basic understanding of how galaxies grow and evolve over time. However, the effect that galaxy mergers have on star formation rates (SFR) is contested, with observations of galaxy mergers showing reduced, enhanced and highly enhanced star formation. We aim to determine the effect of galaxy mergers on the SFR of galaxies using statisticall…
▽ More
Galaxy mergers and interactions are an integral part of our basic understanding of how galaxies grow and evolve over time. However, the effect that galaxy mergers have on star formation rates (SFR) is contested, with observations of galaxy mergers showing reduced, enhanced and highly enhanced star formation. We aim to determine the effect of galaxy mergers on the SFR of galaxies using statistically large samples of galaxies, totalling over 200\,000, over a large redshift range, 0.0 to 4.0. We train and use convolutional neural networks to create binary merger identifications (merger or non-merger) in the SDSS, KiDS and CANDELS imaging surveys. We then compare the galaxy main sequence subtracted SFR of the merging and non-merging galaxies to determine what effect, if any, a galaxy merger has on SFR. We find that the SFR of merging galaxies are not significantly different from the SFR of non-merging systems. The changes in the average SFR seen in the star forming population when a galaxy is merging are small, of the order of a factor of 1.2. However, the higher the SFR above the galaxy main sequence, the higher the fraction of galaxy mergers. Galaxy mergers have little effect on the SFR of the majority of merging galaxies compared to the non-merging galaxies. The typical change in SFR is less than 0.1~dex in either direction. Larger changes in SFR can be seen but are less common. The increase in merger fraction as the distance above the galaxy main sequence increases demonstrates that galaxy mergers can induce starbursts.
△ Less
Submitted 27 August, 2019;
originally announced August 2019.
-
Identifying Galaxy Mergers in Observations and Simulations with Deep Learning
Authors:
W. J. Pearson,
L. Wang,
J. W. Trayford,
C. E. Petrillo,
F. F. S. van der Tak
Abstract:
Mergers are an important aspect of galaxy formation and evolution. We aim to test whether deep learning techniques can be used to reproduce visual classification of observations, physical classification of simulations and highlight any differences between these two classifications. With one of the main difficulties of merger studies being the lack of a truth sample, we can use our method to test b…
▽ More
Mergers are an important aspect of galaxy formation and evolution. We aim to test whether deep learning techniques can be used to reproduce visual classification of observations, physical classification of simulations and highlight any differences between these two classifications. With one of the main difficulties of merger studies being the lack of a truth sample, we can use our method to test biases in visually identified merger catalogues. A convolutional neural network architecture was developed and trained in two ways: one with observations from SDSS and one with simulated galaxies from EAGLE, processed to mimic the SDSS observations. The SDSS images were also classified by the simulation trained network and the EAGLE images classified by the observation trained network. The observationally trained network achieves an accuracy of 91.5% while the simulation trained network achieves 65.2% on the visually classified SDSS and physically classified EAGLE images respectively. Classifying the SDSS images with the simulation trained network was less successful, only achieving an accuracy of 64.6%, while classifying the EAGLE images with the observation network was very poor, achieving an accuracy of only 53.0% with preferential assignment to the non-merger classification. This suggests that most of the simulated mergers do not have conspicuous merger features and visually identified merger catalogues from observations are incomplete and biased towards certain merger types. The networks trained and tested with the same data perform the best, with observations performing better than simulations, a result of the observational sample being biased towards conspicuous mergers. Classifying SDSS observations with the simulation trained network has proven to work, providing tantalizing prospects for using simulation trained networks for galaxy identification in large surveys.
△ Less
Submitted 10 May, 2019; v1 submitted 27 February, 2019;
originally announced February 2019.
-
A multi-wavelength de-blended Herschel view of the statistical properties of dusty star-forming galaxies across cosmic time
Authors:
L. Wang,
W. J. Pearson,
W. Cowley,
J. W. Trayford,
M. Bethermin,
C. Gruppioni,
P. Hurley,
M. J. Michalowski
Abstract:
We aim to study the statistical properties of dusty star-forming galaxies, such as their number counts, luminosity functions (LF) and dust-obscured star-formation rate density (SFRD). We use state-of-the-art de-blended Herschel catalogue in the COSMOS field, generated by combining the Bayesian source extraction tool XID+ and informative prior on the spectral energy distributions, to measure the nu…
▽ More
We aim to study the statistical properties of dusty star-forming galaxies, such as their number counts, luminosity functions (LF) and dust-obscured star-formation rate density (SFRD). We use state-of-the-art de-blended Herschel catalogue in the COSMOS field, generated by combining the Bayesian source extraction tool XID+ and informative prior on the spectral energy distributions, to measure the number counts and LFs at far-infrared (FIR) and sub-millimetre (sub-mm) wavelengths.
Thanks to our de-confusion technique and deep multi-wavelength photometry, we are able to achieve more accurate measurements while probing ten times below the confusion limit. Our number counts at 250 microns agree well with previous Herschel studies. However, our counts at 350 and 500 microns are considerably below previous Herschel results. This is due to previous studies suffering from source confusion which worsens towards longer wavelength. Our number counts at 450 and 870 microns show excellent agreement with previous determinations derived from single dish and interferometric observations. Our measurements of the monochromatic LF and the total IR LF agree well with previous results. The increased dynamic range of our measurements allows us to better measure the faint-end slope of the LF and measure the dust-obscured SFRD out to z~6. We find that the fraction of dust obscured star-formation activity is at its highest around z~1 which then decreases towards both low and high redshift. We do not find a shift of balance between z~3 and z~4 in the cosmic star-formation history from being dominated by unobscured star formation at higher redshift to obscured star formation at lower redshift. However, we do find 3<z<4 to be an interesting transition period as the fraction of the total SFRD that is obscured by dust is significantly lower at higher redshifts.
△ Less
Submitted 25 February, 2019;
originally announced February 2019.
-
Deep Learning for Galaxy Mergers in the Galaxy Main Sequence
Authors:
William J. Pearson,
Lingyu Wang,
James Trayford,
Carlo E. Petrillo,
Floris F. S. van der Tak
Abstract:
Starburst galaxies are often found to be the result of galaxy mergers. As a result, galaxy mergers are often believed to lie above the galaxy main sequence: the tight correlation between stellar mass and star formation rate. Here, we aim to test this claim. Deep learning techniques are applied to images from the Sloan Digital Sky Survey to provide visual-like classifications for over 340 000 objec…
▽ More
Starburst galaxies are often found to be the result of galaxy mergers. As a result, galaxy mergers are often believed to lie above the galaxy main sequence: the tight correlation between stellar mass and star formation rate. Here, we aim to test this claim. Deep learning techniques are applied to images from the Sloan Digital Sky Survey to provide visual-like classifications for over 340 000 objects between redshifts of 0.005 and 0.1. The aim of this classification is to split the galaxy population into merger and non-merger systems and we are currently achieving an accuracy of 91.5%. Stellar masses and star formation rates are also estimated using panchromatic data for the entire galaxy population. With these preliminary data, the mergers are placed onto the full galaxy main sequence, where we find that merging systems lie across the entire star formation rate - stellar mass plane.
△ Less
Submitted 22 January, 2019;
originally announced January 2019.
-
Main sequence of star forming galaxies beyond the Herschel confusion limit
Authors:
W. J. Pearson,
L. Wang,
P. D. Hurley,
K. Małek,
V. Buat,
D. Burgarella,
D. Farrah,
S. J. Oliver,
D. J. B. Smith,
F. F. S. van der Tak
Abstract:
Deep far-infrared (FIR) cosmological surveys are known to be affected by source confusion, causing issues when examining the main sequence (MS) of star forming galaxies. This has typically been partially tackled by the use of stacking. However, stacking only provides the average properties of the objects in the stack. This work aims to trace the MS over $0.2\leq z<6.0$ using the latest de-blended…
▽ More
Deep far-infrared (FIR) cosmological surveys are known to be affected by source confusion, causing issues when examining the main sequence (MS) of star forming galaxies. This has typically been partially tackled by the use of stacking. However, stacking only provides the average properties of the objects in the stack. This work aims to trace the MS over $0.2\leq z<6.0$ using the latest de-blended Herschel photometry, which reaches $\approx10$ times deeper than the 5$σ$ confusion limit in SPIRE. This provides more reliable star formation rates (SFRs), especially for the fainter galaxies, and hence a more reliable MS. We built a pipeline that uses the spectral energy distribution (SED) modelling and fitting tool CIGALE to generate flux density priors in the Herschel SPIRE bands. These priors were then fed into the de-blending tool XID+ to extract flux densities from the SPIRE maps. Multi-wavelength data were combined with the extracted SPIRE flux densities to constrain SEDs and provide stellar mass (M$_{\star}$) and SFRs. These M$_{\star}$ and SFRs were then used to populate the SFR-M$_{\star}$ plane over $0.2\leq z<6.0$. No significant evidence of a high-mass turn-over was found; the best fit is thus a simple two-parameter power law of the form log(SFR)$=α[$log(M$_{\star})-10.5]+β$. The normalisation of the power law increases with redshift, rapidly at $z\lesssim1.8$, from $0.58\pm0.09$ at $z\approx0.37$ to $1.31\pm0.08$ at $z\approx1.8$. The slope is also found to increase with redshift, perhaps with an excess around $1.8\leq z<2.9$. The increasing slope indicates that galaxies become more self-similar as redshift increases, implying that the specific SFR of high-mass galaxies increases over $z=0.2$ to $z=6.0$, becoming closer to that of low-mass galaxies. The excess in the slope at $1.8\leq z<2.9$, if present, coincides with the peak of the cosmic star formation history.
△ Less
Submitted 14 May, 2018; v1 submitted 10 April, 2018;
originally announced April 2018.
-
De-blending Deep Herschel Surveys: A Multi-wavelength Approach
Authors:
W. J. Pearson,
L. Wang,
F. F. S. van der Tak,
P. D. Hurley,
D. Burgarella,
S. J. Oliver
Abstract:
Cosmological surveys in the far infrared are known to suffer from confusion. The Bayesian de-blending tool, XID+, currently provides one of the best ways to de-confuse deep Herschel SPIRE images, using a flat flux density prior. This work is to demonstrate that existing multi-wavelength data sets can be exploited to improve XID+ by providing an informed prior, resulting in more accurate and precis…
▽ More
Cosmological surveys in the far infrared are known to suffer from confusion. The Bayesian de-blending tool, XID+, currently provides one of the best ways to de-confuse deep Herschel SPIRE images, using a flat flux density prior. This work is to demonstrate that existing multi-wavelength data sets can be exploited to improve XID+ by providing an informed prior, resulting in more accurate and precise extracted flux densities. Photometric data for galaxies in the COSMOS field were used to constrain spectral energy distributions (SEDs) using the fitting tool CIGALE. These SEDs were used to create Gaussian prior estimates in the SPIRE bands for XID+. The multi-wavelength photometry and the extracted SPIRE flux densities were run through CIGALE again to allow us to compare the performance of the two priors. Inferred ALMA flux densities (F$^i$), at 870$μ$m and 1250$μ$m, from the best fitting SEDs from the second CIGALE run were compared with measured ALMA flux densities (F$^m$) as an independent performance validation. Similar validations were conducted with the SED modelling and fitting tool MAGPHYS and modified black body functions to test for model dependency. We demonstrate a clear improvement in agreement between the flux densities extracted with XID+ and existing data at other wavelengths when using the new informed Gaussian prior over the original uninformed prior. The residuals between F$^m$ and F$^i$ were calculated. For the Gaussian prior, these residuals, expressed as a multiple of the ALMA error ($σ$), have a smaller standard deviation, 7.95$σ$ for the Gaussian prior compared to 12.21$σ$ for the flat prior, reduced mean, 1.83$σ$ compared to 3.44$σ$, and have reduced skew to positive values, 7.97 compared to 11.50. These results were determined to not be significantly model dependent. This results in statistically more reliable SPIRE flux densities.
△ Less
Submitted 7 April, 2017;
originally announced April 2017.