-
ProPane: Image Warping with Fire
Authors:
A. S. G. Robotham,
R. Tobar,
S. Bellstedt,
S. Casura,
R. H. W. Cook,
J. C. J. D'Silva,
L. J. Davies,
S. P. Driver,
J. Li,
L. P. Garate-Nuñez
Abstract:
In this paper we introduce the software package ProPane, written for the R data analysis language. ProPane combines the full range of wcslib projections with the C++ image manipulation routines provided by the CImg library. ProPane offers routines for image warping and combining (including stacking), and various related tasks such as image alignment tweaking and pixel masking. It can stack an effe…
▽ More
In this paper we introduce the software package ProPane, written for the R data analysis language. ProPane combines the full range of wcslib projections with the C++ image manipulation routines provided by the CImg library. ProPane offers routines for image warping and combining (including stacking), and various related tasks such as image alignment tweaking and pixel masking. It can stack an effectively unlimited number of target frames using multiple parallel cores, and offers threading for many lower level routines. It has been used for a number of current and upcoming large surveys, and we present a range of its capabilities and features. ProPane is already available under a permissive open-source LGPL-3 license at github.com/asgr/ProPane (DOI: 10.5281/zenodo.10057053).
△ Less
Submitted 1 February, 2024; v1 submitted 3 November, 2023;
originally announced November 2023.
-
Quenching massive galaxies across cosmic time with the semi-analytic model SHARK v2.0
Authors:
Claudia D. P. Lagos,
Matias Bravo,
Rodrigo Tobar,
Danail Obreschkow,
Chris Power,
Aaron S. G. Robotham,
Katy L. Proctor,
Samuel Hansen,
Angel Chandro-Gomez,
Julian Carrivick
Abstract:
We introduce version 2.0 of the SHARK semi-analytic model of galaxy formation after many improvements to the physics included. The most significant being: (i) a model describing the exchange of angular momentum (AM) between the interstellar medium and stars; (ii) a new active galactic nuclei feedback model which has two modes, a wind and a jet mode, with the jet mode tied to the jet energy product…
▽ More
We introduce version 2.0 of the SHARK semi-analytic model of galaxy formation after many improvements to the physics included. The most significant being: (i) a model describing the exchange of angular momentum (AM) between the interstellar medium and stars; (ii) a new active galactic nuclei feedback model which has two modes, a wind and a jet mode, with the jet mode tied to the jet energy production; (iii) a model tracking the development of black hole (BH) spins; (iv) more sophisticated modelling of environmental effects on satellite galaxies; and (v) automatic parameter exploration using Particle Swarm Optimisation. We focus on two timely research topics: the structural properties of galaxies and the quenching of massive galaxies. For the former, SHARK v2.0 is capable of producing a more realistic stellar size-mass relation with a plateau marking the transition from disk- to bulge-dominated galaxies, and scaling relations between specific AM and mass that agree well with observations. For the quenching of massive galaxies, SHARK v2.0 produces massive galaxies that are more quenched than the previous version, reproducing well the observed relations between star formation rate (SFR) and stellar mass, and specific SFR and BH mass at $z=0$. SHARK v2.0 produces a number density of massive-quiescent galaxies >1dex higher than the previous version, in good agreement with JWST observations at $z\le 5$; predicts a stellar mass function of passive galaxies in reasonably good agreement with observations at $0.5<z<5$; and environmental quenching to already be effective at $z=5$.
△ Less
Submitted 15 April, 2024; v1 submitted 5 September, 2023;
originally announced September 2023.
-
SimSpin v2.6.0 -- Constructing synthetic spectral IFU cubes for comparison with observational surveys
Authors:
K. E. Harborne,
A. Serene,
E. J. A. Davies,
C. Derkenne,
S. Vaughan,
A. I. Burdon,
C. del P. Lagos,
R. McDermid,
S. O'Toole,
C. Power,
A. S. G. Robotham,
G. Santucci,
R. Tobar
Abstract:
In this work, we present a methodology and a corresponding code-base for constructing mock integral field spectrograph (IFS) observations of simulated galaxies in a consistent and reproducible way. Such methods are necessary to improve the collaboration and comparison of observation and theory results, and accelerate our understanding of how the kinematics of galaxies evolve over time. This code,…
▽ More
In this work, we present a methodology and a corresponding code-base for constructing mock integral field spectrograph (IFS) observations of simulated galaxies in a consistent and reproducible way. Such methods are necessary to improve the collaboration and comparison of observation and theory results, and accelerate our understanding of how the kinematics of galaxies evolve over time. This code, SimSpin, is an open-source package written in R, but also with an API interface such that the code can be interacted with in any coding language. Documentation and individual examples can be found at the open-source website connected to the online repository. SimSpin is already being utilised by international IFS collaborations, including SAMI and MAGPI, for generating comparable data sets from a diverse suite of cosmological hydrodynamical simulations.
△ Less
Submitted 30 August, 2023; v1 submitted 5 July, 2023;
originally announced July 2023.
-
An Empirical Evaluation On the Applicability of the DALiuGE Execution Framework
Authors:
Ying Mei,
Shoulin Wei,
Feng Wang,
Chen Wu,
Rodrigo Tobar,
Mohsim Shaikh,
Hui Deng,
Wei Dai,
Bo Liang,
Andreas Wicenec
Abstract:
The Square Kilometre Array (SKA) project is an international cooperation project to build the largest radio telescope worldwide. Data processing is one of the biggest challenges of building the SKA telescope. As a distributed execution framework, the Data Activated Liu Graph Engine (DALiuGE) was proposed to be one of the candidates for addressing the massive data of the SKA. DALiuGE has many disti…
▽ More
The Square Kilometre Array (SKA) project is an international cooperation project to build the largest radio telescope worldwide. Data processing is one of the biggest challenges of building the SKA telescope. As a distributed execution framework, the Data Activated Liu Graph Engine (DALiuGE) was proposed to be one of the candidates for addressing the massive data of the SKA. DALiuGE has many distinctive features, but its actual ability to handle scientific data is still not evident. In this paper, we perform an objective evaluation of the usability of DALiuGE concerning the execution performance, developer workload, and implementation difficulty of porting the SAGECal to DALiuGE. The evaluation results showed that the DALiuGE enables fast integration of astronomical software, but there are significant differences in the efficiency of different parallel granularities. Even with the deep optimization of the program, there is still a gap between the current DALiuGE and the traditional MPI in execution performance. Therefore, we come to a preliminary conclusion that the DALiuGE has no performance advantage in batch processing of massive data, while it may be more suitable for application scenarios with more customized computational tasks, such as SKA science regional centers.
△ Less
Submitted 24 December, 2021;
originally announced December 2021.
-
Physical properties and evolution of (Sub-)millimeter selected galaxies in the galaxy formation simulation Shark
Authors:
Claudia del P. Lagos,
Elisabete da Cunha,
Aaron S. G. Robotham,
Danail Obreschkow,
Francesco Valentino,
Seiji Fujimoto,
Georgios E. Magdis,
Rodrigo Tobar
Abstract:
We thoroughly explore the properties of (sub)-millimeter (mm) selected galaxies (SMGs) in the Shark semi-analytic model of galaxy formation. Compared to observations, the predicted number counts at wavelengths (lambda) 0.6-2mm and redshift distributions at 0.1-2mm, agree well. At the bright end (>1mJy), Shark galaxies are a mix of mergers and disk instabilities. These galaxies display a stacked FU…
▽ More
We thoroughly explore the properties of (sub)-millimeter (mm) selected galaxies (SMGs) in the Shark semi-analytic model of galaxy formation. Compared to observations, the predicted number counts at wavelengths (lambda) 0.6-2mm and redshift distributions at 0.1-2mm, agree well. At the bright end (>1mJy), Shark galaxies are a mix of mergers and disk instabilities. These galaxies display a stacked FUV-to-FIR spectrum that agrees well with observations. We predict that current optical/NIR surveys are deep enough to detect bright (>1mJy) lambda=0.85-2mm-selected galaxies at z<5, but too shallow to detect counterparts at higher redshift. A James Webb Space Telescope 10,000s survey should detect all counterparts for galaxies with $S_{\rm 0.85mm}>0.01$mJy. We predict SMG's disks contribute significantly (negligibly) to the rest-frame UV (IR). We investigate the 0<z<6 evolution of the intrinsic properties of >1mJy lambda=0.85-2mm-selected galaxies finding their: (i) stellar masses are $>10^{10.2}M_{\odot}$, with the 2mm ones tracing the most massive galaxies ($>10^{11}M_{\odot}$); (ii) specific star formation rates (SFR) are mildly (~3-10x) above the main sequence (MS); (iii) host halo masses are $\gtrsim 10^{12.3}M_{\odot}$, with 2mm galaxies tracing the most massive halos (proto-clusters); (iv) SMGs have lower dust masses ($\approx 10^{8}M_{\odot}$), higher dust temperatures ($\approx 40-45$K) and higher rest-frame V-band attenuation (>1.5) than MS galaxies; (v) sizes decrease with redshift, from 4kpc at z=1 to <1kpc at z=4; (vi) the Carbon Monoxide line spectra of $S_{\rm 0.85mm}>1$mJy sources peak at 4->3. Finally, we study the contribution of SMGs to the molecular gas and cosmic SFR density at 0<z<10, finding that >1mJy sources make a negligible contribution at z>3 and z>5, respectively, suggesting current observations have unveiled the majority of the star formation at 0<z<10.
△ Less
Submitted 13 October, 2020; v1 submitted 19 July, 2020;
originally announced July 2020.
-
From the far-ultraviolet to the far-infrared -- galaxy emission at $0\le z \le 10$ in the Shark semi-analytic model
Authors:
Claudia del P. Lagos,
Aaron S. G. Robotham,
James W. Trayford,
Rodrigo Tobar,
Matías Bravo,
Sabine Bellstedt,
Luke J. M. Davies,
Simon P. Driver,
Pascal J. Elahi,
Danail Obreschkow,
Chris Power
Abstract:
We combine the Shark semi-analytic model of galaxy formation with the ProSpect software tool for spectral energy distribution (SED) generation to study the multi-wavelength emission of galaxies from the far-ultraviolet (FUV) to the far-infrared (FIR) at $0\le z\le 10$. We produce a physical model for the attenuation of galaxies across cosmic time by combining a local Universe empirical relation to…
▽ More
We combine the Shark semi-analytic model of galaxy formation with the ProSpect software tool for spectral energy distribution (SED) generation to study the multi-wavelength emission of galaxies from the far-ultraviolet (FUV) to the far-infrared (FIR) at $0\le z\le 10$. We produce a physical model for the attenuation of galaxies across cosmic time by combining a local Universe empirical relation to compute the dust mass of galaxies from their gas metallicity and mass, attenuation curves derived from radiative transfer calculations of galaxies in the EAGLE hydrodynamic simulation suite, and the properties of Shark galaxies. We are able to produce a wide range of galaxies, from the $z=8$ star-forming galaxies with almost no extinction, $z=2$ submillimeter galaxies, down to the normal star-forming and red sequence galaxies at $z=0$. Quantitatively, we find that Shark reproduces the observed (i) the $z=0$ FUV-to-FIR, (ii) $0\le z\le 3$ rest-frame $K$-band, and (iii) $0\le z\le 10$ rest-frame FUV luminosity functions, (iv) $z\le 8$ UV slopes, (v) the FUV-to-FIR number counts (including the widely disputed 850$μ$m), (vi) redshift distribution of bright $850μ$m galaxies and (vii) the integrated cosmic SED from $z=0$ to $z=1$ to an unprecedented level. This is achieved without the need to invoke changes in the stellar initial mass function, dust-to-metal mass ratio, or metal enrichment timescales. Our model predicts star formation in galaxy disks to dominate in the FUV-to-optical, while bulges dominate at the NIR at all redshifts. The FIR sees a strong evolution in which disks dominate at $z\le 1$ and starbursts (triggered by both galaxy mergers and disk instabilities, in an even mix) dominate at higher redshifts, even out to $z=10$.
△ Less
Submitted 29 August, 2019; v1 submitted 9 August, 2019;
originally announced August 2019.
-
Climbing Halo Merger Trees with TreeFrog
Authors:
Pascal J. Elahi,
Rhys J. J. Poulton,
Rodrigo J. Tobar,
Rodrigo Canas,
Claudia del P. Lagos,
Chris Power,
Aaron S. G. Robotham
Abstract:
We present TreeFrog, a massively parallel halo merger tree builder that is capable comparing different halo catalogues and producing halo merger trees. The code is written in c++11, use the MPI and OpenMP API's for parallelisation, and includes python tools to read/manipulate the data products produced. The code correlates binding energy sorted particle ID lists between halo catalogues, determinin…
▽ More
We present TreeFrog, a massively parallel halo merger tree builder that is capable comparing different halo catalogues and producing halo merger trees. The code is written in c++11, use the MPI and OpenMP API's for parallelisation, and includes python tools to read/manipulate the data products produced. The code correlates binding energy sorted particle ID lists between halo catalogues, determining optimal descendant/progenitor matches using multiple snapshots, a merit function that maximises the number of shared particles using pseudo-radial moments, and a scheme for correcting halo merger tree pathologies. Focusing on VELOCIraptor catalogues for this work, we demonstrate how searching multiple snapshots spanning a dynamical time significantly reduces the number of stranded halos, those lacking a descendant or a progenitor, critically correcting poorly resolved halos. We present a new merit function that improves the distinction between primary and secondary progenitors, reducing tree pathologies. We find FOF accretion rates and merger rates show similar mass ratio dependence. The model merger rates from Poole et al, (2017) agree with the measured net growth of halos through mergers.
△ Less
Submitted 7 May, 2019; v1 submitted 4 February, 2019;
originally announced February 2019.
-
Hunting for Galaxies and Halos in simulations with VELOCIraptor
Authors:
Pascal J. Elahi,
Rodrigo Cañas,
Rhys J. J. Poulton,
Rodrigo J. Tobar,
James S. Willis,
Claudia del P. Lagos,
Chris Power,
Aaron S. G. Robotham
Abstract:
We present VELOCIraptor, a massively parallel galaxy/(sub)halo finder that is also capable of robustly identifying tidally disrupted objects and separate stellar halos from galaxies. The code is written in c++11, use the MPI and OpenMP API's for parallelisation, and includes python tools to read/manipulate the data products produced. We demonstrate the power of the VELOCIraptor (sub)halo finder, s…
▽ More
We present VELOCIraptor, a massively parallel galaxy/(sub)halo finder that is also capable of robustly identifying tidally disrupted objects and separate stellar halos from galaxies. The code is written in c++11, use the MPI and OpenMP API's for parallelisation, and includes python tools to read/manipulate the data products produced. We demonstrate the power of the VELOCIraptor (sub)halo finder, showing how it can identify subhalos deep within the host that have negligible density contrasts to their parent halo. We find a subhalo mass-radial distance dependence: large subhalos with mass ratios of $\gtrsim10^{-2}$ are more common in the central regions that smaller subhalos, a result of dynamical friction and low tidal mass loss rates. This dependence is completely absent in (sub)halo finders in common use, which generally search for substructure in configuration space, yet is present in codes that track particles belonging to halos as they fall into other halos, such as HBT+. VELOCIraptor largely reproduces the dependence seen without tracking, finding a similar radial dependence to HBT+ in well resolved halos from our limited resolution fiducial simulation.
△ Less
Submitted 28 March, 2019; v1 submitted 3 February, 2019;
originally announced February 2019.
-
Shark: introducing an open source, free and flexible semi-analytic model of galaxy formation
Authors:
Claudia del P. Lagos,
Rodrigo J. Tobar,
Aaron S. G. Robotham,
Danail Obreschkow,
Peter D. Mitchell,
Chris Power,
Pascal J. Elahi
Abstract:
We present a new, open source, free semi-analytic model (SAM) of galaxy formation, Shark, designed to be highly flexible and modular, allowing easy exploration of different physical processes and ways of modelling them. We introduce the philosophy behind Shark and provide an overview of the physical processes included in the model. Shark is written in C++11 and has been parallelized with OpenMP. I…
▽ More
We present a new, open source, free semi-analytic model (SAM) of galaxy formation, Shark, designed to be highly flexible and modular, allowing easy exploration of different physical processes and ways of modelling them. We introduce the philosophy behind Shark and provide an overview of the physical processes included in the model. Shark is written in C++11 and has been parallelized with OpenMP. In the released version (v1.1), we implement several different models for gas cooling, active galactic nuclei, stellar and photo-ionisation feedback, and star formation (SF). We demonstrate the basic performance of Shark using the Planck15 cosmology SURFS simulations, by comparing against a large set of observations, including: the stellar mass function (SMF) and stellar-halo mass relation at z=0-4; the cosmic evolution of the star formation rate density (SFRD), stellar mass, atomic and molecular hydrogen; local gas scaling relations; and structural galaxy properties, finding excellent agreement. Significant improvements over previous SAMs are seen in the mass-size relation for disks/bulges, the gas-stellar mass and stellar mass-metallicity relations. To illustrate the power of Shark in exploring the systematic effects of the galaxy formation modelling, we quantify how the scatter of the SF main sequence and the gas scaling relations changes with the adopted SF law, and the effect of the starbursts H$_2$ depletion timescale on the SFRD and $Ω_{\rm H_2}$. We compare Shark with other SAMs and the hydrodynamical simulation EAGLE, and find that SAMs have a much higher halo baryon fractions due to large amounts of intra-halo gas, which in the case of EAGLE is in the intergalactic medium.
△ Less
Submitted 24 September, 2018; v1 submitted 30 July, 2018;
originally announced July 2018.
-
ProFit: Bayesian Profile Fitting of Galaxy Images
Authors:
A. S. G. Robotham,
D. S. Taranu,
R. Tobar,
A. Moffett,
S. P. Driver
Abstract:
We present ProFit, a new code for Bayesian two-dimensional photometric galaxy profile modelling. ProFit consists of a low-level C++ library (libprofit), accessible via a command-line interface and documented API, along with high-level R (ProFit) and Python (PyProFit) interfaces (available at github.com/ICRAR/ libprofit, github.com/ICRAR/ProFit, and github.com/ICRAR/pyprofit respectively). R ProFit…
▽ More
We present ProFit, a new code for Bayesian two-dimensional photometric galaxy profile modelling. ProFit consists of a low-level C++ library (libprofit), accessible via a command-line interface and documented API, along with high-level R (ProFit) and Python (PyProFit) interfaces (available at github.com/ICRAR/ libprofit, github.com/ICRAR/ProFit, and github.com/ICRAR/pyprofit respectively). R ProFit is also available pre-built from CRAN, however this version will be slightly behind the latest GitHub version. libprofit offers fast and accurate two- dimensional integration for a useful number of profiles, including Sersic, Core-Sersic, broken-exponential, Ferrer, Moffat, empirical King, point-source and sky, with a simple mechanism for adding new profiles. We show detailed comparisons between libprofit and GALFIT. libprofit is both faster and more accurate than GALFIT at integrating the ubiquitous Serrsic profile for the most common values of the Serrsic index n (0.5 < n < 8).
The high-level fitting code ProFit is tested on a sample of galaxies with both SDSS and deeper KiDS imaging. We find good agreement in the fit parameters, with larger scatter in best-fit parameters from fitting images from different sources (SDSS vs KiDS) than from using different codes (ProFit vs GALFIT). A large suite of Monte Carlo-simulated images are used to assess prospects for automated bulge-disc decomposition with ProFit on SDSS, KiDS and future LSST imaging. We find that the biggest increases in fit quality come from moving from SDSS- to KiDS-quality data, with less significant gains moving from KiDS to LSST.
△ Less
Submitted 17 January, 2017; v1 submitted 22 November, 2016;
originally announced November 2016.
-
Reddening and metallicity maps of the Milky Way bulge from VVV and 2MASS III. The first global photometric metallicity map of the Galactic bulge
Authors:
O. A. Gonzalez,
M. Rejkuba,
M. Zoccali,
E. Valenti,
D. Minniti,
R. Tobar
Abstract:
We investigate the large scale metallicity distribution in the Galactic bulge, using a large spatial coverage, in order to constrain the bulge formation scenario. We use the VISTA variables in the Via Lactea (VVV) survey data and 2MASS photometry, covering 320 sqdeg of the Galactic bulge, to derive photometric metallicities by interpolating of the (J-Ks)0 colors of individual Red Giant Branch star…
▽ More
We investigate the large scale metallicity distribution in the Galactic bulge, using a large spatial coverage, in order to constrain the bulge formation scenario. We use the VISTA variables in the Via Lactea (VVV) survey data and 2MASS photometry, covering 320 sqdeg of the Galactic bulge, to derive photometric metallicities by interpolating of the (J-Ks)0 colors of individual Red Giant Branch stars based on a set of globular cluster ridge lines. We then use this information to construct the first global metallicity map of the bulge with a resolution of 30'x45'. The metallicity map of the bulge revealed a clear vertical metallicity gradient of ~0.04 dex/deg (~0.28 dex/kpc), with metal-rich stars ([Fe/H]~0) dominating the inner bulge in regions closer to the galactic plane (|b|<5). At larger scale heights, the mean metallicity of the bulge population becomes significantly more metal-poor. This fits in the scenario of a boxy-bulge originated from the vertical inestability of the Galactic bar, formed early via secular evolution of a two component stellar disk. Older, metal-poor stars dominate at higher scale heights due to the non-mixed orbits from the originally hotter thick disk stars.
△ Less
Submitted 1 February, 2013;
originally announced February 2013.
-
Reddening and metallicity maps of the Milky Way bulge from VVV and 2MASS II. The complete high resolution extinction map and implications for Bulge studies
Authors:
O. A. Gonzalez,
M. Rejkuba,
M. Zoccali,
E. Valenti,
D. Minniti,
M. Schultheis,
R. Tobar,
B. Chen
Abstract:
We use the Vista Variables in the Via Lactea (VVV) ESO public survey data to measure extinction values in the complete area of the Galactic bulge covered by the survey at high resolution. We derive reddening values using the method described in Paper I. This is based on measuring the mean (J-Ks) color of red clump giants in small subfields of 2' to 6' in the following bulge area: -10.3<b<+5.1 and…
▽ More
We use the Vista Variables in the Via Lactea (VVV) ESO public survey data to measure extinction values in the complete area of the Galactic bulge covered by the survey at high resolution. We derive reddening values using the method described in Paper I. This is based on measuring the mean (J-Ks) color of red clump giants in small subfields of 2' to 6' in the following bulge area: -10.3<b<+5.1 and -10<l<+10.4. To determine the reddening values E(J-Ks) for each region, we measure the RC color and compare it to the (J-Ks) color of RC stars measured in Baade's window, for which we adopt E(B-V)=0.55. This allows us to construct a reddening map sensitive to small scale variations minimizing the problems arising from differential extinction. The significant reddening variations are clearly observed on spatial scales as small as 2'. We find a good agreement between our extinction measurements and Schlegel maps in the outer bulge, but, as already stated in the literature the Schlegel maps are not reliable for regions within |b| < 6. In the inner regions we compare our results with maps derived from DENIS and Spitzer surveys. While we find good agreement with other studies in the corresponding overlapping regions, our extinction map has better quality due to both higher resolution and a more complete spatial coverage in the Bulge. We investigate the importance of differential reddening and demonstrate the need for high resolution extinction maps for detailed studies of Bulge stellar populations and structure. The extinction variations on scales of up to 2'-6', must be taken into account when analysing the stellar populations of the Bulge.
△ Less
Submitted 18 April, 2012;
originally announced April 2012.