Showing 1–22 of 22 results for author: Zanisi, L

Valid Error Bars for Neural Weather Models using Conformal Prediction

Authors: Vignesh Gopakumar, Joel Oskarrson, Ander Gray, Lorenzo Zanisi, Stanislas Pamela, Daniel Giles, Matt Kusner, Marc Deisenroth

Abstract: Neural weather models have shown immense potential as inexpensive and accurate alternatives to physics-based models. However, most models trained to perform weather forecasting do not quantify the uncertainty associated with their forecasts. This limits the trust in the model and the usefulness of the forecasts. In this work we construct and formalise a conformal prediction framework as a post-pro… ▽ More Neural weather models have shown immense potential as inexpensive and accurate alternatives to physics-based models. However, most models trained to perform weather forecasting do not quantify the uncertainty associated with their forecasts. This limits the trust in the model and the usefulness of the forecasts. In this work we construct and formalise a conformal prediction framework as a post-processing method for estimating this uncertainty. The method is model-agnostic and gives calibrated error bounds for all variables, lead times and spatial locations. No modifications are required to the model and the computational cost is negligible compared to model training. We demonstrate the usefulness of the conformal prediction framework on a limited area neural weather model for the Nordic region. We further explore the advantages of the framework for deterministic and probabilistic models. △ Less

Submitted 20 June, 2024; originally announced June 2024.

Neural-Parareal: Dynamically Training Neural Operators as Coarse Solvers for Time-Parallelisation of Fusion MHD Simulations

Authors: S. J. P. Pamela, N. Carey, J. Brandstetter, R. Akers, L. Zanisi, J. Buchanan, V. Gopakumar, M. Hoelzl, G. Huijsmans, K. Pentland, T. James, G. Antonucci, the JOREK Team

Abstract: The fusion research facility ITER is currently being assembled to demonstrate that fusion can be used for industrial energy production, while several other programmes across the world are also moving forward, such as EU-DEMO, CFETR, SPARC and STEP. The high engineering complexity of a tokamak makes it an extremely challenging device to optimise, and test-based optimisation would be too slow and to… ▽ More The fusion research facility ITER is currently being assembled to demonstrate that fusion can be used for industrial energy production, while several other programmes across the world are also moving forward, such as EU-DEMO, CFETR, SPARC and STEP. The high engineering complexity of a tokamak makes it an extremely challenging device to optimise, and test-based optimisation would be too slow and too costly. Instead, digital design and optimisation must be favored, which requires strongly-coupled suites of High-Performance Computing calculations. In this context, having surrogate models to provide quick estimates with uncertainty quantification is essential to explore and optimise new design options. Furthermore, these surrogates can in turn be used to accelerate simulations in the first place. This is the case of Parareal, a time-parallelisation method that can speed-up large HPC simulations, where the coarse-solver can be replaced by a surrogate. A novel framework, Neural-Parareal, is developed to integrate the training of neural operators dynamically as more data becomes available. For a given input-parameter domain, as more simulations are being run with Parareal, the large amount of data generated by the algorithm is used to train new surrogate models to be used as coarse-solvers for future Parareal simulations, leading to progressively more accurate coarse-solvers, and thus higher speed-up. It is found that such neural network surrogates can be much more effective than traditional coarse-solver in providing a speed-up with Parareal. This study is a demonstration of the convergence of HPC and AI which simply has to become common practice in the world of digital engineering design. △ Less

Submitted 2 May, 2024; originally announced May 2024.

Data efficiency and long term prediction capabilities for neural operator surrogate models of core and edge plasma codes

Authors: N. Carey, L. Zanisi, S. Pamela, V. Gopakumar, J. Omotani, J. Buchanan, J. Brandstetter

Abstract: Simulation-based plasma scenario development, optimization and control are crucial elements towards the successful deployment of next-generation experimental tokamaks and Fusion power plants. Current simulation codes require extremely intensive use of HPC resources that make them unsuitable for iterative or real time applications. Neural network based surrogate models of expensive simulators have… ▽ More Simulation-based plasma scenario development, optimization and control are crucial elements towards the successful deployment of next-generation experimental tokamaks and Fusion power plants. Current simulation codes require extremely intensive use of HPC resources that make them unsuitable for iterative or real time applications. Neural network based surrogate models of expensive simulators have been proposed to speed up such costly workflows. Current efforts in this direction in the Fusion community are mostly limited to point estimates of quantities of interest or simple 1D PDE models, with a few notable exceptions. While the AI literature on methods for neural PDE surrogate models is rich, performance benchmarks for Fusion-relevant 2D fields has so far remained flimited. In this work neural PDE surrogates are trained for the JOREK MHD code and the STORM scrape-off layer code using the PDEArena library (https://github.com/microsoft/pdearena). The performance of these surrogate models is investigated as a function of training set size as well as for long-term predictions. The performance of surrogate models that are trained on either one variable or multiple variables at once is also considered. It is found that surrogates that are trained on more data perform best for both long- and short-term predictions. Additionally, surrogate models trained on multiple variables achieve higher accuracy and more stable performance. Downsampling the training set in time may provide stability in the long term at the expense of the short term predictive capability, but visual inspection of the resulting fields suggests that multiple metrics should be used to evaluate performance. △ Less

Submitted 13 February, 2024; originally announced February 2024.

Comments: IAEA-FEC 2023 manuscript

Plasma Surrogate Modelling using Fourier Neural Operators

Authors: Vignesh Gopakumar, Stanislas Pamela, Lorenzo Zanisi, Zongyi Li, Ander Gray, Daniel Brennand, Nitesh Bhatia, Gregory Stathopoulos, Matt Kusner, Marc Peter Deisenroth, Anima Anandkumar, JOREK Team, MAST Team

Abstract: Predicting plasma evolution within a Tokamak reactor is crucial to realizing the goal of sustainable fusion. Capabilities in forecasting the spatio-temporal evolution of plasma rapidly and accurately allow us to quickly iterate over design and control strategies on current Tokamak devices and future reactors. Modelling plasma evolution using numerical solvers is often expensive, consuming many hou… ▽ More Predicting plasma evolution within a Tokamak reactor is crucial to realizing the goal of sustainable fusion. Capabilities in forecasting the spatio-temporal evolution of plasma rapidly and accurately allow us to quickly iterate over design and control strategies on current Tokamak devices and future reactors. Modelling plasma evolution using numerical solvers is often expensive, consuming many hours on supercomputers, and hence, we need alternative inexpensive surrogate models. We demonstrate accurate predictions of plasma evolution both in simulation and experimental domains using deep learning-based surrogate modelling tools, viz., Fourier Neural Operators (FNO). We show that FNO has a speedup of six orders of magnitude over traditional solvers in predicting the plasma dynamics simulated from magnetohydrodynamic models, while maintaining a high accuracy (MSE in the normalised domain $\approx$ $10^{-5}$). Our modified version of the FNO is capable of solving multi-variable Partial Differential Equations (PDE), and can capture the dependence among the different variables in a single model. FNOs can also predict plasma evolution on real-world experimental data observed by the cameras positioned within the MAST Tokamak, i.e., cameras looking across the central solenoid and the divertor in the Tokamak. We show that FNOs are able to accurately forecast the evolution of plasma and have the potential to be deployed for real-time monitoring. We also illustrate their capability in forecasting the plasma shape, the locations of interactions of the plasma with the central solenoid and the divertor for the full (available) duration of the plasma shot within MAST. The FNO offers a viable alternative for surrogate modelling as it is quick to train and infer, and requires fewer data points, while being able to do zero-shot super-resolution and getting high-fidelity solutions. △ Less

Submitted 18 June, 2024; v1 submitted 10 November, 2023; originally announced November 2023.

Journal ref: Nucl. Fusion 64 056025 (2024)

Efficient training sets for surrogate models of tokamak turbulence with Active Deep Ensembles

Authors: L. Zanisi, A. Ho, T. Madula, J. Barr, J. Citrin, S. Pamela, J. Buchanan, F. Casson, V. Gopakumar, JET contributors

Abstract: Model-based plasma scenario development lies at the heart of the design and operation of future fusion powerplants. Including turbulent transport in integrated models is essential for delivering a successful roadmap towards operation of ITER and the design of DEMO-class devices. Given the highly iterative nature of integrated models, fast machine-learning-based surrogates of turbulent transport ar… ▽ More Model-based plasma scenario development lies at the heart of the design and operation of future fusion powerplants. Including turbulent transport in integrated models is essential for delivering a successful roadmap towards operation of ITER and the design of DEMO-class devices. Given the highly iterative nature of integrated models, fast machine-learning-based surrogates of turbulent transport are fundamental to fulfil the pressing need for faster simulations opening up pulse design, optimization, and flight simulator applications. A significant bottleneck is the generation of suitably large training datasets covering a large volume in parameter space, which can be prohibitively expensive to obtain for higher fidelity codes. In this work, we propose ADEPT (Active Deep Ensembles for Plasma Turbulence), a physics-informed, two-stage Active Learning strategy to ease this challenge. Active Learning queries a given model by means of an acquisition function that identifies regions where additional data would improve the surrogate model. We provide a benchmark study using available data from the literature for the QuaLiKiz quasilinear transport model. We demonstrate quantitatively that the physics-informed nature of the proposed workflow reduces the need to perform simulations in stable regions of the parameter space, resulting in significantly improved data efficiency. We show an up to a factor of 20 reduction in training dataset size needed to achieve the same performance as random sampling. We then validate the surrogates on multichannel integrated modelling of ITG-dominated JET scenarios and demonstrate that they recover the performance of QuaLiKiz to better than 10\%. This matches the performance obtained in previous work, but with two orders of magnitude fewer training data points. △ Less

Submitted 13 October, 2023; originally announced October 2023.

Comments: Submitted to Nuclear Fusion. Comments welcome

Outflows in the Gaseous Discs of Active Galaxies and their impact on Black Hole Scaling Relations

Authors: N. Menci, F. Fiore, F. Shankar, L. Zanisi, C. Feruglio

Abstract: To tackle the still unsolved and fundamental problem of the role of Active Galactic Nuclei (AGN) feedback in shaping galaxies, in this work we implement a new physical treatment of AGN-driven winds into our semi-analytic model of galaxy formation. To each galaxy in our model, we associate solutions for the outflow expansion and the mass outflow rates in different directions, depending on the AGN l… ▽ More To tackle the still unsolved and fundamental problem of the role of Active Galactic Nuclei (AGN) feedback in shaping galaxies, in this work we implement a new physical treatment of AGN-driven winds into our semi-analytic model of galaxy formation. To each galaxy in our model, we associate solutions for the outflow expansion and the mass outflow rates in different directions, depending on the AGN luminosity, on the circular velocity of the host halo, and on gas content of the considered galaxy. To each galaxy we also assign an effective radius derived from energy conservation during merger events, and a stellar velocity dispersion self-consistently computed via Jeans modelling. We derive all the main scaling relations between Black hole (BH) mass and total/bulge stellar mass, velocity dispersion, host halo dark matter mass, and star formation efficiency. We find that our improved AGN feedback mostly controls the dispersion around the relations but plays a subdominant role in shaping slopes and/or normalizations of the scaling relations. Including possible limited-resolution selection biases in the model provides better agreement with the available data. The model does not point to any more fundamental galactic property linked to BH mass, with velocity dispersion playing a similar role with respect to stellar mass, in tension with present data. In line with other independent studies carried out on comprehensive semi-analytic and hydrodynamic galaxy-BH evolution models, our current results signal either an inadequacy of present cosmological models of galaxy formation in fully reproducing the local scaling relations, in terms of both shape and residuals, and/or point to an incompleteness issue affecting the local sample of dynamically-measured BHs. △ Less

Submitted 17 April, 2023; originally announced April 2023.

Comments: 21 pages, accepted for publications in Astronomy & Astrophysics

Fourier Neural Operator for Plasma Modelling

Authors: Vignesh Gopakumar, Stanislas Pamela, Lorenzo Zanisi, Zongyi Li, Anima Anandkumar, MAST Team

Abstract: Predicting plasma evolution within a Tokamak is crucial to building a sustainable fusion reactor. Whether in the simulation space or within the experimental domain, the capability to forecast the spatio-temporal evolution of plasma field variables rapidly and accurately could improve active control methods on current tokamak devices and future fusion reactors. In this work, we demonstrate the util… ▽ More Predicting plasma evolution within a Tokamak is crucial to building a sustainable fusion reactor. Whether in the simulation space or within the experimental domain, the capability to forecast the spatio-temporal evolution of plasma field variables rapidly and accurately could improve active control methods on current tokamak devices and future fusion reactors. In this work, we demonstrate the utility of using Fourier Neural Operator (FNO) to model the plasma evolution in simulations and experiments. Our work shows that the FNO is capable of predicting magnetohydrodynamic models governing the plasma dynamics, 6 orders of magnitude faster than the traditional numerical solver, while maintaining considerable accuracy (NMSE $\sim 10^{-5})$. Our work also benchmarks the performance of the FNO against other standard surrogate models such as Conv-LSTM and U-Net and demonstrate that the FNO takes significantly less time to train, requires less parameters and outperforms other models. We extend the FNO approach to model the plasma evolution observed by the cameras positioned within the MAST spherical tokamak. We illustrate its capability in forecasting the formation of filaments within the plasma as well as the heat deposits. The FNO deployed to model the camera is capable of forecasting the full length of the plasma shot within half the time of the shot duration. △ Less

Submitted 13 February, 2023; originally announced February 2023.

Fourier-RNNs for Modelling Noisy Physics Data

Authors: Vignesh Gopakumar, Stanislas Pamela, Lorenzo Zanisi

Abstract: Classical sequential models employed in time-series prediction rely on learning the mappings from the past to the future instances by way of a hidden state. The Hidden states characterise the historical information and encode the required temporal dependencies. However, most existing sequential models operate within finite-dimensional Euclidean spaces which offer limited functionality when employe… ▽ More Classical sequential models employed in time-series prediction rely on learning the mappings from the past to the future instances by way of a hidden state. The Hidden states characterise the historical information and encode the required temporal dependencies. However, most existing sequential models operate within finite-dimensional Euclidean spaces which offer limited functionality when employed in modelling physics relevant data. Alternatively recent work with neural operator learning within the Fourier space has shown efficient strategies for parameterising Partial Differential Equations (PDE). In this work, we propose a novel sequential model, built to handle Physics relevant data by way of amalgamating the conventional RNN architecture with that of the Fourier Neural Operators (FNO). The Fourier-RNN allows for learning the mappings from the input to the output as well as to the hidden state within the Fourier space associated with the temporal data. While the Fourier-RNN performs identical to the FNO when handling PDE data, it outperforms the FNO and the conventional RNN when deployed in modelling noisy, non-Markovian data. △ Less

Submitted 13 February, 2023; originally announced February 2023.

Testing the key role of the stellar mass-halo mass relation in galaxy merger rates and morphologies via DECODE, a novel Discrete statistical sEmi-empiriCal mODEl

Authors: Hao Fu, Francesco Shankar, Mohammadreza Ayromlou, Max Dickson, Ioanna Koutsouridou, Yetli Rosas-Guevara, Christopher Marsden, Kristina Brocklebank, Mariangela Bernardi, Nikolaos Shiamtanis, Joseph Williams, Lorenzo Zanisi, Viola Allevato, Lumen Boco, Silvia Bonoli, Andrea Cattaneo, Paola Dimauro, Fangzhou Jiang, Andrea Lapi, Nicola Menci, Stefani Petropoulou, Carolin Villforth

Abstract: The relative roles of mergers and star formation in regulating galaxy growth are still a matter of intense debate. We here present our DECODE, a new Discrete statistical sEmi-empiriCal mODEl specifically designed to predict rapidly and efficiently, in a full cosmological context, galaxy assembly and merger histories for any given input stellar mass-halo mass (SMHM) relation. DECODE generates objec… ▽ More The relative roles of mergers and star formation in regulating galaxy growth are still a matter of intense debate. We here present our DECODE, a new Discrete statistical sEmi-empiriCal mODEl specifically designed to predict rapidly and efficiently, in a full cosmological context, galaxy assembly and merger histories for any given input stellar mass-halo mass (SMHM) relation. DECODE generates object-by-object dark matter merger trees (hence discrete) from accurate subhalo mass and infall redshift probability functions (hence statistical) for all subhaloes, including those residing within other subhaloes, with virtually no resolution limits on mass or volume. Merger trees are then converted into galaxy assembly histories via an input, redshift dependent SMHM relation, which is highly sensitive to the significant systematics in the galaxy stellar mass function and on its evolution with cosmic time. DECODE can accurately reproduce the predicted mean galaxy merger rates and assembly histories of hydrodynamic simulations and semi-analytic models, when adopting in input their SMHM relations. In the present work we use DECODE to prove that only SMHM relations implied by stellar mass functions characterized by large abundances of massive galaxies and significant redshift evolution, at least at $M_\star \gtrsim 10^{11} \, M_\odot$, can simultaneously reproduce the local abundances of satellite galaxies, the galaxy (major merger) pairs since $z \sim 3$, and the growth of Brightest Cluster Galaxies. The same models can also reproduce the local fraction of elliptical galaxies, on the assumption that these are strictly formed by major mergers, but not the full bulge-to-disc ratio distributions, which require additional processes. △ Less

Submitted 9 August, 2022; v1 submitted 29 July, 2022; originally announced August 2022.

Comments: MNRAS, accepted, 29 pages, 25 figures

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

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

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

Submitted 13 September, 2021; originally announced September 2021.

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

The size function of massive satellites from the $R_e-R_h$ and $M_{star}-M_h$ relations: constraining the role of environment

Authors: Lorenzo Zanisi, Francesco Shankar, Mariangela Bernardi, Simona Mei, Marc Huertas-Company

Abstract: In previous work we showed that a semi-empirical model in which galaxies in host dark matter haloes are assigned stellar masses via a stellar mass-halo mass (SMHM) relation and sizes ($R_e$) via a linear and tight $R_e-R_h$ relation, can faithfully reproduce the size function of local SDSS central galaxies and the strong size evolution of massive galaxies (MGs, $M_{\rm star}>10^{11.2}M_\odot$). In… ▽ More In previous work we showed that a semi-empirical model in which galaxies in host dark matter haloes are assigned stellar masses via a stellar mass-halo mass (SMHM) relation and sizes ($R_e$) via a linear and tight $R_e-R_h$ relation, can faithfully reproduce the size function of local SDSS central galaxies and the strong size evolution of massive galaxies (MGs, $M_{\rm star}>10^{11.2}M_\odot$). In this third paper of the series, we focus on the population of satellite MGs. We find that without any additional calibration and irrespective of the exact SMHM relation, fraction of quenched galaxies or level of stellar stripping, the same model is able to reproduce the local size function of quiescent satellite MGs in SDSS. In addition, the same model can reproduce the puzzling weak dependence of mean size on host halo mass for both central and satellite galaxies. The model also matches the size function of starforming satellite MGs, after assuming that some of them transform into massive lenticulars in a few Gyr after infalling in the group/cluster environment. However, the vast majority of satellite lenticulars is predicted to form before infall. The $R_e-R_h$ appears to be fundamental to connect galaxies and their host haloes. △ Less

Submitted 27 May, 2021; originally announced May 2021.

Comments: Accepted for publication in MNRAS. 5 pages, 3 figures. Comments and suggestions still welcome

The evolution of compact massive quiescent and starforming galaxies derived from the $R_e-R_h$ and $M_{\rm star}-M_h$ relations

Authors: L. Zanisi, F. Shankar, H. Fu, A. Rodriguez-Puebla, V. Avila-Reese, A. Faisst, E. Daddi, L. Boco, A. Lapi, M. Giavalisco, P. Saracco, F. Buitrago, M. Huertas-Company, A. Puglisi, A. Dekel

Abstract: The mean size ( effective radius $R_e$) of Massive Galaxies (MGs, $M_{\rm star}>10^{11.2}M_\odot$) is observed to increase steadily with cosmic time. It is still unclear whether this trend originates from the size growth of individual galaxies (via, e.g., mergers and/or AGN feedback) or from the inclusion of larger galaxies entering the selection at later epochs (progenitor bias). We here build a… ▽ More The mean size ( effective radius $R_e$) of Massive Galaxies (MGs, $M_{\rm star}>10^{11.2}M_\odot$) is observed to increase steadily with cosmic time. It is still unclear whether this trend originates from the size growth of individual galaxies (via, e.g., mergers and/or AGN feedback) or from the inclusion of larger galaxies entering the selection at later epochs (progenitor bias). We here build a data-driven, flexible theoretical framework to probe the structural evolution of MGs. We assign galaxies to dark matter haloes via stellar mass-halo mass (SMHM) relations with varying high-mass slopes and scatters $σ_{\rm SMHM}$ in stellar mass at fixed halo mass, and assign sizes to galaxies using an empirically-motivated, constant and linear relationship between $R_e$ and the host dark matter halo radius $R_h$. We find that: 1) the fast mean size growth of MGs is well reproduced independently of the shape of the input SMHM relation; 2) the numbers of compact MGs grow steadily until $z\gtrsim2$ and fall off at lower redshifts, suggesting a lesser role of progenitor bias at later epochs; 3) a time-independent scatter $σ_{\rm SMHM}$ is consistent with a scenario in which compact starforming MGs transition into quiescent MGs in a few $10^8$yr with a negligible structural evolution during the compact phase, while a scatter increasing at high redshift implies significant size growth during the starforming phase. A robust measurement of the size function of MGs at high redshift can set strong constraints on the scatter of the SMHM relation and, by extension, on models of galaxy evolution. △ Less

Submitted 20 May, 2021; originally announced May 2021.

Comments: Accepted for publication in MNRAS. Comments still welcome

The MAGPI Survey -- science goals, design, observing strategy, early results and theoretical framework

Authors: C. Foster, J. T. Mendel, C. D. P. Lagos, E. Wisnioski, T. Yuan, F. D'Eugenio, T. M. Barone, K. E. Harborne, S. P. Vaughan, F. Schulze, R. -S. Remus, A. Gupta, F. Collacchioni, D. J. Khim, P. Taylor, R. Bassett, S. M. Croom, R. M. McDermid, A. Poci, A. J. Battisti, J. Bland-Hawthorn, S. Bellstedt, M. Colless, L. J. M. Davies, C. Derkenne , et al. (18 additional authors not shown)

Abstract: We present an overview of the Middle Ages Galaxy Properties with Integral Field Spectroscopy (MAGPI) survey, a Large Program on ESO/VLT. MAGPI is designed to study the physical drivers of galaxy transformation at a lookback time of 3-4 Gyr, during which the dynamical, morphological, and chemical properties of galaxies are predicted to evolve significantly. The survey uses new medium-deep adaptive… ▽ More We present an overview of the Middle Ages Galaxy Properties with Integral Field Spectroscopy (MAGPI) survey, a Large Program on ESO/VLT. MAGPI is designed to study the physical drivers of galaxy transformation at a lookback time of 3-4 Gyr, during which the dynamical, morphological, and chemical properties of galaxies are predicted to evolve significantly. The survey uses new medium-deep adaptive optics aided MUSE observations of fields selected from the GAMA survey, providing a wealth of publicly available ancillary multi-wavelength data. With these data, MAGPI will map the kinematic and chemical properties of stars and ionised gas for a sample of 60 massive (> 7 x 10^10 M_Sun) central galaxies at 0.25 < z <0.35 in a representative range of environments (isolated, groups and clusters). The spatial resolution delivered by MUSE with Ground Layer Adaptive Optics (GLAO, 0.6-0.8 arcsec FWHM) will facilitate a direct comparison with Integral Field Spectroscopy surveys of the nearby Universe, such as SAMI and MaNGA, and at higher redshifts using adaptive optics, e.g. SINS. In addition to the primary (central) galaxy sample, MAGPI will deliver resolved and unresolved spectra for as many as 150 satellite galaxies at 0.25 < z <0.35, as well as hundreds of emission-line sources at z < 6. This paper outlines the science goals, survey design, and observing strategy of MAGPI. We also present a first look at the MAGPI data, and the theoretical framework to which MAGPI data will be compared using the current generation of cosmological hydrodynamical simulations including EAGLE, Magneticum, HORIZON-AGN, and Illustris-TNG. Our results show that cosmological hydrodynamical simulations make discrepant predictions in the spatially resolved properties of galaxies at z ~ 0.3. MAGPI observations will place new constraints and allow for tangible improvements in galaxy formation theory. △ Less

Submitted 14 June, 2021; v1 submitted 27 November, 2020; originally announced November 2020.

Comments: 31 pages, 13 figures, PASA accepted

The relationship between fine galaxy stellar morphology and star formation activity in cosmological simulations: a deep learning view

Authors: Lorenzo Zanisi, Marc Huertas-Company, Francois Lanusse, Connor Bottrell, Annalisa Pillepich, Dylan Nelson, Vicente Rodriguez-Gomez, Francesco Shankar, Lars Hernquist, Avishai Dekel, Berta Margalef-Bentabol, Mark Vogelsberger, Joel Primack

Abstract: Hydrodynamical simulations of galaxy formation and evolution attempt to fully model the physics that shapes galaxies. The agreement between the morphology of simulated and real galaxies, and the way the morphological types are distributed across galaxy scaling relations are important probes of our knowledge of galaxy formation physics. Here we propose an unsupervised deep learning approach to perf… ▽ More Hydrodynamical simulations of galaxy formation and evolution attempt to fully model the physics that shapes galaxies. The agreement between the morphology of simulated and real galaxies, and the way the morphological types are distributed across galaxy scaling relations are important probes of our knowledge of galaxy formation physics. Here we propose an unsupervised deep learning approach to perform a stringent test of the fine morphological structure of galaxies coming from the Illustris and IllustrisTNG (TNG100 and TNG50) simulations against observations from a subsample of the Sloan Digital Sky Survey. Our framework is based on PixelCNN, an autoregressive model for image generation with an explicit likelihood. We adopt a strategy that combines the output of two PixelCNN networks in a metric that isolates the fine morphological details of galaxies from the sky background. We are able to \emph{quantitatively} identify the improvements of IllustrisTNG, particularly in the high-resolution TNG50 run, over the original Illustris. However, we find that the fine details of galaxy structure are still different between observed and simulated galaxies. This difference is driven by small, more spheroidal, and quenched galaxies which are globally less accurate regardless of resolution and which have experienced little improvement between the three simulations explored. We speculate that this disagreement, that is less severe for quenched disky galaxies, may stem from a still too coarse numerical resolution, which struggles to properly capture the inner, dense regions of quenched spheroidal galaxies. △ Less

Submitted 16 December, 2020; v1 submitted 30 June, 2020; originally announced July 2020.

Comments: Accepted for publication in MNRAS. Comments welcome

Probing black hole accretion tracks, scaling relations and radiative efficiencies from stacked X-ray active galactic nuclei

Authors: Francesco Shankar, David H. Weinberg, Christopher Marsden, Philip J. Grylls, Mariangela Bernardi, Guang Yang, Benjamin Moster, Rosamaria Carraro, David M. Alexander, Viola Allevato, Tonima T. Ananna, Angela Bongiorno, Giorgio Calderone, Francesca Civano, Emanuele Daddi, Ivan Delvecchio, Federica Duras, Fabio La Franca, Andrea Lapi, Youjun Lu, Nicola Menci, Mar Mezcua, Federica Ricci, Giulia Rodighiero, Ravi K. Sheth , et al. (3 additional authors not shown)

Abstract: The masses of supermassive black holes at the centres of local galaxies appear to be tightly correlated with the mass and velocity dispersions of their galactic hosts. However, the local Mbh-Mstar relation inferred from dynamically measured inactive black holes is up to an order-of-magnitude higher than some estimates from active black holes, and recent work suggests that this discrepancy arises f… ▽ More The masses of supermassive black holes at the centres of local galaxies appear to be tightly correlated with the mass and velocity dispersions of their galactic hosts. However, the local Mbh-Mstar relation inferred from dynamically measured inactive black holes is up to an order-of-magnitude higher than some estimates from active black holes, and recent work suggests that this discrepancy arises from selection bias on the sample of dynamical black hole mass measurements. In this work we combine X-ray measurements of the mean black hole accretion luminosity as a function of stellar mass and redshift with empirical models of galaxy stellar mass growth, integrating over time to predict the evolving Mbh-Mstar relation. The implied relation is nearly independent of redshift, indicating that stellar and black hole masses grow, on average, at similar rates. Matching the de-biased local Mbh-Mstar relation requires a mean radiative efficiency ~0.15, in line with theoretical expectations for accretion onto spinning black holes. However, matching the "raw" observed relation for inactive black holes requires a mean radiative efficiency around 0.02, far below theoretical expectations. This result provides independent evidence for selection bias in dynamically estimated black hole masses, a conclusion that is robust to uncertainties in bolometric corrections, obscured active black hole fractions, and kinetic accretion efficiency. For our fiducial assumptions, they favour moderate-to-rapid spins of typical supermassive black holes, to achieve a mean radiative efficiency ~0.12-0.20. Our approach has similarities to the classic Soltan analysis, but by using galaxy-based data instead of integrated quantities we are able to focus on regimes where observational uncertainties are minimized. △ Less

Submitted 12 December, 2019; originally announced December 2019.

Comments: 12 pages, 8 Figures. MNRAS, accepted. A discussion around the points raised by arXiv:1909.10821 is included in the Appendix

Galaxy sizes and the galaxy-halo connection -- I: the remarkable tightness of the size distributions

Authors: Lorenzo Zanisi, Francesco Shankar, Andrea Lapi, Nicola Menci, Mariangela Bernardi, Christopher Duckworth, Marc Huertas-Company, Philip Grylls, Paolo Salucci

Abstract: The mass and structural assembly of galaxies is a matter of intense debate. Current theoretical models predict the existence of a linear relationship between galaxy size ($R_e$) and the host dark matter halo virial radius ($R_h$).\\ By making use of semi-empirical models compared to the size distributions of central galaxies from the Sloan Digital Sky Survey, we provide robust constraints on the n… ▽ More The mass and structural assembly of galaxies is a matter of intense debate. Current theoretical models predict the existence of a linear relationship between galaxy size ($R_e$) and the host dark matter halo virial radius ($R_h$).\\ By making use of semi-empirical models compared to the size distributions of central galaxies from the Sloan Digital Sky Survey, we provide robust constraints on the normalization and scatter of the $R_e-R_h$ relation. We explore the parameter space of models in which the $R_e-R_h$ relation is mediated by either the spin parameter or the concentration of the host halo, or a simple constant the nature of which is in principle unknown. We find that the data require extremely tight relations for both early-type and late-type galaxies (ETGs,LTGs), especially for more massive galaxies. These constraints challenge models based solely on angular momentum conservation, which predict significantly wider distributions of galaxy sizes and no trend with stellar mass, if taken at face value. We discuss physically-motivated alterations to the original models that bring the predictions into better agreement with the data. We argue that the measured tight size distributions of SDSS disk galaxies can be reproduced by semi-empirical models in which the $R_e-R_h$ connection is mediated by the \emph{stellar} specific angular momenta $j_{star}.$ We find that current cosmological models of galaxy formation broadly agree with our constraints for LTGs, and justify the strong link between $R_e$ and $j_{star}$ that we propose, however the tightness of the $R_e-R_h$ relation found in such ab-initio theoretical models for ETGs is in tension with our semi-empirical findings. △ Less

Submitted 5 December, 2019; originally announced December 2019.

Comments: 22 pages, 14 figures

Constraining black hole-galaxy scaling relations from the large-scale clustering of Active Galactic Nuclei and implied mean radiative efficiency

Authors: Francesco Shankar, Viola Allevato, Mariangela Bernardi, Christopher Marsden, Andrea Lapi, Nicola Menci, Philip J. Grylls, Mirko Krumpe, Lorenzo Zanisi, Federica Ricci, Fabio La Franca, Ranieri D. Baldi, Jorge Moreno, Ravi K. Sheth

Abstract: A supermassive black hole has been found at the centre of nearly every galaxy observed with sufficient sensitivity. The masses of these black holes are observed to increase with either the total mass or the mean (random) velocity of the stars in their host galaxies. The origin of these correlations remains elusive. Observational systematics and biases severely limit our knowledge of the local demo… ▽ More A supermassive black hole has been found at the centre of nearly every galaxy observed with sufficient sensitivity. The masses of these black holes are observed to increase with either the total mass or the mean (random) velocity of the stars in their host galaxies. The origin of these correlations remains elusive. Observational systematics and biases severely limit our knowledge of the local demography of supermassive black holes thus preventing accurate model comparisons and progress in this field. Here we show that the large-scale spatial distribution of local active galactic nuclei (AGN), believed to be accreting supermassive black holes, can constrain the shape and normalization of the black hole-stellar mass relation thus bypassing resolution-related observational biases. In turn, our results can set more stringent constraints on the so-called "radiative efficiency", a fundamental parameter describing the inner physics of supermassive black holes that is closely linked to their spin, geometry, and ability to release energy. The mean value of the radiative efficiency can be estimated by comparing the average total luminous output of AGN with the relic mass density locked up in quiescent supermassive black holes at galaxy centres today. For currently accepted values of the AGN obscured fractions and bolometric corrections, our newest estimates of the local supermassive black hole mass density favour mean radiative efficiencies of ~10-20\%, suggesting that the vast majority of supermassive black holes are spinning moderately to rapidly. With large-scale AGN surveys coming online, our novel methodology will enable even tighter constraints on the fundamental parameters that regulate the growth of supermassive black holes. △ Less

Submitted 22 October, 2019; originally announced October 2019.

Comments: 31 pages, 7 Figures. Accepted

Predicting fully self-consistent satellite richness, galaxy growth and starformation rates from the STastical sEmi-Empirical modeL STEEL

Authors: Philip J. Grylls, F. Shankar, J. Leja, N. Menci, B. Moster, P. Behroozi, L. Zanisi

Abstract: Observational systematics complicate comparisons with theoretical models limiting understanding of galaxy evolution. In particular, different empirical determinations of the stellar mass function imply distinct mappings between the galaxy and halo masses, leading to diverse galaxy evolutionary tracks. Using our state-of-the-art STatistical sEmi-Empirical modeL, STEEL, we show fully self-consistent… ▽ More Observational systematics complicate comparisons with theoretical models limiting understanding of galaxy evolution. In particular, different empirical determinations of the stellar mass function imply distinct mappings between the galaxy and halo masses, leading to diverse galaxy evolutionary tracks. Using our state-of-the-art STatistical sEmi-Empirical modeL, STEEL, we show fully self-consistent models capable of generating galaxy growth histories that simultaneously and closely agree with the latest data on satellite richness and star-formation rates at multiple redshifts and environments. Central galaxy histories are generated using the central halo mass tracks from state-of-the-art statistical dark matter accretion histories coupled to abundance matching routines. We show that too flat high-mass slopes in the input stellar-mass-halo-mass relations as predicted by previous works, imply non-physical stellar mass growth histories weaker than those implied by satellite accretion alone. Our best-fit models reproduce the satellite distributions at the largest masses and highest redshifts probed, the latest data on star formation rates and its bi-modality in the local Universe, and the correct fraction of ellipticals. Our results are important to predict robust and self-consistent stellar-mass-halo-mass relations and to generate reliable galaxy mock catalogues for the next generations of extra-galactic surveys such as Euclid and LSST. △ Less

Submitted 3 December, 2019; v1 submitted 18 October, 2019; originally announced October 2019.

Comments: 21 Pages

Journal ref: MNRAS, Volume 491, Issue 1, January 2020, Pages 634-654

The Main Sequence at $z \sim 1.3$ contains a sizable fraction of galaxies with compact star formation sizes: a new population of early post-starbursts?

Authors: A. Puglisi, E. Daddi, D. Liu, F. Bournaud, J. D. Silverman, C. Circosta, A. Calabrò, M. Aravena, A. Cibinel, H. Dannerbauer, I. Delvecchio, D. Elbaz, Y. Gao, R. Gobat, S. Jin, E. Le Floc'h, G. E. Magdis, C. Mancini, D. A. Riechers, G. Rodighiero, M. Sargent, F. Valentino, L. Zanisi

Abstract: ALMA measurements for 93 $Herschel$-selected galaxies at $1.1 \leqslant z \leqslant 1.7$ in COSMOS reveal a sizable ($>29$\%) population with compact star formation (SF) sizes, lying on average $> \times 3.6$ below the optical stellar mass ($M_{\star}$)-size relation of disks. This sample widely spans the star-forming Main Sequence (MS), having… ▽ More ALMA measurements for 93 $Herschel$-selected galaxies at $1.1 \leqslant z \leqslant 1.7$ in COSMOS reveal a sizable ($>29$\%) population with compact star formation (SF) sizes, lying on average $> \times 3.6$ below the optical stellar mass ($M_{\star}$)-size relation of disks. This sample widely spans the star-forming Main Sequence (MS), having $10^{8} \leqslant M_{\star} \leqslant 10^{11.5} \ M_{\odot}$ and $20 \leqslant SFR \leqslant 680 \ M_{\odot} \rm yr^{-1}$. The 32 size measurements and 61 upper limits are measured on ALMA images that combine observations of CO(5-4), CO(4-3), CO(2-1) and $λ_{\rm obs} \sim 1.1-1.3 \ \rm mm$ continuum, all tracing the star-forming molecular gas. These compact galaxies have instead normally extended $K_{band}$ sizes, suggesting strong specific $SFR$ gradients. Compact galaxies comprise the $50\pm18 \%$ of MS galaxies at $M_{\star} > 10^{11} M_{\odot}$. This is not expected in standard bi-modal scenarios where MS galaxies are mostly steadily-growing extended disks. We suggest that compact MS objects are early post-starburst galaxies in which the merger-driven boost of SF has subsided. They retain their compact SF size until either further gas accretion restores pre-merger galaxy-wide SF, or until becoming quenched. The fraction of merger-affected SF inside the MS seems thus larger than anticipated and might reach $\sim 50$\% at the highest $M_{\star}$. The presence of large galaxies above the MS demonstrates an overall poor correlation between galaxy SF size and specific $SFR$. △ Less

Submitted 8 May, 2019; originally announced May 2019.

Comments: 7 pages, 5 figures. Accepted for publication in ApJ Letters

A Statistical Semi-Empirical Model: Satellite galaxies in Groups and Clusters

Authors: P. J. Grylls, F. Shankar, L. Zanisi, M. Bernardi

Abstract: We present STEEL a STatistical sEmi-Empirical modeL designed to probe the distribution of satellite galaxies in groups and clusters. Our fast statistical methodology relies on tracing the abundances of central and satellite haloes via their mass functions at all cosmic epochs with virtually no limitation on cosmic volume and mass resolution. From mean halo accretion histories and subhalo mass func… ▽ More We present STEEL a STatistical sEmi-Empirical modeL designed to probe the distribution of satellite galaxies in groups and clusters. Our fast statistical methodology relies on tracing the abundances of central and satellite haloes via their mass functions at all cosmic epochs with virtually no limitation on cosmic volume and mass resolution. From mean halo accretion histories and subhalo mass functions the satellite mass function is progressively built in time via abundance matching techniques constrained by number densities of centrals in the local Universe. By enforcing dynamical merging timescales as predicted by high-resolution N-body simulations, we obtain satellite distributions as a function of stellar mass and halo mass consistent with current data. We show that stellar stripping, star formation, and quenching play all a secondary role in setting the number densities of massive satellites above $M_*\gtrsim 3\times 10^{10}\, M_{\odot}$. We further show that observed star formation rates used in our empirical model over predict low-mass satellites below $M_*\lesssim 3\times 10^{10}\, M_{\odot}$, whereas, star formation rates derived from a continuity equation approach yield the correct abundances similar to previous results for centrals. △ Less

Submitted 30 November, 2018; originally announced December 2018.

Comments: 21 pages, 17 Figures. MNRAS, in press

The Dramatic Size and Kinematic Evolution of Massive Early-Type Galaxies

Authors: A. Lapi, L. Pantoni, L. Zanisi, J. Shi, C. Mancuso, M. Massardi, F. Shankar, A. Bressan, L. Danese

Abstract: [ABRIDGED] We aim to provide a holistic view on the typical size and kinematic evolution of massive early-type galaxies (ETGs), that encompasses their high-$z$ star-forming progenitors, their high-$z$ quiescent counterparts, and their configurations in the local Universe. Our investigation covers the main processes playing a relevant role in the cosmic evolution of ETGs. Specifically, their early… ▽ More [ABRIDGED] We aim to provide a holistic view on the typical size and kinematic evolution of massive early-type galaxies (ETGs), that encompasses their high-$z$ star-forming progenitors, their high-$z$ quiescent counterparts, and their configurations in the local Universe. Our investigation covers the main processes playing a relevant role in the cosmic evolution of ETGs. Specifically, their early fast evolution comprises: biased collapse of the low angular momentum gaseous baryons located in the inner regions of the host dark matter halo; cooling, fragmentation, and infall of the gas down to the radius set by the centrifugal barrier; further rapid compaction via clump/gas migration toward the galaxy center, where strong heavily dust-enshrouded star-formation takes place and most of the stellar mass is accumulated; ejection of substantial gas amount from the inner regions by feedback processes, which causes a dramatic puffing up of the stellar component. In the late slow evolution, passive aging of stellar populations and mass additions by dry merger events occur. We describe these processes relying on prescriptions inspired by basic physical arguments and by numerical simulations, to derive new analytical estimates of the relevant sizes, timescales, and kinematic properties for individual galaxies along their evolution. Then we obtain quantitative results as a function of galaxy mass and redshift, and compare them to recent observational constraints on half-light size $R_e$, on the ratio $v/σ$ between rotation velocity and velocity dispersion (for gas and stars) and on the specific angular momentum $j_\star$ of the stellar component; we find good consistency with the available multi-band data in average values and dispersion, both for local ETGs and for their $z\sim 1-2$ star-forming and quiescent progenitors. △ Less

Submitted 13 March, 2018; originally announced March 2018.

Comments: 28 pages, 10 figures. Accepted by ApJ

Journal ref: Astrophys. J. 857 (2018), 22

Revisiting the bulge-halo conspiracy II: Towards explaining its puzzling dependence on redshift

Authors: Francesco Shankar, Alessandro Sonnenfeld, Philip Grylls, Lorenzo Zanisi, Carlo Nipoti, Kyu-Hyun Chae, Mariangela Bernardi, Carlo Enrico Petrillo, Marc Huertas-Company, Gary A. Mamon, Stewart Buchan

Abstract: We carry out a systematic investigation of the total mass density profile of massive (Mstar~3e11 Msun) early-type galaxies and its dependence on redshift, specifically in the range 0<z<1. We start from a large sample of SDSS early-type galaxies with stellar masses and effective radii measured assuming two different profiles, de Vaucouleurs and Sérsic. We assign dark matter haloes to galaxies via a… ▽ More We carry out a systematic investigation of the total mass density profile of massive (Mstar~3e11 Msun) early-type galaxies and its dependence on redshift, specifically in the range 0<z<1. We start from a large sample of SDSS early-type galaxies with stellar masses and effective radii measured assuming two different profiles, de Vaucouleurs and Sérsic. We assign dark matter haloes to galaxies via abundance matching relations with standard LCDM profiles and concentrations. We then compute the total, mass-weighted density slope at the effective radius gamma', and study its redshift dependence at fixed stellar mass. We find that a necessary condition to induce an increasingly flatter gamma' at higher redshifts, as suggested by current strong lensing data, is to allow the intrinsic stellar profile of massive galaxies to be Sérsic and the input Sérsic index n to vary with redshift approximately as n(z)~(1+z)^(-1). This conclusion holds irrespective of the input Mstar-Mhalo relation, the assumed stellar initial mass function, or even the chosen level of adiabatic contraction in the model. Secondary contributors to the observed redshift evolution of gamma' may come from an increased contribution at higher redshifts of adiabatic contraction and/or bottom-light stellar initial mass functions. The strong lensing selection effects we have simulated seem not to contribute to this effect. A steadily increasing Sérsic index with cosmic time is supported by independent observations, though it is not yet clear whether cosmological hierarchical models (e.g., mergers) are capable of reproducing such a fast and sharp evolution. △ Less

Submitted 21 November, 2017; originally announced November 2017.

Comments: 14 pages, 8 figures. MNRAS, accepted. Main result of the paper in Figure 2