Showing 1–50 of 160 results for author: Feigelson, E D

XUE. Molecular inventory in the inner region of an extremely irradiated Protoplanetary Disk

Authors: María Claudia Ramirez-Tannus, Arjan Bik, Lars Cuijpers, Rens Waters, Christiane Goppl, Thomas Henning, Inga Kamp, Thomas Preibisch, Konstantin V. Getman, Germán Chaparro, Pablo Cuartas-Restrepo, Alex de Koter, Eric D. Feigelson, Sierra L. Grant, Thomas J. Haworth, Sebastián Hernández, Michael A. Kuhn, Giulia Perotti, Matthew S. Povich, Megan Reiter, Veronica Roccatagliata, Elena Sabbi, Benoît Tabone, Andrew J. Winter, Anna F. McLeod , et al. (2 additional authors not shown)

Abstract: We present the first results of the eXtreme UV Environments (XUE) James Webb Space Telescope (JWST) program, that focuses on the characterization of planet forming disks in massive star forming regions. These regions are likely representative of the environment in which most planetary systems formed. Understanding the impact of environment on planet formation is critical in order to gain insights… ▽ More We present the first results of the eXtreme UV Environments (XUE) James Webb Space Telescope (JWST) program, that focuses on the characterization of planet forming disks in massive star forming regions. These regions are likely representative of the environment in which most planetary systems formed. Understanding the impact of environment on planet formation is critical in order to gain insights into the diversity of the observed exoplanet populations. XUE targets 15 disks in three areas of NGC 6357, which hosts numerous massive OB stars, among which some of the most massive stars in our Galaxy. Thanks to JWST we can, for the first time, study the effect of external irradiation on the inner ($< 10$ au), terrestrial-planet forming regions of proto-planetary disks. In this study, we report on the detection of abundant water, CO, CO$_2$, HCN and C$_2$H$_2$ in the inner few au of XUE 1, a highly irradiated disk in NGC 6357. In addition, small, partially crystalline silicate dust is present at the disk surface. The derived column densities, the oxygen-dominated gas-phase chemistry, and the presence of silicate dust are surprisingly similar to those found in inner disks located in nearby, relatively isolated low-mass star-forming regions. Our findings imply that the inner regions of highly irradiated disks can retain similar physical and chemical conditions as disks in low-mass star-forming regions, thus broadening the range of environments with similar conditions for inner disk rocky planet formation to the most extreme star-forming regions in our Galaxy. △ Less

Submitted 18 October, 2023; v1 submitted 17 October, 2023; originally announced October 2023.

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

A study of two periodogram algorithms for improving the detection of small transiting planets

Authors: Yash Gondhalekar, Eric D. Feigelson, Gabriel A. Caceres, Marco Montalto, Snehanshu Saha

Abstract: The sensitivities of two periodograms are compared for weak signal planet detection in transit surveys: the widely used Box-Least Squares (BLS) algorithm following light curve detrending and the Transit Comb Filter (TCF) algorithm following autoregressive ARIMA modeling. Small depth transits are injected into light curves with different simulated noise characteristics. Two measures of spectral pea… ▽ More The sensitivities of two periodograms are compared for weak signal planet detection in transit surveys: the widely used Box-Least Squares (BLS) algorithm following light curve detrending and the Transit Comb Filter (TCF) algorithm following autoregressive ARIMA modeling. Small depth transits are injected into light curves with different simulated noise characteristics. Two measures of spectral peak significance are examined: the periodogram signal-to-noise ratio (SNR) and a False Alarm Probability (FAP) based on the generalized extreme value distribution. The relative performance of the BLS and TCF algorithms for small planet detection is examined for a range of light curve characteristics, including orbital period, transit duration, depth, number of transits, and type of noise. We find that the TCF periodogram applied to ARIMA fit residuals with the SNR detection metric is preferred when short-memory autocorrelation is present in the detrended light curve and even when the light curve noise had white Gaussian noise. BLS is more sensitive to small planets only under limited circumstances with the FAP metric. BLS periodogram characteristics are inferior when autocorrelated noise is present due to heteroscedastic noise and false period detection. Application of these methods to TESS light curves with known small exoplanets confirms our simulation results. The study ends with a decision tree that advises transit survey scientists on procedures to detect small planets most efficiently. The use of ARIMA detrending and TCF periodograms can significantly improve the sensitivity of any transit survey with regularly spaced cadence. △ Less

Submitted 12 December, 2023; v1 submitted 8 August, 2023; originally announced August 2023.

Comments: 30 pages, 13 figures, accepted for publication in ApJL

iid2022: A Workshop on Statistical Methods for Event Data in Astronomy

Authors: Eric D. Feigelson, Massimiliano Bonamente

Abstract: We review the iid2022 workshop on statistical methods for X-ray and $γ$-ray astronomy and high--energy astrophysics event data in astronomy, held in Guntersville, AL, on Nov. 15-18 2022. New methods for faint source detection, spatial point processes, variability and spectral analysis, and machine learning are discussed. Ideas for future developments of advanced methodology are shared. We review the iid2022 workshop on statistical methods for X-ray and $γ$-ray astronomy and high--energy astrophysics event data in astronomy, held in Guntersville, AL, on Nov. 15-18 2022. New methods for faint source detection, spatial point processes, variability and spectral analysis, and machine learning are discussed. Ideas for future developments of advanced methodology are shared. △ Less

Submitted 3 August, 2023; originally announced August 2023.

Comments: 7 pages

Journal ref: Frontiers in Astronomy and Space Sciences, 10:1228508, 2023

An Evenly-Spaced LSST Cadence for Rapidly Variable Stars

Authors: Eric D. Feigelson, Federica B. Bianco, Rosaria Bonito

Abstract: Stars exhibit a bewildering variety of rapidly variable behaviors ranging from explosive magnetic flares to stochastically changing accretion to periodic pulsations or rotation. The principal Rubin Observatory Legacy Survey of Space and Time (LSST) surveys will have cadences too sparse and irregular to capture many of these phenomena. We propose here a LSST micro-survey to observe a single Galacti… ▽ More Stars exhibit a bewildering variety of rapidly variable behaviors ranging from explosive magnetic flares to stochastically changing accretion to periodic pulsations or rotation. The principal Rubin Observatory Legacy Survey of Space and Time (LSST) surveys will have cadences too sparse and irregular to capture many of these phenomena. We propose here a LSST micro-survey to observe a single Galactic field, rich in unobscured stars, in a continuous sequence of 30 second exposures for one long winter night in a single photometric band. The result will be a unique dataset of $\sim 1$ million regularly spaced stellar light curves (LCs). The LCs will constitute a comprehensive collection of late-type stellar flaring, but also other classes like short-period binary systems and cataclysmic variables, young stellar objects and ultra-short period exoplanets. An unknown variety of anomalous Solar System, Galactic and extragalactic variables and transients may also be present. A powerful array of statistical procedures can be applied to individual LCs from the long-standing fields of time series analysis, signal processing and econometrics. Dozens of `features' describing the variability can be extracted and the ensemble of light curves can be subject to advanced machine learning clustering procedures. This will give a unique, authoritative, objective taxonomy of the rapidly variable sky derived from identically cadenced LCs. This micro-survey is best performed early in the Rubin Observatory program, and the results can inform the wider community on the best approaches to variable star identification and classification from the sparse, irregular cadences that dominate the planned surveys. △ Less

Submitted 31 July, 2023; originally announced August 2023.

Comments: Accepted for publication in Astrophysical Journal Supplements. 11 pages

Magnetic Activity-Rotation-Age-Mass Relations in Late Pre-main Sequence Stars

Authors: Konstantin V. Getman, Eric D. Feigelson, Gordon P. Garmire

Abstract: We study the four-dimensional relationships between magnetic activity, rotation, mass and age for solar-type stars in the age range 5-25Myr. This is the late-pre-main sequence (l-PMS) evolutionary phase when rapid changes in star's interior may lead to the changes in magnetic dynamo mechanisms. We carefully derive rotational periods and spot sizes for 471 members of several l-PMS open clusters usi… ▽ More We study the four-dimensional relationships between magnetic activity, rotation, mass and age for solar-type stars in the age range 5-25Myr. This is the late-pre-main sequence (l-PMS) evolutionary phase when rapid changes in star's interior may lead to the changes in magnetic dynamo mechanisms. We carefully derive rotational periods and spot sizes for 471 members of several l-PMS open clusters using photometric light curves from the Zwicky Transient Facility. Magnetic activity was measured in our previous Chandra-based study, and additional rotational data were obtained from other work. Several results emerge. Mass-dependent evolution of rotation through the l-PMS phase agrees with astrophysical models of stellar angular momentum changes, although the data point to a subpopulation of stars with slower initial rotations than commonly assumed. There is a hint of the onset of unsaturated tachoclinal dependency of X-ray activity on rotation, as reported by Argiroffi et al. (2016), but this result is not confidently confirmed. Both X-ray luminosity and star spot area decrease approximately as t^{-1} for solar mass stars suggesting that spot magnetic fields are roughly constant and l-PMS stars follow the universal solar-scaling law between the X-ray luminosity and surface magnetic flux. Assuming convective dynamos are dominant, theoretical magnetic fluxes fail to reveal the universal law for l-PMS stars that enter late Henyey tracks. Altogether we emerge with a few lines of evidence suggesting that the transition from the turbulent to solar-type dynamo occurs at the later stages of l-PMS evolution as stars approach the Zero-Age Main Sequence. △ Less

Submitted 15 May, 2023; originally announced May 2023.

Comments: 29 pages, 11 figures, 4 tables. Accepted for publication in The Astrophysical Journal, May 15, 2022

DIAmante TESS AutoRegressive Planet Search (DTARPS): III. Understanding the DTARPS Candidate Transiting Planet Catalogs

Authors: Elizabeth J. Melton, Eric D. Feigelson, Marco Montalto, Gabriel A. Caceres, Andrew W. Rosenswie, Cullen S. Abelson

Abstract: The DIAmante TESS AutoRegressive Planet Search (DTARPS) project, using novel statistical methods, has identified several hundred candidates for transiting planetary systems obtained from 0.9 million Full Frame Image light curves obtained in the TESS Year 1 southern hemisphere survey (Melton et al. 2022a and 2022b). Several lines of evidence, including limited reconnaissance spectroscopy, indicate… ▽ More The DIAmante TESS AutoRegressive Planet Search (DTARPS) project, using novel statistical methods, has identified several hundred candidates for transiting planetary systems obtained from 0.9 million Full Frame Image light curves obtained in the TESS Year 1 southern hemisphere survey (Melton et al. 2022a and 2022b). Several lines of evidence, including limited reconnaissance spectroscopy, indicate that at least half are true planets rather than False Positives. Here various population properties of these objects are examined. Half of the DTARPS candidates are hot Neptunes, populating the 'Neptune desert' found in Kepler planet samples. The DTARPS samples also identify dozens of Ultra Short Period planets with orbital periods down to 5 hours, high priority systems for atmospheric transimssion spectroscopy, and planets orbiting low-mass M stars. DTARPS methodology is sufficiently well-characterized at each step that preliminary planet occurrence rates can be estimated. Except for the increase in hot Neptunes, DTARPS planet occurrence rates are consistent with Kepler rates. Overall, DTARPS provides one of the largest and most reliable catalog of TESS exoplanet candidates that can be tapped to improve our understanding of various exoplanetary populations and astrophysical processes. △ Less

Submitted 13 February, 2023; originally announced February 2023.

Comments: 29 pages, 16 figures, submitted to the AAS Journals February 13, 2023

DIAmante TESS AutoRegressive Planet Search (DTARPS): II. Hundreds of New TESS Candidate Exoplanets

Authors: Elizabeth J. Melton, Eric D. Feigelson, Marco Montalto, Gabriel A. Caceres, Andrew W. Rosenswie, Cullen S. Abelson

Abstract: The DIAmante TESS AutoRegressive Planet Search (DTARPS) project seeks to identify photometric transiting planets from 976,814 southern hemisphere stars observed in Year 1 of the TESS mission. This paper follows the methodology developed by Melton et al. (Paper I) using light curves extracted and pre-processed by the DIAmante project (Montalto et al. 2020). Paper I emerged with a list of 7,377 ligh… ▽ More The DIAmante TESS AutoRegressive Planet Search (DTARPS) project seeks to identify photometric transiting planets from 976,814 southern hemisphere stars observed in Year 1 of the TESS mission. This paper follows the methodology developed by Melton et al. (Paper I) using light curves extracted and pre-processed by the DIAmante project (Montalto et al. 2020). Paper I emerged with a list of 7,377 light curves with statistical properties characteristic of transiting planets but dominated by False Alarms and False Positives. Here a multistage vetting procedure is applied including: centroid motion and crowding metrics, False Alarm and False Positive reduction, photometric binary elimination, and ephemeris match removal. The vetting produces a catalog of 462 DTARPS Candidates across the southern ecliptic hemisphere and 310 objects in a spatially incomplete Galactic Plane list. Fifty-eight percent were not previously identified as transiting systems. Candidates are flagged for possible blending from nearby stars based on Zwicky Transient Facility data and for possible radial velocity variations based on Gaia satellite data. Orbital periods and planetary radii are refined using astrophysical modeling; the resulting parameters closely match published values for Confirmed Planets. Their properties are discussed in Paper III. △ Less

Submitted 10 April, 2024; v1 submitted 13 February, 2023; originally announced February 2023.

Comments: 34 pages, 11 figures. Published in The Astronomical Journal. Machine Readable Tables for the DTARPS-S Candidate Catalog is available at https://content.cld.iop.org/journals/1538-3881/167/5/203/revision1/ajad29f1t1_mrt.txt with Figure Set at .../ajad29f1figset5.tar.gz. The DTARPS-S Galactic Plane List is available at https://.../ajad29f1t4_mrt.txt with Figure Set .../ajad29f1figset10.tar.gz

Journal ref: Astronomical Journal, 167, 203, 2024

DIAmante TESS AutoRegressive Planet Search (DTARPS): I. Analysis of 0.9 Million Light Curves

Authors: Elizabeth J. Melton, Eric D. Feigelson, Marco Montalto, Gabriel A. Caceres, Andrew W. Rosenswie, Cullen S. Abelson

Abstract: Nearly one million light curves from the TESS Year 1 southern hemisphere extracted from Full Frame Images with the DIAmante pipeline are processed through the AutoRegressive Planet Search statistical procedure. ARIMA models remove trends and lingering autocorrelated noise, the Transit Comb Filter identifies the strongest periodic signal in the light curve, and a Random Forest machine learning clas… ▽ More Nearly one million light curves from the TESS Year 1 southern hemisphere extracted from Full Frame Images with the DIAmante pipeline are processed through the AutoRegressive Planet Search statistical procedure. ARIMA models remove trends and lingering autocorrelated noise, the Transit Comb Filter identifies the strongest periodic signal in the light curve, and a Random Forest machine learning classifier is trained and applied to identify the best potential candidates. Classifier training sets include injections of both planetary transit signals and contaminating eclipsing binaries. The optimized classifier has a True Positive Rate of 92.8% and a False Positive Rate of 0.37% from the labeled training set. The result of this DIAmante TESS autoregressive planet search (DTARPS) analysis is a list of 7,377 potential exoplanet candidates. The classifier has a False Positive Rate of 0.3%, a 64% recall rate for previously confirmed exoplanets, and a 78% negative recall rate for known False Positives. The completeness map of the injected planetary signals shows high recall rates for planets with 8 - 30 R(Earth) radii and periods 0.6-13 days and poor completeness for planets with radii < 2 R(Earth) or periods < 1 day. The list has many False Alarms and False Positives that need to be culled with multifaceted vetting operations (Paper II). △ Less

Submitted 11 April, 2024; v1 submitted 13 February, 2023; originally announced February 2023.

Comments: 34 pages, 24 figures. Published in The Astronomical Journal, 2024 May. A Machine Readable Table for Table 3 is available at https://content.cld.iop.org/journals/1538-3881/167/5/202/revision1/ajad29f0t3_mrt.txt

Journal ref: The Astronomical Journal, 167:203, 2024

Young stellar objects, accretion disks, and their variability with Rubin Observatory LSST

Authors: R. Bonito, L. Venuti, S. Ustamujic, P. Yoachim, R. A. Street, L. Prisinzano, P. Hartigan, M. G. Guarcello, K. G. Stassun, T. Giannini, E. D. Feigelson, A. Caratti o Garatti, S. Orlando, W. I. Clarkson, P. McGehee, E. C. Bellm, J. E. Gizis

Abstract: Vera C. Rubin Observatory, through the Legacy Survey of Space and Time (LSST), will allow us to derive a panchromatic view of variability in young stellar objects (YSOs) across all relevant timescales. Indeed, both short-term variability (on timescales of hours to days) and long-term variability (months to years), predominantly driven by the dynamics of accretion processes in disk-hosting YSOs, ca… ▽ More Vera C. Rubin Observatory, through the Legacy Survey of Space and Time (LSST), will allow us to derive a panchromatic view of variability in young stellar objects (YSOs) across all relevant timescales. Indeed, both short-term variability (on timescales of hours to days) and long-term variability (months to years), predominantly driven by the dynamics of accretion processes in disk-hosting YSOs, can be explored by taking advantage of the multi-band filters option available in Rubin LSST, in particular the $u,g,r,i$ filters that enable us to discriminate between photospheric stellar properties and accretion signatures. The homogeneity and depth of sky coverage that will be achieved with LSST will provide us with a unique opportunity to characterize the time evolution of disk accretion as a function of age and varying environmental conditions (e.g. field crowdedness, massive neighbors, metallicity), by targeting different star-forming regions. In this contribution to the Rubin LSST Survey Strategy Focus Issue, we discuss how implementing a dense observing cadence to explore short-term variability in YSOs represents a key complementary effort to the Wide-Fast-Deep observing mode that will be used to survey the sky over the full duration of the main survey ($\approx$10 years). The combination of these two modes will be vital to investigate the connection between the inner disk dynamics and longer-term eruptive variability behaviors, such as those observed on EXor objects. △ Less

Submitted 2 February, 2023; originally announced February 2023.

Comments: 11 pages, 4 figures, 1 table; accepted for publication in The Astrophysical Journal Supplement Series

DarkMix: Mixture Models for the Detection and Characterization of Dark Matter Halos

Authors: Lluís Hurtado-Gil, Michael A. Kuhn, Pablo Arnalte-Mur, Eric D. Feigelson, Vicent Martínez

Abstract: Dark matter simulations require statistical techniques to properly identify and classify their halos and structures. Nonparametric solutions provide catalogs of these structures but lack the additional learning of a model-based algorithm and might misclassify particles in merging situations. With mixture models, we can simultaneously fit multiple density profiles to the halos that are found in a d… ▽ More Dark matter simulations require statistical techniques to properly identify and classify their halos and structures. Nonparametric solutions provide catalogs of these structures but lack the additional learning of a model-based algorithm and might misclassify particles in merging situations. With mixture models, we can simultaneously fit multiple density profiles to the halos that are found in a dark matter simulation. In this work, we use the Einasto profile (Einasto 1965, 1968, 1969) to model the halos found in a sample of the Bolshoi simulation (Klypin et al. 2011), and we obtain their location, size, shape and mass. Our code is implemented in the R statistical software environment and can be accessed on https://github.com/LluisHGil/darkmix. △ Less

Submitted 8 August, 2022; originally announced August 2022.

Comments: 25 pages, 22 figures, 5 tables

Time domain methods for X-ray and gamma-ray astronomy

Authors: Eric D. Feigelson, Vinay L. Kashyap, Aneta Siemiginowska

Abstract: A variety of statistical methods for understanding variability in the time domain for low count rate X-ray and gamma-ray sources are explored. Variability can be detected using nonparametric (Anderson-Darling and overdispersion tests) and parametric (sequential likelihood-based tests) tools. Once detected, variability can be characterized by nonparametric (autocorrelation function, structure funct… ▽ More A variety of statistical methods for understanding variability in the time domain for low count rate X-ray and gamma-ray sources are explored. Variability can be detected using nonparametric (Anderson-Darling and overdispersion tests) and parametric (sequential likelihood-based tests) tools. Once detected, variability can be characterized by nonparametric (autocorrelation function, structure function,wavelet analysis) and parametric (multiple change point model such as Bayesian Blocks, integer autoregressive models, C-statistic and Poisson regression) methods. New multidimensional variability detection approaches are outlined. Software packages designed for high energy data analysis are deficient but tools are available in the R statistical software environment. Most of the methods presented here are not commonly used in high energy astronomy. △ Less

Submitted 16 March, 2022; originally announced March 2022.

Comments: 26 pages, 9 figures. To appear in Handbook for X-ray and Gamma-Ray Astrophysics, Volume 4: Analysis techniques, Section XVIII: Timing Analysis (Belloni & Bhattacharya, eds., Springer)

The Effect of Molecular Cloud Properties on the Kinematics of Stars Formed in the Trifid Region

Authors: Michael A. Kuhn, Lynne A. Hillenbrand, Eric D. Feigelson, Ian Fowler, Konstantin V. Getman, Patrick S. Broos, Matthew S. Povich, Mariusz Gromadzki

Abstract: The dynamical states of molecular clouds may affect the properties of the stars they form. In the vicinity of the Trifid Nebula ($d=1180\pm25$ pc), the main star cluster (Trifid Main) lies within an expanding section of the molecular cloud; however, ~0.3 deg to the north (Trifid North), the cloud's velocity structure is more tranquil. We acquired a Chandra X-ray observation to identify pre-main-se… ▽ More The dynamical states of molecular clouds may affect the properties of the stars they form. In the vicinity of the Trifid Nebula ($d=1180\pm25$ pc), the main star cluster (Trifid Main) lies within an expanding section of the molecular cloud; however, ~0.3 deg to the north (Trifid North), the cloud's velocity structure is more tranquil. We acquired a Chandra X-ray observation to identify pre-main-sequence stars in Trifid North, complementing a previous observation of Trifid Main. In Trifid North, we identified 51 candidate pre-main-sequence stars, of which 13 are high-confidence Trifid members based on Gaia EDR3 parallaxes and proper motions. We also re-analyzed membership of Trifid Main and separated out multiple background stellar associations. Trifid North represents a stellar population ~10% as rich as Trifid Main that formed in a separate part of the cloud. The 1D stellar velocity dispersion in Trifid North ($0.6\pm0.2$ km/s) is three times lower than in Trifid Main ($1.9\pm0.2$ km/s). Furthermore, in Trifid Main, proper motions indicate that the portion of the star cluster superimposed on the optical nebula is expanding. Expansion of the HII region around the O-star HD 164492A, and the resulting gas expulsion, can explain both the motions of the stars and gas in Trifid Main. Contrary to previous studies, we find no evidence that a cloud-cloud collision triggered star formation in the region. △ Less

Submitted 15 September, 2022; v1 submitted 14 March, 2022; originally announced March 2022.

Comments: Accepted for publication in ApJ. 26 pages, 18 figures, and 4 tables. Minor style changes have been applied to match the published version

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

Authors: Konstantin V. Getman, Eric D. Feigelson, Gordon P. Garmire, Patrick S. Broos, Michael A. Kuhn, Thomas Preibisch, Vladimir S. Airapetian

Abstract: 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… ▽ More 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. △ Less

Submitted 11 June, 2022; v1 submitted 3 March, 2022; originally announced March 2022.

Comments: 39 pages, 12 figures, 8 tables. Accepted for publication in The Astrophysical Journal, June 10, 2022

Incorporating Measurement Error in Astronomical Object Classification

Authors: Sarah Shy, Hyungsuk Tak, Eric D. Feigelson, John D. Timlin, G. Jogesh Babu

Abstract: 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… ▽ More 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. △ Less

Submitted 2 May, 2022; v1 submitted 13 December, 2021; originally announced December 2021.

Journal ref: The Astronomical Journal, 164 (1), 6, 2022

A hot Mars-sized exoplanet transiting an M dwarf

Authors: Caleb I. Cañas, Suvrath Mahadevan, William D. Cochran, Chad F. Bender, Eric D. Feigelson, C. E. Harman, Ravi Kumar Kopparapu, Gabriel A. Caceres, Scott A. Diddams, Michael Endl, Eric B. Ford, Samuel Halverson, Fred Hearty, Sinclaire Jones, Shubham Kanodia, Andrea S. J. Lin, Andrew J. Metcalf, Andrew Monson, Joe P. Ninan, Lawrence W. Ramsey, Paul Robertson, Arpita Roy, Christian Schwab, Guðmundur Stefánsson

Abstract: We validate the planetary nature of an ultra-short period planet orbiting the M dwarf KOI-4777. We use a combination of space-based photometry from Kepler, high-precision, near-infrared Doppler spectroscopy from the Habitable-zone Planet Finder, and adaptive optics imaging to characterize this system. KOI-4777.01 is a Mars-sized exoplanet ($\mathrm{R}_{p}=0.51 \pm 0.03R_{\oplus}$) orbiting the hos… ▽ More We validate the planetary nature of an ultra-short period planet orbiting the M dwarf KOI-4777. We use a combination of space-based photometry from Kepler, high-precision, near-infrared Doppler spectroscopy from the Habitable-zone Planet Finder, and adaptive optics imaging to characterize this system. KOI-4777.01 is a Mars-sized exoplanet ($\mathrm{R}_{p}=0.51 \pm 0.03R_{\oplus}$) orbiting the host star every 0.412-days ($\sim9.9$-hours). This is the smallest validated ultra-short period planet known and we see no evidence for additional massive companions using our HPF RVs. We constrain the upper $3σ$ mass to $M_{p}<0.34~\mathrm{M_\oplus}$ by assuming the planet is less dense than iron. Obtaining a mass measurement for KOI-4777.01 is beyond current instrumental capabilities. △ Less

Submitted 7 December, 2021; originally announced December 2021.

Comments: 25 pages, 8 figures, accepted for publication in AJ

Journal ref: AJ, 163, 3 (2022)

X-ray Super-Flares From Pre-Main Sequence Stars: Flare Modeling

Authors: Konstantin V. Getman, Eric D. Feigelson, Gordon P. Garmire

Abstract: Getman et al. (2021) reports the discovery, energetics, frequencies, and effects on environs of $>1000$ X-ray super-flares with X-ray energies $E_X \sim 10^{34}-10^{38}$~erg from pre-main sequence (PMS) stars identified in the $Chandra$ MYStIX and SFiNCs surveys. Here we perform detailed plasma evolution modeling of $55$ bright MYStIX/SFiNCs super-flares from these events. They constitute a large… ▽ More Getman et al. (2021) reports the discovery, energetics, frequencies, and effects on environs of $>1000$ X-ray super-flares with X-ray energies $E_X \sim 10^{34}-10^{38}$~erg from pre-main sequence (PMS) stars identified in the $Chandra$ MYStIX and SFiNCs surveys. Here we perform detailed plasma evolution modeling of $55$ bright MYStIX/SFiNCs super-flares from these events. They constitute a large sample of the most powerful stellar flares analyzed in a uniform fashion. They are compared with published X-ray super-flares from young stars in the Orion Nebula Cluster, older active stars, and the Sun. Several results emerge. First, the properties of PMS X-ray super-flares are independent of the presence or absence of protoplanetary disks inferred from infrared photometry, supporting the solar-type model of PMS flaring magnetic loops with both footpoints anchored in the stellar surface. Second, most PMS super-flares resemble solar long duration events (LDEs) that are associated with coronal mass ejections. Slow rise PMS super-flares are an interesting exception. Third, strong correlations of super-flare peak emission measure and plasma temperature with the stellar mass are similar to established correlations for the PMS X-ray emission composed of numerous smaller flares. Fourth, a new correlation of loop geometry is linked to stellar mass; more massive stars appear to have thicker flaring loops. Finally, the slope of a long-standing relationship between the X-ray luminosity and magnetic flux of various solar-stellar magnetic elements appears steeper in PMS super-flares than for solar events. △ Less

Submitted 20 July, 2021; v1 submitted 15 June, 2021; originally announced June 2021.

Comments: 30 pages, 11 figures, 2 tables. Accepted for publication in The Astrophysical Journal, July 2021

X-ray Super-Flares From Pre-Main Sequence Stars: Flare Energetics And Frequency

Authors: Konstantin V. Getman, Eric D. Feigelson

Abstract: Solar-type stars exhibit their highest levels of magnetic activity during their early convective pre-main sequence (PMS) phase of evolution. The most powerful PMS flares, super-flares and mega-flares, have peak X-ray luminosities of $\log(L_X)=30.5-34.0$~erg~s$^{-1}$ and total energies $\log(E_X)=34-38$~erg. Among $>24,000$ X-ray selected young ($t<5$~Myr) members of 40 nearby star-forming regions… ▽ More Solar-type stars exhibit their highest levels of magnetic activity during their early convective pre-main sequence (PMS) phase of evolution. The most powerful PMS flares, super-flares and mega-flares, have peak X-ray luminosities of $\log(L_X)=30.5-34.0$~erg~s$^{-1}$ and total energies $\log(E_X)=34-38$~erg. Among $>24,000$ X-ray selected young ($t<5$~Myr) members of 40 nearby star-forming regions from our earlier $Chandra$ MYStIX and SFiNCs surveys, we identify and analyze a well-defined sample of 1,086 X-ray super-flares and mega-flares, the largest sample ever studied. Most are considerably more powerful than optical/X-ray super-flares detected on main sequence stars. This study presents energy estimates of these X-ray flares and the properties of their host stars. These events are produced by young stars of all masses over evolutionary stages ranging from protostars to diskless stars, with the occurrence rate positively correlated with stellar mass. Flare properties are indistinguishable for disk-bearing and diskless stars indicating star-disk magnetic fields are not involved. A slope $α\simeq2$ in the flare energy distributions $dN/dE_X \propto E_X^{-α}$ is consistent with those of optical/X-ray flaring from older stars and the Sun. Mega-flares ($\log(E_X) > 36.2$~erg) from solar-mass stars have occurrence rate of $1.7_{-0.6}^{+1.0}$ flares/star/year and contribute at least $10-20$\% to the total PMS X-ray energetics. These explosive events may have important astrophysical effects on protoplanetary disk photoevaporation, ionization of disk gas, production of spallogenic radionuclides in disk solids, and hydrodynamic escape of young planetary atmospheres. Our following paper details plasma and magnetic loop modeling of the $>50$ brightest X-ray mega-flares. △ Less

Submitted 10 May, 2021; originally announced May 2021.

Comments: 36 pages, 12 figures, 4 tables. Accepted for publication in The Astrophysical Journal

21st Century Statistical and Computational Challenges in Astrophysics

Authors: Eric D. Feigelson, Rafael S. de Souza, Emille E. O. Ishida, Gutti Jogesh Babu

Abstract: Modern astronomy has been rapidly increasing our ability to see deeper into the universe, acquiring enormous samples of cosmic populations. Gaining astrophysical insights from these datasets requires a wide range of sophisticated statistical and machine learning methods. Long-standing problems in cosmology include characterization of galaxy clustering and estimation of galaxy distances from photom… ▽ More Modern astronomy has been rapidly increasing our ability to see deeper into the universe, acquiring enormous samples of cosmic populations. Gaining astrophysical insights from these datasets requires a wide range of sophisticated statistical and machine learning methods. Long-standing problems in cosmology include characterization of galaxy clustering and estimation of galaxy distances from photometric colors. Bayesian inference, central to linking astronomical data to nonlinear astrophysical models, addresses problems in solar physics, properties of star clusters, and exoplanet systems. Likelihood-free methods are growing in importance. Detection of faint signals in complicated noise is needed to find periodic behaviors in stars and detect explosive gravitational wave events. Open issues concern treatment of heteroscedastic measurement errors and understanding probability distributions characterizing astrophysical systems. The field of astrostatistics needs increased collaboration with statisticians in the design and analysis stages of research projects, and to jointly develop new statistical methodologies. Together, they will draw more astrophysical insights into astronomical populations and the cosmos itself. △ Less

Submitted 26 May, 2020; originally announced May 2020.

Comments: Accepted for publication in volume 8 of Annual Reviews of Statistics and Its Application. 26 pages, 7 figures

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.

Comments: arXiv admin note: substantial text overlap with arXiv:1905.05116

Report number: FERMILAB-FN-1093-A-AE-SCD

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

Authors: Konstantin V. Getman, Eric D. Feigelson, Michael A. Kuhn, Patrick S. Broos, Gordon P. Garmire

Abstract: 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 c… ▽ More 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. △ Less

Submitted 15 October, 2019; v1 submitted 10 October, 2019; originally announced October 2019.

Comments: 33 pages, 11 figures, 3 machine readable tables. To appear in The Astronomical Journal

Advanced Astrophysics Discovery Technology in the Era of Data Driven Astronomy

Authors: Richard K. Barry, Jogesh G. Babu, John G. Baker, Eric D. Feigelson, Amanpreet Kaur, Alan J. Kogut, Steven B. Kraemer, James P. Mason, Piyush Mehrotra, Gregory Olmschenk, Jeremy D. Schnittman, Amalie Stokholm, Eric R. Switzer, Brian A. Thomas, Raymond J. Walker

Abstract: Experience suggests that structural issues in how institutional Astrophysics approaches data-driven science and the development of discovery technology may be hampering the community's ability to respond effectively to a rapidly changing environment in which increasingly complex, heterogeneous datasets are challenging our existing information infrastructure and traditional approaches to analysis.… ▽ More Experience suggests that structural issues in how institutional Astrophysics approaches data-driven science and the development of discovery technology may be hampering the community's ability to respond effectively to a rapidly changing environment in which increasingly complex, heterogeneous datasets are challenging our existing information infrastructure and traditional approaches to analysis. We stand at the confluence of a new epoch of multimessenger science, remote co-location of data and processing power and new observing strategies based on miniaturized spacecraft. Significant effort will be required by the community to adapt to this rapidly evolving range of possible discovery moduses. In the suggested creation of a new Astrophysics element, Advanced Astrophysics Discovery Technology, we offer an affirmative solution that places the visibility of discovery technologies at a level that we suggest is fully commensurate with their importance to the future of the field. △ Less

Submitted 24 July, 2019; originally announced July 2019.

Comments: White paper submitted to the ASTRO2020 decadal survey

Gaia Stellar Kinematics in the Head of the Orion A Cloud: Runaway Stellar Groups and Gravitational Infall

Authors: K. V. Getman, E. D. Feigelson, M. A. Kuhn, G. P. Garmire

Abstract: This work extends previous kinematic studies of young stars in the Head of the Orion A cloud (OMC-1/2/3/4/5). It is based on large samples of infrared, optical, and X-ray selected pre-main sequence stars with reliable radial velocities and Gaia-derived parallaxes and proper motions. Stellar kinematic groups are identified assuming they mimic the motion of their parental gas. Several groups are fou… ▽ More This work extends previous kinematic studies of young stars in the Head of the Orion A cloud (OMC-1/2/3/4/5). It is based on large samples of infrared, optical, and X-ray selected pre-main sequence stars with reliable radial velocities and Gaia-derived parallaxes and proper motions. Stellar kinematic groups are identified assuming they mimic the motion of their parental gas. Several groups are found to have peculiar kinematics: the NGC 1977 cluster and two stellar groups in the Extended Orion Nebula (EON) cavity are caught in the act of departing their birthplaces. The abnormal motion of NGC 1977 may have been caused by a global hierarchical cloud collapse, feedback by massive Ori OB1ab stars, supersonic turbulence, cloud-cloud collision, and/or slingshot effect; the former two models are favored by us. EON groups might have inherited anomalous motions of their parental cloudlets due to small-scale `rocket effects' from nearby OB stars. We also identify sparse stellar groups to the east and west of Orion A that are drifting from the central region, possibly a slowly expanding halo of the Orion Nebula Cluster. We confirm previously reported findings of varying line-of-sight distances to different parts of the cloud's Head with associated differences in gas velocity. Three-dimensional movies of star kinematics show contraction of the groups of stars in OMC-1 and global contraction of OMC-123 stars. Overall, the Head of Orion A region exhibits complex motions consistent with theoretical models involving hierarchical gravitational collapse in (possibly turbulent) clouds with OB stellar feedback. △ Less

Submitted 24 May, 2019; originally announced May 2019.

Comments: Accepted for publication in MNRAS. 26 pages, 13 figures. The two 3-D stellar kinematic movies, aimed as Supplementary Materials, can be found on YouTube at: https://youtu.be/B4GHCVvCYfo (`restricted' sample) and https://youtu.be/6fUu8sP0QFI (`full' sample)

AutoRegressive Planet Search: Application to the Kepler Mission

Authors: Gabriel A. Caceres, Eric D. Feigelson, G. Jogesh Babu, Natalia Bahamonde, Alejandra Christen, Karine Bertin, Cristian Meza, Michel Curé

Abstract: The 4-year light curves of 156,717 stars observed with NASA's Kepler mission are analyzed using the AutoRegressive Planet Search (ARPS) methodology described by Caceres et al. (2019). The three stages of processing are: maximum likelihood ARIMA modeling of the light curves to reduce stellar brightness variations; constructing the Transit Comb Filter periodogram to identify transit-like periodic di… ▽ More The 4-year light curves of 156,717 stars observed with NASA's Kepler mission are analyzed using the AutoRegressive Planet Search (ARPS) methodology described by Caceres et al. (2019). The three stages of processing are: maximum likelihood ARIMA modeling of the light curves to reduce stellar brightness variations; constructing the Transit Comb Filter periodogram to identify transit-like periodic dips in the ARIMA residuals; Random Forest classification trained on Kepler Team confirmed planets using several dozen features from the analysis. Orbital periods between 0.2 and 100 days are examined. The result is a recovery of 76% of confirmed planets, 97% when period and transit depth constraints are added. The classifier is then applied to the full Kepler dataset; 1,004 previously noticed and 97 new stars have light curve criteria consistent with the confirmed planets, after subjective vetting removes clear False Alarms and False Positive cases. The 97 Kepler ARPS Candidate Transits mostly have periods $P<10$ days; many are UltraShort Period hot planets with radii $<1$% of the host star. Extensive tabular and graphical output from the ARPS time series analysis is provided to assist in other research relating to the Kepler sample. △ Less

Submitted 23 May, 2019; originally announced May 2019.

Comments: 66 pages with 26 figures and 5 tables to appear in the Astronomical Journal. A version with high-resolution graphics, machine readable tables and FigureSet for the 97KACTs is available at https://drive.google.com/drive/folders/107xZIAj3C0HHqsW66Xc-8FA6O30EncVt?usp=sharing

AutoRegressive Planet Search: Feasibility Study for Irregular Time Series

Authors: Andrew M. Stuhr, Eric D. Feigelson, Gabriel A. Caceres, Joel D. Hartman

Abstract: Sensitive signal processing methods are needed to detect transiting planets from ground-based photometric surveys. Caceres et al. (2019) show that the AutoRegressive Planet Search (ARPS) method --- a combination of autoregressive integrated moving average (ARIMA) parametric modeling, a new Transit Comb Filter (TCF) periodogram, and machine learning classification --- is effective when applied to e… ▽ More Sensitive signal processing methods are needed to detect transiting planets from ground-based photometric surveys. Caceres et al. (2019) show that the AutoRegressive Planet Search (ARPS) method --- a combination of autoregressive integrated moving average (ARIMA) parametric modeling, a new Transit Comb Filter (TCF) periodogram, and machine learning classification --- is effective when applied to evenly spaced light curves from space-based missions. We investigate here whether ARIMA and TCF will be effective for ground-based survey light curves that are often sparsely sampled with high noise levels from atmospheric and instrumental conditions. The ARPS procedure is applied to selected light curves with strong planetary signals from the Kepler mission that have been altered to simulate the conditions of ground-based exoplanet surveys. Typical irregular cadence patterns are used from the HATSouth survey. We also evaluate recovery of known planets from HATSouth. Simulations test transit signal recovery as a function of cadence pattern and duration, stellar magnitude, planet orbital period and transit depth. Detection rates improve for shorter periods and deeper transits. The study predicts that the ARPS methodology will detect planets with $\gtrsim 0.1$\% transit depth and periods $\lesssim 40$ days in HATSouth stars brighter than $\sim$15 mag. ARPS methodology is therefore promising for planet discovery from ground-based exoplanet surveys with sufficiently dense cadence patterns. △ Less

Submitted 9 May, 2019; originally announced May 2019.

Comments: 26 pages, 7 figures, 1 machine readable table. Accepted for publication in the Astronomical Journal

Characterizing Variable Stars in a Single Night with LSST

Authors: Eric D. Feigelson, Frederica Bianco, Sara Bonito

Abstract: Stars exhibit a bewildering variety of variable behaviors ranging from explosive magnetic flares to stochastically changing accretion to periodic pulsations or rotations. The principal LSST surveys will have cadences too sparse and irregular to capture most of these phenomena. A novel idea is proposed here to observe a single Galactic field, rich in unobscured stars, in a continuous sequence of… ▽ More Stars exhibit a bewildering variety of variable behaviors ranging from explosive magnetic flares to stochastically changing accretion to periodic pulsations or rotations. The principal LSST surveys will have cadences too sparse and irregular to capture most of these phenomena. A novel idea is proposed here to observe a single Galactic field, rich in unobscured stars, in a continuous sequence of $\sim 15$ second exposures for one long winter night in a single photometric band. The result will be a unique dataset of $\sim 1$ million regularly spaced stellar lightcurves. The lightcurves will gives a particularly comprehensive collection of dM star variability. A powerful array of statistical procedures can be applied to the ensemble of lightcurves from the long-standing fields of time series analysis, signal processing and econometrics. Dozens of `features' describing the variability can be extracted and subject to machine learning classification, giving a unique authoritative objective classification of rapidly variable stars. The most effective features can then inform the wider LSST community on the best approaches to variable star identification and classification from the sparse, irregular cadences that dominate the LSST project. △ Less

Submitted 23 January, 2019; originally announced January 2019.

Comments: White Paper for LSST Call for Observing Strategies for DeepDrilling and Minisurveys

Autoregressive Times Series Methods for Time Domain Astronomy

Authors: Eric D. Feigelson, G. Jogesh Babu, Gabriel A. Caceres

Abstract: Celestial objects exhibit a wide range of variability in brightness at different wavebands. Surprisingly, the most common methods for characterizing time series in statistics -- parametric autoregressive modeling -- is rarely used to interpret astronomical light curves. We review standard ARMA, ARIMA and ARFIMA (autoregressive moving average fractionally integrated) models that treat short-memory… ▽ More Celestial objects exhibit a wide range of variability in brightness at different wavebands. Surprisingly, the most common methods for characterizing time series in statistics -- parametric autoregressive modeling -- is rarely used to interpret astronomical light curves. We review standard ARMA, ARIMA and ARFIMA (autoregressive moving average fractionally integrated) models that treat short-memory autocorrelation, long-memory $1/f^α$ `red noise', and nonstationary trends. Though designed for evenly spaced time series, moderately irregular cadences can be treated as evenly-spaced time series with missing data. Fitting algorithms are efficient and software implementations are widely available. We apply ARIMA models to light curves of four variable stars, discussing their effectiveness for different temporal characteristics. A variety of extensions to ARIMA are outlined, with emphasis on recently developed continuous-time models like CARMA and CARFIMA designed for irregularly spaced time series. Strengths and weakness of ARIMA-type modeling for astronomical data analysis and astrophysical insights are reviewed. △ Less

Submitted 23 January, 2019; originally announced January 2019.

Comments: 17 pages, 4 figures, published in 'Frontiers of Physics', vol 6, id. 80 (2018)

AutoRegressive Planet Search: Methodology

Authors: Gabriel A. Caceres, Eric D. Feigelson, G. Jogesh Babu, Natalia Bahamonde, Alejandra Christen, Karine Bertin, Cristian Meza, Michel Curé

Abstract: The detection of periodic signals from transiting exoplanets is often impeded by extraneous aperiodic photometric variability, either intrinsic to the star or arising from the measurement process. Frequently, these variations are autocorrelated wherein later flux values are correlated with previous ones. In this work, we present the methodology of the Autoregessive Planet Search (ARPS) project whi… ▽ More The detection of periodic signals from transiting exoplanets is often impeded by extraneous aperiodic photometric variability, either intrinsic to the star or arising from the measurement process. Frequently, these variations are autocorrelated wherein later flux values are correlated with previous ones. In this work, we present the methodology of the Autoregessive Planet Search (ARPS) project which uses Autoregressive Integrated Moving Average (ARIMA) and related statistical models that treat a wide variety of stochastic processes, as well as nonstationarity, to improve detection of new planetary transits. Providing a time series is evenly spaced or can be placed on an evenly spaced grid with missing values, these low-dimensional parametric models can prove very effective. We introduce a planet-search algorithm to detect periodic transits in the residuals after the application of ARIMA models. Our matched-filter algorithm, the Transit Comb Filter (TCF), is closely related to the traditional Box-fitting Least Squares and provides an analogous periodogram. Finally, if a previously identified or simulated sample of planets is available, selected scalar features from different stages of the analysis -- the original light curves, ARIMA fits, TCF periodograms, and folded light curves -- can be collectively used with a multivariate classifier to identify promising candidates while efficiently rejecting false alarms. We use Random Forests for this task, in conjunction with Receiver Operating Characteristic (ROC) curves, to define discovery criteria for new, high fidelity planetary candidates. The ARPS methodology can be applied to both evenly spaced satellite light curves and densely cadenced ground-based photometric surveys. △ Less

Submitted 14 May, 2019; v1 submitted 15 January, 2019; originally announced January 2019.

Comments: 40 pages, 12 figures, to appear in the Astronomical Journal

Kinematics in Young Star Clusters and Associations with Gaia DR2

Authors: Michael A. Kuhn, Lynne A. Hillenbrand, Alison Sills, Eric D. Feigelson, Konstantin V. Getman

Abstract: The Gaia mission has opened a new window into the internal kinematics of young star clusters at the sub-km/s level, with implications for our understanding of how star clusters form and evolve. We use a sample of 28 clusters and associations with ages from 1-5 Myr, where lists of members are available from previous X-ray, optical, and infrared studies. Proper motions from Gaia DR2 reveals that at… ▽ More The Gaia mission has opened a new window into the internal kinematics of young star clusters at the sub-km/s level, with implications for our understanding of how star clusters form and evolve. We use a sample of 28 clusters and associations with ages from 1-5 Myr, where lists of members are available from previous X-ray, optical, and infrared studies. Proper motions from Gaia DR2 reveals that at least 75% of these systems are expanding; however, rotation is only detected in one system. Typical expansion velocities are on the order of ~0.5 km/s, and, in several systems, there is a positive radial gradient in expansion velocity. Systems that are still embedded in molecular clouds are less likely to be expanding than those that are partially or fully revealed. One-dimensional velocity dispersions, which range from 1 to 3 km/s, imply that most of the stellar systems in our sample are supervirial and that some are unbound. In star-forming regions that contain multiple clusters or subclusters, we find no evidence that these groups are coalescing, implying that hierarchical cluster assembly, if it occurs, must happen rapidly during the embedded stage. △ Less

Submitted 11 November, 2018; v1 submitted 5 July, 2018; originally announced July 2018.

Comments: Accepted for publication in the Astrophysical Journal. Supplementary material is available at https://drive.google.com/drive/folders/1Cis_YIbItQ4sbItLUGz_BayOkkCVVkfM?usp=sharing. Changes: A typo in the distance to the Carina Nebula on page 5 was corrected

The APOGEE-2 Survey of the Orion Star Forming Complex II: Six-dimensional structure

Authors: Marina Kounkel, Kevin Covey, Genaro Suárez, Carlos Román-Zúñiga, Jesus Hernandez, Keivan Stassun, Karl O Jaehnig, Eric D. Feigelson, Karla Peña Ramírez, Alexandre Roman-Lopes, Nicola Da Rio, Guy S. Stringfellow, J. Serena Kim, Jura Borissova, José G. Fernández-Trincado, Adam Burgasser, D. A. García-Hernández, Olga Zamora, Kaike Pan, Christian Nitschelm

Abstract: We present an analysis of spectrosopic and astrometric data from APOGEE-2 and Gaia DR2 to identify structures towards the Orion Complex. By applying a hierarchical clustering algorithm to the 6-dimensional stellar data, we identify spatially and/or kinematically distinct groups of young stellar objects with ages ranging from 1 to 12 Myr. We also investigate the star forming history within the Orio… ▽ More We present an analysis of spectrosopic and astrometric data from APOGEE-2 and Gaia DR2 to identify structures towards the Orion Complex. By applying a hierarchical clustering algorithm to the 6-dimensional stellar data, we identify spatially and/or kinematically distinct groups of young stellar objects with ages ranging from 1 to 12 Myr. We also investigate the star forming history within the Orion Complex, and identify peculiar sub-clusters. With this method we reconstruct the older populations in the regions that are presently largely devoid of molecular gas, such as Orion C (which includes the $σ$ Ori cluster), and Orion D (the population that traces Ori OB1a, OB1b, and Orion X). We report on the distances, kinematics, and ages of the groups within the Complex. The Orion D groups is in the process of expanding. On the other hand, Orion B is still in the process of contraction. In $λ$ Ori the proper motions are consistent with a radial expansion due to an explosion from a supernova; the traceback age from the expansion exceeds the age of the youngest stars formed near the outer edges of the region, and their formation would have been triggered when they were half-way from the cluster center to their current positions. We also present a comparison between the parallax and proper motion solutions obtained by Gaia DR2, and those obtained towards star-forming regions by Very Long Baseline Array. △ Less

Submitted 25 June, 2018; v1 submitted 12 May, 2018; originally announced May 2018.

Comments: Accepted to AJ. Interactive figure 10 will be in the online version of the journal and is currently available at http://mkounkel.com/ori3d/

Intracluster Age Gradients In Numerous Young Stellar Clusters

Authors: Konstantin V. Getman, Eric D. Feigelson, Michael A. Kuhn, Matthew R. Bate, Patrick S. Broos, Gordon P. Garmire

Abstract: The pace and pattern of star formation leading to rich young stellar clusters is quite uncertain. In this context, we analyze the spatial distribution of ages within 19 young (median t<3 Myr on the Siess et al. (2000) timescale), morphologically simple, isolated, and relatively rich stellar clusters. Our analysis is based on young stellar object samples from the MYStIX and SFiNCs surveys, and a ne… ▽ More The pace and pattern of star formation leading to rich young stellar clusters is quite uncertain. In this context, we analyze the spatial distribution of ages within 19 young (median t<3 Myr on the Siess et al. (2000) timescale), morphologically simple, isolated, and relatively rich stellar clusters. Our analysis is based on young stellar object samples from the MYStIX and SFiNCs surveys, and a new estimator of pre-main sequence (PMS) stellar ages, AgeJX, derived from X-ray and near-infrared photometric data. Median cluster ages are computed within four annular subregions of the clusters. We confirm and extend the earlier result of Getman et al. (2014): 80% percent of the clusters show age trends where stars in cluster cores are younger than in outer regions. Our cluster stacking analyses establish the existence of an age gradient to high statistical significance in several ways. Time scales vary with the choice of PMS evolutionary model; the inferred median age gradient across the studied clusters ranges from 0.75 Myr/pc to 1.5 Myr/pc. The empirical finding reported in the present study -- late or continuing formation of stars in the cores of star clusters with older stars dispersed in the outer regions -- has a strong foundation with other observational studies and with the astrophysical models like the global hierarchical collapse model of Vazquez-Semadeni et al. (2017). △ Less

Submitted 13 April, 2018; originally announced April 2018.

Comments: Recently published in MNRAS; 11 pages, 5 figures, and 4 tables

Journal ref: MNRAS 476 (2018) 1213G

Circumstellar Disk Lifetimes In Numerous Galactic Young Stellar Clusters

Authors: Alexander J. W. Richert, Konstantin V. Getman, Eric D. Feigelson, Michael A. Kuhn, Patrick S. Broos, Matthew S. Povich, Matthew R. Bate, Gordon P. Garmire

Abstract: Photometric detections of dust circumstellar disks around pre-main sequence (PMS) stars, coupled with estimates of stellar ages, provide constraints on the time available for planet formation. Most previous studies on disk longevity, starting with Haisch, Lada & Lada (2001), use star samples from PMS clusters but do not consider datasets with homogeneous photometric sensitivities and/or ages place… ▽ More Photometric detections of dust circumstellar disks around pre-main sequence (PMS) stars, coupled with estimates of stellar ages, provide constraints on the time available for planet formation. Most previous studies on disk longevity, starting with Haisch, Lada & Lada (2001), use star samples from PMS clusters but do not consider datasets with homogeneous photometric sensitivities and/or ages placed on a uniform timescale. Here we conduct the largest study to date of the longevity of inner dust disks using X-ray and 1--8 micrometre infrared photometry from the MYStIX and SFiNCs projects for 69 young clusters in 32 nearby star-forming regions with ages t<=5 Myr. Cluster ages are derived by combining the empirical AgeJX method with PMS evolutionary models, which treat dynamo-generated magnetic fields in different ways. Leveraging X-ray data to identify disk-free objects, we impose similar stellar mass sensitivity limits for disk-bearing and disk-free YSOs while extending the analysis to stellar masses as low as M~0.1 Mo. We find that the disk longevity estimates are strongly affected by the choice of PMS evolutionary model. Assuming a disk fraction of 100% at zero age, the inferred disk half-life changes significantly, from t1/2 ~ 1.3--2 Myr to t1/2 ~ 3.5 Myr when switching from non-magnetic to magnetic PMS models. In addition, we find no statistically significant evidence that disk fraction varies with stellar mass within the first few Myr of life for stars with masses <2 Mo, but our samples may not be complete for more massive stars. The effects of initial disk fraction and star-forming environment are also explored. △ Less

Submitted 13 April, 2018; originally announced April 2018.

Comments: Accepted for publication in MNRAS (on April 11, 2018); 18 pages, 9 figures, and 4 tables

Young Star Clusters In Nearby Molecular Clouds

Authors: Konstantin V. Getman, Michael A. Kuhn, Eric D. Feigelson, Patrick S. Broos, Matthew R. Bate, Gordon P. Garmire

Abstract: The SFiNCs (Star Formation in Nearby Clouds) project is an X-ray/infrared study of the young stellar populations in 22 star forming regions with distances <=1 kpc designed to extend our earlier MYStIX survey of more distant clusters. Our central goal is to give empirical constraints on cluster formation mechanisms. Using parametric mixture models applied homogeneously to the catalog of SFiNCs youn… ▽ More The SFiNCs (Star Formation in Nearby Clouds) project is an X-ray/infrared study of the young stellar populations in 22 star forming regions with distances <=1 kpc designed to extend our earlier MYStIX survey of more distant clusters. Our central goal is to give empirical constraints on cluster formation mechanisms. Using parametric mixture models applied homogeneously to the catalog of SFiNCs young stars, we identify 52 SFiNCs clusters and 19 unclustered stellar structures. The procedure gives cluster properties including location, population, morphology, association to molecular clouds, absorption, age (AgeJX), and infrared spectral energy distribution (SED) slope. Absorption, SED slope, and AgeJX are age indicators. SFiNCs clusters are examined individually, and collectively with MYStIX clusters, to give the following results. (1) SFiNCs is dominated by smaller, younger, and more heavily obscured clusters than MYStIX. (2) SFiNCs cloud-associated clusters have the high ellipticities aligned with their host molecular filaments indicating morphology inherited from their parental clouds. (3) The effect of cluster expansion is evident from the radius-age, radius-absorption, and radius-SED correlations. Core radii increase dramatically from ~0.08 to ~0.9 pc over the age range 1--3.5 Myr. Inferred gas removal timescales are longer than 1 Myr. (4) Rich, spatially distributed stellar populations are present in SFiNCs clouds representing early generations of star formation. An Appendix compares the performance of the mixture models and nonparametric Minimum Spanning Tree to identify clusters. This work is a foundation for future SFiNCs/MYStIX studies including disk longevity, age gradients, and dynamical modeling. △ Less

Submitted 13 April, 2018; originally announced April 2018.

Comments: Accepted for publication in MNRAS (on February 19, 2018); 30 pages, 14 figures, and 8 tables

The Massive Star-Forming Regions Omnibus X-Ray Catalog, Second Installment

Authors: Leisa K. Townsley, Patrick S. Broos, Gordon P. Garmire, Gemma E. Anderson, Eric D. Feigelson, Tim Naylor, Matthew S. Povich

Abstract: We present the second installment of the Massive Star-forming Regions (MSFRs) Omnibus X-ray Catalog (MOXC2), a compilation of X-ray point sources detected in Chandra/ACIS observations of 16 Galactic MSFRs and surrounding fields. MOXC2 includes 13 ACIS mosaics, three containing a pair of unrelated MSFRs at different distances, with a total catalog of 18,396 point sources. The MSFRs sampled range ov… ▽ More We present the second installment of the Massive Star-forming Regions (MSFRs) Omnibus X-ray Catalog (MOXC2), a compilation of X-ray point sources detected in Chandra/ACIS observations of 16 Galactic MSFRs and surrounding fields. MOXC2 includes 13 ACIS mosaics, three containing a pair of unrelated MSFRs at different distances, with a total catalog of 18,396 point sources. The MSFRs sampled range over distances of 1.3 kpc to 6 kpc and populations varying from single massive protostars to the most massive Young Massive Cluster known in the Galaxy. By carefully detecting and removing X-ray point sources down to the faintest statistically-significant limit, we facilitate the study of the remaining unresolved X-ray emission. Through comparison with mid-infrared images that trace photon-dominated regions and ionization fronts, we see that the unresolved X-ray emission is due primarily to hot plasmas threading these MSFRs, the result of feedback from the winds and supernovae of massive stars. The 16 MSFRs studied in MOXC2 more than double the MOXC1 sample, broadening the parameter space of ACIS MSFR explorations and expanding Chandra's substantial contribution to contemporary star formation science. △ Less

Submitted 13 February, 2018; originally announced February 2018.

Comments: 60 pages, 24 figures. Accepted by ApJS. Figures have been significantly downgraded to meet arXiv requirements. See https://doi.org/10.5281/zenodo.1067748 for higher-resolution versions of this paper and for supplemental data products

Mixture Models in Astronomy

Authors: Michael A. Kuhn, Eric D. Feigelson

Abstract: Mixture models combine multiple components into a single probability density function. They are a natural statistical model for many situations in astronomy, such as surveys containing multiple types of objects, cluster analysis in various data spaces, and complicated distribution functions. This chapter in the CRC Handbook of Mixture Analysis is concerned with astronomical applications of mixture… ▽ More Mixture models combine multiple components into a single probability density function. They are a natural statistical model for many situations in astronomy, such as surveys containing multiple types of objects, cluster analysis in various data spaces, and complicated distribution functions. This chapter in the CRC Handbook of Mixture Analysis is concerned with astronomical applications of mixture models for cluster analysis, classification, and semi-parametric density estimation. We present several classification examples from the literature, including identification of a new class, analysis of contaminants, and overlapping populations. In most cases, mixtures of normal (Gaussian) distributions are used, but it is sometimes necessary to use different distribution functions derived from astrophysical experience. We also address the use of mixture models for the analysis of spatial distributions of objects, like galaxies in redshift surveys or young stars in star-forming regions. In the case of galaxy clustering, mixture models may not be the optimal choice for understanding the homogeneous and isotropic structure of voids and filaments. However, we show that mixture models, using astrophysical models for star clusters, may provide a natural solution to the problem of subdividing a young stellar population into subclusters. Finally, we explore how mixture models can be used for mathematically advanced modeling of data with heteroscedastic uncertainties or missing values, providing two example algorithms, the measurement error regression model of Kelly (2007) and the Extreme Deconvolution model of Bovy et al. (2011). The challenges presented by astronomical science, aided by the public availability of catalogs from major surveys and missions, are a rich area for collaboration between statisticians and astronomers. △ Less

Submitted 28 January, 2019; v1 submitted 29 November, 2017; originally announced November 2017.

Comments: This is a revised preprint of a chapter in the Handbook of Mixture Analysis, edited by S. Früwirth-Schnatter, G. Celeux, and C. P. Robert, published in the "Handbooks of Modern Statistical Methods" series by Chapman & Hall/CRC, 2018; revisions include several additional references and several changes to mathematical notation

The Structure of the Young Star Cluster NGC 6231. II. Structure, Formation, and Fate

Authors: Michael A. Kuhn, Konstantin V. Getman, Eric D. Feigelson, Alison Sills, Mariusz Gromadzki, Nicolás Medina, Jordanka Borissova, Radostin Kurtev

Abstract: The young cluster NGC 6231 (stellar ages ~2-7 Myr) is observed shortly after star-formation activity has ceased. Using the catalog of 2148 probable cluster members obtained from Chandra, VVV, and optical surveys (Paper I), we examine the cluster's spatial structure and dynamical state. The spatial distribution of stars is remarkably well fit by an isothermal sphere with moderate elongation, while… ▽ More The young cluster NGC 6231 (stellar ages ~2-7 Myr) is observed shortly after star-formation activity has ceased. Using the catalog of 2148 probable cluster members obtained from Chandra, VVV, and optical surveys (Paper I), we examine the cluster's spatial structure and dynamical state. The spatial distribution of stars is remarkably well fit by an isothermal sphere with moderate elongation, while other commonly used models like Plummer spheres, multivariate normal distributions, or power-law models are poor fits. The cluster has a core radius of $1.2\pm0.1$ pc and a central density of ~200 stars pc$^{-3}$. The distribution of stars is mildly mass segregated. However, there is no radial stratification of the stars by age. Although most of the stars belong to a single cluster, a small subcluster of stars is found superimposed on the main cluster, and there are clumpy non-isotropic distributions of stars outside ~4 core radii. When the size, mass, and age of NGC 6231 are compared to other young star clusters and subclusters in nearby active star-forming regions, it lies at the high-mass end of the distribution but along the same trend line. This could result from similar formation processes, possibly hierarchical cluster assembly. We argue that NGC 6231 has expanded from its initial size but that it remains gravitationally bound. △ Less

Submitted 4 October, 2017; originally announced October 2017.

Comments: Accepted for publication in The Astronomical Journal; 27 pages, 2 tables, and 14 figures

The Structure of the Young Star Cluster NGC 6231. I. Stellar Population

Authors: Michael A. Kuhn, Nicolás Medina, Konstantin V. Getman, Eric D. Feigelson, Mariusz Gromadzki, Jordanka Borissova, Radostin Kurtev

Abstract: NGC 6231 is a young cluster (age ~2-7 Myr) dominating the Sco OB1 association (distance ~1.59 kpc) with ~100 O and B stars and a large pre-main-sequence stellar population. We combine a reanalysis of archival Chandra X-ray data with multi-epoch NIR photometry from the VVV survey and published optical catalogs to obtain a catalog of 2148 probable cluster members. This catalog is 70% larger than pre… ▽ More NGC 6231 is a young cluster (age ~2-7 Myr) dominating the Sco OB1 association (distance ~1.59 kpc) with ~100 O and B stars and a large pre-main-sequence stellar population. We combine a reanalysis of archival Chandra X-ray data with multi-epoch NIR photometry from the VVV survey and published optical catalogs to obtain a catalog of 2148 probable cluster members. This catalog is 70% larger than previous censuses of probable cluster members in NGC 6231, and it includes many low-mass stars detected in the NIR but not in the optical and some B-stars without previously noted X-ray counterparts. In addition, we identify 295 NIR variables, about half of which are expected to be pre-main-sequence stars. With the more-complete sample, we estimate a total population in the Chandra field of 5700-7500 cluster members down to 0.08 $M_\odot$ (assuming a universal initial mass function) with a completeness limit at 0.5 $M_\odot$. A decrease in stellar X-ray luminosities is noted relative to other younger clusters. However, within the cluster, there is little variation in the distribution of X-ray luminosities for ages less than 5 Myr. X-ray spectral hardness for B stars may be useful for distinguishing between early-B stars with X-rays generated in stellar winds and B-star systems with X-rays from a pre-main-sequence companions (>35% of B stars). A small fraction of pre-main-sequence stars have unusually high X-ray median energies or reddened near-infrared colors, which might be explained by absorption from thick or edge-on disks or being a background field star. △ Less

Submitted 31 May, 2017; originally announced June 2017.

Comments: Accepted for publication in The Astronomical Journal; 40 pages, 5 tables, and 17 figures; full versions of Tables 1, 3, and 4 are available at http://astro.psu.edu/mystix

Multiwavelength Studies of Young OB Associations

Authors: Eric D. Feigelson

Abstract: We discuss how contemporary multiwavelength observations of young OB-dominated clusters address long-standing astrophysical questions: Do clusters form rapidly or slowly with an age spread? When do clusters expand and disperse to constitute the field star population? Do rich clusters form by amalgamation of smaller subclusters? What is the pattern and duration of cluster formation in massive star… ▽ More We discuss how contemporary multiwavelength observations of young OB-dominated clusters address long-standing astrophysical questions: Do clusters form rapidly or slowly with an age spread? When do clusters expand and disperse to constitute the field star population? Do rich clusters form by amalgamation of smaller subclusters? What is the pattern and duration of cluster formation in massive star forming regions (MSFRs)? Past observational difficulties in obtaining good stellar censuses of MSFRs have been alleviated in recent studies that combine X-ray and infrared surveys to obtain rich, though still incomplete, censuses of young stars in MSFRs. We describe here one of these efforts, the MYStIX project, that produced a catalog of 31,784 probable members of 20 MSFRs. We find that age spread within clusters are real in the sense that the stars in the core formed after the cluster halo. Cluster expansion is seen in the ensemble of (sub)clusters, and older dispersing populations are found across MSFRs. Direct evidence for subcluster merging is still unconvincing. Long-lived, asynchronous star formation is pervasive across MSFRs. △ Less

Submitted 26 April, 2017; originally announced April 2017.

Comments: 22 pages, 9 figures. To appear in "The Origin of Stellar Clusters", edited by Steven Stahler, Springer, 2017, in press

A probabilistic approach to emission-line galaxy classification

Authors: R. S. de Souza, M. L. L. Dantas, M. V. Costa-Duarte, E. D. Feigelson, M. Killedar, P. -Y. Lablanche, R. Vilalta, A. Krone-Martins, R. Beck, F. Gieseke

Abstract: We invoke a Gaussian mixture model (GMM) to jointly analyse two traditional emission-line classification schemes of galaxy ionization sources: the Baldwin-Phillips-Terlevich (BPT) and $\rm W_{Hα}$ vs. [NII]/H$α$ (WHAN) diagrams, using spectroscopic data from the Sloan Digital Sky Survey Data Release 7 and SEAGal/STARLIGHT datasets. We apply a GMM to empirically define classes of galaxies in a thre… ▽ More We invoke a Gaussian mixture model (GMM) to jointly analyse two traditional emission-line classification schemes of galaxy ionization sources: the Baldwin-Phillips-Terlevich (BPT) and $\rm W_{Hα}$ vs. [NII]/H$α$ (WHAN) diagrams, using spectroscopic data from the Sloan Digital Sky Survey Data Release 7 and SEAGal/STARLIGHT datasets. We apply a GMM to empirically define classes of galaxies in a three-dimensional space spanned by the $\log$ [OIII]/H$β$, $\log$ [NII]/H$α$, and $\log$ EW(H$α$), optical parameters. The best-fit GMM based on several statistical criteria suggests a solution around four Gaussian components (GCs), which are capable to explain up to 97 per cent of the data variance. Using elements of information theory, we compare each GC to their respective astronomical counterpart. GC1 and GC4 are associated with star-forming galaxies, suggesting the need to define a new starburst subgroup. GC2 is associated with BPT's Active Galaxy Nuclei (AGN) class and WHAN's weak AGN class. GC3 is associated with BPT's composite class and WHAN's strong AGN class. Conversely, there is no statistical evidence -- based on four GCs -- for the existence of a Seyfert/LINER dichotomy in our sample. Notwithstanding, the inclusion of an additional GC5 unravels it. The GC5 appears associated to the LINER and Passive galaxies on the BPT and WHAN diagrams respectively. Subtleties aside, we demonstrate the potential of our methodology to recover/unravel different objects inside the wilderness of astronomical datasets, without lacking the ability to convey physically interpretable results. The probabilistic classifications from the GMM analysis are publicly available within the COINtoolbox (https://cointoolbox.github.io/GMM\_Catalogue/). △ Less

Submitted 18 August, 2017; v1 submitted 22 March, 2017; originally announced March 2017.

Comments: Accepted for publication in MNRAS

Candidate X-ray-Emitting OB Stars in the MYStIX Massive Star-Forming Regions

Authors: Matthew S. Povich, Heather A. Busk, Eric D. Feigelson, Leisa K. Townsley, Michael A. Kuhn

Abstract: Massive, O and early B-type (OB) stars remain incompletely catalogued in the nearby Galaxy due to high extinction, bright visible and infrared nebular emission in H II regions, and high field star contamination. These difficulties are alleviated by restricting the search to stars with X-ray emission. Using the X-ray point sources from the Massive Young star-forming complex Study in Infrared and X-… ▽ More Massive, O and early B-type (OB) stars remain incompletely catalogued in the nearby Galaxy due to high extinction, bright visible and infrared nebular emission in H II regions, and high field star contamination. These difficulties are alleviated by restricting the search to stars with X-ray emission. Using the X-ray point sources from the Massive Young star-forming complex Study in Infrared and X-rays (MYStIX) survey of OB-dominated regions, we identify 98 MYStIX candidate OB (MOBc) stars by fitting their 1-8 micron spectral energy distributions (SEDs) with reddened stellar atmosphere models. We identify 27 additional MOBc stars based on JHK photometry of X-ray stars lacking SED fitting. These candidate OB stars indicate that the current census of stars earlier than B1, taken across the 18 MYStIX regions studied, is less than 50% complete. We also fit the SEDs of 239 previously-published OB stars to measure interstellar extinction and bolometric luminosities, revealing six candidate massive binary systems and five candidate O-type (super)giants. As expected, candidate OB stars have systematically higher extinction than previously-published OB stars. Notable results for individual regions include: identification of the OB population of a recently discovered massive cluster in NGC 6357; an older OB association in the M17 complex; and new massive luminous O stars near the Trifid Nebula. In several relatively poorly-studied regions (RCW 38, NGC 6334, NGC 6357, Trifid, and NGC 3576), the OB populations may increase by factors of >2. △ Less

Submitted 23 January, 2017; originally announced January 2017.

Comments: 49 pages, 4 figure sets, 5 tables, accepted for publication in ApJ

The changing landscape of astrostatistics and astroinformatics

Authors: Eric D. Feigelson

Abstract: The history and current status of the cross-disciplinary fields of astrostatistics and astroinformatics are reviewed. Astronomers need a wide range of statistical methods for both data reduction and science analysis. With the proliferation of high-throughput telescopes, efficient large scale computational methods are also becoming essential. However, astronomers receive only weak training in these… ▽ More The history and current status of the cross-disciplinary fields of astrostatistics and astroinformatics are reviewed. Astronomers need a wide range of statistical methods for both data reduction and science analysis. With the proliferation of high-throughput telescopes, efficient large scale computational methods are also becoming essential. However, astronomers receive only weak training in these fields during their formal education. Interest in the fields is rapidly growing with conferences organized by scholarly societies, textbooks and tutorial workshops, and research studies pushing the frontiers of methodology. R, the premier language of statistical computing, can provide an important software environment for the incorporation of advanced statistical and computational methodology into the astronomical community. △ Less

Submitted 19 December, 2016; originally announced December 2016.

Comments: 7 pages, to appear in 'Astroinformatics', IAU Symposium #325, M. Brescia et al. (eds.), Cambridge University Press (2017)

Star Formation In Nearby Clouds (SFiNCs): X-ray And Infrared Source Catalogs And Membership

Authors: Konstantin V. Getman, Patrick S. Broos, Michael A. Kuhn, Eric D. Feigelson, Alexander J. W. Richert, Yosuke Ota, Matthew R. Bate, Gordon P. Garmire

Abstract: The Star Formation in Nearby Clouds (SFiNCs) project is aimed at providing detailed study of the young stellar populations and star cluster formation in nearby 22 star forming regions (SFRs) for comparison with our earlier MYStIX survey of richer, more distant clusters. As a foundation for the SFiNCs science studies, here, homogeneous data analyses of the Chandra X-ray and Spitzer mid-infrared arc… ▽ More The Star Formation in Nearby Clouds (SFiNCs) project is aimed at providing detailed study of the young stellar populations and star cluster formation in nearby 22 star forming regions (SFRs) for comparison with our earlier MYStIX survey of richer, more distant clusters. As a foundation for the SFiNCs science studies, here, homogeneous data analyses of the Chandra X-ray and Spitzer mid-infrared archival SFiNCs data are described, and the resulting catalogs of over 15300 X-ray and over 1630000 mid-infrared point sources are presented. On the basis of their X-ray/infrared properties and spatial distributions, nearly 8500 point sources have been identified as probable young stellar members of the SFiNCs regions. Compared to the existing X-ray/mid-infrared publications, the SFiNCs member list increases the census of YSO members by 6-200% for individual SFRs and by 40% for the merged sample of all 22 SFiNCs SFRs. △ Less

Submitted 15 December, 2016; originally announced December 2016.

Comments: Accepted for publication in ApJS; 21 text pages, 11 tables, 22 figures (including 12 figure sets and SPCM Atlas). Full versions of all figure sets and tables can be found on-line at: https://drive.google.com/drive/folders/0B4lwbriAuoXfVDNHckRnLV92Rzg?usp=sharing . The 5GB SPCM Atlas can be found on-line at: https://drive.google.com/drive/folders/0B4lwbriAuoXfSXJtbzB6V2tyNTA?usp=sharing

Young Stellar Populations in MYStIX Star Forming Regions: Candidate Protostars

Authors: Gregory Romine, Eric D. Feigelson, Konstantin V. Getman, Michael A. Kuhn, Matthew S. Povich

Abstract: 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 contami… ▽ More 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. △ Less

Submitted 26 September, 2016; v1 submitted 21 September, 2016; originally announced September 2016.

Comments: 35 pages, 1 figure, 3 tables. Accepted for publication in the American Astronomical Society Journals. A full version with figure set and electronic tables is available (with other MYStIX papers) at http://www.astro.psu.edu/mystix

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

Authors: Alexander J. W. Richert, Eric D. Feigelson, Konstantin V. Getman, Michael A. Kuhn

Abstract: 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 protoplanet… ▽ More 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. △ Less

Submitted 6 August, 2015; originally announced August 2015.

Comments: 24 pages, 32 figures. ApJ, accepted

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

Authors: Michael A. Kuhn, Eric D. Feigelson, Konstantin V. Getman, Alison Sills, Matthew R. Bate, Jordanka Borissova

Abstract: 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, an… ▽ More 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. △ Less

Submitted 20 July, 2015; originally announced July 2015.

Comments: Accepted for publication in The Astrophysical Journal ; 48 pages, 13 figures, and 6 tables

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

Authors: Michael A. Kuhn, Konstantin V. Getman, Eric D. Feigelson

Abstract: 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 eff… ▽ More 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. △ Less

Submitted 21 January, 2015; originally announced January 2015.

Comments: Accepted for publication in The Astrophysical Journal; 45 pages, 9 figures, and 3 tables. Supplemental materials from the online-version of this article, including a machine-readable table and "Data behind the Figure" (*.fits), are available at the public MYStIX website (http://astro.psu.edu/mystix/)

The Spatial Structure of Young Stellar Clusters. I. Subclusters

Authors: Michael A. Kuhn, Eric D. Feigelson, Konstantin V. Getman, Adrian J. Baddeley, Patrick S. Broos, Alison Sills, Matthew R. Bate, Matthew S. Povich, Kevin L. Luhman, Heather A. Busk, Tim Naylor, Robert R. King

Abstract: 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… ▽ More 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. △ Less

Submitted 17 March, 2014; originally announced March 2014.

Comments: Accepted for publication in The Astrophysical Journal; 67 pages, 9 figures, and 4 tables. Supplemental materials from the online-version of this article, including machine-readable tables and source code in R are available at the public MYStIX website (http://astro.psu.edu/mystix/)

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

Authors: Konstantin V. Getman, Eric D. Feigelson, Michael A. Kuhn

Abstract: 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 l… ▽ More 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. △ Less

Submitted 11 March, 2014; originally announced March 2014.

Comments: Accepted for publication in ApJ; 42 pages, 5 figures, and 1 table; High quality version is at http://astro.psu.edu/mystix

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

Authors: Konstantin V. Getman, Eric D. Feigelson, Michael A. Kuhn, Patrick S. Broos, Leisa K. Townsley, Tim Naylor, Matthew S. Povich, Kevin L. Luhman, Gordon P. Garmire

Abstract: 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… ▽ More 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. △ Less

Submitted 11 March, 2014; originally announced March 2014.

Comments: Accepted for publication in ApJ; 89 pages, 23 figures, 2 Tables; High quality version is at http://astro.psu.edu/mystix

The Massive Star-forming Regions Omnibus X-ray Catalog

Authors: Leisa K. Townsley, Patrick S. Broos, Gordon P. Garmire, Jeroen Bouwman, Matthew S. Povich, Eric D. Feigelson, Konstantin V. Getman, Michael A. Kuhn

Abstract: 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 unreso… ▽ More 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. △ Less

Submitted 11 March, 2014; originally announced March 2014.

Comments: Accepted to ApJS, March 3, 2014. 51 pages, 25 figures

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

Authors: Patrick S. Broos, Konstantin V. Getman, Matthew S. Povich, Eric D. Feigelson, Leisa K. Townsley, Tim Naylor, Michael A. Kuhn, R. R. King, Heather A. Busk

Abstract: 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 publ… ▽ More 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 △ Less

Submitted 17 September, 2013; originally announced September 2013.

Comments: 47 pages, 7 figures. Accepted for publication of Astrophysical Journal Supplements (with 6 other MYStIX papers). Complete versions of all MYStIX papers are available at http://www.astro.psu.edu/mystix