sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Kuhn, Michael A.

Normalized to: Kuhn, M.

46 article(s) in total. 279 co-authors, from 1 to 28 common article(s). Median position in authors list is 3,0.

[1]  [pdf] - 2127879
The Formation of a Stellar Association in the NGC 7000/IC 5070 Complex: Results from Kinematic Analysis of Stars and Gas
Comments: ApJ (in press). 22 figures and 7 tables. Full versions of Tables 2 and 7 and the extra panels for Figures 12, 16, and 17 are available at
Submitted: 2020-06-15, last modified: 2020-07-06
We examine the clustering and kinematics of young stellar objects (YSOs) in the North America/Pelican Nebulae, as revealed by Gaia astrometry, in relation to the structure and motions of the molecular gas, as indicated in molecular line maps. The Gaia parallaxes and proper motions allow us to significantly refine previously published lists of YSOs, demonstrating that many of the objects previously thought to form a distributed population turn out to be non-members. The members are subdivided into at least 6 spatio-kinematic groups, each of which is associated with its own molecular cloud component or components. Three of the groups are expanding, with velocity gradients of 0.3-0.5 km/s/pc, up to maximum velocities of ~8 km/s away from the groups' centers. The two known O-type stars associated with the region, 2MASS J20555125+4352246 and HD 199579, are rapidly escaping one of these groups, following the same position-velocity relation as the low-mass stars. We calculate that a combination of gas expulsion and tidal forces from the clumpy distribution of molecular gas could impart the observed velocity gradients within the groups. However, on a global scale, the relative motions of the groups do not appear either divergent or convergent. The velocity dispersion of the whole system is consistent with the kinetic energy gained due to gravitational collapse of the complex. Most of the stellar population has ages similar to the free-fall timescales for the natal clouds. Thus, we suggest the nearly free-fall collapse of a turbulent molecular cloud as the most likely scenario for star formation in this complex.
[2]  [pdf] - 2025985
IRAS 09002-4732: A Laboratory for the Formation of Rich Stellar Clusters
Comments: 33 pages, 11 figures, 3 machine readable tables. To appear in The Astronomical Journal
Submitted: 2019-10-10, last modified: 2019-10-15
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.
[3]  [pdf] - 1971291
A Comparison of the X-ray Properties of FU Ori-type Stars to Generic Young Stellar Objects
Comments: 17 pages, 6 figures, 3 tables. Accepted for publication in ApJ
Submitted: 2019-07-07, last modified: 2019-08-12
Like other young stellar objects (YSOs), FU Ori-type stars have been detected as strong X-ray emitters. However, little is known about how the outbursts of these stars affect their X-ray properties. We assemble available X-ray data from XMM Newton and Chandra observations of 16 FU Ori stars, including a new XMM Newton observation of Gaia 17bpi during its optical rise phase. Of these stars, six were detected at least once, while 10 were non-detections, for which we calculate upper limits on intrinsic X-ray luminosity ($L_X$) as a function of plasma temperature ($kT$) and column density ($N_H$). The detected FU Ori stars tend to be more X-ray luminous than typical for non-outbursting YSOs, based on comparison to a sample of low-mass stars in the Orion Nebula Cluster. FU Ori stars with high $L_X$ have been observed both at the onset of their outbursts and decades later. We use the Kaplan-Meier estimator to investigate whether the higher X-ray luminosities for FU Ori stars is characteristic or a result of selection effects, and we find the difference to be statistically significant ($p<0.01$) even when non-detections are taken into account. The additional X-ray luminosity of FU Ori stars relative to non-outbursting YSOs cannot be explained by accretion shocks, given the high observed plasma temperatures. This suggests that, for many FU Ori stars, either 1) the outburst leads to a restructuring of the magnetosphere in a way that enhances X-ray emission, or 2) FU Ori outbursts are more likely to occur among YSOs with the highest quiescent X-ray luminosity.
[4]  [pdf] - 1894388
Gaia Stellar Kinematics in the Head of the Orion A Cloud: Runaway Stellar Groups and Gravitational Infall
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: (`restricted' sample) and (`full' sample)
Submitted: 2019-05-24
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.
[5]  [pdf] - 1905743
An initial overview of the extent and structure of recent star formation within the Serpens Molecular Cloud using Gaia Data Release 2
Comments: 25 pages, accepted by ApJ
Submitted: 2019-04-08
The dense clusters within the Serpens Molecular Cloud are among the most active regions of nearby star formation. In this paper, we use Gaia DR2 parallaxes and proper motions to statistically measure 1167 kinematic members of Serpens, few of which were previously identified, to evaluate the star formation history of the complex. The optical members of Serpens are concentrated in three distinct groups located at 380 to 480 pc; the densest clusters are still highly obscured by optically-thick dust and have few optical members. The total population of young stars and protostars in Serpens is at least 2000 stars, including past surveys that were most sensitive to protostars and disks, and may be far higher. Distances to dark clouds measured from deficits in star counts are consistent with the distances to the optical star clusters. The Serpens Molecular Cloud is seen in the foreground of the Aquila Rift, dark clouds located at 600 to 700 pc, and behind patchy extinction, here called the Serpens Cirrus, located at ~250 pc. Based on the lack of a distributed population of older stars, the star formation rate throughout the Serpens Molecular Cloud increased by at least a factor of 20 within the past ~5 Myr. The optically bright stars in Serpens Northeast are visible because their natal molecular cloud has been eroded and not because they were flung outwards from a central factory of star formation. The separation between subclusters of 20 to 100 pc and the absence of an older population leads to speculation that an external forcing was needed to trigger the active star formation.
[6]  [pdf] - 1850859
Astro2020 Science White Paper: The Next Decade of Astroinformatics and Astrostatistics
Comments: Submitted to the Astro2020 Decadal Survey call for science white papers
Submitted: 2019-03-15
Over the past century, major advances in astronomy and astrophysics have been largely driven by improvements in instrumentation and data collection. With the amassing of high quality data from new telescopes, and especially with the advent of deep and large astronomical surveys, it is becoming clear that future advances will also rely heavily on how those data are analyzed and interpreted. New methodologies derived from advances in statistics, computer science, and machine learning are beginning to be employed in sophisticated investigations that are not only bringing forth new discoveries, but are placing them on a solid footing. Progress in wide-field sky surveys, interferometric imaging, precision cosmology, exoplanet detection and characterization, and many subfields of stellar, Galactic and extragalactic astronomy, has resulted in complex data analysis challenges that must be solved to perform scientific inference. Research in astrostatistics and astroinformatics will be necessary to develop the state-of-the-art methodology needed in astronomy. Overcoming these challenges requires dedicated, interdisciplinary research. We recommend: (1) increasing funding for interdisciplinary projects in astrostatistics and astroinformatics; (2) dedicating space and time at conferences for interdisciplinary research and promotion; (3) developing sustainable funding for long-term astrostatisics appointments; and (4) funding infrastructure development for data archives and archive support, state-of-the-art algorithms, and efficient computing.
[7]  [pdf] - 1890352
The Zwicky Transient Facility: Science Objectives
Graham, Matthew J.; Kulkarni, S. R.; Bellm, Eric C.; Adams, Scott M.; Barbarino, Cristina; Blagorodnova, Nadejda; Bodewits, Dennis; Bolin, Bryce; Brady, Patrick R.; Cenko, S. Bradley; Chang, Chan-Kao; Coughlin, Michael W.; De, Kishalay; Eadie, Gwendolyn; Farnham, Tony L.; Feindt, Ulrich; Franckowiak, Anna; Fremling, Christoffer; Gal-yam, Avishay; Gezari, Suvi; Ghosh, Shaon; Goldstein, Daniel A.; Golkhou, V. Zach; Goobar, Ariel; Ho, Anna Y. Q.; Huppenkothen, Daniela; Ivezic, Zeljko; Jones, R. Lynne; Juric, Mario; Kaplan, David L.; Kasliwal, Mansi M.; Kelley, Michael S. P.; Kupfer, Thomas; Lee, Chien-De; Lin, Hsing Wen; Lunnan, Ragnhild; Mahabal, Ashish A.; Miller, Adam A.; Ngeow, Chow-Choong; Nugent, Peter; Ofek, Eran O.; Prince, Thomas A.; Rauch, Ludwig; van Roestel, Jan; Schulze, Steve; Singer, Leo P.; Sollerman, Jesper; Taddia, Francesco; Yan, Lin; Ye, Quan-Zhi; Yu, Po-Chieh; Andreoni, Igor; Barlow, Tom; Bauer, James; Beck, Ron; Belicki, Justin; Biswas, Rahul; Brinnel, Valery; Brooke, Tim; Bue, Brian; Bulla, Mattia; Burdge, Kevin; Burruss, Rick; Connolly, Andrew; Cromer, John; Cunningham, Virginia; Dekany, Richard; Delacroix, Alex; Desai, Vandana; Duev, Dmitry A.; Hacopians, Eugean; Hale, David; Helou, George; Henning, John; Hover, David; Hillenbrand, Lynne A.; Howell, Justin; Hung, Tiara; Imel, David; Ip, Wing-Huen; Jackson, Edward; Kaspi, Shai; Kaye, Stephen; Kowalski, Marek; Kramer, Emily; Kuhn, Michael; Landry, Walter; Laher, Russ R.; Mao, Peter; Masci, Frank J.; Monkewitz, Serge; Murphy, Patrick; Nordin, Jakob; Patterson, Maria T.; Penprase, Bryan; Porter, Michael; Rebbapragada, Umaa; Reiley, Dan; Riddle, Reed; Rigault, Mickael; Rodriguez, Hector; Rusholme, Ben; van Santen, Jakob; Shupe, David L.; Smith, Roger M.; Soumagnac, Maayane T.; Stein, Robert; Surace, Jason; Szkody, Paula; Terek, Scott; van Sistine, Angela; van Velzen, Sjoert; Vestrand, W. Thomas; Walters, Richard; Ward, Charlotte; Zhang, Chaoran; Zolkower, Jeffry
Comments: 26 pages, 7 figures, Published in PASP Focus Issue on the Zwicky Transient Facility
Submitted: 2019-02-05
The Zwicky Transient Facility (ZTF), a public-private enterprise, is a new time domain survey employing a dedicated camera on the Palomar 48-inch Schmidt telescope with a 47 deg$^2$ field of view and 8 second readout time. It is well positioned in the development of time domain astronomy, offering operations at 10% of the scale and style of the Large Synoptic Survey Telescope (LSST) with a single 1-m class survey telescope. The public surveys will cover the observable northern sky every three nights in g and r filters and the visible Galactic plane every night in g and r. Alerts generated by these surveys are sent in real time to brokers. A consortium of universities which provided funding ("partnership") are undertaking several boutique surveys. The combination of these surveys producing one million alerts per night allows for exploration of transient and variable astrophysical phenomena brighter than r $\sim$ 20.5 on timescales of minutes to years. We describe the primary science objectives driving ZTF including the physics of supernovae and relativistic explosions, multi-messenger astrophysics, supernova cosmology, active galactic nuclei and tidal disruption events, stellar variability, and Solar System objects.
[8]  [pdf] - 1822733
Mixture Models in Astronomy
Comments: This is a revised preprint of a chapter in the Handbook of Mixture Analysis, edited by S. Fr\"uwirth-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
Submitted: 2017-11-29, last modified: 2019-01-28
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.
[9]  [pdf] - 1806290
Gaia 17bpi: An FU Ori Type Outburst
Comments: Accepted to ApJ
Submitted: 2018-12-17
We report on the source Gaia 17bpi and identify it as a new, ongoing FU Ori type outburst, associated with a young stellar object. The optical lightcurve from Gaia exhibited a 3.5 mag rise with the source appearing to plateau in mid/late 2018. Mid-infrared observations from NEOWISE also show a $>$3 mag rise that occurred in two stages, with the second one coincident with the optical brightening, and the first one preceding the optical brightening by $\sim$1.5 years. We model the outburst as having started between October and December of 2014. This wavelength-dependent aspect of young star accretion-driven outbursts has never been documented before. Both the mid-infrared and the optical colors of the object become bluer as the outburst proceeds. Optical spectroscopic characteristics in the outburst phase include: a GK-type absorption spectrum, strong wind/outflow in e.g. Mgb, NaD, H$\alpha$, KI, OI, and CaII profiles, and detection of LiI 6707 \AA. The infrared spectrum in the outburst phase is similar to that of an M-type spectrum, notably exhibiting prominent $H_2O$ and $^{12}$CO (2-0) bandhead absorption in the K-band, and likely HeI wind in the Y-band. The new FU Ori source Gaia 17bpi is associated with a little-studied dark cloud in the galactic plane, located at a distance of 1.27 kpc.
[10]  [pdf] - 1811070
Kinematics in Young Star Clusters and Associations with Gaia DR2
Comments: Accepted for publication in the Astrophysical Journal. Supplementary material is available at Changes: A typo in the distance to the Carina Nebula on page 5 was corrected
Submitted: 2018-07-05, last modified: 2018-11-11
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.
[11]  [pdf] - 1779765
New Galactic Star Clusters Discovered in the Disk Area of the VVVX Survey
Comments: 19 pages, 10 figures, published in the MNRAS
Submitted: 2018-11-05
The "VISTA Variables in the Via Lactea eXtended (VVVX)" ESO Public Survey is a near-infrared photometric sky survey that covers nearly 1700 sq. deg towards the Galactic disk and bulge. It is well-suited to search for new open clusters, hidden behind dust and gas. The pipeline processed and calibrated Ks-band tile images of 40% of the disk area covered by VVVX was visually inspected for stellar over-densities. Then, we identified cluster candidates by examination of the composite JHKs color images. The color-magnitude diagrams of the cluster candidates are constructed. Whenever possible the Gaia DR2 parameters are used to calculate the mean proper motions, radial velocities, reddening and distances. We report the discovery of 120 new infrared clusters and stellar groups. Approximately, half of them (47%) are faint, compact, highly reddened, and they seem to be associated with other indicators of recent star formation, such as nearby Young Stellar Objects, Masers, H II regions or bubbles. The preliminary distance determinations allow us to trace the clusters up to 4.5 kpc, but most of the cluster candidates are centered at 2.2 kpc. The mean proper motions of the clusters, show that in general, they follow the disk motion of the Galaxy.
[12]  [pdf] - 1743752
An Automated tool to detect variable sources in the Vista Variables in the V\'{\i}a L\'actea Survey. The VVV Variables (V$^{4}$) catalog of tiles d001 and d002
Comments: 27 pages, 19 figures
Submitted: 2018-06-11
Time-varying phenomena are one of the most substantial sources of astrophysical information and their study has led to many fundamental discoveries in modern astronomy. We have developed an automated tool to search and analyze variable sources in the near infrared $\rm K_{s}$-band, using the data from the Vista Variables in the V\'ia L\'actea (VVV) ESO Public Large Survey. This process relies on the characterization of variable sources using different variability indices, calculated from time series generated with Point Spread Function photometry of sources under analysis. In particular, we used two main indices: the total amplitude $\rm \Delta K_s$ and the eta index, $\eta$, to identify variable sources. Once variable objects are identified, periods are determined with Generalized Lomb-Scargle periodograms, and the Information Potential Metric. Variability classes are assigned according to a compromise between comparisons with VVV Templates and the period of the variability. The automated tool is applied on VVV tiles d001 and d002 and led to discovery of 200 variable sources. We detected 70 irregular variable sources and 130 periodic ones. In addition nine open cluster candidates projected in the region are analyzed, the infrared variable candidates found around these clusters are further scrutinized by cross-matching their locations against emission star candidates from VPHAS+ survey $\rm H_{\alpha}$ color cuts.
[13]  [pdf] - 1666208
Intracluster Age Gradients In Numerous Young Stellar Clusters
Comments: Recently published in MNRAS; 11 pages, 5 figures, and 4 tables
Submitted: 2018-04-13
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).
[14]  [pdf] - 1675030
Young Star Clusters In Nearby Molecular Clouds
Comments: Accepted for publication in MNRAS (on February 19, 2018); 30 pages, 14 figures, and 8 tables
Submitted: 2018-04-13
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.
[15]  [pdf] - 1675031
Circumstellar Disk Lifetimes In Numerous Galactic Young Stellar Clusters
Comments: Accepted for publication in MNRAS (on April 11, 2018); 18 pages, 9 figures, and 4 tables
Submitted: 2018-04-13
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.
[16]  [pdf] - 1608533
VIRAC: The VVV Infrared Astrometric Catalogue
Comments: Accepted for publication in MNRAS. See
Submitted: 2017-10-24
We present VIRAC version 1, a near-infrared proper motion and parallax catalogue of the VISTA VVV survey for 312,587,642 unique sources averaged across all overlapping pawprint and tile images covering 560 deg$^2$ of the bulge of the Milky Way and southern disk. The catalogue includes 119 million high quality proper motion measurements, of which 47 million have statistical uncertainties below 1 mas yr$^{-1}$. In the 11$<K_s<$14 magnitude range the high quality motions have a median uncertainty of 0.67 mas yr$^{-1}$. The catalogue also includes 6,935 sources with quality-controlled 5 $\sigma$ parallaxes with a median uncertainty of 1.1 mas. The parallaxes show reasonable agreement with the TYCHO-Gaia Astrometric Solution (TGAS), though caution is advised for data with modest significance. The SQL database housing the data is made available via the web. We give example applications for studies of Galactic structure, nearby objects (low mass stars and brown dwarfs, subdwarfs, white dwarfs) and kinematic distance measurements of YSOs. Nearby objects discovered include LTT 7251 B, an L7 benchmark companion to a G dwarf with over 20 published elemental abundances, a bright L sub-dwarf, VVV 1256-6202, with extremely blue colours and nine new members of the 25 pc sample. We also demonstrate why this catalogue remains useful in the era of Gaia. Future versions will be based on profile fitting photometry, use the Gaia absolute reference frame and incorporate the longer time baseline of the VVV extended survey (VVVX).
[17]  [pdf] - 1589294
The Structure of the Young Star Cluster NGC 6231. II. Structure, Formation, and Fate
Comments: Accepted for publication in The Astronomical Journal; 27 pages, 2 tables, and 14 figures
Submitted: 2017-10-04
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.
[18]  [pdf] - 1586512
C$_{60}^+$ and the Diffuse Interstellar Bands: An Independent Laboratory Check
Submitted: 2017-07-28
In 2015, Campbell et al. (Nature 523, 322) presented spectroscopic laboratory gas phase data for the fullerene cation, C$_{60}^+$, that coincide with reported astronomical spectra of two diffuse interstellar band (DIB) features at 9633 and 9578 \AA. In the following year additional laboratory spectra were linked to three other and weaker DIBs at 9428, 9366, and 9349 \AA. The laboratory data were obtained using wavelength-dependent photodissociation spectroscopy of small (up to three) He-tagged C$_{60}^+-$He$_n$ ion complexes, yielding rest wavelengths for the bare C$_{60}^+$ cation by correcting for the He-induced wavelength shifts. Here we present an alternative approach to derive the rest wavelengths of the four most prominent C$_{60}^+$ absorption features, using high resolution laser dissociation spectroscopy of C$_{60}^+$ embedded in ultracold He droplets. Accurate wavelengths of the bare fullerene cation are derived based on linear wavelength shifts recorded for He$_n$C$_{60}^+$ species with $n$ up to 32. A careful analysis of all available data results in precise rest wavelengths (in air) for the four most prominent C$_{60}^+$ bands: 9631.9(1) \AA, 9576.7(1) \AA, 9427.5(1) \AA, and 9364.9(1) \AA. The corresponding band widths have been derived and the relative band intensity ratios are discussed.
[19]  [pdf] - 1584108
The Structure of the Young Star Cluster NGC 6231. I. Stellar Population
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
Submitted: 2017-05-31
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.
[20]  [pdf] - 1564021
Candidate X-ray-Emitting OB Stars in the MYStIX Massive Star-Forming Regions
Comments: 49 pages, 4 figure sets, 5 tables, accepted for publication in ApJ
Submitted: 2017-01-23
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.
[21]  [pdf] - 1557793
Star Formation In Nearby Clouds (SFiNCs): X-ray And Infrared Source Catalogs And Membership
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: . The 5GB SPCM Atlas can be found on-line at:
Submitted: 2016-12-15
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.
[22]  [pdf] - 1503548
Young Stellar Clusters Containing Massive Young Stellar Objects in the VVV Survey
Comments: 45 pages; 17 figures; published in AJ
Submitted: 2016-10-26
The purpose of this research is to study the connection of global properties of eight young stellar clusters projected in the Vista Variables in the Via Lactea (VVV) ESO Large Public Survey disk area and their young stellar object population. The analysis in based on the combination of spectroscopic parallax-based reddening and distance determinations with main sequence and pre-main sequence ishochrone fitting to determine the basic parameters (reddening, age, distance) of the sample clusters. The lower mass limit estimations show that all clusters are low or intermediate mass (between 110 and 1800 Mo), the slope Gamma of the obtained present-day mass functions of the clusters is close to the Kroupa initial mass function. On the other hand, the young stellar objects in the surrounding cluster's fields are classified by low resolution spectra, spectral energy distribution fit with theoretical predictions, and variability, taking advantage of multi-epoch VVV observations. All spectroscopically confirmed young stellar objects (except one) are found to be massive (more than 8 Mo). Using VVV and GLIMPSE color-color cuts we have selected a large number of new young stellar object candidates, which are checked for variability and 57% are found to show at least low-amplitude variations. In few cases it was possible to distinguish between YSO and AGB classification on the basis of the light curves.
[23]  [pdf] - 1531685
Young Stellar Populations in MYStIX Star Forming Regions: Candidate Protostars
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
Submitted: 2016-09-21, last modified: 2016-09-26
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.
[24]  [pdf] - 1576377
VVV High proper motion stars I. The catalogue of bright Ks < 13.5 stars
Comments: 15 pages, 7 figures, accepted in the MNRAS
Submitted: 2016-09-12
Knowledge of the stellar content near the Sun is important for a broad range of topics ranging from the search for planets to the study of Milky Way structure. The most powerful method for identifying potentially nearby stars is proper motion (PM) surveys. All old optical surveys avoid, or are at least substantially incomplete, near the Galactic plane. The depth and breadth of the "Vista Variables in Via Lactea" (VVV) near-IR survey significantly improves this situation. Taking advantage of the VVV survey database, we have measured PMs in the densest regions of the MW bulge and southern plane in order to complete the census of nearby objects. We have developed a custom PM pipeline based on VVV catalogues from the Cambridge Astronomy Survey Unit (CASU), by comparing the first epoch of JHKs with the multi-epoch Ks-bands acquired later. Taking advantage of the large time baseline between the 2MASS and the VVV observations, we also obtained 2MASS-VVV PMs. We present a near-IR proper motion catalogue for the whole area of the VVV survey, which includes 3003 moving stellar sources. All of these have been visually inspected and are real PM objects. Our catalogue is in very good agreement with the proper motion data supplied in IR catalogues outside the densest zone of the MW. The majority of the PM objects in our catalogue are nearby M-dwarfs, as expected. This new database allow us to identify 57 common proper motion binary candidates, among which are two new systems within ~30pc of the Sun.
[25]  [pdf] - 1457503
Near-infrared photometry and spectroscopy of the low Galactic latitude globular cluster 2MASS-GC03
Comments: Accepted for publication, MNRAS
Submitted: 2016-07-18
We present deep near-infrared photometry and spectroscopy of the globular cluster 2MASS-GC03 projected in the Galactic disk using MMIRS on the Clay telescope (Las Campanas Observatory) and VISTA Variables in the Via Lactea survey (VVV) data. Most probable cluster member candidates were identified from near-infrared photometry. Out of ten candidates that were followed-up spectroscopically, five have properties of cluster members, from which we calculate <[Fe/H]> = -0.9 +- 0.2 and a radial velocity of v_r > = -78 +- 12km/s. A distance of 10.8kpc is estimated from 3 likely RRLyrae members. Given that the cluster is currently at a distance of 4.2kpc from the Galactic center, the cluster's long survival time of an estimated 11.3 +- 1.2Gyr strengthens the case for its globular-cluster nature. The cluster has a hint of elongation in the direction of the Galactic center.
[26]  [pdf] - 1365701
Ultra-Low-Temperature Reactions of Carbon Atoms with Hydrogen Molecules
Submitted: 2016-02-29
The reactions of carbon atoms with dihydrogen have been investigated in liquid helium droplets at $T$ = 0.37 K. A calorimetric technique was applied to monitor the energy released in the reaction. The barrierless reaction between a single carbon atom and a single dihydrogen molecule was detected. Reactions between dihydrogen clusters and carbon atoms have been studied by high-resolution mass spectrometry. The formation of hydrocarbon cations of the type C$_m$H$_n^+$, with $m$ = 1-4 and $n$ = 1-15 was observed. With enhanced concentration of dihydrogen, the mass spectra demonstrated the main "magic" peak assigned to the CH$_5^+$ cation. A simple formation pathway and the high stability of this cation suggest its high abundance in the interstellar medium.
[27]  [pdf] - 1374349
Massive open star clusters using the VVV survey V. Young clusters with an OB stellar population
Comments: Accepted for publication, A&A
Submitted: 2016-02-15, last modified: 2016-02-18
The ESO public survey VISTA Variables in the V\'ia L\'actea (VVV) has contributed with deep multi-epoch photometry of the Galactic bulge and the adjacent part of the disk over 526 square degrees. More than a hundred cluster candidates have been reported thanks to this survey. We present the fifth article in a series of papers focused on young and massive clusters discovered in the VVV survey. In this paper, we present the physical characterization of five clusters with a spectroscopically confirmed OB-type stellar population. To characterize the clusters, we used near-infrared photometry ($J$, $H,$ and $K_S$) from the VVV survey and near-infrared $K$-band spectroscopy from ISAAC at VLT, following the methodology presented in the previous articles of the series. All clusters in our sample are very young (ages between 1-20 Myr), and their total mass are between $(1.07^{+0.40}_{-0.30})\cdot10^2$ $M_{\odot}$ and $(4.17^{+4.15}_{-2.08})\cdot10^3$ $M_{\odot}$. We observed a relation between the clusters total mass $M_{ecl}$ and the mass of their most massive stellar member $m_{max}$, even for clusters with an age $<$ 10 Myr.
[28]  [pdf] - 1351603
Observational constraints on star cluster formation theory - I. The mass-radius relation
Comments: 13 pages, 7 figures, accepted by A&A
Submitted: 2015-12-01
Stars form predominantly in groups usually denoted as clusters or associations. The observed stellar groups display a broad spectrum of masses, sizes and other properties, so it is often assumed that there is no underlying structure in this diversity. Here we show that the assumption of an unstructured multitude of cluster or association types might be misleading. Current data compilations of clusters show correlations between cluster mass, size, age, maximum stellar mass etc. In this first paper we take a closer look at the correlation of cluster mass and radius. We use literature data to explore relations in cluster and molecular core properties in the solar neighborhood. We show that for embedded clusters in the solar neighborhood there exists a clear correlation between cluster mass and half-mass radius of the form $M_c = C R_c^{\gamma}$ with gamma = 1.7 +/-0.2. This correlation holds for infra red K band data as well as X-ray sources and for clusters containing a hundred stars up to those consisting of a few tens of thousands of stars. The correlation is difficult to verify for clusters containing <30 stars due to low-number statistics. Dense clumps of gas are the progenitors of the embedded clusters. We find a similar slope for the mass-size relation of dense, massive clumps as for the embedded star clusters. This might point at a direct translation from gas to stellar mass: however, it is difficult to relate size measurements for clusters (stars) to those for gas profiles. Taking into account multiple paths for clump mass into cluster mass, we obtain an average star-formation efficiency of 18%{+9.3}{-5.7} for the embedded clusters in the solar neighborhood. The derived mass-radius relation gives constraints for the theory of clustered star formation. Analytical models and simulations of clustered star formation have to reproduce this relation in order to be realistic (abridged)
[29]  [pdf] - 1312078
Massive open star clusters using the VVV survey IV. WR 62-2, a new very massive star in the core of the VVV CL041 cluster
Comments: 8 pages, 7 figures, 2 tables, Accepted to A&A
Submitted: 2015-10-08
Context The ESO Public Survey VISTA Variables in the V\'ia L\'actea (VVV) provides deep multi-epoch infrared observations for an unprecedented 562 sq. degrees of the Galactic bulge and adjacent regions of the disk. Nearly 150 new open clusters and cluster candidates have been discovered in this survey. Aims We present the fourth article in a series of papers focussed on young and massive clusters discovered in the VVV survey. This article is dedicated to the cluster VVV CL041, which contains a new very massive star candidate, WR 62-2. Methods Following the methodology presented in the first paper of the series, wide-field, deep JHKs VVV observations, combined with new infrared spectroscopy, are employed to constrain fundamental parameters (distance, reddening, mass, age) of VVV CL041. Results We confirm that the cluster VVV CL041 is a young (less than 4 Myrs) and massive (3 +/- 2 x 10^3 Msol) cluster, and not a simple asterism. It is located at a distance of 4.2 +/- 0.9 kpc, and its reddening is A_V = 8.0 +/- 0.2 mag, which is slightly lower than the average for the young clusters towards the centre of the Galaxy. Spectral analysis shows that the most luminous star of the cluster, of the WN8h spectral type, is a candidate to have an initial mass larger than 100 Msol.
[30]  [pdf] - 1289903
A not so massive cluster hosting a very massive star
Comments: 1 page, 1 figure, to be published in the "International Workshop on Wolf-Rayet Stars conference proceedings" by Universit\"atsverlag Potsdam (editors: W.-R. Hamann, A. Sander, and H. Todt)
Submitted: 2015-10-07
We present the first physical characterization of the young open cluster VVV CL041. We spectroscopically observed the cluster main-sequence stellar population and a very-massive star candidate: WR62-2. CMFGEN modeling to our near-infrared spectra indicates that WR62-2 is a very luminous (10$^{6.4\pm0.2} L_{\odot}$) and massive ($\sim80 M_{\odot}$) star.
[31]  [pdf] - 1280965
No evidence for protoplanetary disk destruction by OB stars in the MYStIX sample
Comments: 24 pages, 32 figures. ApJ, accepted
Submitted: 2015-08-06
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.
[32]  [pdf] - 1300238
The Spatial Structure of Young Stellar Clusters. III. Physical Properties and Evolutionary States
Comments: Accepted for publication in The Astrophysical Journal ; 48 pages, 13 figures, and 6 tables
Submitted: 2015-07-20
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.
[33]  [pdf] - 1224114
The Spatial Structure of Young Stellar Clusters. II. Total Young Stellar Populations
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 (
Submitted: 2015-01-21
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.
[34]  [pdf] - 1208430
The Spatial Structure of Young Stellar Clusters. I. Subclusters
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 (
Submitted: 2014-03-17
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.
[35]  [pdf] - 1208320
Age Gradients in the Stellar Populations of Massive Star Forming Regions Based on a New Stellar Chronometer
Comments: Accepted for publication in ApJ; 89 pages, 23 figures, 2 Tables; High quality version is at
Submitted: 2014-03-11
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.
[36]  [pdf] - 1208321
Core-Halo Age Gradients and Star Formation in the Orion Nebula and NGC~2024 Young Stellar Clusters
Comments: Accepted for publication in ApJ; 42 pages, 5 figures, and 1 table; High quality version is at
Submitted: 2014-03-11
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.
[37]  [pdf] - 1208305
The Massive Star-forming Regions Omnibus X-ray Catalog
Comments: Accepted to ApJS, March 3, 2014. 51 pages, 25 figures
Submitted: 2014-03-11
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.
[38]  [pdf] - 1521922
A Massive Young Star-Forming Complex Study in Infrared and X-ray: X-ray Sources in Ten Star Forming Regions
Comments: 32 pages, 6 figures. Accepted for publication in Astrophysical Journal Supplements (with 6 other MYStIX papers). Complete versions of all MYStIX papers are available at
Submitted: 2013-09-17
The Massive Young star-forming complex Study in Infrared and X-ray (MYStIX) uses data from the Chandra X-ray Observatory to identify and characterize the young stellar populations in twenty Galactic (d < 4 kpc) massive star-forming regions. Here, the X-ray analysis for Chandra ACIS-I observations of ten of the MYStIX fields is described, and a catalog of >10,000 X-ray sources is presented. In comparison to other published Chandra source lists for the same regions, the number of MYStIX detected faint X-ray sources in a region is often doubled. While the higher catalog sensitivity increases the chance of false detections, it also increases the number of matches to infrared stars. X-ray emitting contaminants include foreground stars, background stars, and extragalactic sources. The X-ray properties of sources in these classes are discussed.
[39]  [pdf] - 1179328
Identifying Young Stars in Massive Star-Forming Regions for the MYStIX Project
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
Submitted: 2013-09-17
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
[40]  [pdf] - 1179327
The MYStIX InfraRed-Excess Source Catalog
Comments: 67 pages, 11 figures. Accepted for publication in Astrophysical Journal Supplements (with 6 other MYStIX papers). Complete versions of all MYStIX papers are available at
Submitted: 2013-09-17
The MYStIX project (Massive Young Star-Forming Complex Study in Infrared and X-rays) provides a comparative study of 20 Galactic massive star-forming complexes (d = 0.4 to 3.6 kpc). Probable stellar members in each target complex are identified using X-ray and/or infrared data via two pathways: (1) X-ray detections of young/massive stars with coronal activity/strong winds; or (2) infrared excess (IRE) selection of young stellar objects (YSOs) with circumstellar disks and/or protostellar envelopes. We present the methodology for the second pathway, using Spitzer/IRAC, 2MASS, and UKIRT imaging and photometry. Although IRE selection of YSOs is welltrodden territory, MYStIX presents unique challenges. We combine IR spectral energy distribution (SED) fitting with IR color cuts and spatial clustering analysis to identify IRE sources and isolate probable YSO members in each MYStIX target field from the myriad types of contaminating sources that can resemble YSOs: extragalactic sources, evolved stars, nebular knots, and even unassociated foreground/background YSOs. Applying our methodology consistently across 18 of the target complexes, we produce the MYStIX IRE Source (MIRES) Catalog comprising 20,719 sources, including 8686 probable stellar members of the MYStIX target complexes. We also classify the SEDs of 9365 IR counterparts to MYStIX X-ray sources to assist the first pathway, the identification of X-ray detected stellar members. The MIRES catalog provides a foundation for follow-up studies of diverse phenomena related to massive star cluster formation, including protostellar outflows, circumstellar disks, and sequential star formation triggered by massive star feedback processes.
[41]  [pdf] - 1179319
Overview of the Massive Young Star-Forming Complex Study in Infrared and X-ray (MYStIX) Project
Comments: 74 pages, 10 figures. To appear in Astrophysical Journal Supplements (with 6 other MYStIX papers). A version with complete figures is available at
Submitted: 2013-09-17
MYStIX (Massive Young Star-Forming Complex Study in Infrared and Xray) seeks to characterize 20 OB-dominated young clusters and their environs at distances d < 4 kpc using imaging detectors on the Chandra X-ray Observatory, Spitzer Space Telescope, and the United Kingdom InfraRed Telescope. The observational goals are to construct catalogs of star-forming complex stellar members with well-defined criteria, and maps of nebular gas (particularly of hot X-ray emitting plasma) and dust. A catalog of MYStIX Probable Complex Members (MPCMs) with several hundred OB stars and > 30, 000 low mass premain sequence is assembled. This sample and related data products will be used to seek new empirical constraints on theoretical models of cluster formation and dynamics, mass segregation, OB star formation, star formation triggering on the periphery of HII regions, the survivability of protoplanetary disks in HII regions. This paper give an introduction and overview of the project, covering the data analysis methodology and application to two star forming regions, NGC 2264 and the Trifid Nebula.
[42]  [pdf] - 1521923
A Massive Young Star-Forming Complex Study in Infrared and X-ray: Mid-Infrared Observations and Catalogs
Comments: 34 pages, 10 figures. Accepted for Astrophysical Journal Supplements (with 6 other MYStIX papers). Complete versions of this and other MYStIX papers are available at
Submitted: 2013-09-17
Spitzer IRAC observations and stellar photometric catalogs are presented for the Massive Young Star-Forming Complex Study in the Infrared and X-ray (MYStIX). MYStIX is a multiwavelength census of young stellar members of twenty nearby (d < 4 kpc), Galactic, star-forming regions that contain at least one O star. All regions have data available from the Spitzer Space Telescope, consisting of GLIMPSE or other published catalogs for eleven regions and results of our own photometric analysis of archival data for the remaining nine regions. This paper seeks to construct deep and reliable catalogs of sources from the Spitzer images. Mid-infrared study of these regions faces challenges of crowding and high nebulosity. Our new catalogs typically contain fainter sources than existing Spitzer studies, which improves the match rate to Chandra X-ray sources that are likely to be young stars, but increases the possibility of spurious point-source detections, especially peaks in the nebulosity. IRAC color-color diagrams help distinguish spurious detections of nebular PAH emission from the infrared excess associated with dusty disks around young stars. The distributions of sources on the mid-infrared color-magnitude and color-color diagrams reflect differences between MYStIX regions, including astrophysical effects such as stellar ages and disk evolution.
[43]  [pdf] - 550563
MYStIX First Results: Spatial Structures of Massive Young Stellar Clusters
Comments: 4 pages, 1 figure; to appear in The Labyrinth of Star Formation, (eds.) D. Stamatellos, S. Goodwin, and D. Ward-Thompson, Springer, in press
Submitted: 2012-08-16
Observations of the spatial distributions of young stars in star-forming regions can be linked to the theory of clustered star formation using spatial statistical methods. The MYStIX project provides rich samples of young stars from the nearest high-mass star-forming regions. Maps of stellar surface density reveal diverse structure and subclustering. Young stellar clusters and subclusters are fit with isothermal spheres and ellipsoids using the Bayesian Information Criterion to estimate the number of subclusters. Clustering is also investigated using Cartwright and Whitworth's Q statistic and the inhomogeneous two-point correlation function. Mass segregation is detected in several cases, in both centrally concentrated and fractally structured star clusters, but a few clusters are not mass segregated.
[44]  [pdf] - 1150574
The Elephant Trunk Nebula and the Trumpler 37 cluster: Contribution of triggered star formation to the total population of an HII region
Comments: Accepted for publication in MNRAS; 28 pages, 18 figures
Submitted: 2012-08-07
Rich young stellar clusters produce HII regions whose expansion into the nearby molecular cloud is thought to trigger the formation of new stars. However, the importance of this mode of star formation is uncertain. This investigation seeks to quantify triggered star formation (TSF) in IC 1396A (a.k.a., the Elephant Trunk Nebula), a bright rimmed cloud (BRC) on the periphery of the nearby giant HII region IC 1396 produced by the Trumpler 37 cluster. X-ray selection of young stars from Chandra X-ray Observatory data is combined with existing optical and infrared surveys to give a more complete census of the TSF population. Over 250 young stars in and around IC 1396A are identified; this doubles the previously known population. A spatio-temporal gradient of stars from the IC 1396A cloud toward the primary ionizing star HD 206267 is found. We argue that the TSF mechanism in IC 1396A is the radiation-driven implosion process persisting over several million years. Analysis of the X-ray luminosity and initial mass functions indicates that >140 stars down to 0.1 Msun were formed by TSF. Considering other BRCs in the IC 1396 HII region, we estimate the TSF contribution for the entire HII region exceeds 14-25% today, and may be higher over the lifetime of the HII region. Such triggering on the periphery of HII regions may be a significant mode of star formation in the Galaxy.
[45]  [pdf] - 1041548
A Chandra Observation of the Obscured Star-Forming Complex W40
Comments: Accepted for publication in The Astrophysical Journal. 60 pages, 16 figures
Submitted: 2010-10-26
The young stellar cluster illuminating the W40 H II region, one of the nearest massive star forming regions, has been observed with the ACIS detector on board the Chandra X-ray Observatory. Due to its high obscuration, this is a poorly-studied stellar cluster with only a handful of bright stars visible in the optical band, including three OB stars identified as primary excitation sources. We detect 225 X-ray sources, of which 85% are confidently identified as young stellar members of the region. Two potential distances of the cluster, 260 pc and 600 pc, are used in the paper. Supposing the X-ray luminosity function to be universal, it supports a 600 pc distance as a lower limit for W40 and a total population of at least 600 stars down to 0.1 Mo under the assumption of a coeval population with a uniform obscuration. In fact, there is strong spatial variation in Ks-band-excess disk fraction and non-uniform obscuration due to a dust lane that is identified in absorption in optical, infrared and X-ray. The dust lane is likely part of a ring of material which includes the molecular core within W40. In contrast to the likely ongoing star formation in the dust lane, the molecular core is inactive. The star cluster has a spherical morphology, an isothermal sphere density profile, and mass segregation down to 1.5 Mo. However, other cluster properties, including a \leq{1} Myr age estimate and ongoing star formation, indicate that the cluster is not dynamically relaxed. X-ray diffuse emission and a powerful flare from a young stellar object are also reported.
[46]  [pdf] - 10493
Chandra spectroscopy of the hot star beta Crucis and the discovery of a pre-main-sequence companion
Comments: Accepted for publication in MNRAS; 19 pages, 15 figures, some in color; version with higher-resolution figures available at
Submitted: 2008-02-27
In order to test the O star wind-shock scenario for X-ray production in less luminous stars with weaker winds, we made a pointed 74 ks observation of the nearby early B giant, beta Cru (B0.5 III), with the Chandra HETGS. We find that the X-ray spectrum is quite soft, with a dominant thermal component near 3 million K, and that the emission lines are resolved but quite narrow, with half-widths of 150 km/s. The forbidden-to-intercombination line ratios of Ne IX and Mg XI indicate that the hot plasma is distributed in the wind, rather than confined near the photosphere. It is difficult to understand the X-ray data in the context of the standard wind-shock paradigm for OB stars, primarily because of the narrow lines, but also because of the high X-ray production efficiency. A scenario in which the bulk of the outer wind is shock heated is broadly consistent with the data, but not very well motivated theoretically. It is possible that magnetic channeling could explain the X-ray properties, although no field has been detected on beta Cru. We detected periodic variability in the hard (hnu > 1 keV) X-rays, modulated on the known optical period of 4.58 hours, which is the period of the primary beta Cep pulsation mode for this star. We also have detected, for the first time, an apparent companion to beta Cru at a projected separation of 4 arcsec. This companion was likely never seen in optical images because of the presumed very high contrast between it and beta Cru in the optical. However, the brightness contrast in the X-ray is only 3:1, which is consistent with the companion being an X-ray active low-mass pre-main-sequence star. The companion's X-ray spectrum is relatively hard and variable, as would be expected from a post T Tauri star.