CIGALEMC: GALAXY PARAMETER ESTIMATION USING A MARKOV CHAIN MONTE CARLO APPROACH WITH CIGALE

Serra, Paolo; Amblard, Alexandre; Temi, Pasquale; Im, Stephen; Burgarella, Denis; Giovannoli, Elodie; Buat, Veronique; Noll, Stefan

doi:10.1088/0004-637X/740/1/22; COUNTRY OF INPUT: INTERNATIONAL ATOMIC ENERGY AGENCY (IAEA)

Title: CIGALEMC: GALAXY PARAMETER ESTIMATION USING A MARKOV CHAIN MONTE CARLO APPROACH WITH CIGALE

Full Record
Other Related Research

Abstract

We introduce a fast Markov Chain Monte Carlo (MCMC) exploration of the astrophysical parameter space using a modified version of the publicly available code Code Investigating GALaxy Emission (CIGALE). The original CIGALE builds a grid of theoretical spectral energy distribution (SED) models and fits to photometric fluxes from ultraviolet to infrared to put constraints on parameters related to both formation and evolution of galaxies. Such a grid-based method can lead to a long and challenging parameter extraction since the computation time increases exponentially with the number of parameters considered and results can be dependent on the density of sampling points, which must be chosen in advance for each parameter. MCMC methods, on the other hand, scale approximately linearly with the number of parameters, allowing a faster and more accurate exploration of the parameter space by using a smaller number of efficiently chosen samples. We test our MCMC version of the code CIGALE (called CIGALEMC) with simulated data. After checking the ability of the code to retrieve the input parameters used to build the mock sample, we fit theoretical SEDs to real data from the well-known and -studied Spitzer Infrared Nearby Galaxy Survey sample. We discuss constraints on the parameters andmore »« less

Authors:

Serra, Paolo; Amblard, Alexandre; Temi, Pasquale; Im, Stephen ^[1]; Burgarella, Denis; Giovannoli, Elodie; Buat, Veronique ^[2]; Noll, Stefan ^[3]

Astrophysics Branch, NASA/Ames Research Center, MS 245-6, Moffett Field, CA 94035 (United States)
Observatoire Astronomique de Marseille-Provence, 38 rue Frederic Joliot-Curie, 13388 Marseille Cedex 13 (France)
Institut fuer Astro- und Teilchenphysik, Universitaet Innsbruck, Technikerstr.25/8, 6020 Innsbruck (Austria)

Publication Date:: Mon Oct 10 00:00:00 EDT 2011

OSTI Identifier:: 21587427

Resource Type:: Journal Article

Journal Name:: Astrophysical Journal

Additional Journal Information:: Journal Volume: 740; Journal Issue: 1; Other Information: DOI: 10.1088/0004-637X/740/1/22; Country of input: International Atomic Energy Agency (IAEA); Journal ID: ISSN 0004-637X

Country of Publication:: United States

Language:: English

Subject:: 79 ASTROPHYSICS, COSMOLOGY AND ASTRONOMY; C CODES; DATA ANALYSIS; ENERGY SPECTRA; GALACTIC EVOLUTION; GALAXIES; MONTE CARLO METHOD; SIMULATION; CALCULATION METHODS; COMPUTER CODES; EVOLUTION; SPECTRA

Citation Formats


                    Serra, Paolo, Amblard, Alexandre, Temi, Pasquale, Im, Stephen, Burgarella, Denis, Giovannoli, Elodie, Buat, Veronique, and Noll, Stefan. CIGALEMC: GALAXY PARAMETER ESTIMATION USING A MARKOV CHAIN MONTE CARLO APPROACH WITH CIGALE.  United States: N. p., 2011. 
        Web.  doi:10.1088/0004-637X/740/1/22; COUNTRY OF INPUT: INTERNATIONAL ATOMIC ENERGY AGENCY (IAEA).

Copy to clipboard


                    Serra, Paolo, Amblard, Alexandre, Temi, Pasquale, Im, Stephen, Burgarella, Denis, Giovannoli, Elodie, Buat, Veronique, & Noll, Stefan. CIGALEMC: GALAXY PARAMETER ESTIMATION USING A MARKOV CHAIN MONTE CARLO APPROACH WITH CIGALE.  United States.  https://doi.org/10.1088/0004-637X/740/1/22; COUNTRY OF INPUT: INTERNATIONAL ATOMIC ENERGY AGENCY (IAEA)

Copy to clipboard


                    Serra, Paolo, Amblard, Alexandre, Temi, Pasquale, Im, Stephen, Burgarella, Denis, Giovannoli, Elodie, Buat, Veronique, and Noll, Stefan. 2011.  
        "CIGALEMC: GALAXY PARAMETER ESTIMATION USING A MARKOV CHAIN MONTE CARLO APPROACH WITH CIGALE".  United States.  https://doi.org/10.1088/0004-637X/740/1/22; COUNTRY OF INPUT: INTERNATIONAL ATOMIC ENERGY AGENCY (IAEA).

Copy to clipboard


                    
@article{osti_21587427,

  title        = {CIGALEMC: GALAXY PARAMETER ESTIMATION USING A MARKOV CHAIN MONTE CARLO APPROACH WITH CIGALE},

  author       = {Serra, Paolo and Amblard, Alexandre and Temi, Pasquale and Im, Stephen and Burgarella, Denis and Giovannoli, Elodie and Buat, Veronique and Noll, Stefan},

  abstractNote = {We introduce a fast Markov Chain Monte Carlo (MCMC) exploration of the astrophysical parameter space using a modified version of the publicly available code Code Investigating GALaxy Emission (CIGALE). The original CIGALE builds a grid of theoretical spectral energy distribution (SED) models and fits to photometric fluxes from ultraviolet to infrared to put constraints on parameters related to both formation and evolution of galaxies. Such a grid-based method can lead to a long and challenging parameter extraction since the computation time increases exponentially with the number of parameters considered and results can be dependent on the density of sampling points, which must be chosen in advance for each parameter. MCMC methods, on the other hand, scale approximately linearly with the number of parameters, allowing a faster and more accurate exploration of the parameter space by using a smaller number of efficiently chosen samples. We test our MCMC version of the code CIGALE (called CIGALEMC) with simulated data. After checking the ability of the code to retrieve the input parameters used to build the mock sample, we fit theoretical SEDs to real data from the well-known and -studied Spitzer Infrared Nearby Galaxy Survey sample. We discuss constraints on the parameters and show the advantages of our MCMC sampling method in terms of accuracy of the results and optimization of CPU time.},

  doi          = {10.1088/0004-637X/740/1/22; COUNTRY OF INPUT: INTERNATIONAL ATOMIC ENERGY AGENCY (IAEA)},

  url          = {https://www.osti.gov/biblio/21587427},
  journal      = {Astrophysical Journal},

  issn         = {0004-637X},

  number       = 1,

  volume       = 740,

  place        = {United States},

  year         = {Mon Oct 10 00:00:00 EDT 2011},

  month        = {Mon Oct 10 00:00:00 EDT 2011}

}

Copy to clipboard

Journal Article:

https://doi.org/10.1088/0004-637X/740/1/22; COUNTRY OF INPUT: INTERNATIONAL ATOMIC ENERGY AGENCY (IAEA)

Other availability

Find in Google Scholar

Search WorldCat to find libraries that may hold this journal

Save / Share:

Export Metadata

Save to My Library

Similar records in OSTI.GOV collections:

SPECTRAL ENERGY DISTRIBUTION FITTING WITH MARKOV CHAIN MONTE CARLO: METHODOLOGY AND APPLICATION TO z = 3.1 Ly{alpha}-EMITTING GALAXIES

Journal Article Acquaviva, Viviana; Gawiser, Eric; Guaita, Lucia - Astrophysical Journal

We present GalMC, a Markov Chain Monte Carlo (MCMC) algorithm designed to fit the spectral energy distributions (SEDs) of galaxies to infer physical properties such as age, stellar mass, dust reddening, metallicity, redshift, and star formation rate. We describe the features of the code and the extensive tests conducted to ensure that our procedure leads to unbiased parameter estimation and accurate evaluation of uncertainties. We compare its performance to grid-based algorithms, showing that the efficiency in CPU time is {approx}100 times better for MCMC for a three-dimensional parameter space and increasing with the number of dimensions. We use GalMC tomore »« less
https://doi.org/10.1088/0004-637X/737/2/47
DIGS: deep inference of galaxy spectra with neural posterior estimation

Journal Article Khullar, Gourav; Nord, Brian; Ćiprijanović, Aleksandra; ... - Machine Learning: Science and Technology

Abstract With the advent of billion-galaxy surveys with complex data, the need of the hour is to efficiently model galaxy spectral energy distributions (SEDs) with robust uncertainty quantification. The combination of simulation-based inference (SBI) and amortized neural posterior estimation (NPE) has been successfully used to analyse simulated and real galaxy photometry both precisely and efficiently. In this work, we utilise this combination and build on existing literature to analyse simulated noisy galaxy spectra. Here, we demonstrate a proof-of-concept study of spectra that is (a) an efficient analysis of galaxy SEDs and inference of galaxy parameters with physically interpretable uncertainties; andmore »« less
https://doi.org/10.1088/2632-2153/ac98f4
A critical assessment of stellar mass measurement methods

Journal Article Mobasher, Bahram; Dahlen, Tomas; Ferguson, Henry; ... - Astrophysical Journal

This is the second paper in a series aimed at investigating the main sources of uncertainty in measuring the observable parameters in galaxies from their spectral energy distributions (SEDs). In the first paper we presented a detailed account of the photometric redshift measurements and an error analysis of this process. In this paper we perform a comprehensive study of the main sources of random and systematic error in stellar mass estimates for galaxies, and their relative contributions to the associated error budget. Since there is no prior knowledge of the stellar mass of galaxies (unlike their photometric redshifts), we usemore »« less
https://doi.org/10.1088/0004-637X/808/1/101
Cosmological Evidence Modelling: a new simulation-based approach to constrain cosmology on non-linear scales

Journal Article Lange, Johannes; van den Bosch, Frank; Zentner, Andrew; ... - Monthly Notices of the Royal Astronomical Society

Extracting accurate cosmological information from galaxy–galaxy and galaxy–matter correlation functions on non-linear scales ($${\lesssim } 10 \, h^{-1}{\rm {Mpc}}$$) requires cosmological simulations. Additionally, one has to marginalize over several nuisance parameters of the galaxy–halo connection. However, the computational cost of such simulations prohibits naive implementations of stochastic posterior sampling methods like Markov chain Monte Carlo (MCMC) that would require of order $$\mathcal {O}(10^6)$$ samples in cosmological parameter space. Several groups have proposed surrogate models as a solution: a so-called emulator is trained to reproduce observables for a limited number of realizations in parameter space. Afterwards, this emulator is used asmore »« less
Cited by 19
https://doi.org/10.1093/mnras/stz2664

Full Text Available
A LINK TO THE PAST: USING MARKOV CHAIN MONTE CARLO FITTING TO CONSTRAIN FUNDAMENTAL PARAMETERS OF HIGH-REDSHIFT GALAXIES

Journal Article Pirzkal, N; Rothberg, B; Koekemoer, Anton; ... - Astrophysical Journal

We have developed a new method for fitting spectral energy distributions (SEDs) to identify and constrain the physical properties of high-redshift (4 < z < 8) galaxies. Our approach uses an implementation of Bayesian based Markov Chain Monte Carlo that we have dubbed '{pi}MC{sup 2}'. It allows us to compare observations to arbitrarily complex models and to compute 95% credible intervals that provide robust constraints for the model parameters. The work is presented in two sections. In the first, we test {pi}MC{sup 2} using simulated SEDs to not only confirm the recovery of the known inputs but to assess themore »« less
https://doi.org/10.1088/0004-637X/748/2/122

Similar Records