Source author record

Eric D. Feigelson

Eric D. Feigelson appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

39works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

39 published item(s)

preprint2022arXiv

Evolution of X-ray Activity in <25 Myr Old Pre-Main Sequence Stars

Measuring the evolution of X-ray emission from pre-main sequence (PMS) stars gives insight into two issues: the response of magnetic dynamo processes to changes in interior structure and the effects of high-energy radiation on protoplanetary disks and primordial planetary atmospheres. We present a sample of 6,003 stars with ages 7-25Myr in ten nearby open clusters from Chandra X-ray and Gaia-EDR3 surveys. Combined with previous results in large samples of younger (<5Myr) stars in MYStIX and SFiNCs star forming regions, mass-stratified activity-age relations are derived for early phases of stellar evolution. X-ray luminosity (Lx) is constant during the first few Myr, possibly due to the presence of extended X-ray coronas insensitive to temporal changes in stellar size. Lx then decays during the 7-25Myr period, more rapidly as stellar mass increases. This decay is interpreted as decreasing efficiency of the alpha^2 dynamo as radiative cores grow and a solar-type alpha-Omega dynamo emerges. For more massive 3.5-7Mo fully radiative stars, the X-ray emission plummets indicating lack of an effective magnetic dynamo. The findings provide improved measurements of high energy radiation effects on circumstellar material, first the protoplanetary disk and then the atmospheres of young planets. The observed X-ray luminosities can be so high that an inner Earth-mass rocky, unmagnetized planet around a solar-mass PMS star might lose its primary and secondary atmospheres within a few-several million years. PMS X-ray emission may thus have a significant impact on evolution of early planetary atmospheres and the conditions promoting the rise of habitability.

preprint2022arXiv

Incorporating Measurement Error in Astronomical Object Classification

Most general-purpose classification methods, such as support-vector machine (SVM) and random forest (RF), fail to account for an unusual characteristic of astronomical data: known measurement error uncertainties. In astronomical data, this information is often given in the data but discarded because popular machine learning classifiers cannot incorporate it. We propose a simulation-based approach that incorporates heteroscedastic measurement error into existing classification method to better quantify uncertainty in classification. The proposed method first simulates perturbed realizations of the data from a Bayesian posterior predictive distribution of a Gaussian measurement error model. Then, a chosen classifier is fit to each simulation. The variation across the simulations naturally reflects the uncertainty propagated from the measurement errors in both labeled and unlabeled data sets. We demonstrate the use of this approach via two numerical studies. The first is a thorough simulation study applying the proposed procedure to SVM and RF, which are well-known hard and soft classifiers, respectively. The second study is a realistic classification problem of identifying high-$z$ $(2.9 \leq z \leq 5.1)$ quasar candidates from photometric data. The data are from merged catalogs of the Sloan Digital Sky Survey, the $Spitzer$ IRAC Equatorial Survey, and the $Spitzer$-HETDEX Exploratory Large-Area Survey. The proposed approach reveals that out of 11,847 high-$z$ quasar candidates identified by a random forest without incorporating measurement error, 3,146 are potential misclassifications with measurement error. Additionally, out of $1.85$ million objects not identified as high-$z$ quasars without measurement error, 936 can be considered new candidates with measurement error.

preprint2019arXiv

IRAS 09002-4732: A Laboratory for the Formation of Rich Stellar Clusters

IRAS 09002-4732 is a poorly studied embedded cluster of stars in the Vela Molecular Ridge at a distance of 1.7kpc. Deep observations with the Chandra X-ray Observatory, combined with existing optical and infrared surveys, produce a catalog of 441 probable pre-main sequence members of the region. The stellar spatial distribution has two components: most stars reside in a rich, compact, elliptical cluster, but a minority reside within a molecular filament several parsecs long that straddles the cluster. The filament has active distributed star formation with dozens of unclustered protostars. The cluster pre-main sequence population is $\leq 0.8$ Myr old and deeply embedded; its most massive member is extremely young producing an ultracompact H II region. The cluster total population deduced from the X-ray luminosity function is surprisingly rich, twice that of the Orion Nebula Cluster. The cluster core is remarkably dense where strong N-body interactions should be occurring; its Initial Mass Function may be deficient in massive stars. We infer that IRAS 09002-4732 is a rare case where a rich cluster is forming today in a molecular filament, consistent with astrophysical models of cluster formation in clouds that involve the hierarchical formation and merging of groups in molecular filaments.

preprint2016arXiv

Young Stellar Populations in MYStIX Star Forming Regions: Candidate Protostars

The Massive Young Star Forming Complex in Infrared and X-ray (MYStIX) project provides a new census on stellar members of massive star forming regions within 4 kpc. Here the MYStIX Infrared Excess catalog (MIRES) and Chandra-based X-ray photometric catalogs are mined to obtain high-quality samples of Class I protostars using criteria designed to reduce extragalactic and Galactic field star contamination. A total of 1,109 MYStIX Candidate Protostars (MCPs) are found in 14 star forming regions. Most are selected from protoplanetary disk infrared excess emission, but 20% are found from their ultrahard X-ray spectra from heavily absorbed magnetospheric flare emission. Two-thirds of the MCP sample is newly reported here. The resulting samples are strongly spatially associated with molecular cores and filaments on Herschel far-infrared maps. This spatial agreement and other evidence indicate that the MCP sample has high reliability with relatively few 'false positives' from contaminating populations. But the limited sensitivity and sparse overlap among the infrared and X-ray subsamples indicate that the sample is very incomplete with many 'false negatives'. Maps, tables, and source descriptions are provided to guide further study of star formation in these regions. In particular, the nature of ultrahard X-ray protostellar candidates without known infrared counterparts needs to be eludicated.

preprint2015arXiv

No evidence for protoplanetary disk destruction by OB stars in the MYStIX sample

HST images of proplyds in the Orion Nebula, as well as submillimeter/radio measurements, show that the dominant O7 star Theta1 Ori C photoevaporates nearby disks around pre-main sequence stars. Theory predicts that massive stars photoevaporate disks within distances of order 0.1 pc. These findings suggest that young, OB-dominated massive H II regions are inhospitable to the survival of protoplanetary disks, and subsequently to the formation and evolution of planets. In the current work, we test this hypothesis using large samples of pre-main sequence stars in 20 massive star-forming regions selected with X-ray and infrared photometry in the MYStIX survey. Complete disk destruction would lead to a deficit of cluster members with excess in JHKs and Spitzer/IRAC bands in the vicinity of O stars. In four MYStIX regions containing O stars and a sufficient surface density of disk-bearing sources to reliably test for spatial avoidance, we find no evidence for the depletion of inner disks around pre-main sequence stars in the vicinity of O-type stars, even very luminous O2-O5 stars. These results suggest that massive star-forming regions are not very hostile to the survival of protoplanetary disks and, presumably, to the formation of planets.

preprint2015arXiv

The Spatial Structure of Young Stellar Clusters. II. Total Young Stellar Populations

We investigate the intrinsic stellar populations (estimated total numbers of OB and pre-main-sequence stars down to 0.1 Mo) that are present in 17 massive star-forming regions (MSFRs) surveyed by the MYStIX project. The study is based on the catalog of >31,000 MYStIX Probable Complex Members with both disk-bearing and disk-free populations, compensating for extinction, nebulosity, and crowding effects. Correction for observational sensitivities is made using the X-ray Luminosity Function (XLF) and the near-infrared Initial Mass Function (IMF)--a correction that is often not made by infrared surveys of young stars. The resulting maps of the projected structure of the young stellar populations, in units of intrinsic stellar surface density, allow direct comparison between different regions. Several regions have multiple dense clumps, similar in size and density to the Orion Nebula Cluster. The highest projected density of ~34,000 stars/pc^2 is found in the core of the RCW38 cluster. Histograms of surface density show different ranges of values in different regions, supporting the conclusion of Bressert et al. (2010, B10) that no universal surface-density threshold can distinguish between clustered and distributed star-formation. However, a large component of the young stellar population of MSFRs resides in dense environments of 200-10,000 stars/pc^2 (including within the nearby Orion molecular clouds), and we find that there is no evidence for the B10 conclusion that such dense regions form an extreme "tail" of the distribution. Tables of intrinsic populations for these regions are used in our companion study of young cluster properties and evolution.

preprint2015arXiv

The Spatial Structure of Young Stellar Clusters. III. Physical Properties and Evolutionary States

We analyze the physical properties of stellar clusters that are detected in massive star-forming regions in the MYStIX project--a comparative, multiwavelength study of young stellar clusters within 3.6 kpc that contain at least one O-type star. Tabulated properties of subclusters in these regions include physical sizes and shapes, intrinsic numbers of stars, absorptions by the molecular clouds, and median subcluster ages. Physical signs of dynamical evolution are present in the relations of these properties, including statistically significant correlations between subcluster size, central density, and age, which are likely the result of cluster expansion after gas removal. We argue that many of the subclusters identified in Paper I are gravitationally bound because their radii are significantly less than what would be expected from freely expanding clumps of stars with a typical initial stellar velocity dispersion of ~3 km/s for star-forming regions. We explore a model for cluster formation in which structurally simpler clusters are built up hierarchically through the mergers of subclusters--subcluster mergers are indicated by an inverse relation between the numbers of stars in a subcluster and their central densities (also seen as a density vs. radius relation that is less steep than would be expected from pure expansion). We discuss implications of these effects for the dynamical relaxation of young stellar clusters.

preprint2014arXiv

Age Gradients in the Stellar Populations of Massive Star Forming Regions Based on a New Stellar Chronometer

A major impediment to understanding star formation in massive star forming regions (MSFRs) is the absence of a reliable stellar chronometer to unravel their complex star formation histories. We present a new estimation of stellar ages using a new method that employs near-infrared (NIR) and X-ray photometry, AgeJX. Stellar masses are derived from X-ray luminosities using the Lx - Mass relation from the Taurus cloud. J-band luminosities are compared to mass-dependent pre-main-sequence evolutionary models to estimate ages. AgeJX is sensitive to a wide range of evolutionary stages, from disk-bearing stars embedded in a cloud to widely dispersed older pre-main sequence stars. The MYStIX (Massive Young Star-Forming Complex Study in Infrared and X-ray) project characterizes 20 OB-dominated MSFRs using X-ray, mid-infrared, and NIR catalogs. The AgeJX method has been applied to 5525 out of 31,784 MYStIX Probable Complex Members. We provide a homogeneous set of median ages for over a hundred subclusters in 15 MSFRs; median subcluster ages range between 0.5 Myr and 5 Myr. The important science result is the discovery of age gradients across MYStIX regions. The wide MSFR age distribution appears as spatially segregated structures with different ages. The AgeJX ages are youngest in obscured locations in molecular clouds, intermediate in revealed stellar clusters, and oldest in distributed populations. The NIR color index J-H, a surrogate measure of extinction, can serve as an approximate age predictor for young embedded clusters.

preprint2014arXiv

Core-Halo Age Gradients and Star Formation in the Orion Nebula and NGC~2024 Young Stellar Clusters

We analyze age distributions of two nearby rich stellar clusters, the NGC 2024 (Flame Nebula) and Orion Nebula Cluster (ONC) in the Orion molecular cloud complex. Our analysis is based on samples from the MYStIX survey and a new estimator of pre-main sequence (PMS) stellar ages, AgeJX, derived from X-ray and near-infrared photometric data. To overcome the problem of uncertain individual ages and large spreads of age distributions for entire clusters, we compute median ages and their confidence intervals of stellar samples within annular subregions of the clusters. We find core-halo age gradients in both the NGC 2024 cluster and ONC: PMS stars in cluster cores appear younger and thus were formed later than PMS stars in cluster peripheries. These findings are further supported by the spatial gradients in the disk fraction and K-band excess frequency. Our age analysis is based on AgeJX estimates for PMS stars, and is independent of any consideration of OB stars. The result has important implications for the formation of young stellar clusters. One basic implication is that clusters form slowly and the apparent age spreads in young stellar clusters, which are often controversial, are (at least in part) real. The result further implies that simple models where clusters form inside-out are incorrect, and more complex models are needed. We provide several star formation scenarios that alone or in combination may lead to the observed core-halo age gradients.

preprint2014arXiv

The Massive Star-forming Regions Omnibus X-ray Catalog

We present the Massive Star-forming Regions (MSFRs) Omnibus X-ray Catalog (MOXC), a compendium of X-ray point sources from {\em Chandra}/ACIS observations of a selection of MSFRs across the Galaxy, plus 30 Doradus in the Large Magellanic Cloud. MOXC consists of 20,623 X-ray point sources from 12 MSFRs with distances ranging from 1.7 kpc to 50 kpc. Additionally, we show the morphology of the unresolved X-ray emission that remains after the catalogued X-ray point sources are excised from the ACIS data, in the context of \Spitzer\ and {\em WISE} observations that trace the bubbles, ionization fronts, and photon-dominated regions that characterize MSFRs. In previous work, we have found that this unresolved X-ray emission is dominated by hot plasma from massive star wind shocks. This diffuse X-ray emission is found in every MOXC MSFR, clearly demonstrating that massive star feedback (and the several-million-degree plasmas that it generates) is an integral component of MSFR physics.

preprint2014arXiv

The Spatial Structure of Young Stellar Clusters. I. Subclusters

The clusters of young stars in massive star-forming regions show a wide range of sizes, morphologies, and numbers of stars. Their highly subclustered structures are revealed by the MYStIX project's sample of 31,754 young stars in nearby sites of star formation (regions at distances <3.6 kpc that contain at least one O-type star.) In 17 of the regions surveyed by MYStIX, we identify subclusters of young stars using finite mixture models -- collections of isothermal ellipsoids that model individual subclusters. Maximum likelihood estimation is used to estimate the model parameters, and the Akaike Information Criterion is used to determine the number of subclusters. This procedure often successfully finds famous subclusters, such as the BN/KL complex behind the Orion Nebula Cluster and the KW-object complex in M17. A catalog of 142 subclusters is presented, with 1 to 20 subclusters per region. The subcluster core radius distribution for this sample is peaked at 0.17 pc with a standard deviation of 0.43 dex, and subcluster core radius is negatively correlated with gas/dust absorption of the stars -- a possible age effect. Based on the morphological arrangements of subclusters, we identify four classes of spatial structure: long chains of subclusters, clumpy structures, isolated clusters with a core-halo structure, and isolated clusters well fit by a single isothermal ellipsoid.

preprint2013arXiv

A Massive Young Star-Forming Complex Study in Infrared and X-ray: Mid-Infrared Observations and Catalogs

Spitzer IRAC observations and stellar photometric catalogs are presented for the Massive Young Star-Forming Complex Study in the Infrared and X-ray (MYStIX). MYStIX is a multiwavelength census of young stellar members of twenty nearby (d < 4 kpc), Galactic, star-forming regions that contain at least one O star. All regions have data available from the Spitzer Space Telescope, consisting of GLIMPSE or other published catalogs for eleven regions and results of our own photometric analysis of archival data for the remaining nine regions. This paper seeks to construct deep and reliable catalogs of sources from the Spitzer images. Mid-infrared study of these regions faces challenges of crowding and high nebulosity. Our new catalogs typically contain fainter sources than existing Spitzer studies, which improves the match rate to Chandra X-ray sources that are likely to be young stars, but increases the possibility of spurious point-source detections, especially peaks in the nebulosity. IRAC color-color diagrams help distinguish spurious detections of nebular PAH emission from the infrared excess associated with dusty disks around young stars. The distributions of sources on the mid-infrared color-magnitude and color-color diagrams reflect differences between MYStIX regions, including astrophysical effects such as stellar ages and disk evolution.

preprint2013arXiv

A Massive Young Star-Forming Complex Study in Infrared and X-ray: X-ray Sources in Ten Star Forming Regions

The Massive Young star-forming complex Study in Infrared and X-ray (MYStIX) uses data from the Chandra X-ray Observatory to identify and characterize the young stellar populations in twenty Galactic (d < 4 kpc) massive star-forming regions. Here, the X-ray analysis for Chandra ACIS-I observations of ten of the MYStIX fields is described, and a catalog of >10,000 X-ray sources is presented. In comparison to other published Chandra source lists for the same regions, the number of MYStIX detected faint X-ray sources in a region is often doubled. While the higher catalog sensitivity increases the chance of false detections, it also increases the number of matches to infrared stars. X-ray emitting contaminants include foreground stars, background stars, and extragalactic sources. The X-ray properties of sources in these classes are discussed.

preprint2013arXiv

Bayesian Matching for X-ray and Infrared Sources in the MYStIX Project

Identifying the infrared counterparts of X-ray sources in Galactic Plane fields such as those of the MYStIX project presents particular difficulties due to the high density of infrared sources. This high stellar density makes it inevitable that a large fraction of X-ray positions will have a faint field star close to them, which standard matching techniques may incorrectly take to be the counterpart. Instead we use the infrared data to create a model of both the field star and counterpart magnitude distributions, which we then combine with a Bayesian technique to yield a probability that any star is the counterpart of an X-ray source. In our more crowded fields, between 10 and 20% of counterparts that would be identified on the grounds of being the closest star to X-ray position within a 99% confidence error circle are instead identified by the Bayesian technique as field stars. These stars are preferentially concentrated at faint magnitudes. Equally importantly the technique also gives a probability that the true counterpart to the X-ray source falls beneath the magnitude limit of the infrared catalog. In deriving our method, we place it in the context of other procedures for matching astronomical catalogs.

preprint2013arXiv

Identifying Young Stars in Massive Star-Forming Regions for the MYStIX Project

The Massive Young star-forming Complex Study in Infrared and X-rays (MYStIX) project requires samples of young stars that are likely members of 20 nearby Galactic massive star-forming regions. Membership is inferred from statistical classification of X-ray sources, from detection of a robust infrared excess that is best explained by circumstellar dust in a disk or infalling envelope, and from published spectral types that are unlikely to be found among field stars. We present the MYStIX membership lists here, and describe in detail the statistical classification of X-ray sources via a \Naive Bayes Classi

preprint2013arXiv

New Organizations to Support Astroinformatics and Astrostatistics

In the past two years, the environment within which astronomers conduct their data analysis and management has rapidly changed. Working Groups associated with international societies and Big Data projects have emerged to support and stimulate the new fields of astroinformatics and astrostatistics. Sponsoring societies include the Intenational Statistical Institute, International Astronomical Union, American Astronomical Society, and Large Synoptic Survey Telescope project. They enthusiastically support cross-disciplinary activities where the advanced capabilities of computer science, statistics and related fields of applied mathematics are applied to advance research on planets, stars, galaxies and the Universe. The ADASS community is encouraged to join these organizations and to explore and engage in their public communication Web site, the Astrostatistics and Astroinformatics Portal (http://asaip.psu.edu).

preprint2013arXiv

Overview of the Massive Young Star-Forming Complex Study in Infrared and X-ray (MYStIX) Project

MYStIX (Massive Young Star-Forming Complex Study in Infrared and Xray) seeks to characterize 20 OB-dominated young clusters and their environs at distances d < 4 kpc using imaging detectors on the Chandra X-ray Observatory, Spitzer Space Telescope, and the United Kingdom InfraRed Telescope. The observational goals are to construct catalogs of star-forming complex stellar members with well-defined criteria, and maps of nebular gas (particularly of hot X-ray emitting plasma) and dust. A catalog of MYStIX Probable Complex Members (MPCMs) with several hundred OB stars and > 30, 000 low mass premain sequence is assembled. This sample and related data products will be used to seek new empirical constraints on theoretical models of cluster formation and dynamics, mass segregation, OB star formation, star formation triggering on the periphery of HII regions, the survivability of protoplanetary disks in HII regions. This paper give an introduction and overview of the project, covering the data analysis methodology and application to two star forming regions, NGC 2264 and the Trifid Nebula.

preprint2013arXiv

The MYStIX InfraRed-Excess Source Catalog

The MYStIX project (Massive Young Star-Forming Complex Study in Infrared and X-rays) provides a comparative study of 20 Galactic massive star-forming complexes (d = 0.4 to 3.6 kpc). Probable stellar members in each target complex are identified using X-ray and/or infrared data via two pathways: (1) X-ray detections of young/massive stars with coronal activity/strong winds; or (2) infrared excess (IRE) selection of young stellar objects (YSOs) with circumstellar disks and/or protostellar envelopes. We present the methodology for the second pathway, using Spitzer/IRAC, 2MASS, and UKIRT imaging and photometry. Although IRE selection of YSOs is welltrodden territory, MYStIX presents unique challenges. We combine IR spectral energy distribution (SED) fitting with IR color cuts and spatial clustering analysis to identify IRE sources and isolate probable YSO members in each MYStIX target field from the myriad types of contaminating sources that can resemble YSOs: extragalactic sources, evolved stars, nebular knots, and even unassociated foreground/background YSOs. Applying our methodology consistently across 18 of the target complexes, we produce the MYStIX IRE Source (MIRES) Catalog comprising 20,719 sources, including 8686 probable stellar members of the MYStIX target complexes. We also classify the SEDs of 9365 IR counterparts to MYStIX X-ray sources to assist the first pathway, the identification of X-ray detected stellar members. The MIRES catalog provides a foundation for follow-up studies of diverse phenomena related to massive star cluster formation, including protostellar outflows, circumstellar disks, and sequential star formation triggered by massive star feedback processes.

preprint2013arXiv

The Progression of Star Formation in the Rosette Molecular Cloud

Using Spitzer Space Telescope and Chandra X-ray Observatory data, we identify YSOs in the Rosette Molecular Cloud (RMC). By being able to select cluster members and classify them into YSO types, we are able to track the progression of star formation locally within the cluster environments and globally within the cloud. We employ nearest neighbor method (NNM) analysis to explore the density structure of the clusters and YSO ratio mapping to study age progressions in the cloud. We find a relationship between the YSO ratios and extinction which suggests star formation occurs preferentially in the densest parts of the cloud and that the column density of gas rapidly decreases as the region evolves. This suggests rapid removal of gas may account for the low star formation efficiencies observed in molecular clouds. We find that the overall age spread across the RMC is small. Our analysis suggests that star formation started throughout the complex around the same time. Age gradients in the cloud appear to be localized and any effect the HII region has on the star formation history is secondary to that of the primordial collapse of the cloud.

preprint2013arXiv

The Wide-Field Near Infrared Data: Optimal Photometry in Crowded Fields

We present JHK infrared data from the UK Infrared Telescope for a subset of the regions of the MYStIX (Massive Young Star-Forming Complex Study in Infrared and X-ray) survey. Some of the data were obtained specifically for the MYStIX project, and some as part of the UKIRT Infrared Deep Sky Survey's Galactic Plane Survey. In most of these fields crowding is a significant issue for aperture photometry, and so we have re-extracted the photometry from the processed images using an optimal extraction technique, and we describe how we adapt the optimal technique to mitigate the effects of crowding.

preprint2013arXiv

VOStat: A Statistical Web Service for Astronomers

VOStat is a Web service providing interactive statistical analysis of astronomical tabular datasets. It is integrated into the suite of analysis and visualization tools associated with the international Virtual Observatory (VO) through the SAMP communication system. A user supplies VOStat with a dataset extracted from the VO, or otherwise acquired, and chooses among $\sim 60$ statistical functions. These include data transformations, plots and summaries, density estimation, one- and two-sample hypothesis tests, global and local regressions, multivariate analysis and clustering, spatial analysis, directional statistics, survival analysis (for censored data like upper limits), and time series analysis. The statistical operations are performed using the public domain {\bf R} statistical software environment, including a small fraction of its $>4000$ {\bf CRAN} add-on packages. The purpose of VOStat is to facilitate a wider range of statistical analyses than are commonly used in astronomy, and to promote use of more advanced methodology in {\bf R} and {\bf CRAN}.

preprint2012arXiv

MYStIX First Results: Spatial Structures of Massive Young Stellar Clusters

Observations of the spatial distributions of young stars in star-forming regions can be linked to the theory of clustered star formation using spatial statistical methods. The MYStIX project provides rich samples of young stars from the nearest high-mass star-forming regions. Maps of stellar surface density reveal diverse structure and subclustering. Young stellar clusters and subclusters are fit with isothermal spheres and ellipsoids using the Bayesian Information Criterion to estimate the number of subclusters. Clustering is also investigated using Cartwright and Whitworth's Q statistic and the inhomogeneous two-point correlation function. Mass segregation is detected in several cases, in both centrally concentrated and fractally structured star clusters, but a few clusters are not mass segregated.

preprint2012arXiv

Statistical Methods for Astronomy

This review outlines concepts of mathematical statistics, elements of probability theory, hypothesis tests and point estimation for use in the analysis of modern astronomical data. Least squares, maximum likelihood, and Bayesian approaches to statistical inference are treated. Resampling methods, particularly the bootstrap, provide valuable procedures when distributions functions of statistics are not known. Several approaches to model selection and good- ness of fit are considered. Applied statistics relevant to astronomical research are briefly discussed: nonparametric methods for use when little is known about the behavior of the astronomical populations or processes; data smoothing with kernel density estimation and nonparametric regression; unsupervised clustering and supervised classification procedures for multivariate problems; survival analysis for astronomical datasets with nondetections; time- and frequency-domain times series analysis for light curves; and spatial statistics to interpret the spatial distributions of points in low dimensions. Two types of resources are presented: about 40 recommended texts and monographs in various fields of statistics, and the public domain R software system for statistical analysis. Together with its \sim 3500 (and growing) add-on CRAN packages, R implements a vast range of statistical procedures in a coherent high-level language with advanced graphics.

preprint2012arXiv

The Elephant Trunk Nebula and the Trumpler 37 cluster: Contribution of triggered star formation to the total population of an HII region

Rich young stellar clusters produce HII regions whose expansion into the nearby molecular cloud is thought to trigger the formation of new stars. However, the importance of this mode of star formation is uncertain. This investigation seeks to quantify triggered star formation (TSF) in IC 1396A (a.k.a., the Elephant Trunk Nebula), a bright rimmed cloud (BRC) on the periphery of the nearby giant HII region IC 1396 produced by the Trumpler 37 cluster. X-ray selection of young stars from Chandra X-ray Observatory data is combined with existing optical and infrared surveys to give a more complete census of the TSF population. Over 250 young stars in and around IC 1396A are identified; this doubles the previously known population. A spatio-temporal gradient of stars from the IC 1396A cloud toward the primary ionizing star HD 206267 is found. We argue that the TSF mechanism in IC 1396A is the radiation-driven implosion process persisting over several million years. Analysis of the X-ray luminosity and initial mass functions indicates that >140 stars down to 0.1 Msun were formed by TSF. Considering other BRCs in the IC 1396 HII region, we estimate the TSF contribution for the entire HII region exceeds 14-25% today, and may be higher over the lifetime of the HII region. Such triggering on the periphery of HII regions may be a significant mode of star formation in the Galaxy.

preprint2011arXiv

A Catalog of Chandra X-ray Sources in the Carina Nebula

We present a catalog of ~14,000 X-ray sources observed by the ACIS instrument on the Chandra X-ray Observatory within a 1.42 square degree survey of the Great Nebula in Carina, known as the Chandra Carina Complex Project (CCCP). This study appears in a Special Issue of the ApJS devoted to the CCCP. Here, we describe the data reduction and analysis procedures performed on the X-ray observations, including calibration and cleaning of the X-ray event data, point source detection, and source extraction. The catalog appears to be complete across most of the field to an absorption-corrected total-band luminosity of ~10^{30.7} erg/s for a typical low-mass pre-main sequence star. Counterparts to the X-ray sources are identified in a variety of visual, near-infrared, and mid-infrared surveys. The X-ray and infrared source properties presented here form the basis of many CCCP studies of the young stellar populations in Carina.

preprint2011arXiv

A Chandra ACIS Study of the Young Star Cluster Trumpler 15 in Carina and Correlation with Near-infrared Sources

Using the highest-resolution X-ray observation of the Trumpler 15 star cluster taken by the Chandra X-ray Observatory, we estimate the total size of its stellar population by comparing the X-ray luminosity function of the detected sources to a calibrator cluster, and identify for the first time a significant fraction (~14%) of its individual members. The highest-resolution near-IR observation of Trumpler 15 (taken by the HAWK-I instrument on the VLT) was found to detect most of our X-ray selected sample of cluster members, with a K-excess disk frequency of 3.8+-0.7%. The near-IR data, X-ray luminosity function, and published spectral types of the brightest members support a cluster age estimate (5-10 Myr) that is older than those for the nearby Trumpler 14 and Trumpler 16 clusters, and suggest that high-mass members may have already exploded as supernovae. The morphology of the inner ~0.7 pc core of the cluster is found to be spherical. However, the outer regions (beyond 2 pc) are elongated, forming an `envelope' of stars that, in projection, appears to connect Trumpler 15 to Trumpler 14; this morphology supports the view that these clusters are physically associated. Clear evidence of mass segregation is seen. This study appears in a Special Issue of the ApJS devoted to the Chandra Carina Complex Project (CCCP), a 1.42 square degree Chandra X-ray survey of the Great Nebula in Carina.

preprint2011arXiv

A Naive Bayes Source Classifier for X-ray Sources

The Chandra Carina Complex Project (CCCP) provides a sensitive X-ray survey of a nearby starburst region over >1 square degree in extent. Thousands of faint X-ray sources are found, many concentrated into rich young stellar clusters. However, significant contamination from unrelated Galactic and extragalactic sources is present in the X-ray catalog. We describe the use of a naive Bayes classifier to assign membership probabilities to individual sources, based on source location, X-ray properties, and visual/infrared properties. For the particular membership decision rule adopted, 75% of CCCP sources are classified as members, 11% are classified as contaminants, and 14% remain unclassified. The resulting sample of stars likely to be Carina members is used in several other studies, which appear in a Special Issue of the ApJS devoted to the CCCP.

preprint2011arXiv

Source Contamination in X-ray Studies of Star-Forming Regions: Application to the Chandra Carina Complex Project

We describe detailed simulations of X-ray-emitting populations to evaluate the levels of contamination by both Galactic and extragalactic X-ray sources unrelated to a star-forming region under study. For Galactic contaminations, we consider contribution from main-sequence stars and giants (not including cataclysmic variables and other classes of accretion-driven X-ray binary systems) as they make the dominant contribution at the position of the Carina Nebula. The simulations take into consideration a variety of technical factors involving a Galactic population synthesis model, stellar X-ray luminosity functions, Chandra telescope response, source detection methodology, and possible spatial variations in the X-ray background and absorption through molecular clouds. When applied to the 1.42 square-degree field of the Chandra Carina Complex Project (CCCP), the simulations predict ~5000 contaminating sources (1 source per square arcminute of the survey), evenly distributed across the field. The results of the simulations are further employed in a companion CCCP study to assign membership probabilities to individual sources.

preprint2011arXiv

The Chandra Carina Complex Project View of Trumpler 16

Trumpler 16 is a well--known rich star cluster containing the eruptive supergiant $η$ Carinæ and located in the Carina star-forming complex. In the context of the Chandra Carina Complex Project, we study Trumpler 16 using new and archival X-ray data. A revised X-ray source list of the Trumpler 16 region contains 1232 X-ray sources including 1187 likely Carina members. These are matched to 1047 near-infrared counterparts detected by the HAWK-I instrument at the VLT allowing for better selection of cluster members. The cluster is irregular in shape. Although it is roughly circular, there is a high degree of sub-clustering, no noticeable central concentration and an extension to the southeast. The high--mass stars show neither evidence of mass segregation nor evidence of strong differential extinction. The derived power-law slope of the X-ray luminosity function for Trumpler 16 reveals a much steeper function than the Orion Nebula Cluster implying different ratio of solar- to higher-mass stars. We estimate the total Trumpler 16 pre-main sequence population to be > 6500 Class II and Class III X-ray sources. An overall K-excess disk frequency of ~ 8.9% is derived using the X-ray selected sample, although there is some variation among the sub-clusters, especially in the Southeastern extension. X-ray emission is detected from 29 high--mass stars with spectral types between B2 and O3.

preprint2011arXiv

X-ray Star Clusters in the Carina Complex

The distribution of young stars found in the Chandra Carina Complex Project (CCCP) is examined for clustering structure. X-ray surveys are advantageous for identifying young stellar populations compared to optical and infrared surveys in suffering less contamination from nebular emission and Galactic field stars. The analysis is based on smoothed maps of a spatially complete subsample of about 3000 brighter X-ray sources classified as Carina members, and about 10,000 stars from the full CCCP sample. The principal known clusters are recovered, and some additional smaller groups are identified. No rich embedded clusters are present, although a number of sparse groups are found. The CCCP reveals considerable complexity in clustering properties. The Trumpler 14 and 15 clusters have rich stellar populations in unimodal, centrally concentrated structures several parsecs across. Non-spherical internal structure is seen, and large-scale low surface density distributions surround these rich clusters. Trumpler 16, in contrast, is comprised of several smaller clusters within a circular boundary. Collinder 228 is a third type of cluster which extends over tens of parsecs with many sparse compact groups likely arising from triggered star formation processes. A widely dispersed, but highly populous, distribution of X-ray stars across the about 50 pc CCCP mosaic supports a model of past generations of star formation in the region. Collinder 234, a group of massive stars without an associated cluster of pre-main sequence stars, may be part of this dispersed population.

preprint2010arXiv

A Chandra Observation of the Obscured Star-Forming Complex W40

The young stellar cluster illuminating the W40 H II region, one of the nearest massive star forming regions, has been observed with the ACIS detector on board the Chandra X-ray Observatory. Due to its high obscuration, this is a poorly-studied stellar cluster with only a handful of bright stars visible in the optical band, including three OB stars identified as primary excitation sources. We detect 225 X-ray sources, of which 85% are confidently identified as young stellar members of the region. Two potential distances of the cluster, 260 pc and 600 pc, are used in the paper. Supposing the X-ray luminosity function to be universal, it supports a 600 pc distance as a lower limit for W40 and a total population of at least 600 stars down to 0.1 Mo under the assumption of a coeval population with a uniform obscuration. In fact, there is strong spatial variation in Ks-band-excess disk fraction and non-uniform obscuration due to a dust lane that is identified in absorption in optical, infrared and X-ray. The dust lane is likely part of a ring of material which includes the molecular core within W40. In contrast to the likely ongoing star formation in the dust lane, the molecular core is inactive. The star cluster has a spherical morphology, an isothermal sphere density profile, and mass segregation down to 1.5 Mo. However, other cluster properties, including a \leq{1} Myr age estimate and ongoing star formation, indicate that the cluster is not dynamically relaxed. X-ray diffuse emission and a powerful flare from a young stellar object are also reported.

preprint2010arXiv

A Chandra Study of the Rosette Star-Forming Complex. III. The NGC 2237 Cluster and the Region's Star Formation History

We present Chandra X-ray images of the NGC 2237 young star cluster on the periphery of the Rosette Nebula. We detect 168 X-ray sources, 80% of which have counterparts in USNO, 2MASS, and deep FLAMINGOS images. These constitute the first census of the cluster members with 0.2<~M<~2 Msun. Star locations in near-infrared color-magnitude diagrams indicate a cluster age around 2 Myr with a visual extinction of 1<Av<3 at 1.4 kpc, the distance of the Rosette Nebula's main cluster NGC 2244. We derive the K-band luminosity function and the X-ray luminosity function of the cluster, which indicate a population ~400-600 stars. The X-ray-selected sample shows a K-excess disk frequency of 13%. The young Class II counterparts are aligned in an arc ~3 pc long suggestive of a triggered formation process induced by the O stars in NGC 2244. The diskless Class III sources are more dispersed. Several X-ray emitting stars are located inside the molecular cloud and around gaseous pillars projecting from the cloud. These stars, together with a previously unreported optical outflow originating inside the cloud, indicate that star formation is continuing at a low level and the cluster is still growing. This X-ray view of young stars on the western side of the Rosette Nebula complements our earlier studies of the central cluster NGC 2244 and the embedded clusters on the eastern side of the Nebula. The large scale distribution of the clusters and molecular material is consistent with a scenario in which the rich central NGC 2244 cluster formed first, and its expanding HII region triggered the formation of the now-unobscured clusters RMC XA and NGC 2237. A large swept-up shell material around the HII region is now in a second phase of collect-and-collapse fragmentation, leading to the recent formation of subclusters. Other clusters deeper in the molecular cloud appear unaffected by the Nebula expansion.

preprint2010arXiv

Cross-disciplinary research in astronomy

In the distant past, astronomy was often intertwined with religion into a unified cosmos. As science became a distinct cultural enterprise, astronomy has witnessed a variety of rich interactions with other fields. Mathematical statistics was stimulated in the 19th century by astronomical problems, and today astrostatistics is a small but growing cross-disciplinary field advancing methodology to address challenges in astronomical data analysis. Throughout the 20th century, astronomy became closely allied with physics such that astronomy and astrophysics are now profoundly intertwined. Physical chemistry played a major role in the identification of molecules in the Milky Way Galaxy, and astrochemistry is now an active subfield giving insights into cosmic molecular processes. The importance of cross-disciplinary interactions with engineering (for instrumentation), Earth sciences (for planetary studies), computer science (for astroinformatics) and life sciences (for astrobiology) is also growing. Cross-disciplinary research has been essential both for crucial discoveries in astronomy and for improving the quality of astronomical research. It should be fostered with increased flexibility in the training of young astronomers and with sufficient funding to nurture these fields.

preprint2010arXiv

Innovations in the Analysis of Chandra-ACIS Observations

As members of the instrument team for the Advanced CCD Imaging Spectrometer (ACIS) on NASA's Chandra X-ray Observatory and as Chandra General Observers, we have developed a wide variety of data analysis methods that we believe are useful to the Chandra community, and have constructed a significant body of publicly-available software (the ACIS Extract package) addressing important ACIS data and science analysis tasks. This paper seeks to describe these data analysis methods for two purposes: to document the data analysis work performed in our own science projects, and to help other ACIS observers judge whether these methods may be useful in their own projects (regardless of what tools and procedures they choose to implement those methods). The ACIS data analysis recommendations we offer here address much of the workflow in a typical ACIS project, including data preparation, point source detection via both wavelet decomposition and image reconstruction, masking point sources, identification of diffuse structures, event extraction for both point and diffuse sources, merging extractions from multiple observations, nonparametric broad-band photometry, analysis of low-count spectra, and automation of these tasks. Many of the innovations presented here arise from several, often interwoven, complications that are found in many Chandra projects: large numbers of point sources (hundreds to several thousand), faint point sources, misaligned multiple observations of an astronomical field, point source crowding, and scientifically relevant diffuse emission.

preprint2009arXiv

A Chandra Study of the Rosette Star-Forming Complex. II. Clusters in the Rosette Molecular Cloud

We explore here the young stellar populations in the Rosette Molecular Cloud (RMC) region with high spatial resolution X-ray images from the Chandra X-ray Observatory, which are effective in locating weak-lined T Tauri stars as well as disk-bearing young stars. A total of 395 X-ray point sources are detected, 299 of which (76%) have an optical or near-infrared (NIR) counterpart identified from deep FLAMINGOS images. From X-ray and mass sensitivity limits, we infer a total population of about 1700 young stars in the survey region. Based on smoothed stellar surface density maps, we investigate the spatial distribution of the X-ray sources and define three distinctive structures and substructures within them. Structures B and C are associated with previously known embedded IR clusters, while structure A is a new X-ray-identified unobscured cluster. A high mass protostar RMCX #89 = IRAS 06306+0437 and its associated sparse cluster is studied. The different subregions are not coeval but do not show a simple spatial-age pattern. Disk fractions vary between subregions and are generally 20% of the total stellar population inferred from the X-ray survey. The data are consistent with speculations that triggered star formation around the HII region is present in the RMC, but do not support a simple sequential triggering process through the cloud interior. While a significant fraction of young stars are located in a distributed population throughout the RMC region, it is not clear they originated in clustered environments.

preprint2009arXiv

Methods for Estimating Fluxes and Absorptions of Faint X-ray Sources

X-ray sources with very few counts can be identified with low-noise X-ray detectors such as ACIS onboard the Chandra X-ray Observatory. These sources are often too faint for parametric spectral modeling using well-established methods such as fitting with XSPEC. We discuss the estimation of apparent and intrinsic broad-band X-ray fluxes and soft X-ray absorption from gas along the line of sight to these sources, using nonparametric methods. Apparent flux is estimated from the ratio of the source count rate to the instrumental effective area averaged over the chosen band. Absorption, intrinsic flux, and errors on these quantities are estimated from comparison of source photometric quantities with those of high S/N spectra that were simulated using spectral models characteristic of the class of astrophysical sources under study. The concept of this method is similar to the long-standing use of color-magnitude diagrams in optical and infrared astronomy, with X-ray median energy replacing color index and X-ray source counts replacing magnitude. Our nonparametric method is tested against the apparent spectra of 2000 faint sources in the Chandra observation of the rich young stellar cluster in the M17 HII region. We show that the intrinsic X-ray properties can be determined with little bias and reasonable accuracy using these observable photometric quantities without employing often uncertain and time-consuming methods of non-linear parametric spectral modeling. Our method is calibrated for thermal spectra characteristic of stars in young stellar clusters, but recalibration should be possible for some other classes of faint X-ray sources such as extragalactic AGN.

preprint2009arXiv

Stellar Clusters in the NGC 6334 Star Forming Complex

The full stellar population of NGC 6334, one of the most spectacular regions of massive star formation in the nearby Galaxy, have not been well-sampled in past studies. We analyze here a mosaic of two Chandra X-ray Observatory images of the region using sensitive data analysis methods, giving a list of 1607 faint X-ray sources with arcsecond positions and approximate line-of-sight absorption. About 95 percent of these are expected to be cluster members, most lower mass pre-main sequence stars. Extrapolating to low X-ray levels, the total stellar population is estimated to be 20-30,000 pre-main sequence stars. The X-ray sources show a complicated spatial pattern with about 10 distinct star clusters. The heavily-obscured clusters are mostly associated with previously known far-infrared sources and radio HII regions. The lightly-obscured clusters are mostly newly identified in the X-ray images. Dozens of likely OB stars are found, both in clusters and dispersed throughout the region, suggesting that star formation in the complex has proceeded over millions of years. A number of extraordinarily heavily absorbed X-ray sources are associated with the active regions of star formation.

preprint2005arXiv

Global X-ray properties of the Orion Nebula region

Based on the Chandra Orion Ultradeep Project (COUP) observation, we establish the global X-ray properties of the stellar population associated with the Orion Nebula. Three components contribute roughly equally to the integrated COUP luminosity in the hard (2-8 keV) X-ray band: several OB stars, 822 lightly obscured cool stars in the Orion Nebula Cluster (ONC), and 559 heavily obscured stars. ONC stars 0.5-2 pc from the center show a spatial asymmetry consistent with violent relaxation in the stellar dynamics. The obscured COUP sources concentrate around both OMC-1 molecular cores; these small-scale structures indicate ages t < 0.1 Myr. The X-ray luminosity function (XLF) of the lightly obscured sample is roughly lognormal in shape. The obscured population is deficient in lower-luminosity stars, perhaps due to localized circumstellar material. Mass-stratified XLFs show that one-third of the Orion Nebula region hard-band emission is produced by the bright O6 star theta-1 Ori C and half is produced by lower mass pre-main sequence stars with masses 0.3 < M < 3 Mo. Very low mass stars contribute little to the cluster X-ray emission. Using the hard band emission, we show that young stellar clusters like the ONC can be readily detected and resolved with Chandra across the Galactic disk, even in the presence of heavy obscuration. The Orion Nebula sample is a valuable template for studies of distant clusters. For example, the peak of the XLF shape can serve as a standard candle for a new distance measure to distant young stellar clusters, and the presence of a neon emission line complex around 1 keV can serve as a diagnostic for young stars.

preprint2005arXiv

The Origin of T Tauri X-ray Emission: New Insights from the Chandra Orion Ultradeep Project

We use the data of the Chandra Orion Ultradeep Project (COUP) to study the nearly 600 X-ray sources that can be reliably identified with optically well characterized T Tauri stars (TTS) in the Orion Nebula Cluster. We detect X-ray emission from more than 97% of the optically visible late-type (spectral types F to M) cluster stars. This proofs that there is no ``X-ray quiet'' population of late-type stars with suppressed magnetic activity. All TTS with known rotation periods lie in the saturated or super-saturated regime of the relation between activity and Rossby numbers seen for main-sequence (MS) stars, but the TTS show a much larger scatter in X-ray activity than seen for the MS stars. Strong near-linear relations between X-ray luminosities, bolometric luminosities and mass are present. We also find that the fractional X-ray luminosity rises slowly with mass over the 0.1 - 2 M_sun range. The plasma temperatures determined from the X-ray spectra of the TTS are much hotter than in MS stars, but seem to follow a general solar-stellar correlation between plasma temperature and activity level. The large scatter about the relations between X-ray activity and stellar parameters seems to be related to the influence of accretion on the X-ray emission. While the X-ray activity of the non-accreting TTS is consistent with that of rapidly rotating MS stars, the accreting stars are less X-ray active (by a factor of ~2-3 on average) and produce much less well defined correlations than the non-accretors. We discuss possible reasons for the suppression of X-ray emission by accretion and the implications of our findings on long-standing questions related to the origin of the X-ray emission from young stars.