-
The Chandra Source Catalog Release 2 Series
Authors:
Ian N. Evans,
Janet D. Evans,
J. Rafael Martínez-Galarza,
Joseph B. Miller,
Francis A. Primini,
Mojegan Azadi,
Douglas J. Burke,
Francesca M. Civano,
Raffaele D'Abrusco,
Giuseppina Fabbiano,
Dale E. Graessle,
John D. Grier,
John C. Houck,
Jennifer Lauer,
Michael L. McCollough,
Michael A. Nowak,
David A. Plummer,
Arnold H. Rots,
Aneta Siemiginowska,
Michael S. Tibbetts
Abstract:
The Chandra Source Catalog (CSC) is a virtual X-ray astrophysics facility that enables both detailed individual source studies and statistical studies of large samples of X-ray sources detected in ACIS and HRC-I imaging observations obtained by the Chandra X-ray Observatory. The catalog provides carefully-curated, high-quality, and uniformly calibrated and analyzed tabulated positional, spatial, p…
▽ More
The Chandra Source Catalog (CSC) is a virtual X-ray astrophysics facility that enables both detailed individual source studies and statistical studies of large samples of X-ray sources detected in ACIS and HRC-I imaging observations obtained by the Chandra X-ray Observatory. The catalog provides carefully-curated, high-quality, and uniformly calibrated and analyzed tabulated positional, spatial, photometric, spectral, and temporal source properties, as well as science-ready X-ray data products. The latter includes multiple types of source- and field-based FITS format products that can be used as a basis for further research, significantly simplifying followup analysis of scientifically meaningful source samples. We discuss in detail the algorithms used for the CSC Release 2 Series, including CSC 2.0, which includes 317,167 unique X-ray sources on the sky identified in observations released publicly through the end of 2014, and CSC 2.1, which adds Chandra data released through the end of 2021 and expands the catalog to 407,806 sources. Besides adding more recent observations, the CSC Release 2 Series includes multiple algorithmic enhancements that provide significant improvements over earlier releases. The compact source sensitivity limit for most observations is ~5 photons over most of the field of view, which is ~2x fainter than Release 1, achieved by co-adding observations and using an optimized source detection approach. A Bayesian X-ray aperture photometry code produces robust fluxes even in crowded fields and for low count sources. The current release, CSC 2.1, is tied to the Gaia-CRF3 astrometric reference frame for the best sky positions for catalog sources.
△ Less
Submitted 15 July, 2024;
originally announced July 2024.
-
Unsupervised Machine Learning for the Classification of Astrophysical X-ray Sources
Authors:
Víctor Samuel Pérez-Díaz,
Juan Rafael Martínez-Galarza,
Alexander Caicedo,
Raffaele D'Abrusco
Abstract:
The automatic classification of X-ray detections is a necessary step in extracting astrophysical information from compiled catalogs of astrophysical sources. Classification is useful for the study of individual objects, statistics for population studies, as well as for anomaly detection, i.e., the identification of new unexplored phenomena, including transients and spectrally extreme sources. Desp…
▽ More
The automatic classification of X-ray detections is a necessary step in extracting astrophysical information from compiled catalogs of astrophysical sources. Classification is useful for the study of individual objects, statistics for population studies, as well as for anomaly detection, i.e., the identification of new unexplored phenomena, including transients and spectrally extreme sources. Despite the importance of this task, classification remains challenging in X-ray astronomy due to the lack of optical counterparts and representative training sets. We develop an alternative methodology that employs an unsupervised machine learning approach to provide probabilistic classes to Chandra Source Catalog sources with a limited number of labeled sources, and without ancillary information from optical and infrared catalogs. We provide a catalog of probabilistic classes for 8,756 sources, comprising a total of 14,507 detections, and demonstrate the success of the method at identifying emission from young stellar objects, as well as distinguishing between small-scale and large-scale compact accretors with a significant level of confidence. We investigate the consistency between the distribution of features among classified objects and well-established astrophysical hypotheses such as the unified AGN model. This provides interpretability to the probabilistic classifier. Code and tables are available publicly through GitHub. We provide a web playground for readers to explore our final classification at https://umlcaxs-playground.streamlit.app.
△ Less
Submitted 22 January, 2024;
originally announced January 2024.
-
EWOCS-I: The catalog of X-ray sources in Westerlund 1 from the Extended Westerlund 1 and 2 Open Clusters Survey
Authors:
M. G. Guarcello,
E. Flaccomio,
J. F. Albacete-Colombo,
V. Almendros-Abad,
K. Anastasopoulou,
M. Andersen,
C. Argiroffi,
A. Bayo,
E. S. Bartlett,
N. Bastian,
M. De Becker,
W. Best,
R. Bonito,
A. Borghese,
D. Calzetti,
R. Castellanos,
C. Cecchi-Pestellini,
S. Clark,
C. J. Clarke,
F. Coti Zelati,
F. Damiani,
J. J. Drake,
M. Gennaro,
A. Ginsburg,
E. K. Grebel
, et al. (26 additional authors not shown)
Abstract:
Context. With a mass exceeding several 10^4 solar masses and a rich and dense population of massive stars, supermassive young star clusters represent the most massive star-forming environment that is dominated by the feedback from massive stars and gravitational interactions among stars. Aims. In this paper we present the "Extended Westerlund 1 and 2 Open Clusters Survey" (EWOCS) project, which ai…
▽ More
Context. With a mass exceeding several 10^4 solar masses and a rich and dense population of massive stars, supermassive young star clusters represent the most massive star-forming environment that is dominated by the feedback from massive stars and gravitational interactions among stars. Aims. In this paper we present the "Extended Westerlund 1 and 2 Open Clusters Survey" (EWOCS) project, which aims to investigate the influence of the starburst environment on the formation of stars and planets, and on the evolution of both low and high mass stars. The primary targets of this project are Westerlund 1 and 2, the closest supermassive star clusters to the Sun. Methods. The project is based primarily on recent observations conducted with the Chandra and JWST observatories. Specifically, the Chandra survey of Westerlund 1 consists of 36 new ACIS-I observations, nearly co-pointed, for a total exposure time of 1 Msec. Additionally, we included 8 archival Chandra/ACIS-S observations. This paper presents the resulting catalog of X-ray sources within and around Westerlund 1. Sources were detected by combining various existing methods, and photon extraction and source validation were carried out using the ACIS-Extract software. Results. The EWOCS X-ray catalog comprises 5963 validated sources out of the 9420 initially provided to ACIS-Extract, reaching a photon flux threshold of approximately 2x10^-8 photons/cm^2/s. The X-ray sources exhibit a highly concentrated spatial distribution, with 1075 sources located within the central 1 arcminute. We have successfully detected X-ray emissions from 126 out of the 166 known massive stars of the cluster, and we have collected over 71000 photons from the magnetar CXO J164710.20-455217
△ Less
Submitted 15 December, 2023; v1 submitted 14 December, 2023;
originally announced December 2023.
-
The Future of Astronomical Data Infrastructure: Meeting Report
Authors:
Michael R. Blanton,
Janet D. Evans,
Dara Norman,
William O'Mullane,
Adrian Price-Whelan,
Luca Rizzi,
Alberto Accomazzi,
Megan Ansdell,
Stephen Bailey,
Paul Barrett,
Steven Berukoff,
Adam Bolton,
Julian Borrill,
Kelle Cruz,
Julianne Dalcanton,
Vandana Desai,
Gregory P. Dubois-Felsmann,
Frossie Economou,
Henry Ferguson,
Bryan Field,
Dan Foreman-Mackey,
Jaime Forero-Romero,
Niall Gaffney,
Kim Gillies,
Matthew J. Graham
, et al. (47 additional authors not shown)
Abstract:
The astronomical community is grappling with the increasing volume and complexity of data produced by modern telescopes, due to difficulties in reducing, accessing, analyzing, and combining archives of data. To address this challenge, we propose the establishment of a coordinating body, an "entity," with the specific mission of enhancing the interoperability, archiving, distribution, and productio…
▽ More
The astronomical community is grappling with the increasing volume and complexity of data produced by modern telescopes, due to difficulties in reducing, accessing, analyzing, and combining archives of data. To address this challenge, we propose the establishment of a coordinating body, an "entity," with the specific mission of enhancing the interoperability, archiving, distribution, and production of both astronomical data and software. This report is the culmination of a workshop held in February 2023 on the Future of Astronomical Data Infrastructure. Attended by 70 scientists and software professionals from ground-based and space-based missions and archives spanning the entire spectrum of astronomical research, the group deliberated on the prevailing state of software and data infrastructure in astronomy, identified pressing issues, and explored potential solutions. In this report, we describe the ecosystem of astronomical data, its existing flaws, and the many gaps, duplication, inconsistencies, barriers to access, drags on productivity, missed opportunities, and risks to the long-term integrity of essential data sets. We also highlight the successes and failures in a set of deep dives into several different illustrative components of the ecosystem, included as an appendix.
△ Less
Submitted 7 November, 2023;
originally announced November 2023.
-
Linking Anomalous Behaviour with Stellar Properties: An Unsupervised Exploration of TESS Light Curves
Authors:
Dennis A. Crake,
Juan Rafael Martínez-Galarza
Abstract:
With the upcoming plethora of astronomical time-domain datasets and surveys, anomaly detection as a way to discover new types of variable stars and transients has inspired a new wave of research. Yet, the fundamental definition of what constitutes an anomaly and how this depends on the overall properties of the population of light curves studied remains a discussed issue. Building on a previous st…
▽ More
With the upcoming plethora of astronomical time-domain datasets and surveys, anomaly detection as a way to discover new types of variable stars and transients has inspired a new wave of research. Yet, the fundamental definition of what constitutes an anomaly and how this depends on the overall properties of the population of light curves studied remains a discussed issue. Building on a previous study focused on Kepler light curves, we present an analysis that uses the Unsupervised Random Forest to search for anomalies in TESS light curves. We provide a catalogue of anomalous light curves, classify them according to their variability characteristics and associate their anomalous nature to any particular evolutionary stage or astrophysical configuration. For anomalies belonging to known classes (e.g. eclipsing binaries), we have investigated which physical parameters drive the anomaly score. We find a combination of unclassified anomalies and objects of a known class with outlying physical configurations, such as rapid pulsators, deep eclipsing binaries of long periods, and irregular light curves due to obscuration in YSOs. Remarkably, we find that the set of anomalous types differ between the Kepler and TESS datasets, indicating that the overall properties of the parent population are an important driver of anomalous behaviour.
△ Less
Submitted 24 January, 2023;
originally announced January 2023.
-
Wavelength Calibration and Resolving Power of the JWST MIRI Medium Resolution Spectrometer
Authors:
A. Labiano,
I. Argyriou,
J. Alvarez-Marquez,
A. Glasse,
A. Glauser,
P. Patapis,
D. Law,
B. R. Brandl,
K. Justtanont,
F. Lahuis,
J. R. Martinez-Galarza,
M. Mueller,
A. Noriega-Crespo,
P. Royer,
B. Shaughnessy,
B. Vandenbussche
Abstract:
The Mid-Infrared Instrument (MIRI) on-board JWST will provide imaging, coronagraphy, low-resolution spectroscopy and medium-resolution spectroscopy at unprecedented sensitivity levels in the mid-infrared wavelength range. The Medium-Resolution Spectrometer (MRS) of MIRI is an integral field spectrograph that provides diffraction-limited spectroscopy between 4.9 and 28.3 um, within a FOV varying fr…
▽ More
The Mid-Infrared Instrument (MIRI) on-board JWST will provide imaging, coronagraphy, low-resolution spectroscopy and medium-resolution spectroscopy at unprecedented sensitivity levels in the mid-infrared wavelength range. The Medium-Resolution Spectrometer (MRS) of MIRI is an integral field spectrograph that provides diffraction-limited spectroscopy between 4.9 and 28.3 um, within a FOV varying from 13 to 56" square. From ground testing, we calculate the physical parameters essential to general observers and calibrating the wavelength solution and resolving power of the MRS is critical for maximising the scientific performance of the instrument. We have used ground-based observations of discrete spectral features in combination with Fabry-Perot etalon spectra to characterize the wavelength solution and spectral resolving power of the MRS. We present the methodology used to derive the MRS spectral characterisation, which includes the precise wavelength coverage of each MRS sub-band, computation of the resolving power as a function of wavelength, and measuring slice-dependent spectral distortions. The resolving power varies from R3500 in channel 1 to R1500 in channel 4. Based on the ground test data, the wavelength calibration accuracy is estimated to be below one tenth of a pixel, with small systematic shifts due to the target position within a slice for unresolved sources, that have a maximum amplitude of about 0.25 spectral resolution elements. Based on ground test data, the MRS complies with the spectral requirements for both the R and wavelength accuracy for which it was designed. We also present the commissioning strategies and targets that will be followed to update the spectral characterisation of the MRS.
△ Less
Submitted 8 September, 2021;
originally announced September 2021.
-
Dust-Enshrouded AGN can Dominate Host-Galaxy-Scale Cold-Dust Emission
Authors:
Jed McKinney,
Christopher C. Hayward,
Lee J. Rosenthal,
Juan Rafael Martinez-Galarza,
Alexandra Pope,
Anna Sajina,
Howard A. Smith
Abstract:
It is widely assumed that long-wavelength infrared (IR) emission from cold dust (T~20-40K) is a reliable tracer of star formation even in the presence of a bright active galactic nucleus (AGN). Based on radiative transfer (RT) models of clumpy AGN tori, hot dust emission from the torus contributes negligibly to the galaxy spectral energy distribution (SED) at $λ\ga100$ \micron. However, these mode…
▽ More
It is widely assumed that long-wavelength infrared (IR) emission from cold dust (T~20-40K) is a reliable tracer of star formation even in the presence of a bright active galactic nucleus (AGN). Based on radiative transfer (RT) models of clumpy AGN tori, hot dust emission from the torus contributes negligibly to the galaxy spectral energy distribution (SED) at $λ\ga100$ \micron. However, these models do not include AGN heating of host-galaxy-scale diffuse dust, which may have far-IR (FIR) colors comparable to cold diffuse dust heated by stars. To quantify the contribution of AGN heating to host-galaxy-scale cold dust emission at $λ\ga100$ \micron, we perform dust RT calculations on a simulated galaxy merger both including and excluding the bright AGN that it hosts. By differencing the SEDs yielded by RT calculations with and without AGN that are otherwise identical, we quantify the FIR cold dust emission arising solely from re-processed AGN photons. In extreme cases, AGN-heated host-galaxy-scale dust can increase galaxy-integrated FIR flux densities by factors of 2-4; star formation rates calculated from the FIR luminosity assuming no AGN contribution can overestimate the true value by comparable factors. Because the FIR colors of such systems are similar to those of purely star-forming galaxies and redder than torus models, broadband SED decomposition may be insufficient for disentangling the contributions of stars and heavily dust-enshrouded AGN in the most IR-luminous galaxies. We demonstrate how kpc-scale resolved observations can be used to identify deeply dust-enshrouded AGN with cool FIR colors when spectroscopic and/or X-ray detection methods are unavailable.
△ Less
Submitted 26 July, 2021; v1 submitted 23 March, 2021;
originally announced March 2021.
-
Results of the Photometric LSST Astronomical Time-series Classification Challenge (PLAsTiCC)
Authors:
R. Hložek,
K. A. Ponder,
A. I. Malz,
M. Dai,
G. Narayan,
E. E. O. Ishida,
T. Allam Jr,
A. Bahmanyar,
R. Biswas,
L. Galbany,
S. W. Jha,
D. O. Jones,
R. Kessler,
M. Lochner,
A. A. Mahabal,
K. S. Mandel,
J. R. Martínez-Galarza,
J. D. McEwen,
D. Muthukrishna,
H. V. Peiris,
C. M. Peters,
C. N. Setzer
Abstract:
Next-generation surveys like the Legacy Survey of Space and Time (LSST) on the Vera C. Rubin Observatory will generate orders of magnitude more discoveries of transients and variable stars than previous surveys. To prepare for this data deluge, we developed the Photometric LSST Astronomical Time-series Classification Challenge (PLAsTiCC), a competition which aimed to catalyze the development of ro…
▽ More
Next-generation surveys like the Legacy Survey of Space and Time (LSST) on the Vera C. Rubin Observatory will generate orders of magnitude more discoveries of transients and variable stars than previous surveys. To prepare for this data deluge, we developed the Photometric LSST Astronomical Time-series Classification Challenge (PLAsTiCC), a competition which aimed to catalyze the development of robust classifiers under LSST-like conditions of a non-representative training set for a large photometric test set of imbalanced classes. Over 1,000 teams participated in PLAsTiCC, which was hosted in the Kaggle data science competition platform between Sep 28, 2018 and Dec 17, 2018, ultimately identifying three winners in February 2019. Participants produced classifiers employing a diverse set of machine learning techniques including hybrid combinations and ensemble averages of a range of approaches, among them boosted decision trees, neural networks, and multi-layer perceptrons. The strong performance of the top three classifiers on Type Ia supernovae and kilonovae represent a major improvement over the current state-of-the-art within astronomy. This paper summarizes the most promising methods and evaluates their results in detail, highlighting future directions both for classifier development and simulation needs for a next generation PLAsTiCC data set.
△ Less
Submitted 22 December, 2020;
originally announced December 2020.
-
Disentangling the AGN and Star-Formation Contributions to the Radio-X-ray Emission of Radio-Loud Quasars at 1<z<2
Authors:
Mojegan Azadi,
Belinda Wilkes,
Joanna Kuraszkiewicz,
Jonathan McDowell,
Ralf Siebenmorgen,
Matthew Ashby,
Mark Birkinshaw,
Diana Worrall,
Natasha Abrams,
Peter Barthel,
Giovanni Fazio,
Martin Haas,
Sóley Hyman,
Rafael Martínez-Galarza,
Eileen Meyer
Abstract:
We constrain the emission mechanisms responsible for the prodigious electromagnetic output generated by active galactic nuclei (AGN) and their host galaxies with a novel state-of-the-art AGN radio- to-X-ray spectral energy distribution model fitting code (ARXSED). ARXSED combines multiple components to fit the spectral energy distributions (SEDs) of AGN and their host galaxies. Emission components…
▽ More
We constrain the emission mechanisms responsible for the prodigious electromagnetic output generated by active galactic nuclei (AGN) and their host galaxies with a novel state-of-the-art AGN radio- to-X-ray spectral energy distribution model fitting code (ARXSED). ARXSED combines multiple components to fit the spectral energy distributions (SEDs) of AGN and their host galaxies. Emission components include radio structures such as lobes and jets, infrared emission from the AGN torus, visible-to-X-ray emission from the accretion disk, and radio-to-ultraviolet emission from the host galaxy. Applying ARXSED to the radio SEDs of 20 3CRR quasars at 1 < z < 2 verifies the need for more than a simple power law when compact radio structures are present. The non-thermal emission contributes 91%-57% of the observed-frame 1.25mm to 850μm flux, and this component must be accounted for when using these wavelengths to estimate star-formation properties. We predict the presence of strong radio-linked X-ray emission in more than half the sample sources. ARXSED estimates median (and the associated first and third quartile ranges) BH mass of $2.9_{1.7}^{6.0} \times 10^9~\rm M_{\odot}$, logarithm of Eddington ratio of $ -1.0_{-1.2}^{-0.6} $, and spin of $ 0.98_{0.94}^{0.99} $ for our sample. The inferred AGN torus and accretion disk parameters agree with those estimated from spectroscopic analyses of similar samples in the literature. We present the median intrinsic SED of the luminous radio-loud quasars at 1 < z < 2 ; this SED represents a significant improvement in the way each component is modeled.
△ Less
Submitted 24 February, 2023; v1 submitted 5 November, 2020;
originally announced November 2020.
-
A method for finding anomalous astronomical light curves and their analogs
Authors:
Juan Rafael Martínez-Galarza,
Federica Bianco,
Dennis Crake,
Kushal Tirumala,
Ashish A. Mahabal,
Matthew J. Graham,
Daniel Giles
Abstract:
Our understanding of the Universe has profited from deliberate, targeted studies of known phenomena, as well as from serendipitous, unexpected discoveries, such as the discovery of a complex variability pattern in the direction of KIC 8462852 (Boyajian's star). Upcoming surveys, such as the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST), will explore the parameter space of astrop…
▽ More
Our understanding of the Universe has profited from deliberate, targeted studies of known phenomena, as well as from serendipitous, unexpected discoveries, such as the discovery of a complex variability pattern in the direction of KIC 8462852 (Boyajian's star). Upcoming surveys, such as the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST), will explore the parameter space of astrophysical transients at all time scales, and offer the opportunity to discover even more extreme examples of unexpected phenomena. We investigate strategies to identify novel objects and to contextualize them within large time-series data sets in order to facilitate the discovery of new classes of objects, as well as the physical interpretation of their anomalous nature. We develop a method that combines tree-based and manifold-learning algorithms for anomaly detection in order to perform two tasks: 1) identify and rank anomalous objects in a time-domain dataset; and 2) group those anomalies according to their similarity in order to identify analogs. We achieve the latter by combining an anomaly score from a tree-based method with a dimensionality manifold-learning reduction strategy. Clustering in the reduced space allows for the successful identification of anomalies and analogs. We also assess the impact of pre-processing and feature engineering schemes and investigate the astrophysical nature of the objects that our models identify as anomalous by augmenting the Kepler data with Gaia color and luminosity information. We find that multiple models, used in combination, are a promising strategy to identify novel light curves and light curve families.
△ Less
Submitted 15 September, 2021; v1 submitted 14 September, 2020;
originally announced September 2020.
-
The AGN contribution to the UV-FIR luminosities of interacting galaxies and its role in identifying the Main Sequence
Authors:
Andrés F. Ramos P.,
M. L. N. Ashby,
Howard A. Smith,
Juan R. Martínez-Galarza,
Aliza G. Beverage,
Jeremy Dietrich,
Mario-A. Higuera-G.,
Aaron S. Weiner
Abstract:
Emission from active galactic nuclei (AGNs) is known to play an important role in the evolution of many galaxies including luminous and ultraluminous systems (U/LIRGs), as well as merging systems. However, the extent, duration, and exact effects of its influence are still imperfectly understood. To assess the impact of AGNs on interacting systems, we present a Spectral Energy Distribution (SED) an…
▽ More
Emission from active galactic nuclei (AGNs) is known to play an important role in the evolution of many galaxies including luminous and ultraluminous systems (U/LIRGs), as well as merging systems. However, the extent, duration, and exact effects of its influence are still imperfectly understood. To assess the impact of AGNs on interacting systems, we present a Spectral Energy Distribution (SED) analysis of a sample of 189 nearby galaxies. We gather and systematically re-reduce archival broad-band imaging mosaics from the ultraviolet to the far-infrared using data from GALEX, SDSS, 2MASS, IRAS, WISE, Spitzer and Herschel. We use spectroscopy from Spitzer/IRS to obtain fluxes from fine-structure lines that trace star formation and AGN activity. Utilizing the SED modelling and fitting tool CIGALE, we derive the physical conditions of the ISM, both in star-forming regions and in nuclear regions dominated by the AGN in these galaxies. We investigate how the star formation rates (SFRs) and the fractional AGN contributions ($f_{\rm{AGN}}$) depend on stellar mass, galaxy type, and merger stage. We find that luminous galaxies more massive than about $10^{10} \rm{M}_{*}$ are likely to deviate significantly from the conventional galaxy main-sequence relation. Interestingly, infrared AGN luminosity and stellar mass in this set of objects are much tighter than SFR and stellar mass. We find that buried AGNs may occupy a locus between bright starbursts and pure AGNs in the $f_{\rm{AGN}}$-[Ne V]/[Ne II] plane. We identify a modest correlation between $f_{\rm{AGN}}$ and mergers in their later stages.
△ Less
Submitted 11 September, 2020;
originally announced September 2020.
-
Algorithms and Statistical Models for Scientific Discovery in the Petabyte Era
Authors:
Brian Nord,
Andrew J. Connolly,
Jamie Kinney,
Jeremy Kubica,
Gautaum Narayan,
Joshua E. G. Peek,
Chad Schafer,
Erik J. Tollerud,
Camille Avestruz,
G. Jogesh Babu,
Simon Birrer,
Douglas Burke,
João Caldeira,
Douglas A. Caldwell,
Joleen K. Carlberg,
Yen-Chi Chen,
Chuanfei Dong,
Eric D. Feigelson,
V. Zach Golkhou,
Vinay Kashyap,
T. S. Li,
Thomas Loredo,
Luisa Lucie-Smith,
Kaisey S. Mandel,
J. R. Martínez-Galarza
, et al. (13 additional authors not shown)
Abstract:
The field of astronomy has arrived at a turning point in terms of size and complexity of both datasets and scientific collaboration. Commensurately, algorithms and statistical models have begun to adapt --- e.g., via the onset of artificial intelligence --- which itself presents new challenges and opportunities for growth. This white paper aims to offer guidance and ideas for how we can evolve our…
▽ More
The field of astronomy has arrived at a turning point in terms of size and complexity of both datasets and scientific collaboration. Commensurately, algorithms and statistical models have begun to adapt --- e.g., via the onset of artificial intelligence --- which itself presents new challenges and opportunities for growth. This white paper aims to offer guidance and ideas for how we can evolve our technical and collaborative frameworks to promote efficient algorithmic development and take advantage of opportunities for scientific discovery in the petabyte era. We discuss challenges for discovery in large and complex data sets; challenges and requirements for the next stage of development of statistical methodologies and algorithmic tool sets; how we might change our paradigms of collaboration and education; and the ethical implications of scientists' contributions to widely applicable algorithms and computational modeling. We start with six distinct recommendations that are supported by the commentary following them. This white paper is related to a larger corpus of effort that has taken place within and around the Petabytes to Science Workshops (https://petabytestoscience.github.io/).
△ Less
Submitted 4 November, 2019;
originally announced November 2019.
-
The Photometric LSST Astronomical Time-series Classification Challenge (PLAsTiCC): Data set
Authors:
The PLAsTiCC team,
Tarek Allam Jr.,
Anita Bahmanyar,
Rahul Biswas,
Mi Dai,
Lluís Galbany,
Renée Hložek,
Emille E. O. Ishida,
Saurabh W. Jha,
David O. Jones,
Richard Kessler,
Michelle Lochner,
Ashish A. Mahabal,
Alex I. Malz,
Kaisey S. Mandel,
Juan Rafael Martínez-Galarza,
Jason D. McEwen,
Daniel Muthukrishna,
Gautham Narayan,
Hiranya Peiris,
Christina M. Peters,
Kara Ponder,
Christian N. Setzer,
The LSST Dark Energy Science Collaboration,
The LSST Transients
, et al. (1 additional authors not shown)
Abstract:
The Photometric LSST Astronomical Time Series Classification Challenge (PLAsTiCC) is an open data challenge to classify simulated astronomical time-series data in preparation for observations from the Large Synoptic Survey Telescope (LSST), which will achieve first light in 2019 and commence its 10-year main survey in 2022. LSST will revolutionize our understanding of the changing sky, discovering…
▽ More
The Photometric LSST Astronomical Time Series Classification Challenge (PLAsTiCC) is an open data challenge to classify simulated astronomical time-series data in preparation for observations from the Large Synoptic Survey Telescope (LSST), which will achieve first light in 2019 and commence its 10-year main survey in 2022. LSST will revolutionize our understanding of the changing sky, discovering and measuring millions of time-varying objects.
In this challenge, we pose the question: how well can we classify objects in the sky that vary in brightness from simulated LSST time-series data, with all its challenges of non-representativity? In this note we explain the need for a data challenge to help classify such astronomical sources and describe the PLAsTiCC data set and Kaggle data challenge, noting that while the references are provided for context, they are not needed to participate in the challenge.
△ Less
Submitted 28 September, 2018;
originally announced October 2018.
-
The Photometric LSST Astronomical Time-series Classification Challenge (PLAsTiCC): Selection of a performance metric for classification probabilities balancing diverse science goals
Authors:
A. I. Malz,
R. Hložek,
T. Allam Jr,
A. Bahmanyar,
R. Biswas,
M. Dai,
L. Galbany,
E. E. O. Ishida,
S. W. Jha,
D. O. Jones,
R. Kessler,
M. Lochner,
A. A. Mahabal,
K. S. Mandel,
J. R. Martínez-Galarza,
J. D. McEwen,
D. Muthukrishna,
G. Narayan,
H. Peiris,
C. M. Peters,
K. A. Ponder,
C. N. Setzer,
The LSST Dark Energy Science Collaboration,
The LSST Transients,
Variable Stars Science Collaboration
Abstract:
Classification of transient and variable light curves is an essential step in using astronomical observations to develop an understanding of their underlying physical processes. However, upcoming deep photometric surveys, including the Large Synoptic Survey Telescope (LSST), will produce a deluge of low signal-to-noise data for which traditional labeling procedures are inappropriate. Probabilistic…
▽ More
Classification of transient and variable light curves is an essential step in using astronomical observations to develop an understanding of their underlying physical processes. However, upcoming deep photometric surveys, including the Large Synoptic Survey Telescope (LSST), will produce a deluge of low signal-to-noise data for which traditional labeling procedures are inappropriate. Probabilistic classification is more appropriate for the data but are incompatible with the traditional metrics used on deterministic classifications. Furthermore, large survey collaborations intend to use these classification probabilities for diverse science objectives, indicating a need for a metric that balances a variety of goals. We describe the process used to develop an optimal performance metric for an open classification challenge that seeks probabilistic classifications and must serve many scientific interests. The Photometric LSST Astronomical Time-series Classification Challenge (PLAsTiCC) is an open competition aiming to identify promising techniques for obtaining classification probabilities of transient and variable objects by engaging a broader community both within and outside astronomy. Using mock classification probability submissions emulating archetypes of those anticipated of PLAsTiCC, we compare the sensitivity of metrics of classification probabilities under various weighting schemes, finding that they yield qualitatively consistent results. We choose as a metric for PLAsTiCC a weighted modification of the cross-entropy because it can be meaningfully interpreted. Finally, we propose extensions of our methodology to ever more complex challenge goals and suggest some guiding principles for approaching the choice of a metric of probabilistic classifications.
△ Less
Submitted 31 July, 2021; v1 submitted 28 September, 2018;
originally announced September 2018.
-
Unraveling the Spectral Energy Distributions of Clustered YSOs
Authors:
Juan R. Martínez-Galarza,
Pavlos Protopapas,
Howard A. Smith,
Esteban F. E. Morales
Abstract:
Stars form in clustered environments, but how they form when the available resources are shared is still not well understood. A related question is whether the IMF is in fact universal across galactic environments, a galactic initial mass function (IGIMF), or whether it is an average of local IMFs. One of the long-standing problems in resolving this question and in the study of young clusters is o…
▽ More
Stars form in clustered environments, but how they form when the available resources are shared is still not well understood. A related question is whether the IMF is in fact universal across galactic environments, a galactic initial mass function (IGIMF), or whether it is an average of local IMFs. One of the long-standing problems in resolving this question and in the study of young clusters is observational: the emission from multiple sources is frequently seen as blended because at different wavelengths or with different telescopes the beam sizes are different. The confusion hinders our ability to fully characterize clustered star formation. Here we present a new method that uses a genetic algorithm and Bayesian inference to fit the blended SEDs and images of individual YSOs in confused clusters. We apply this method to the infrared photometry of a sample comprising 70 Spitzer-selected, low-mass ($M_{\rm{cl}}<100~\rm{M}_{\odot}$) young clusters in the galactic plane, and use the derived physical parameters to investigate the distributions of masses and evolutionary stages of clustered YSOs, and the implications of those distributions for studies of the IMF and the different models of star formation. We find that for low-mass clusters composed of class I and class II YSOs, there exists a non-trivial relationship between the total stellar mass of the cluster ($M_{\rm{cl}}$) and the mass of its most massive member ($m_{\rm{max}}$). The properties of the derived correlation are most compatible with the random sampling of a Kroupa IMF, with a fundamental high-mass limit of $150~\rm{M}_{\odot}$. Our results are also compatible with SPH models that predict a dynamical termination of the accretion in protostars, with massive stars undergoing this stopping at later times in their evolution.
△ Less
Submitted 28 March, 2018;
originally announced March 2018.
-
The AGN Luminosity Fraction in Merging Galaxies
Authors:
Jeremy Dietrich,
Aaron S. Weiner,
Matthew L. N. Ashby,
Christopher C. Hayward,
Juan Rafael Martínez-Galarza,
Andrés F. Ramos Padilla,
Lee Rosenthal,
Howard A. Smith,
S. P. Willner,
Andreas Zezas
Abstract:
Galaxy mergers are key events in galaxy evolution, often causing massive starbursts and fueling active galactic nuclei (AGN). In these highly dynamic systems, it is not yet precisely known how much starbursts and AGN respectively contribute to the total luminosity, at what interaction stages they occur, and how long they persist. Here we estimate the fraction of the bolometric infrared (IR) lumino…
▽ More
Galaxy mergers are key events in galaxy evolution, often causing massive starbursts and fueling active galactic nuclei (AGN). In these highly dynamic systems, it is not yet precisely known how much starbursts and AGN respectively contribute to the total luminosity, at what interaction stages they occur, and how long they persist. Here we estimate the fraction of the bolometric infrared (IR) luminosity that can be attributed to AGN by measuring and modeling the full ultraviolet to far-infrared spectral energy distributions (SEDs) in up to 33 broad bands for 24 merging galaxies with the Code for Investigating Galaxy Emission. In addition to a sample of 12 confirmed AGN in late-stage mergers, found in the $Infrared$ $Array$ $Satellite$ Revised Bright Galaxy Sample or Faint Source Catalog, our sample includes a comparison sample of 12 galaxy mergers from the $Spitzer$ Interacting Galaxies Survey, mostly early-stage. We perform identical SED modeling of simulated mergers to validate our methods, and we supplement the SED data with mid-IR spectra of diagnostic lines obtained with $Spitzer$ InfraRed Spectrograph. The estimated AGN contributions to the IR luminosities vary from system to system from 0% up to 91% but are significantly greater in the later-stage, more luminous mergers, consistent with what is known about galaxy evolution and AGN triggering.
△ Less
Submitted 29 June, 2018; v1 submitted 12 January, 2018;
originally announced January 2018.
-
Young Stellar Objects in the Massive Star-Forming Regions W51 and W43
Authors:
G. Saral,
J. L. Hora,
M. Audard,
X. P. Koenig,
J. R. Martínez-Galarza,
F. Motte,
Q. Nguyen-Luong,
A. T. Saygac,
H. A. Smith
Abstract:
We present the results of our investigation of the star-forming complexes W51 and W43, two of the brightest in the first Galactic quadrant. In order to determine the young stellar object (YSO) populations in W51 and W43 we used color-magnitude relations based on Spitzer mid-infrared and 2MASS/UKIDSS near-infrared data. We identified 302 Class I YSOs and 1178 Class II/transition disk candidates in…
▽ More
We present the results of our investigation of the star-forming complexes W51 and W43, two of the brightest in the first Galactic quadrant. In order to determine the young stellar object (YSO) populations in W51 and W43 we used color-magnitude relations based on Spitzer mid-infrared and 2MASS/UKIDSS near-infrared data. We identified 302 Class I YSOs and 1178 Class II/transition disk candidates in W51, and 917 Class I YSOs and 5187 Class II/transition disk candidates in W43. We also identified tens of groups of YSOs in both regions using the Minimal Spanning Tree (MST) method. We found similar cluster densities in both regions even though Spitzer was not able to probe the densest part of W43. By using the Class II/I ratios, we traced the relative ages within the regions and based on the morphology of the clusters we argue that several sites of star formation are independent of one another in terms of their ages and physical conditions. We used spectral energy distribution (SED)-fitting to identify the massive YSO (MYSO) candidates since they play a vital role in the star formation process and then examined them to see if they are related to any massive star formation tracers such as UCH II regions, masers and dense fragments. We identified 17 MYSO candidates in W51, and 14 in W43, respectively and found that groups of YSOs hosting MYSO candidates are positionally associated with H II regions in W51, though we do not see any MYSO candidates associated with previously identified massive dense fragments in W43.
△ Less
Submitted 22 April, 2017; v1 submitted 8 March, 2017;
originally announced March 2017.
-
Merger Signatures in the Dynamics of Star-forming Gas
Authors:
Chao-Ling Hung,
Christopher C. Hayward,
Howard A. Smith,
Matthew L. N. Ashby,
Lauranne Lanz,
Juan R. Martínez-Galarza,
D. B. Sanders,
Andreas Zezas
Abstract:
Spatially resolved kinematics have been used to determine the dynamical status of star-forming galaxies with ambiguous morphologies, and constrain the importance of galaxy interactions during the assembly of galaxies. However, measuring the importance of interactions or galaxy merger rates requires knowledge of the systematics in kinematic diagnostics and the visible time with merger indicators. W…
▽ More
Spatially resolved kinematics have been used to determine the dynamical status of star-forming galaxies with ambiguous morphologies, and constrain the importance of galaxy interactions during the assembly of galaxies. However, measuring the importance of interactions or galaxy merger rates requires knowledge of the systematics in kinematic diagnostics and the visible time with merger indicators. We analyze the dynamics of star-forming gas in a set of binary merger hydrodynamic simulations with stellar mass ratios of 1:1 and 1:4. We find that the evolution of kinematic asymmetries traced by star-forming gas mirrors morphological asymmetries derived from mock optical images, in which both merger indicators show the largest deviation from isolated disks during strong interaction phases. Based on a series of simulations with various initial disk orientations, orbital parameters, gas fractions, and mass ratios, we find that the merger signatures are visible for ~0.2-0.4 Gyr with kinematic merger indicators but can be approximately twice as long for equal-mass mergers of massive gas-rich disk galaxies designed to be analogs of z~2-3 submillimeter galaxies. Merger signatures are most apparent after the second passage and before the black holes coalescence, but in some cases they persist up to several hundred Myr after coalescence. About 20-60% of the simulated galaxies are not identified as mergers during the strong interaction phase, implying that galaxies undergoing violent merging process do not necessarily exhibit highly asymmetric kinematics in their star-forming gas. The lack of identifiable merger signatures in this population can lead to an underestimation of merger abundances in star-forming galaxies, and including them in samples of star-forming disks may bias the measurements of disk properties such as intrinsic velocity dispersion.
△ Less
Submitted 26 November, 2015;
originally announced November 2015.
-
The Mid-Infrared Instrument for the James Webb Space Telescope, VI: The Medium Resolution Spectrometer
Authors:
Martyn Wells,
J. -W. Pel,
Alistair Glasse,
G. S. Wright,
Gabby Aitink-Kroes,
Ruyman Azzollini,
Steven Beard,
B. R. Brandl,
Angus Gallie,
V. C. Geers,
A. M. Glauser,
Peter Hastings,
Th. Henning,
Rieks Jager,
K. Justtanont,
Bob Kruizinga,
Fred Lahuis,
David Lee,
I. Martinez-Delgado,
J. R. Martinez-Galarza,
M. Meijers,
Jane E. Morrison,
Friedrich Mueller,
Thodori Nakos,
Brian O'Sullivan
, et al. (13 additional authors not shown)
Abstract:
We describe the design and performance of the Medium Resolution Spectrometer (MRS) for the JWST-MIRI instrument. The MRS incorporates four coaxial spectral channels in a compact opto-mechanical layout that generates spectral images over fields of view up to 7.7 X 7.7 arcseconds in extent and at spectral resolving powers ranging from 1,300 to 3,700. Each channel includes an all-reflective integral…
▽ More
We describe the design and performance of the Medium Resolution Spectrometer (MRS) for the JWST-MIRI instrument. The MRS incorporates four coaxial spectral channels in a compact opto-mechanical layout that generates spectral images over fields of view up to 7.7 X 7.7 arcseconds in extent and at spectral resolving powers ranging from 1,300 to 3,700. Each channel includes an all-reflective integral field unit (IFU): an 'image slicer' that reformats the input field for presentation to a grating spectrometer. Two 1024 X 1024 focal plane arrays record the output spectral images with an instantaneous spectral coverage of approximately one third of the full wavelength range of each channel. The full 5 to 28.5 micron spectrum is then obtained by making three exposures using gratings and pass-band-determining filters that are selected using just two three-position mechanisms. The expected on-orbit optical performance is presented, based on testing of the MIRI Flight Model and including spectral and spatial coverage and resolution. The point spread function of the reconstructed images is shown to be diffraction limited and the optical transmission is shown to be consistent with the design expectations.
△ Less
Submitted 12 August, 2015;
originally announced August 2015.
-
A WISE Census of Young Stellar Objects in Perseus OB2 Association
Authors:
Mohaddesseh Azimlu,
J. Rafael Martinez-Galarza,
August A. Muench
Abstract:
We have performed a WISE (Wide-Field Infrared Survey Explorer) based study to identify and characterize young stellar objects (YSOs) in 12x12 degree Perseus OB2 association. Spectral energy distribution (SED) slope in range of 3.4-12 micron and a 5sigma selection criteria were used to select our initial sample. Further manual inspection reduced our final catalog to 156 known and 119 YSO candidate.…
▽ More
We have performed a WISE (Wide-Field Infrared Survey Explorer) based study to identify and characterize young stellar objects (YSOs) in 12x12 degree Perseus OB2 association. Spectral energy distribution (SED) slope in range of 3.4-12 micron and a 5sigma selection criteria were used to select our initial sample. Further manual inspection reduced our final catalog to 156 known and 119 YSO candidate. The spatial distribution of newly found YSOs all over the field shows an older generation of star formation which most of its massive members have evolved into main sequence stars. In contrast, the majority of younger members lie within the Perseus molecular cloud and currently active star forming clusters such as NGC1333 and IC348. We also identified additional 66 point sources which passed YSO selection criteria but are likely AGB stars. However their spatial distribution suggests that they may contain a fraction of the YSOs. Comparing our results with the commonly used color-color selections, we found that while color selection method fails in picking up bright but evolved weak disks, our SED fitting method can identify such sources, including transitional disks. In addition we have less contamination with background sources such as galaxies, but in a price of loosing fainter (Jmag > 12) YSOs. Finally we employed a Bayesian Monte Carlo SED fitting method to determine the characteristics of each YSO candidate. Distribution of SED slopes and model driven age and mass confirms separated YSO populations with suggested three age groups of younger than 1 Myr old, 1-5 Myr old, and older than 5 Myrs which agrees with the age of Per OB2 association and currently star forming sites within the cloud.
△ Less
Submitted 10 July, 2015;
originally announced July 2015.
-
Variations of the ISM Compactness Across the Main Sequence of Star-Forming Galaxies: Observations and Simulations
Authors:
Juan R. Martínez-Galarza,
Howard A. Smith,
Lauranne Lanz,
Christopher C. Hayward,
Andreas Zezas,
Lee Rosenthal,
Aaron Weiner,
Chao-Ling Hung,
Matthew L. N. Ashby,
Brent Groves
Abstract:
(abridged) The majority of star-forming galaxies follow a simple empirical correlation in the star formation rate (SFR) versus stellar mass ($M_*$) plane, usually referred to as the star formation Main Sequence (MS). Here we combine a set of hydro-dynamical simulations of interacting galactic disks with state-of-the-art radiative transfer codes to analyze how the evolution of mergers is reflecte…
▽ More
(abridged) The majority of star-forming galaxies follow a simple empirical correlation in the star formation rate (SFR) versus stellar mass ($M_*$) plane, usually referred to as the star formation Main Sequence (MS). Here we combine a set of hydro-dynamical simulations of interacting galactic disks with state-of-the-art radiative transfer codes to analyze how the evolution of mergers is reflected upon the properties of the MS. We present \textsc{Chiburst}, a Markov Chain Monte Carlo (MCMC) Spectral Energy Distribution (SED) code that fits the multi-wavelength, broad-band photometry of galaxies and derives stellar masses, star formation rates, and geometrical properties of the dust distribution. We apply this tool to the SEDs of simulated mergers and compare the derived results with the reference output from the simulations. Our results indicate that changes in the SEDs of mergers as they approach coalescence and depart from the MS are related to an evolution of dust geometry in scales larger than a few hundred parsecs. This is reflected in a correlation between the specific star formation rate (sSFR), and the compactness parameter $\mathcal{C}$, that parametrizes this geometry and hence the evolution of dust temperature ($T_{\rm{dust}}$) with time. As mergers approach coalescence, they depart from the MS and increase their compactness, which implies that moderate outliers of the MS are consistent with late-type mergers. By further applying our method to real observations of Luminous Infrared Galaxies (LIRGs), we show that the merger scenario is unable to explain these extreme outliers of the MS. Only by significantly increasing the gas fraction in the simulations are we able to reproduce the SEDs of LIRGs.
△ Less
Submitted 29 December, 2015; v1 submitted 8 December, 2014;
originally announced December 2014.
-
Optical spectroscopic observations of gamma-ray blazars candidates I: preliminary results
Authors:
A. Paggi,
D. Milisavljevic,
N. Masetti,
E. Jiménez-Bailón,
V. Chavushyan,
R. D'Abrusco,
F. Massaro,
M. Giroletti,
H. A. Smith,
R. Margutti,
G. Tosti,
J. R. Martínez-Galarza,
H. Otí-Floranes,
M. Landoni,
J. E. Grindlay,
S. Funk
Abstract:
A significant fraction ($\sim 30$ %) of the gamma-ray sources listed in the second $\textit{Fermi}$ LAT (2FGL) catalog is still of unknown origin, being not yet associated with counterparts at lower energies. Using the available information at lower energies and optical spectroscopy on the selected counterparts of these gamma-ray objects we can pinpoint their exact nature. Here we present a pilot…
▽ More
A significant fraction ($\sim 30$ %) of the gamma-ray sources listed in the second $\textit{Fermi}$ LAT (2FGL) catalog is still of unknown origin, being not yet associated with counterparts at lower energies. Using the available information at lower energies and optical spectroscopy on the selected counterparts of these gamma-ray objects we can pinpoint their exact nature. Here we present a pilot project pointing to assess the effectiveness of the several classification methods developed to select gamma-ray blazar candidates. To this end, we report optical spectroscopic observations of a sample of 5 gamma-ray blazar candidates selected on the basis of their infrared WISE colors or of their low-frequency radio properties. Blazars come in two main classes: BL Lacs and FSRQs, showing similar optical spectra except for the stronger emission lines of the latter. For three of our sources the almost featureless optical spectra obtained confirm their BL Lac nature, while for the source WISEJ022051.24+250927.6 we observe emission lines with equivalent width $EW\sim 31$ $Å$, identifying it as a FSRQ with $z = 0.48$. The source WISEJ064459.38+603131.7, although not featuring a clear radio counterpart, shows a blazar-like spectrum with weak emission lines with $EW \sim 7$ $Å$, yielding a redshift estimate of $z=0.36$. In addition we report optical spectroscopic observations of 4 WISE sources associated with known gamma-ray blazars without a firm classification or redshift estimate. For all of these latter sources we confirm a BL Lac classification, with a tentative redshift estimate for the source WISEJ100800.81+062121.2 of $z = 0.65$.
△ Less
Submitted 21 April, 2014;
originally announced April 2014.
-
The total infrared luminosity may significantly overestimate the star formation rate of recently quenched galaxies
Authors:
Christopher C. Hayward,
Lauranne Lanz,
Matthew L. N. Ashby,
Giovanni Fazio,
Lars Hernquist,
Juan Rafael Martínez-Galarza,
Kai Noeske,
Howard A. Smith,
Stijn Wuyts,
Andreas Zezas
Abstract:
The total infrared (IR) luminosity is very useful for estimating the star formation rate (SFR) of galaxies, but converting the IR luminosity into an SFR relies on assumptions that do not hold for all galaxies. We test the effectiveness of the IR luminosity as an SFR indicator by applying it to synthetic spectral energy distributions generated from three-dimensional hydrodynamical simulations of is…
▽ More
The total infrared (IR) luminosity is very useful for estimating the star formation rate (SFR) of galaxies, but converting the IR luminosity into an SFR relies on assumptions that do not hold for all galaxies. We test the effectiveness of the IR luminosity as an SFR indicator by applying it to synthetic spectral energy distributions generated from three-dimensional hydrodynamical simulations of isolated disc galaxies and galaxy mergers. In general, the SFR inferred from the IR luminosity agrees well with the true instantaneous SFR of the simulated galaxies. However, for the major mergers in which a strong starburst is induced, the SFR inferred from the IR luminosity can overestimate the instantaneous SFR during the post-starburst phase by greater than two orders of magnitude. Even though the instantaneous SFR decreases rapidly after the starburst, the stars that were formed in the starburst remain dust-obscured and thus produce significant IR luminosity. Consequently, use of the IR luminosity as an SFR indicator may cause one to conclude that post-starburst galaxies are still star-forming, whereas in reality, star formation was recently quenched.
△ Less
Submitted 19 September, 2014; v1 submitted 31 January, 2014;
originally announced February 2014.
-
Ongoing Massive Star Formation in NGC 604
Authors:
Juan R. Martínez-Galarza,
Deidre Hunter,
Brent Groves,
Bernhard Brandl
Abstract:
NGC 604 is the second most massive H II region in the Local Group, thus an important laboratory for massive star formation. Using a combination of observational and analytical tools that include Spitzer spectroscopy, Herschel photometry, Chandra imaging, and Bayesian Spectral Energy Distribution fitting, we investigate the physical conditions in NGC 604, and quantify the amount of massive star for…
▽ More
NGC 604 is the second most massive H II region in the Local Group, thus an important laboratory for massive star formation. Using a combination of observational and analytical tools that include Spitzer spectroscopy, Herschel photometry, Chandra imaging, and Bayesian Spectral Energy Distribution fitting, we investigate the physical conditions in NGC 604, and quantify the amount of massive star formation currently taking place. We derive an average age of 4 +/- 1 Myr and a total stellar mass of 1.6 (+1.6)(-1.0) x 10^5 M_sun for the entire region, in agreement with previous optical studies. Across the region we find an effect of the X-ray field on both the abundance of aromatic molecules and the [Si II] emission. Within NGC 604 we identify several individual bright infrared sources with diameters of about 15 pc and luminosity weighted masses between 10^3 M_sun and 10^4 M_sun. Their spectral properties indicate that some of these sources are embedded clusters in process of formation, which together account for ~8% of the total stellar mass in the NGC 604 system. The variations of the radiation field strength across NGC 604 are consistent with a sequential star formation scenario, with at least two bursts in the last few million years. Our results indicate that massive star formation in NGC 604 is still ongoing, likely triggered by the earlier bursts.
△ Less
Submitted 24 October, 2012; v1 submitted 19 October, 2012;
originally announced October 2012.
-
The Physical Conditions in Starbursts Derived from Bayesian Fitting of Mid-IR SEDS: 30 Doradus as a Template
Authors:
J. R. Martínez-Galarza,
B. Groves,
B. Brandl,
G. E. de Messieres,
R. Indebetouw,
M. A. Dopita
Abstract:
To understand and interpret the observed Spectral Energy Distributions (SEDs) of starbursts, theoretical or semi-empirical SED models are necessary. Yet, while they are well-founded in theory, independent verification and calibration of these models, including the exploration of possible degeneracies between their parameters, are rarely made. As a consequence, a robust fitting method that leads to…
▽ More
To understand and interpret the observed Spectral Energy Distributions (SEDs) of starbursts, theoretical or semi-empirical SED models are necessary. Yet, while they are well-founded in theory, independent verification and calibration of these models, including the exploration of possible degeneracies between their parameters, are rarely made. As a consequence, a robust fitting method that leads to unique and reproducible results has been lacking. Here we introduce a novel approach based on Bayesian analysis to fit the Spitzer-IRS spectra of starbursts using the SED models proposed by Groves et al. (2008). We demonstrate its capabilities and verify the agreement between the derived best fit parameters and actual physical conditions by modelling the nearby, well-studied, giant HII region 30 Dor in the LMC. The derived physical parameters, such as cluster mass, cluster age, ISM pressure and covering fraction of photodissociation regions, are representative of the 30 Dor region. The inclusion of the emission lines in the modelling is crucial to break degeneracies. We investigate the limitations and uncertainties by modelling sub-regions, which are dominated by single components, within 30 Dor. A remarkable result for 30 Doradus in particular is a considerable contribution to its mid-infrared spectrum from hot ({\simeq} 300K) dust. The demonstrated success of our approach will allow us to derive the physical conditions in more distant, spatially unresolved starbursts.
△ Less
Submitted 29 June, 2011;
originally announced June 2011.
-
Wavelength calibration of the JWST-MIRI medium resolution spectrometer
Authors:
J. R. Martinez-Galarza,
A. M. Glauser,
A. Hernan-Caballero,
R. Azzollini,
A. Glasse,
S. Kendrew,
B. Brandl,
F. Lahuis
Abstract:
We present the wavelength and spectral resolution characterisation of the Integral Field Unit (IFU) Medium Resolution Spectrometer for the Mid-InfraRed Instrument (MIRI), to fly onboard the James Webb Space Telescope in 2014. We use data collected using the Verification Model of the instrument and develop an empirical method to calibrate properties such as wavelength range and resolving power in a…
▽ More
We present the wavelength and spectral resolution characterisation of the Integral Field Unit (IFU) Medium Resolution Spectrometer for the Mid-InfraRed Instrument (MIRI), to fly onboard the James Webb Space Telescope in 2014. We use data collected using the Verification Model of the instrument and develop an empirical method to calibrate properties such as wavelength range and resolving power in a portion of the spectrometer's full spectral range (5-28 microns). We test our results against optical models to verify the system requirements and combine them with a study of the fringing pattern in the instrument's detector to provide a more accurate calibration. We show that MIRI's IFU spectrometer will be able to produce spectra with a resolving power above R=2800 in the wavelength range 6.46-7.70 microns, and that the unresolved spectral lines are well fitted by a Gaussian profile.
△ Less
Submitted 27 June, 2010;
originally announced June 2010.
-
Infrared Emission by Dust Around lambda Bootis Stars: Debris Disks or Thermally Emitting Nebulae?
Authors:
J. R. Martínez-Galarza,
I. Kamp,
K. Y. L. Su,
A. Gáspár,
G. Rieke,
E. E. Mamajek
Abstract:
We present a model that describes stellar infrared excesses due to heating of the interstellar (IS) dust by a hot star passing through a diffuse IS cloud. This model is applied to six lambda Bootis stars with infrared excesses. Plausible values for the IS medium (ISM) density and relative velocity between the cloud and the star yield fits to the excess emission. This result is consistent with th…
▽ More
We present a model that describes stellar infrared excesses due to heating of the interstellar (IS) dust by a hot star passing through a diffuse IS cloud. This model is applied to six lambda Bootis stars with infrared excesses. Plausible values for the IS medium (ISM) density and relative velocity between the cloud and the star yield fits to the excess emission. This result is consistent with the diffusion/accretion hypothesis that lambda Bootis stars (A- to F-type stars with large underabundances of Fe-peak elements) owe their characteristics to interactions with the ISM. This proposal invokes radiation pressure from the star to repel the IS dust and excavate a paraboloidal dust cavity in the IS cloud, while the metal-poor gas is accreted onto the stellar photosphere. However, the measurements of the infrared excesses can also be fit by planetary debris disk models. A more detailed consideration of the conditions to produce lambda Bootis characteristics indicates that the majority of infrared-excess stars within the Local Bubble probably have debris disks. Nevertheless, more distant stars may often have excesses due to heating of interstellar material such as in our model.
△ Less
Submitted 11 December, 2008;
originally announced December 2008.
-
Modeling the Infrared Bow Shock at delta Velorum: Implications for Studies of Debris Disks and lambda Bootis Stars
Authors:
A. Gáspár,
K. Y. L. Su,
G. H. Rieke,
Z. Balog,
I. Kamp,
J. R. Martínez-Galarza,
K. Stapelfeldt
Abstract:
We have discovered a bow shock shaped mid-infrared excess region in front of delta Velorum using 24 micron observations obtained with the Multiband Imaging Photometer for Spitzer (MIPS). The excess has been classified as a debris disk from previous infrared observations. Although the bow shock morphology was only detected in the 24 micron observations, its excess was also resolved at 70 micron.…
▽ More
We have discovered a bow shock shaped mid-infrared excess region in front of delta Velorum using 24 micron observations obtained with the Multiband Imaging Photometer for Spitzer (MIPS). The excess has been classified as a debris disk from previous infrared observations. Although the bow shock morphology was only detected in the 24 micron observations, its excess was also resolved at 70 micron. We show that the stellar heating of an ambient interstellar medium (ISM) cloud can produce the measured flux. Since delta Velorum was classified as a debris disk star previously, our discovery may call into question the same classification of other stars. We model the interaction of the star and ISM, producing images that show the same geometry and surface brightness as is observed. The modeled ISM is 15 times overdense relative to the average Local Bubble value, which is surprising considering the close proximity (24 pc) of delta Velorum.
The abundance anomalies of lambda Bootis stars have been previously explained as arising from the same type of interaction of stars with the ISM. Low resolution optical spectra of delta Velorum show that it does not belong to this stellar class. The star therefore is an interesting testbed for the ISM accretion theory of the lambda Bootis phenomenon.
△ Less
Submitted 26 September, 2007;
originally announced September 2007.