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.

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/).

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.

This Astro2020 State of the Profession Consideration White Paper highlights the growth of astrostatistics and astroinformatics in astronomy, identifies key issues hampering the maturation of these new subfields, and makes recommendations for structural improvements at different levels that, if acted upon, will make significant positive impacts across astronomy.

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.

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.

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.

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.

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.

We write in response to the call from the 2020 Decadal Survey to submit white papers illustrating the most pressing scientific questions in astrophysics for the coming decade. We propose exploration as the central question for the Decadal Committee's discussions.The history of astronomy shows that paradigm changing discoveries are not driven by well formulated scientific questions, based on the knowledge of the time. They were instead the result of the increase in discovery space fostered by new telescopes and instruments. An additional tool for increasing the discovery space is provided by the analysis and mining of the increasingly larger amount of archival data available to astronomers. Revolutionary observing facilities, and the state of the art astronomy archives needed to support these facilities, will open up the universe to new discovery. Here we focus on exploration for compact objects and multi messenger science. This white paper includes science examples of the power of the discovery approach, encompassing all the areas of astrophysics covered by the 2020 Decadal Survey.

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.

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.

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^\alpha$ `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.

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.

One of the aims of LSST is to perform a systematic survey of star clusters and star forming regions (SFRs) in our Galaxy. In particular, the observations obtained with LSST will make a big difference in Galactic regions that have been poorly studied in the past, such as the anticenter and the disk beyond the Galactic center, and they will have a strong impact in discovering new distant SFRs. These results can be achieved by exploiting the exquisite depth that will be attained if the wide-fast-deep (WFD) observing strategy of the main survey is also adopted for the Galactic plane, in the g, r, and i filters.

Young stars exhibit short-term photometric variability caused by mass accretion events from circumstellar disks, the presence of dusty warps within the inner disks, starspots that rotate across the stellar surfaces, and flares. Long-term variability also occurs owing to starspot longevity and cycles, and from changes in stellar angular momenta and activity as the stars age. We propose to observe the Carina star-forming region in different bands with a cadence of 30 minutes every night for one week per year to clarify the nature of both the short-term and long-term variability of the thousands of young stars in this region. By obtaining well-sampled multicolor lightcurves of this dense young cluster, LSST would acquire the first statistically significant data on how these objects vary on both short and long timescales. This information will allow us to relate the observed variability to stellar properties such as mass, age, binarity, and to environmental properties such as location within or exterior to the H II region, and to the presence or absence of a circumstellar disk.

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.

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 $\sigma$ 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 $\lambda$ 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.

The Orion Star Forming Complex (OSFC) is a central target for the APOGEE-2 Young Cluster Survey. Existing membership catalogs span limited portions of the OSFC, reflecting the difficulty of selecting targets homogeneously across this extended, highly structured region. We have used data from wide field photometric surveys to produce a less biased parent sample of young stellar objects (YSOs) with infrared (IR) excesses indicative of warm circumstellar material or photometric variability at optical wavelengths across the full 420 square degrees extent of the OSFC. When restricted to YSO candidates with H < 12.4, to ensure S/N ~100 for a six visit source, this uniformly selected sample includes 1307 IR excess sources selected using criteria vetted by Koenig & Liesawitz and 990 optical variables identified in the Pan-STARRS1 3$\pi$ survey: 319 sources exhibit both optical variability and evidence of circumstellar disks through IR excess. Objects from this uniformly selected sample received the highest priority for targeting, but required fewer than half of the fibers on each APOGEE-2 plate. We fill the remaining fibers with previously confirmed and new color-magnitude selected candidate OSFC members. Radial velocity measurements from APOGEE-1 and new APOGEE-2 observations taken in the survey's first year indicate that ~90% of the uniformly selected targets have radial velocities consistent with Orion membership.The APOGEE-2 Orion survey will include >1100 bona fide YSOs whose uniform selection function will provide a robust sample for comparative analyses of the stellar populations and properties across all sub-regions of Orion.

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).

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.

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.

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.

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.

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\alpha}$ vs. [NII]/H$\alpha$ (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$\beta$, $\log$ [NII]/H$\alpha$, and $\log$ EW(H${\alpha}$), 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/).

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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

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.

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

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.

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.

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

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.

Stars over a wide range of masses and evolutionary stages are nowadays known to emit X-rays. This X-ray emission is a unique probe of the most energetic phenomena occurring in the circumstellar environment of these stars, and provides precious insight on magnetic phenomena or hydrodynamic shocks. Owing to its large collecting area, Athena+ will open up an entirely new window on these phenomena. Indeed, Athena+ will not only allow us to study many more objects with an unprecedented spectral resolution, but will also pioneer the study of the dynamics of these objects via time-resolved high-resolution spectroscopy. In this way, Athena+ will be a unique tool to study accretion processes in TTauri stars, flaring activity in young stars, dynamos in ultra-cool dwarfs, small and large-scale structures in the winds of single massive stars, wind interactions in massive binary systems, hot bubbles in planetary nebula... All these studies will lead to a deeper understanding of yet poorly understood processes which have profound impact in star and planetary system formation as well as in feedback processes on Galactic scale.

This White Paper, submitted to the recent ESA call for science themes to define its future large missions, advocates the need for a transformational leap in our understanding of two key questions in astrophysics: 1) How does ordinary matter assemble into the large scale structures that we see today? 2) How do black holes grow and shape the Universe? Hot gas in clusters, groups and the intergalactic medium dominates the baryonic content of the local Universe. To understand the astrophysical processes responsible for the formation and assembly of these large structures, it is necessary to measure their physical properties and evolution. This requires spatially resolved X-ray spectroscopy with a factor 10 increase in both telescope throughput and spatial resolving power compared to currently planned facilities. Feedback from supermassive black holes is an essential ingredient in this process and in most galaxy evolution models, but it is not well understood. X-ray observations can uniquely reveal the mechanisms launching winds close to black holes and determine the coupling of the energy and matter flows on larger scales. Due to the effects of feedback, a complete understanding of galaxy evolution requires knowledge of the obscured growth of supermassive black holes through cosmic time, out to the redshifts where the first galaxies form. X-ray emission is the most reliable way to reveal accreting black holes, but deep survey speed must improve by a factor ~100 over current facilities to perform a full census into the early Universe. The Advanced Telescope for High Energy Astrophysics (Athena+) mission provides the necessary performance (e.g. angular resolution, spectral resolution, survey grasp) to address these questions and revolutionize our understanding of the Hot and Energetic Universe. These capabilities will also provide a powerful observatory to be used in all areas of astrophysics.

The Hubble Space Telescope (HST) Treasury Program on the Orion Nebula Cluster has used 104 orbits of HST time to image the Great Orion Nebula region with the Advanced Camera for Surveys (ACS), the Wide-Field/Planetary Camera 2 (WFPC2) and the Near Infrared Camera and Multi Object Spectrograph (NICMOS) instruments in 11 filters ranging from the U-band to the H-band equivalent of HST. The program has been intended to perform the definitive study of the stellar component of the ONC at visible wavelengths, addressing key questions like the cluster IMF, age spread, mass accretion, binarity and cirumstellar disk evolution. The scanning pattern allowed to cover a contiguous field of approximately 600 square arcminutes with both ACS and WFPC2, with a typical exposure time of approximately 11 minutes per ACS filter, corresponding to a point source depth AB(F435W) = 25.8 and AB(F775W)=25.2 with 0.2 magnitudes of photometric error. We describe the observations, data reduction and data products, including images, source catalogs and tools for quick look preview. In particular, we provide ACS photometry for 3399 stars, most of them detected at multiple epochs, WFPC2 photometry for 1643 stars, 1021 of them detected in the U-band, and NICMOS JH photometry for 2116 stars. We summarize the early science results that have been presented in a number of papers. The final set of images and the photometric catalogs are publicly available through the archive as High Level Science Products at the STScI Multimission Archive hosted by the Space Telescope Science Institute.

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

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

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

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.

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.

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

The VAO (Virtual Astronomical Observatory) Science Council (VAO-SC) met on July 27-28, 2011 at the Harvard-Smithsonian Center for Astrophysics in Cambridge MA, to review the VAO performance during its first year of operations. In this meeting the VAO demonstrated the new tools for astronomers that are being released in September 2011 and presented plans for the second year of activities, resulting from studies conducted during the first year. This document contains the recommendations of the VAO-SC for the second year of activity of the VAO.

We present a deep Spitzer/IRAC survey of the OB association IC 1795 carried out to investigate the evolution of protoplanetary disks in regions of massive star formation. Combining Spitzer/IRAC data with Chandra/ACIS observations, we find 289 cluster members. An additional 340 sources with an infrared excess, but without X-ray counterpart, are classified as cluster member candidates. Both surveys are complete down to stellar masses of about 1 Msun. We present pre-main sequence isochrones computed for the first time in the Spitzer/IRAC colors. The age of the cluster, determined via the location of the Class III sources in the [3.6]-[4.5]/[3.6] color-magnitude diagram, is in the range of 3 - 5 Myr. As theoretically expected, we do not find any systematic variation in the spatial distribution of disks within 0.6 pc of either O-type star in the association. However, the disk fraction in IC 1795 does depend on the stellar mass: sources with masses >2 Msun have a disk fraction of ~20%, while lower mass objects (2-0.8 Msun) have a disk fraction of ~50%. This implies that disks around massive stars have a shorter dissipation timescale.

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

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

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

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

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

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

The Great Nebula in Carina provides an exceptional view into the violent massive star formation and feedback that typifies giant HII regions and starburst galaxies. We have mapped the Carina star-forming complex in X-rays, using archival Chandra data and a mosaic of 20 new 60ks pointings using the Chandra X-ray Observatory's Advanced CCD Imaging Spectrometer, as a testbed for understanding recent and ongoing star formation and to probe Carina's regions of bright diffuse X-ray emission. This study has yielded a catalog of properties of >14,000 X-ray point sources; >9800 of them have multiwavelength counterparts. Using Chandra's unsurpassed X-ray spatial resolution, we have separated these point sources from the extensive, spatially-complex diffuse emission that pervades the region; X-ray properties of this diffuse emission suggest that it traces feedback from Carina's massive stars. In this introductory paper, we motivate the survey design, describe the Chandra observations, and present some simple results, providing a foundation for the 15 papers that follow in this Special Issue and that present detailed catalogs, methods, and science results.

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.

Recommendations of the VAO-Science Council following the meeting of March 26-27, 2010. Refer to the webpage http://www.aui.edu/vao.php

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

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

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

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

V2129 Oph is a 1.35 solar mass classical T Tauri star, known to possess a strong and complex magnetic field. By extrapolating from an observationally derived magnetic surface map, obtained through Zeeman-Doppler imaging, models of V2129 Oph's corona have been constructed, and used to make predictions regarding the global X-ray emission measure, the amount of modulation of X-ray emission, and the density of accretion shocks. In late June 2009 we will under take an ambitious multi-wavelength, multi-observing site, and near contemporaneous campaign, combining spectroscopic optical, nIR, UV, X-ray, spectropolarimetric and photometric monitoring. This will allow the validity of the 3D field topologies derived via field extrapolation to be determined.

We present IRS spectra and revised MIPS photometry for the 18 members of the Eta Chamaeleontis cluster. Aged 8 Myr, the Eta Cha cluster is one of the few nearby regions within the 5-10 Myr age range, during which the disk fraction decreases dramatically and giant planet formation must come to an end. For the 15 low-mass members, we measure a disk fraction ~50%, high for their 8 Myr age, and 4 of the 8 disks lack near-IR excesses, consistent with the empirical definition of "transition'' disks. Most of the disks are comparable to geometrically flat disks. The comparison with regions of different ages suggests that at least some of the "transition" disks may represent the normal type of disk around low-mass stars. Therefore, their flattened structure and inner holes may be related to other factors (initial masses of the disk and the star, environment, binarity), rather than to pure time evolution. We analyze the silicate dust in the disk atmosphere, finding moderate crystalline fractions (~10-30%) and typical grain sizes ~1-3 micron, without any characteristic trend in the composition. These results are common to other regions of different ages, suggesting that the initial grain processing occurs very early in the disk lifetime (<1 Myr). Large grain sizes in the disk atmosphere cannot be used as a proxy for age, but are likely related to higher disk turbulence. The dust mineralogy varies between the 8-12micron and the 20-30 micron features, suggesting high temperature dust processing and little radial mixing. Finally, the analysis of IR and optical data on the B9 star Eta Cha reveals that it is probably surrounded by a young debris disk with a large inner hole, instead of being a classical Be star.

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

The Cepheus B (CepB) molecular cloud and a portion of the nearby CepOB3b OB association, one of the most active regions of star formation within 1 kpc, have been observed with the IRAC detector on board the Spitzer Space Telescope. The goals are to study protoplanetary disk evolution and processes of sequential triggered star formation in the region. Out of ~400 pre-main sequence (PMS) stars selected with an earlier Chandra X-ray Observatory observation, 95% are identified with mid-infrared sources and most of these are classified as diskless or disk-bearing stars. The discovery of the additional >200 IR-excess low-mass members gives a combined Chandra+Spitzer PMS sample complete down to 0.5 Mo outside of the cloud, and somewhat above 1 Mo in the cloud. Analyses of the nearly disk-unbiased combined Chandra+Spitzer selected stellar sample give several results. Our major finding is a spatio-temporal gradient of young stars from the hot molecular core towards the primary ionizing O star HD 217086. This strongly supports the radiation driven implosion (RDI) model of triggered star formation in the region. The empirical estimate for the shock velocity of 1 km/s is very similar to theoretical models of RDI in shocked molecular clouds...ABRIDGED... Other results include: 1. agreement of the disk fractions, their mass dependency, and fractions of transition disks with other clusters; 2. confirmation of the youthfulness of the embedded CepB cluster; 3. confirmation of the effect of suppression of time-integrated X-ray emission in disk-bearing versus diskless systems.

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

While it may seem counterintuitive that X-ray astronomy should give any insights into low-temperature planetary systems, planets orbit stars whose magnetized surfaces divert a small fraction of the stellar energy into high energy products: coronal UV and X-rays, flare X-rays and energetic particles, and a high-velocity stellar wind. In our Solar System, X-ray emission gives unique insights into the solar activity, planetary atmospheres, cometary comae, charge exchange physics, and space weather across the Solar System. The stellar activity of young stars is greatly elevated and can substantially affect protoplanetary disks and planet formation processes. We highlight six studies achievable with the planned International X-ray Observatory which address in unique ways issues in planetary sciences: probing X-ray irradiation of protoplanetary disks with the iron fluorescent line and its effects on disk turbulence; study the complex charge-exchange X-ray emission from Jupiter and the Martian exosphere; elucidate charge-exchange processes in cometary comae; understanding heliospheric charge-exchange emission and the interpretation of the soft X-ray background; and examining the magnetic activity of planet-hosting stars and its evaporation of planetary atmospheres.

Most stars are born in rich young stellar clusters (YSCs) embedded in giant molecular clouds. The most massive stars live out their short lives there, profoundly influencing their natal environments by ionizing HII regions, inflating wind-blown bubbles, and soon exploding as supernovae. Thousands of lower-mass pre-main sequence stars accompany the massive stars, and the expanding HII regions paradoxically trigger new star formation as they destroy their natal clouds. While this schematic picture is established, our understanding of the complex astrophysical processes involved in clustered star formation have only just begun to be elucidated. The technologies are challenging, requiring both high spatial resolution and wide fields at wavelengths that penetrate obscuring molecular material and remove contaminating Galactic field stars. We outline several important projects for the coming decade: the IMFs and structures of YSCs; triggered star formation around YSC; the fate of OB winds; the stellar populations of Infrared Dark Clouds; the most massive star clusters in the Galaxy; tracing star formation throughout the Galactic Disk; the Galactic Center region and YSCs in the Magellanic Clouds. Programmatic recommendations include: developing a 30m-class adaptive optics infrared telescope; support for high-resolution and wide field X-ray telescopes; large-aperture sub-millimeter and far-infrared telescopes; multi-object infrared spectrographs; and both numerical and analytical theory.

Perhaps more than other physical sciences, astronomy is frequently statistical in nature. The objects under study are inaccessible to direct manipulation in the laboratory, so the astronomer is restricted to observing a few external characteristics and inferring underlying properties and physics. Astronomy played a profound role in the historical development of statistics from the ancient Greeks through the 19th century. But the fields drifted apart in the 20th century as astronomy turned towards astrophysics and statistics towards human affairs. Today we see a resurgence in astrostatistical activity with the proliferation of survey mega-datasets and the need to link complicated data to nonlinear astrophysical models. Several contemporary astrostatistical challenges are outlined: heteroscedastic measurement errors, censoring and truncation in multivariate databases; time series analysis of variable objects including dynamical models of extrasolar planetary systems; treatments of faint sources and other Poisson processes; the anisotropic spatial point process of galaxy clustering; and model fitting and selection for the cosmic microwave background.

While the exceptional sensitivity of Chandra and XMM-Newton has resulted in revolutionary studies of the Galactic neighborhood in the soft (<10 keV) X-ray band, there are many open questions. We discuss these issues and how they would be addressed by very wide-area (> 100 sq. deg.) X-ray surveys.

Pre-main sequence (PMS) stars are known to produce powerful X-ray flares which resemble magnetic reconnection solar flares scaled by factors up to 10^4. However, numerous puzzles are present including the structure of X-ray emitting coronae and magnetospheres, effects of protoplanetary disks, and effects of stellar rotation. To investigate these issues in detail, we examine 216 of the brightest flares from 161 PMS stars observed in the Chandra Orion Ultradeep Project (COUP). These constitute the largest homogeneous dataset of PMS, or indeed stellar flares at any stellar age, ever acquired. Our effort is based on a new flare spectral analysis technique that avoids nonlinear parametric modeling. It can be applied to much weaker flares and is more sensitive than standard methods. We provide a catalog with >30 derived flare properties and an electronic atlas for this unique collection of stellar X-ray flares. The current study (Paper I) examines the flare morphologies, and provides general comparison of COUP flare characteristics with those of other active X-ray stars and the Sun. Paper II will concentrate on relationships between flare behavior, protoplanetary disks, and other stellar properties. Several results are obtained. First, the COUP flares studied here are among the most powerful, longest, and hottest stellar X-ray flares ever studied. Second, no significant statistical differences in peak flare luminosity or temperature distributions are found among different morphological flare classes, suggesting a common underlying mechanism for all flares. Third, comparison with the general solar-scaling laws indicates that COUP flares may not fit adequately proposed power-temperature and duration-temperature solar-stellar fits. Fourth, COUP super-hot flares are found to be brighter but shorter than ... ABRIDGED

We study the properties of powerful X-ray flares from 161 pre-main sequence (PMS) stars observed with the Chandra X-ray Observatory in the Orion Nebula region. Relationships between flare properties, protoplanetary disks and accretion are examined in detail to test models of star-disk interactions at the inner edge of the accretion disks. Previous studies had found no differences in flaring between diskfree and accreting systems other than a small overall diminution of X-ray luminosity in accreting systems. The most important finding is that X-ray coronal extents in fast-rotating diskfree stars can significantly exceed the Keplerian corotation radius, whereas X-ray loop sizes in disky and accreting systems do not exceed the corotation radius. This is consistent with models of star-disk magnetic interaction where the inner disk truncates and confines the PMS stellar magnetosphere. We also find two differences between flares in accreting and diskfree PMS stars. First, a subclass of super-hot flares with peak plasma temperatures exceeding 100 MK are preferentially present in accreting systems. Second, we tentatively find that accreting stars produce flares with shorter durations. Both results may be consequences of the distortion and destabilization of the stellar magnetosphere by the interacting disk. Finally, we find no evidence that any flare types, even slow-rise flat-top flares are produced in star-disk magnetic loops. All are consistent with enhanced solar long-duration events with both footprints anchored in the stellar surface.

The origin and evolution of the X-rays in very young stellar objects (YSOs) are not yet well understood since it is very hard to observe YSOs in the protostellar phase. We study the X-ray properties of Class 0-I objects in the Orion Nebula Cluster (ONC) and compare them with those of the more evolved Class II and III members. Using Chandra Orion Ultradeep Project (COUP) data, we study the X-ray properties of stars in different evolutionary classes: luminosities, NH, temperatures and time variability are compared in order to understand if the interaction between the circumstellar material and the central object can influence the X-ray emission. We have assembled the deepest and most complete photometric catalog of objects in the ONC region from the UV to 8 microns using data from HST, WFI@2.2m ESO and ISPI@4m CTIO telescopes, and Spitzer IRAC. We select high probability candidate Class 0-I protostars, distinguishing between those having a spectral energy distribution which rises from K up to 8 microns (Class 0-Ia) from those where the SED rises from K up to 4.5 microns and decreasing afterwards (Class 0-Ib). We select a sample of bona fide Class II stars and a set of Class III stars with IR emission consistent with normal photospheres. Our principal result is that Class 0-Ia objects are significantly less luminous in X-rays, both in the total and hard bands, than the more evolved Class II stars with mass larger than 0.5 Msun; these latter show X-ray luminosities similar to those of Class 0-Ib stars. This result supports the hypothesis that the onset of X-ray emission occurs at a very early stage of star formation. Temporal variability and spectral properties of Class 0-I stars are similar to those of the more evolved Class II and III objects, except for a larger absorption likely due to gas in the circumstellar material.

We present the first high spatial resolution Chandra X-ray study of NGC 2244, the 2 Myr old stellar cluster immersed in the Rosette Nebula. Over 900 X-ray sources are detected; 77% have optical or FLAMINGOS near-infrared (NIR) stellar counterparts and are mostly previously uncatalogued young stellar cluster members. All known OB stars with spectral type earlier than B1 are detected and the X-ray selected stellar population is estimated to be nearly complete between 0.5 and 3 Msun. The X-ray luminosity function (XLF) ranges from 29.4<logLx<32 ergs/s in the hard (2-8keV) band. By comparing the NGC 2244 and Orion Nebula Cluster XLFs, we estimate a total population of 2000 stars in NGC 2244. A number of further results emerge from our analysis: The XLF and the associated K-band luminosity function indicate a normal Salpeter initial mass function (IMF) for NGC 2244. This is inconsistent with the top-heavy IMF reported from earlier optical studies that lacked a good census of <4Msun stars. The spatial distribution of X-ray stars is strongly concentrated around the central O5 star, HD 46150. The other early O star, HD 46223, has few companions. The cluster's stellar radial density profile shows two distinctive structures. This double structure, combined with the absence of mass segregation, indicates that this cluster is not in dynamical equilibrium. The spatial distribution of X-ray selected K-excess disk stars and embedded stars is asymmetric with an apparent deficit towards the north. The fraction of X-ray-selected cluster members with K-band excesses caused by inner protoplanetary disks is 6%, slightly lower than the 10% disk fraction estimated from the FLAMINGOS study based on the NIR-selected sample. This is due to the high efficiency of X-ray surveys in locating disk-free T Tauri stars.[Abridged]

We report on accretion- and outflow-related X-rays from T Tauri stars, based on results from the "XMM-Newton Extended Survey of the Taurus Molecular Cloud." X-rays potentially form in shocks of accretion streams near the stellar surface, although we hypothesize that direct interactions between the streams and magnetic coronae may occur as well. We report on the discovery of a "soft excess" in accreting T Tauri stars supporting these scenarios. We further discuss a new type of X-ray source in jet-driving T Tauri stars. It shows a strongly absorbed coronal component and a very soft, weakly absorbed component probably related to shocks in microjets. The excessive coronal absorption points to dust-depletion in the accretion streams.

An 800 sq-arcmin mosaic image of the W3 star forming complex obtained with the Chandra X-ray Observatory gives a valuable new view of the spatial structure of its young stellar populations. The Chandra image reveals about 1300 faint X-ray sources, most of which are PMS stars in the cloud. Some, but not all, of the high-mass stars producing hypercompact and ultracompact H II (UCHII) regions are also seen, as reported in a previous study. The Chandra images reveal three dramatically different embedded stellar populations. The W3 Main cluster extends over 7 pc with about 900 X-ray stars in a nearly-spherical distribution centered on the well-studied UCHII regions and high-mass protostars. The cluster surrounding the prototypical UCHII region W3(OH) shows a much smaller (<0.6 pc), asymmetrical, and clumpy distribution of about 50 PMS stars. The massive star ionizing the W3 North H II region is completely isolated without any accompanying PMS stars. In W3 Main, the inferred ages of the widely distributed PMS stars are significantly older than the inferred ages of the central OB stars illuminating the UCHIIs. We suggest that different formation mechanisms are necessary to explain the diversity of the W3 stellar populations: cluster-wide gravitational collapse with delayed OB star formation in W3 Main, collect-and-collapse triggering by shock fronts in W3(OH), and a runaway O star or isolated massive star formation in W3 North.

Context. X-ray flares are common phenomena in pre-main sequence stars. Their analysis gives insights into the physics at work in young stellar coronae. The Orion Nebula Cluster offers a unique opportunity to study large samples of young low mass stars. This work is part of the Chandra Orion Ultradeep project (COUP), an ~10 day long X-ray observation of the Orion Nebula Cluster (ONC). Aims. Our main goal is to statistically characterize the flare-like variability of 165 low mass (0.1-0.3 M_sun) ONC members in order to test and constrain the physical scenario in which flares explain all the observed emission. Methods. We adopt a maximum likelihood piece-wise representation of the observed X-ray light curves and detect flares by taking into account both the amplitude and time derivative of the count-rate. We then derive the frequency and energy distribution of the flares. Results. The high energy tail of the energy distribution of flares is well described by a power-law with index 2.2. We test the hypothesis that light curves are built entirely by overlapping flares with a single power law energy distribution. We constrain the parameters of this simple model for every single light curve. The analysis of synthetic light curves obtained from the model indicates a good agreement with the observed data. Comparing low mass stars with stars in the mass interval (0.9-1.2M_sun), we establish that, at ~1 Myr, low mass and solar mass stars of similar X-ray luminosity have very similar flare frequencies. Conclusions. Our observational results are consistent with the following model/scenario: the light curves are entirely built by over- lapping flares with a power-law intensity distribution; the intense flares are individually detected, while the weak ones merge and form a pseudo-quiescent level, which we indicate as the characteristic level.

The mysterious high galactic latitude cometary globule CG12 has been observed with the ACIS detector on board the Chandra X-ray Observatory. We detect 128 X-ray sources; half are likely young stars formed within the globule's head. This new population of >50 T-Tauri stars and one new embedded protostar is far larger than the previously reported few intermediate-mass and two protostellar members of the cloud. Most of the newly discovered stars have masses 0.2-0.7Mo, and 9-15% have K-band excesses from inner protoplanetary disks. X-ray properties provide an independent distance estimate consistent with CG12's unusual location >200pc above the Galactic plane. The star formation efficiency in CG12 appears to be 15-35%, far above that seen in other triggered molecular globules. The median photometric age found for the T-Tauri population is ~4Myr with a large spread of <1-20Myr and ongoing star formation in the molecular cores. The stellar age and spatial distributions are inconsistent with a simple radiation driven implosion (RDI) model, and suggest either that CG12 is an atypically large shocked globule, or it has been subject to several distinct episodes of triggering and ablation. We report a previously unnoticed group of B-type stars northwest of CG12 which may be the remnants of an OB association which produced multiple supernova explosions that could have shocked and ablated the cloud over a 15-30Myr period. HD120958 (B3e), the most luminous member of the group, may be currently driving an RDI shock into the CG12 cloud.

We present the results of a high-resolution X-ray imaging study of the stellar population in the Galactic massive star-forming region RCW49 and its central OB association Westerlund 2. We obtained a 40 ks X-ray image of a 17'x17' field using the Chandra X-ray Observatory and deep NIR images using the Infrared Survey Facility in a concentric 8'3x8'3 region. We detected 468 X-ray sources and identified optical, NIR, and Spitzer Space Telescope MIR counterparts for 379 of them. The unprecedented spatial resolution and sensitivity of the X-ray image, enhanced by optical and infrared imaging data, yielded the following results: (1) The central OB association Westerlund 2 is resolved for the first time in the X-ray band. X-ray emission is detected from all spectroscopically-identified early-type stars in this region. (2) Most (86%) X-ray sources with optical or infrared identifications are cluster members in comparison with a control field in the Galactic Plane. (3) A loose constraint (2--5 kpc) for the distance to RCW49 is derived from the mean X-ray luminosity of T Tauri stars. (4) The cluster X-ray population consists of low-mass pre--main-sequence and early-type stars as obtained from X-ray and NIR photometry. About 30 new OB star candidates are identified. (5) We estimate a cluster radius of 6'--7' based on the X-ray surface number density profiles. (6) A large fraction (90%) of cluster members are identified individually using complimentary X-ray and MIR excess emission. (7) The brightest five X-ray sources, two Wolf-Rayet stars and three O stars, have hard thermal spectra.

Following the Chandra Orion Ultradeep Project (COUP) observation, we have studied the chemical composition of the hot plasma in a sample of 146 X-ray bright pre-main sequence stars in the Orion Nebula Cluster. We report measurements of individual element abundances for a subsample of 86 slightly-absorbed and bright X-ray sources, using low resolution X-ray spectra obtained from the Chandra ACIS instrument. The X-ray emission originates from a plasma with temperatures and elemental abundances very similar to those of active coronae in older stars. A clear pattern of abundances vs. First Ionization Potential (FIP) is evident if solar photospheric abundances are assumed as reference. The results are validated by extensive simulations. The observed abundance distributions are compatible with a single pattern of abundances for all stars, although a weak dependence on flare loop size may be present. The abundance of calcium is the only one which appears to vary substantially between stars, but this quantity is affected by relatively large uncertainties. The ensemble properties of the X-ray bright COUP sources confirm that the iron in the emitting plasma is underabundant with respect to both the solar composition and to the average stellar photospheric values. Comparison of the present plasma abundances with those of the stellar photospheres and those of the gaseous component of the nebula, indicates a good agreement for all the other elements with available measurements, and in particular for the high-FIP elements (Ne, Ar, O, and S) and for the low-FIP element Si. We conclude that there is evidence of a significant chemical fractionation effect only for iron, which appears to be depleted by a factor 1.5--3 with respect to the stellar composition.

We present a statistical analysis of simultaneous optical and X-ray light curves, spanning 600 ks, for 814 pre-main-sequence (PMS) stars in the Orion Nebula Cluster. The aim of this study is to establish the relationship, if any, between the sites of optical and X-ray variability, and thereby to elucidate the origins of X-ray production in PMS stars. In a previous paper we showed that optical and X-ray variability in PMS stars are very rarely time-correlated. Here, using time-averaged variability indicators to examine the joint occurrences of optical and X-ray variability, we confirm that the two forms of variability are not directly causally related. However, a strong and highly statistically significant correlation is found between optical variability and X-ray luminosity. As this correlation is found to be independent of accretion activity, we argue that X-ray production in PMS stars must instead be intimately connected with the presence and strength of optically variable, magnetically active surface regions (i.e. spots) on these stars. Moreover, because X-ray variability and optical variability are rarely time-correlated, we conclude that the sites of X-ray production are not exclusively co-spatial with these regions. We argue that solar-analog coronae, heated by topologically complex fields, can explain these findings.

We report here results from a Chandra ACIS observation of the stellar populations in and around the Messier 17 H II region. The field reveals 886 sources, 771 of which have stellar counterparts in infrared images. In addition to comprehensive tables of X-ray source properties, several results are presented: * The X-ray Luminosity Function is calibrated to that of the Orion Nebula Cluster population to infer a total population of roughly 8000--10,000 stars * About 40% of the ACIS sources are heavily obscured with A_V > 10 mag. Some are concentrated around well-studied star-forming regions but most are distributed across the field. X-ray emission is detected from 64 of the hundreds of Class I protostar candidates that can be identified by near- and mid-infrared colors. These constitute the most likely protostar candidates known in M17. * The spatial distribution of X-ray stars is complex: we find a new embedded cluster, a 2 pc-long arc of young stars along the southwest edge of the M17 H II region, and 0.1 pc substructure within various populations. These structures may indicate that the populations are dynamically young. * All (14/14) of the known O stars but only about half (19/34) of the known B0--B3 stars in the M17 field are detected. Six of these stars exhibit very hard thermal plasma components (kT>4 keV) that may be due to colliding wind binaries. More than 100 candidate new OB stars are found, including 28 X-ray detected intermediate- and high-mass protostar candidates with infrared excesses. * Only a small fraction (perhaps 10%) of X-ray selected high- and intermediate- mass stars exhibit K-band emitting protoplanetary disks, providing further evidence that inner disks evolve very rapidly around more massive stars.

Stable, hydrogen-burning, M dwarf stars comprise about 75% of all stars in the Galaxy. They are extremely long-lived and because they are much smaller in mass than the Sun (between 0.5 and 0.08 MSun), their temperature and stellar luminosity are low and peaked in the red. We have re-examined what is known at present about the potential for a terrestrial planet forming within, or migrating into, the classic liquid-surface-water habitable zone close to an M dwarf star. Observations of protoplanetary disks suggest that planet-building materials are common around M dwarfs, but N-body simulations differ in their estimations of the likelihood of potentially-habitable, wet planets residing within their habitable zones, which are only ~ 1/5 to 1/50 of the width of that for a G star. Particularly in light of the claimed detection of the planets with masses as small as 5.5 and 7.5 MEarth orbiting M stars, there seems no reason to exclude the possibility of terrestrial planets. Tidally locked synchronous rotation within the narrow habitable zone doesn't necessarily lead to atmospheric collapse, and active stellar flaring may not be as much of an evolutionarily disadvantageous factor as has previously been supposed. We conclude that M dwarf stars may indeed be viable hosts for planets on which the origin and evolution of life can occur. A number of planetary processes such as cessation of geothermal activity, or thermal and non-thermal atmospheric loss processes may limit the duration of planetary habitability to periods far shorter than the extreme lifetime of the M dwarf star. Nevertheless, it makes sense to include M dwarf stars in programs that seek to find habitable worlds and evidence of life.

We report the Chandra detection of hard X-ray emission from the Welch ring in W49A, an organized structure of ultra-compact (UC) HII regions containing a dozen nascent early-type stars. Two UC HII regions are associated with hard X-ray emission in a deep Advanced CCD Imaging Spectrometer image exposed for 96.7 ks. One of the two X-ray sources has no near-infrared counterpart and is extended by ~5 arcsec, or ~0.3 pc, at a distance of ~11.4 kpc, which is spatially aligned with the cometary radio continuum emission associated with the UC HII region. The X-ray spectrum of the emission, when fit with a thermal model, indicates a heavily absorbed plasma with extinction of \~5x10^{23}/cm^{2}, temperature of ~7 keV, and X-ray luminosity in the 3.0-8.0 keV band of ~3x10^{33} ergs/s. Both the luminosity and the size of the emission resemble the extended hard emission found in UC HII regions in Sagittarius B2, yet they are smaller by an order of magnitude than the emission found in massive star clusters such as NGC 3603. Three possibilities are discussed for the cause of the hard extended emission in the Welch ring: an ensemble of unresolved point sources, shocked interacting winds of the young O stars, and a wind-blown bubble interacting with ambient cold matter.

The formation of planets is directly linked to the evolution of the circumstellar (CS) disk from which they are born. The dissipation timescales of CS disks are, therefore, of direct astrophysical importance in evaluating the time available for planet formation. We employ Spitzer Space Telescope spectra to complete the CS disk census for the late-type members of the ~8 Myr-old eta Chamaeleontis star cluster. Of the 15 K- and M-type members, eight show excess emission. We find that the presence of a CS disk is anti-correlated with binarity, with all but one disk associated with single stars. With nine single stars in total, about 80% retain a CS disk. Of the six known or suspected close binaries the only CS disk is associated with the primary of RECX 9. No circumbinary disks have been detected. We also find that stars with disks are slow rotators with surface values of specific angular momentum j = 2-15 j_sun. All high specific angular momentum systems with j = 20-30 j_sun are confined to the primary stars of binaries. This provides novel empirical evidence for rotational disk locking and again demonstrates the much shorter disk lifetimes in close binary systems compared to single star systems. We estimate the characteristic mean disk dissipation timescale to be ~5 Myr and ~9 Myr for the binary and single star systems, respectively.

The IC 1396N cometary globule within the large nearby HII region IC 1396 has been observed with the ACIS detector on board the Chandra X-ray Observatory. We detect 117 X-ray sources, of which ~50-60 are likely members of the young open cluster Trumpler~37 dispersed throughout the HII region, and 25 are associated with young stars formed within the globule. Spitzer/2MASS photometry shows the X-ray population is very young: 3 older Class III stars, 16 classical T Tauri stars, 6 protostars including a Class 0/I system. We infer a total T Tauri population of ~30 stars in the globule, including the undetected population, with a star formation efficiency of 1-4%. An elongated source spatial distribution with an age gradient oriented towards the exciting star is discovered in the X-ray population of IC 1396N, supporting similar findings in other cometary globules. The geometric and age distribution is consistent with the RDI model for triggered star formation in CGs by HII region shocks. The inferred velocity of the shock front propagating into the globule is ~0.6km/s. The large number of X-ray-luminous protostars in the globule suggests either an unusually high ratio of Class I/0 vs. Class II/III stars, or a non-standard IMF favoring higher mass stars by the triggering process. The Chandra source associated with the luminous Class 0/I protostar IRAS 21391+5802 is one of the youngest stars ever detected in the X-ray band. We also establish for the first time that the X-ray absorption in protostars arises from the local infalling envelopes rather than ambient molecular cloud material.

We present the first high spatial resolution X-ray study of the massive star forming region NGC 6357, obtained in a 38 ks Chandra/ACIS observation. Inside the brightest constituent of this large HII region complex is the massive open cluster Pismis 24. It contains two of the brightest and bluest stars known, yet remains poorly studied; only a handful of optically bright stellar members have been identified. We investigate the cluster extent and Initial Mass Function and detect ~800 X-ray sources with a limiting sensitivity of 10^{30} ergs s^{-1}; this provides the first reliable probe of the rich intermediate-mass and low-mass population of this massive cluster, increasing the number of known members from optical study by a factor of ~50. The high luminosity end (log L_h[2-8 keV]\ge 30.3 ergs s^{-1}) of the observed X-ray luminosity function in NGC 6357 is clearly consistent with a power law relation as seen in the Orion Nebula Cluster and Cepheus B, yielding the first estimate of NGC 6357's total cluster population, a few times the known Orion population. We investigate the structure of the cluster, finding small-scale substructures superposed on a spherical cluster with 6 pc extent, and discuss its relationship to the nebular morphology. The long-standing Lx - 10^{-7}L_{bol} correlation for O stars is confirmed. Twenty-four candidate O stars and one possible new obscured massive YSO or Wolf-Rayet star are presented. Many cluster members are estimated to be intermediate-mass stars from available infrared photometry (assuming an age of 1 Myr), but only a few exhibit K-band excess. We report the first detection of X-ray emission from an Evaporating Gaseous Globule at the tip of a molecular pillar; this source is likely a B0-B2 protostar.

(abridged:) The XMM-Newton Extended Survey of the Taurus Molecular Cloud (XEST) surveys the most populated ~5 square degrees of the Taurus star formation region, using the XMM-Newton X-ray observatory to study the thermal structure, variability, and long-term evolution of hot plasma, to investigate the magnetic dynamo, and to search for new potential members of the association. Many targets are also studied in the optical, and high-resolution X-ray grating spectroscopy has been obtained for selected bright sources. The X-ray spectra have been coherently analyzed with two different thermal models (2-component thermal model, and a continuous emission measure distribution model). We present overall correlations with fundamental stellar parameters that were derived from the previous literature. A few detections from Chandra observations have been added. The present overview paper introduces the project and provides the basic results from the X-ray analysis of all sources detected in the XEST survey.Comprehensive tables summarize the stellar properties of all targets surveyed. The survey goes deeper than previous X-ray surveys of Taurus by about an order of magnitude and for the first time systematically accesses very faint and strongly absorbed TMC objects. We find a detection rate of 85% and 98% for classical and weak-line T Tau stars (CTTS resp. WTTS), and identify about half of the surveyed protostars and brown dwarfs. Overall, 136 out of 169 surveyed stellar systems are detected. We describe an X-ray luminosity vs. mass correlation, discuss the distribution of X-ray-to-bolometric luminosity ratios, and show evidence for lower X-ray luminosities in CTTS compared to WTTS. Detailed analysis (e.g., variability, rotation-activity relations, influence of accretion on X-rays) will be discussed in a series of accompanying papers.

We present a database of $BVRI$ time-series photometry of the Orion Nebula Cluster obtained with two ground-based telescopes at different longitudes to provide simultaneous coverage with the 13-d Chandra observation of the cluster. The resulting database of simultaneous optical and X-ray light curves for some 800 pre--main-sequence (PMS) stars represents, by a factor of hundreds, the largest synoptic, multi--wavelength-regime, time-series study of young stars to date. This database will permit detailed analyses of the relationship between optical and X-ray variability among a statistically significant ensemble of PMS stars, with the goal of elucidating the origins of PMS X-ray production. In this first paper, we present the optical observations, describe the combined X-ray/optical database, and perform an analysis of time-correlated variability in the optical and X-ray light curves. We identify 40 stars (representing 5% of our study sample) with possible time-correlated optical and X-ray variability. Examples of both positive and negative time-correlations are found, possibly representing X-ray flares and persistent coronal features associated with both cool and hot surface spots (i.e. magnetically active regions and accretion shocks). We also find two possible examples of "white-light" flares coincident with X-ray flares; these may correspond to the impulsive heating phase in solar-analog flares. However, though interesting, these represent unusual cases. More generally, we find very little evidence to suggest a direct causal link between the sources of optical and X-ray variability in PMS stars. The conclusion that accretion is a primary driver of X-ray production in PMS stars is not supported by our findings.

Although the environments of star and planet formation are thermodynamically cold, substantial X-ray emission from 10-100 MK plasmas is present. In low mass pre-main sequence stars, X-rays are produced by violent magnetic reconnection flares. In high mass O stars, they are produced by wind shocks on both stellar and parsec scales. The recent Chandra Orion Ultradeep Project, XMM-Newton Extended Survey of Taurus, and Chandra studies of more distant high-mass star forming regions reveal a wealth of X-ray phenomenology and astrophysics. X-ray flares mostly resemble solar-like magnetic activity from multipolar surface fields, although extreme flares may arise in field lines extending to the protoplanetary disk. Accretion plays a secondary role. Fluorescent iron line emission and absorption in inclined disks demonstrate that X-rays can efficiently illuminate disk material. The consequent ionization of disk gas and irradiation of disk solids addresses a variety of important astrophysical issues of disk dynamics, planet formation, and meteoritics. New observations of massive star forming environments such as M 17, the Carina Nebula and 30 Doradus show remarkably complex X-ray morphologies including the low-mass stellar population, diffuse X-ray flows from blister HII regions, and inhomogeneous superbubbles. X-ray astronomy is thus providing qualitatively new insights into star and planet formation.

We report the detection during the Chandra Orion Ultradeep Project (COUP) of two soft, constant, and faint X-ray sources associated with the Herbig-Haro object HH210. HH210 is located at the tip of the NNE finger of the emission line system bursting out of the BN-KL complex, northwest of the Trapezium cluster in the OMC-1 molecular cloud. Using a recent Halpha image obtained with the ACS imager on board HST, and taking into account the known proper motions of HH210 emission knots, we show that the position of the brightest X-ray source, COUP703, coincides with the emission knot 154-040a of HH210, which is the emission knot of HH210 having the highest tangential velocity (425 km/s). The second X-ray source, COUP704, is located on the complicated emission tail of HH210 close to an emission line filament and has no obvious optical/infrared counterpart. Spectral fitting indicates for both sources a plasma temperature of ~0.8 MK and absorption-corrected X-ray luminosities of about 1E30 erg/s (0.5-2.0 keV). These X-ray sources are well explained by a model invoking a fast-moving, radiative bow shock in a neutral medium with a density of ~12000 cm^{-3}. The X-ray detection of COUP704 therefore reveals, in the complicated HH210 region, an energetic shock not yet identified at other wavelengths.

Cepheus B (Cep B) molecular cloud and a portion of the nearby Cep OB3b OB association, one of the most active regions of star formation within 1 kpc, has been observed with the ACIS detector on board the Chandra X-ray Observatory. We detect 431 X-ray sources, of which 89% are confidently identified as clustered pre-main sequence stars. Two main results are obtained. First, we provide the best census to date for the stellar population of the region. We identify many members of two rich stellar clusters: the lightly obscured Cep OB3b association, and the deeply embedded cluster in Cep B whose existence was previously traced only by a handful of radio sources and T Tauri stars. Second, we find a discrepancy between the X-ray Luminosity Functions of the Cep OB3b and the Orion Nebula Cluster. This may be due to different Initial Mass Functions of two regions (excess of ~0.3 solar mass stars), or different age distributions. Several other results are obtained. A diffuse X-ray component seen in the field is attributed to the integrated emission of unresolved low mass PMS stars. The X-ray emission from HD 217086 (O7n), the principle ionizing source of the region, follows the standard model involving many small shocks in an unmagnetized radiatively accelerated wind. The X-ray source #294 joins a number of similar superflare PMS stars where long magnetic structures may connect the protoplanetary disk to the stellar surface.

We present an X-ray tour of diffuse emission in the 30 Doradus star-forming complex in the Large Magellanic Cloud using high-spatial-resolution X-ray images and spatially-resolved spectra obtained with the Advanced CCD Imaging Spectrometer aboard the Chandra X-ray Observatory. The dominant X-ray feature of the 30 Doradus nebula is the intricate network of diffuse emission generated by interacting stellar winds and supernovae working together to create vast superbubbles filled with hot plasma. We construct maps of the region showing variations in plasma temperature (T = 3--9 million degrees), absorption (N_H = 1--6 x 10^{21} cm^{-2}), and absorption-corrected X-ray surface brightness (S_X = 3--126 x 10^{31} ergs s^{-1} pc^{-2}). Enhanced images reveal the pulsar wind nebula in the composite supernova remnant N157B and the Chandra data show spectral evolution from non-thermal synchrotron emission in the N157B core to a thermal plasma in its outer regions. In a companion paper we show that R136, the central massive star cluster, is resolved at the arcsecond level into almost 100 X-ray sources. Through X-ray studies of 30 Doradus the complete life cycle of such a massive stellar cluster can be revealed.

We have studied the X-ray point source population of the 30 Doradus star-forming complex in the Large Magellanic Cloud using high-spatial-resolution X-ray images and spatially-resolved spectra obtained with the Advanced CCD Imaging Spectrometer (ACIS) aboard the Chandra X-ray Observatory. Here we describe the X-ray sources in a 17' x 17' field centered on R136, the massive star cluster at the center of the main 30 Dor nebula. We detect 20 of the 32 Wolf-Rayet stars in the ACIS field. R136 is resolved at the subarcsecond level into almost 100 X-ray sources, including many typical O3--O5 stars as well as a few bright X-ray sources previously reported. Over two orders of magnitude of scatter in L_X is seen among R136 O stars, suggesting that X-ray emission in the most massive stars depends critically on the details of wind properties and binarity of each system, rather than reflecting the widely-reported characteristic value L_X/L_bol ~ 10^{-7}. Such a canonical ratio may exist for single massive stars in R136, but our data are too shallow to confirm this relationship. Through this and future X-ray studies of 30 Doradus, the complete life cycle of a massive stellar cluster can be revealed.

Following the suggestion of Zuckerman et al. (2001, ApJ, 562, L87), we consider the evidence that 51 Eri (spectral type F0) and GJ 3305 (M0), historically classified as unrelated main sequence stars in the solar neighborhood, are instead a wide physical binary system and members of the young beta Pic moving group (BPMG). The BPMG is the nearest (d < 50 pc) of several groups of young stars with ages around 10 Myr that are kinematically convergent with the Oph-Sco-Cen Association (OSCA), the nearest OB star association. Combining SAAO optical photometry, Hobby-Eberly Telescope high-resolution spectroscopy, Chandra X-ray data, and UCAC2 catalog kinematics, we confirm with high confidence that the system is indeed extremely young. GJ 3305 itself exhibits very strong magnetic activity but has rapidly depleted most of its lithium. The 51 Eri/GJ 3305 system is the westernmost known member of the OSCA, lying 110 pc from the main subgroups. The system is similar to the BPMG wide binary HD 172555/CD -64d1208 and the HD 104237 quintet, suggesting that dynamically fragile multiple systems can survive the turbulent environments of their natal giant molecular cloud complexes, while still being imparted high dispersion velocities. Nearby young systems such as these are excellent targets for evolved circumstellar disk and planetary studies, having stellar ages comparable to that of the late phases of planet formation.

In January 2003, the Chandra Orion Ultradeep Project (COUP) detected about 1400 young stars during a 13.2 day observation of the Orion Nebula Cluster (ONC). This paper studies a well-defined sample of 28 solar-mass COUP sources to characterize the magnetic activity of analogs of the young Sun and thereby to improve understanding of the effects of solar X-rays on the solar nebula during the era of planet formation. We find that active young Suns spend 70% of their time in a characteristic state with relatively constant flux and magnetically confined plasma with temperatures kT_2 = 2.1 * kT_1. During characteristic periods, the 0.5-8 keV X-ray luminosity is about 0.03% of the bolometric luminosity. One or two powerful flares per week with peak luminosities logL_x ~ 30-32 erg/s are typically superposed on this characteristic emission accompanied by heating of the hot plasma component from ~2.4 keV to ~7 keV at the flare peak. The energy distribution of flares superposed on the characteristic emission level follows the relationship dN/dE ~ E^-1.7. The flare rates are consistent with the production of sufficiently energetic protons to spawn a spallogenic origin of some important short-lived radionuclides found in ancient meteorites. The X-rays can ionize gas in the circumstellar disk at a rate of 6 10^-9 ionizations per second at 1 AU from the central star, orders of magnitude above cosmic ray ionization rates. The estimated energetic particle fluences are sufficient to a account for many isotopic anomalies observed in meteoritic inclusions.

We investigate the X-ray and near-infrared emission properties of a sample of pre-main sequence (PMS) stellar systems in the Orion Nebula Cluster (ONC) that display evidence for circumstellar disks (``proplyds'') and optical jets in Hubble Space Telescope (HST) imaging. Approximately 70% of ~140 proplyds were detected as X-ray sources in the Chandra Orion Ultradeep Program observation of the ONC, including ~25% of proplyds that do not display central stars in HST imaging. Many proplyds display near-infrared excesses, suggesting disk accretion is ongoing onto the central, PMS stars. Among circumstellar disks that are detected in absorption in HST imaging, we find that X-ray absorbing column and apparent disk inclination are well correlated, providing insight into the disk scale heights and metal abundances of UV- and X-ray-irradiated protoplanetary disks. Jet sources display some of the largest near-infrared excesses among the X-ray-detected proplyds, suggesting that the origin of the jets is closely related to ongoing, PMS stellar accretion. Some jet sources display double-peaked X-ray spectral energy distributions, with soft components indicative of strong shocks in the jet collimation region. These results support models in which circumstellar disks collimate and/or launch jets from young stellar objects and, furthermore, demonstrate that star-disk-jet interactions may contribute to PMS X-ray emission.

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

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

The Chandra X-ray Observatory is providing remarkable new views of massive star-forming regions, revealing all stages in the life cycle of high-mass stars and their effects on their surroundings. We present a Chandra tour of several high-mass star-forming regions, highlighting physical processes that characterize the life of a cluster of high-mass stars, from deeply-embedded cores too young to have established an HII region to superbubbles so large that they shape our views of galaxies. Along the way we see that X-ray observations reveal hundreds of stellar sources powering great HII region complexes, suffused by both hard and soft diffuse X-ray structures caused by fast O-star winds thermalized in wind-wind collisions or by termination shocks against the surrounding media. Finally, we examine the effects of the deaths of high-mass stars that remained close to their birthplaces, exploding as supernovae within the superbubbles that these clusters created. We present new X-ray results on W51 IRS2E and 30 Doradus and we introduce new data on Trumpler 14 in Carina and the W3 HII region complexes W3 Main and W3(OH).

Astronomical observations of flares from analogs of the early Sun have the potential to give critical insights into the high energy irradiation environment of protoplanetary disks. Solar-mass young stellar objects are significantly more X-ray luminous than the typical low-mass T Tauri star. They undergo frequent strong flaring on a several day time scale. Very powerful flares also occur, but on a longer time frame. The hard X-ray spectrum of these stars become even harder during flaring. The X-rays from these sun-like young stellar objects have the potential to ionize circumstellar material out to large distances. Three specific illustrations are given of the effects of the X-rays: The physics and chemistry of the atmospheres of the inner accretion disks; the ionization level at the disk midplane, important for the viability of the magnetorotational instability; and the nuclear fluence in the irradiation zone just interior to the inner edge of the disk, important in local irradiation scenarios for producing the short-lived radionuclides found in meteorites.

We investigate the spatial distribution of X-ray emitting plasma in a sample of young Orion Nebula Cluster stars by modulation of their X-ray light-curves due to stellar rotation. The study, part of the Chandra Orion Ultradeep Project (COUP), is made possible by the exceptional length of the observation: 10 days of ACIS integration during a time span of 13 days, yielding a total of 1616 detected sources in the 17x17 arcmin field of view. We here focus on a subsample of 233 X-ray-bright stars with known rotational periods. We search for X-ray modulation using the Lomb Normalized Periodogram method. X-ray modulation related to the rotation period is detected in at least 23 stars with periods between 2 and 12 days and relative amplitudes ranging from 20% to 70%. In 16 cases, the X-ray modulation period is similar to the stellar rotation period while in seven cases it is about half that value, possibly due to the presence of X-ray emitting structures at opposite stellar longitudes. These results constitute the largest sample of low mass stars in which X-ray rotational modulation has been observed. The detection of rotational modulation indicates that the X-ray emitting regions are distributed inhomogeneneously in longitude and do not extend to distances significantly larger than the stellar radius. Modulation is observed in stars with saturated activity levels (L_X/L_bol ~ 10^(-3)) showing that saturation is not due to the filling of the stellar surface with X-ray emitting regions.

We have analyzed a number of intense X-ray flares observed in the Chandra Orion Ultradeep Project (COUP), a 13 days observation of the Orion Nebula Cluster (ONC). Analysis of the flare decay allows to determine the size, peak density and magnetic field of the flaring structure. A total of 32 events (the most powerful 1% of COUP flares), have sufficient statistics for the analysis. A broad range of decay times (from 10 to 400 ks) are present in the sample. Peak flare temperatures are often very high, with half of the flares in the sample showing temperatures in excess of 100 MK. Significant sustained heating is present in the majority of the flares. The magnetic structures which are found, are in a number of cases very long, with semi-lengths up to 10^12 cm, implying the presence of magnetic fields of hundreds of G extending to comparable distance from the stellar photosphere. These very large sizes for the flaring structures ($ >> R_*) are not found in more evolved stars, where, almost invariably, the same type of analysis results in structures with L <= R_*. As the majority of young stars in the ONC are surrounded by disks, we speculate that the large magnetic structures which confine the flaring plasma are actually the same type of structures which channel the plasma in the magnetospheric accretion paradigm, connecting the star's photosphere with the accretion disk.

We study the relation between age and magnetic activity in late-type pre-main sequence (PMS) stars, for the first time using mass-stratified subsamples. The effort is based on the Chandra Orion Ultradeep Project (COUP) which provides very sensitive and homogenous X-ray data on a uniquely large sample of 481 optically well-characterized low-extinction low-mass members of the Orion Nebula Cluster, for which individual stellar masses and ages could be determined. More than 98 percent of the stars in this sample are detected as X-ray sources. Within the PMS phase for stellar ages in the range $\sim 0.1-10$ Myr, we establish a mild decay in activity with stellar age $\tau$ roughly as $L_{\rm X} \propto \tau^{-1/3}$. On longer timescales, when the Orion stars are compared to main sequence stars, the X-ray luminosity decay law for stars in the $0.5 < M < 1.2$ M$_\odot$ mass range is more rapid with $L_{\rm X} \propto \tau^{-0.75}$ over the wide range of ages $5 < \log \tau < 9.5$ yr. The magnetic activity history for M stars with masses $0.1 < M < 0.4 M_\odot$ is distinctly different. Only a mild decrease in X-ray luminosity, and even a mild increase in $L_{\rm X}/L_{\rm bol}$ and $F_{\rm X}$, is seen over the 1-100 Myr range, though the X-ray emission does decay over long timescales on the main sequence. Together with COUP results on the absence of a rotation-activity relation in Orion stars, we find that the activity-age decay is strong across the entire history of solar-type stars but is not attributable to rotational deceleration during the early epochs. A combination of tachocline and distributed convective dynamos may be operative in young solar-type stars. The results for the lowest mass stars are most easily understood by the dominance of convective dynamos during both the PMS and main sequence phases.

We use the sensitive X-ray data from the Chandra Orion Ultradeep Project (COUP) to study the X-ray properties of 34 spectroscopically-identified brown dwarfs with near-infrared spectral types between M6 and M9 in the core of the Orion Nebula Cluster. Nine of the 34 objects are clearly detected as X-ray sources. The apparently low detection rate is in many cases related to the substantial extinction of these brown dwarfs; considering only the BDs with $A_V \leq 5$ mag, nearly half of the objects (7 out of 16) are detected in X-rays. Our 10-day long X-ray lightcurves of these objects exhibit strong variability, including numerous flares. While one of the objects was only detected during a short flare, a statistical analysis of the lightcurves provides evidence for continuous (`quiescent') emission in addition to flares for all other objects. Of these, the $\sim$ M9 brown dwarf COUP 1255 = HC 212 is one of the coolest known objects with a clear detection of quiescent X-ray emission. The X-ray properties (spectra, fractional X-ray luminosities, flare rates) of these young brown dwarfs are similar to those of the low-mass stars in the ONC, and thus there is no evidence for changes in the magnetic activity around the stellar/substellar boundary, which lies at $\sim$ M6 for ONC sources. Since the X-ray properties of the young brown dwarfs are also similar to those of M6--M9 field stars, the key to the magnetic activity in very cool objects seems to be the effective temperature, which determines the degree of ionization in the atmosphere.

The X-ray properties of twenty ~1 Myr old O, B, and A stars of the Orion Trapezium are examined with data from the Chandra Orion Ultradeep Project (COUP). On the basis of simple theories for X-ray emission, we define two classes separated at spectral type B4: hotter stars have strong winds that may give rise to X-ray emission in small- or large-scale wind shocks, and cooler stars that should be X-ray dark due to their weaker winds and absence of outer convection zones where dynamos can generate magnetic fields. Only two of the massive stars show exclusively the constant soft-spectrum emission expected from the standard model for X-ray emission from hot stars involving many small shocks in an unmagnetized radiatively accelerated wind. Most of the other massive O7-B3 stars exhibit some combination of soft-spectrum wind emission, hard-spectrum flaring, and/or rotational modulation indicating large-scale inhomogeneity, suggesting magnetic confinement of winds with large-scale shocks. Most of the stars in the weak-wind class exhibit X-ray flares and have luminosities that are consistent with magnetic activity from known or unseen low-mass companions. All non-detections belong to the weak-wind class.

The Chandra Orion Ultradeep project (COUP) observation described in a companion paper by Getman et al. provides an exceptionally deep X-ray survey of the Orion Nebula Cluster and associated embedded young stellar objects. Membership of the region is important for studies of the stellar IMF, cluster dynamics, and star formation. The COUP study detected 1616 X-ray sources. In this study we confirm cloud membership for 1315 stars, identify 16 probable foreground field stars having optical counterparts with discrepant proper motions, and classify the remaining 285 X-ray sources, of which 51 are lightly and 234 heavily obscured. The 51 lightly obscured sources without known counterparts fall into three groups. (i) Sixteen are likely new members of the Orion Nebula Cluster. (ii) Two with unusually soft and non-flaring X-ray emission appear to be associated with nebular shocks, and may be new examples of X-rays produced at the bow shocks of Herbig-Haro outflows. (iii) The remaining thirty three are very weak uncertain sources, possibly spurious. Out of 234 heavily absorbed sources without optical or near-infrared counterparts 75 COUP sources are likely new embedded cloud members (with membership for 42 confirmed by powerful X-ray flares), and the remaining 159 are likely extragalactic AGN seen through the molecular cloud, as argued by a careful simulation of the extragalactic background population. Finally, a few new binary companions to Orion stars may have been found, but most cases of proximate COUP sources can be attributed to chance superpositions in this crowded field.

We present a study of the X-ray sources detected in the vicinity of two density peaks in the Orion Molecular Cloud 1 (OMC-1) behind the Orion Nebula Cluster (ONC), as seen in the exceptionally deep (~10 days) exposure of the Chandra Orion Ultradeep Project (COUP). We focus on a 40"*50" region around the Becklin-Neugebauer object and Kleinmann-Low nebula (collectively BN-KL) and a 60"*75" region around OMC-1S, a secondary star-forming peak some 90" south of BN-KL. Forty-three and sixty X-ray sources were detected in BN-KL and OMC-1S, respectively, of which half and one-third, respectively, were found to be foreground members of the ONC, while the remaining sources are identified as obscured X-ray sources with column densities 22<log(NH)/cm^{-2})<24. All but 5 and 18 of these sources have near-infrared stellar counterparts, and 22 of these appear to be young stellar objects (YSOs). X-ray sources are found close to four luminous mid-IR sources namely BN, IRc3-i2, IRc2-C, and Source n; their X-ray variability and spectral properties are typical of coronal activity in low-mass stars rather than wind emission from massive stars, suggesting that the X-ray emission may be arising from companions. The X-ray light curve of the X-ray source close to BN shows a periodicity of ~8.3 days and from an X-ray image deconvolution of the region around BN, we conclude that either BN itself or a low mass companion with a projected separation of ~200 AU was detected. On the other hand, no emission is seen from the bright radio Source I, held by some to be the main source of luminosity in BN-KL. In OMC-1S, Chandra unveils a new subcluster of seven YSOs without infrared counterparts (abridged).

The study of magnetic activity of pre-main sequence (PMS) stars through their X-ray emission is entering a mature phase. We review here two recent developments. First, we present some early findings from the Chandra Orion Ultradeep Project (COUP) relating to the age-activity-rotation relations of PMS stars. COUP is a nearly continuous exposure of the Orion Nebula Cluster over 13 days and gives the most comprehensive view to date of X-ray emission from the PMS population. We find that the activity-rotation connection that dominates main sequence activity behavior is absent in PMS stars, although the activity-age relation is present as well as links to stellar mass and volume. Second, we discuss investigations of effects of PMS X-rays on their cold protoplanetary disks. It is possible that the the ionization, heating and spallogenic effects of high energy photons and particles produced in PMS stellar flares have significant effects on disk dynamics, heating and chemistry. X-ray ionization may induce MHD turbulence in the disk which, in turn, may regulate disk evolution and planet formation processes such as Jovian planet migration.

Sensitive imaging X-ray observations of young stellar clusters (YSCs, ages less than about 10 Myr) are valuable tools for the acquisition of an unbiased census of cluster members needed for Initial Mass Function (IMF) studies. Several dozen YSCs, both nearby and across the Galactic disk, have been observed with the Chandra and XMM-Newton satellites, detecting >10,000 low-mass cluster members. Many of these samples should be nearly complete down to about 1 Mo. An important additional benefit is that the YSC X-ray luminosity function appears to be universal with a lognormal shape, providing a new standard candle for measurement of YSC distances.

We present the result of a systematic search for the iron Kalpha fluorescent line at 6.4 keV among 1616 X-ray sources detected by ultra-deep Chandra observations of the Orion Nebula Cluster and the obscured Orion Molecular Cloud 1 population as part of the Chandra Orion Ultra-deep Project (COUP). Seven sources are identified to have an excess emission at 6.4 keV among 127 control sample sources with significant counts in the 6.0-9.0 keV band. These seven sources are young stellar objects (YSOs) characterized by intense flare-like flux variations, thermal spectra, and near-infrared (NIR) counterparts. The observed equivalent widths of the line cannot be attributed to the fluorescence by interstellar or circumstellar matter along the line of sight. The X-ray spectral fits and NIR colors of the 6.4 keV sources show that these sources have X-ray absorption of > 1x10^22 cm^(-2) and NIR excess emission, which is not expected when the fluorescence occurs at the stellar photosphere. We therefore conclude that the iron fluorescent line of YSOs arises from reflection off of circumstellar disks, which are irradiated by the hard X-ray continuum emission of magnetic reconnection flares.

We present a description of the data reduction methods and the derived catalog of more than 1600 X-ray point sources from the exceptionally deep January 2003 Chandra X-ray Observatory observation of the Orion Nebula Cluster and embedded populations around OMC-1. The observation was obtained with Chandra's Advanced CCD Imaging Spectrometer (ACIS) and has been nicknamed the Chandra Orion Ultradeep Project (COUP). With an 838 ks exposure made over a continuous period of 13.2 days, the COUP observation provides the most uniform and comprehensive dataset on the X-ray emission of normal stars ever obtained in the history of X-ray astronomy.

We examine the X-ray spectra and variability of the sample of X-ray sources with L_X = 10^{31}-10^{33} erg s^{-1} identified within the inner 9' of the Galaxy. Very few of the sources exhibit intra-day or inter-month variations. We find that the spectra of the point sources near the Galactic center are very hard between 2--8 keV, even after accounting for absorption. When modeled as power laws the median photon index is Gamma=0.7, while when modeled as thermal plasma we can only obtain lower limits to the temperature of kT>8 keV. The combined spectra of the point sources is similarly hard, with a photon index of Gamma=0.8. Strong line emission is observed from low-ionization, He-like, and H-like Fe, both in the average spectra and in the brightest individual sources. The line ratios of the highly-ionized Fe in the average spectra are consistent with emission from a plasma in thermal equilibrium. This line emission is observed whether average spectra are examined as a function of the count rate from the source, or as a function of the hardness ratios of individual sources. This suggests that the hardness of the spectra may in fact to due local absorption that partially-covers the X-ray emitting regions in the Galactic center systems. We suggest that most of these sources are intermediate polars, which (1) often exhibit hard spectra with prominent Fe lines, (2) rarely exhibit either flares on short time scales or changes in their mean X-ray flux on long time scales, and (3) are the most numerous hard X-ray sources with comparable luminosities in the Galaxy.

Sensitive X-ray imaging surveys provide a new and effective tool to establish the census of pre-main sequence (PMS) stars in nearby young stellar clusters. We report here a deep Chandra X-ray Observatory observation of PMS stars in the Chamaeleon I North cloud, achieving a limiting luminosity of log L_t ~ 27 erg/s (0.5-8 keV band) in a 0.8x0.8 pc region. Of the 107 X-ray sources, 37 are associated with Galactic stars of which 27 are previously recognized cloud members. These include three PMS brown dwarfs; the protostellar brown dwarf ISO 192 has a particularly high level of magnetic activity. Followup optical photometry and spectroscopy establishes that 9-10 of the Chandra sources are probably magnetically active background stars. Several previously proposed cloud members are also inferred to be interlopers due to the absence of X-ray emission at the level expected from the L_t - K correlation. No new X-ray discovered stars were confidently found despite the high sensitivity. From these findings, we argue that the sample of 27 PMS cloud members in the Chandra field is uncontaminated and complete down to K = 12 or M ~ 0.1 Mo. The initial mass function (IMF) derived from our sample is deficient in 0.1-0.3 Mo stars compared to the IMF of the rich Orion Nebula Cluster and other Galactic populations. We can not discriminate whether this is due to different star formation processes, mass segregation, or dynamical ejection of lower mass stars.

We present Chandra/ACIS images of several high-mass star-forming regions. The massive stellar clusters powering these HII regions are resolved at the arcsecond level into hundreds of stellar sources, similar to those seen in closer young stellar clusters. However, we also detect diffuse X-ray emission on parsec scales that is spatially and spectrally distinct from the point source population. For nearby regions (e.g. M17 and Rosette) the emission is soft, with plasma temperatures less than 10 million degrees, in contrast to what is seen in more distant complexes (e.g. RCW49, W51). This extended emission most likely arises from the fast O-star winds thermalized either by wind-wind collisions or by a termination shock against the surrounding media. We have established that only a small portion of the wind energy and mass appears in the observed diffuse X-ray plasma; in the blister HII regions, we suspect that most of it flows without cooling into the low-density interstellar medium through blow-outs or fissures in the surrounding neutral material. These data provide compelling observational evidence that strong wind shocks are present in HII regions.

We examine the spectrum of diffuse emission detected in the 17' by 17' field around Sgr A* during 625 ks of Chandra observations. The spectrum exhibits He-like and H-like lines from Si, S, Ar, Ca, and Fe, that are consistent with originating in a two-temperature plasma, as well as a prominent low-ionization Fe line. The cooler, kT=0.8 keV plasma differs in surface brightness across the image by a factor of 9. This soft plasma is probably heated by supernovae. The radiative cooling rate of the plasma within the inner 20 pc of the Galaxy could be balanced by 1% of the kinetic energy of one supernova every 300,000 y. The hotter, kT=8 keV component is more spatially uniform, ranging over a factor of 2 in surface brightness. The intensity of the hard plasma is correlated with that of the soft, but they are probably only indirectly related, because supernova remnants are not observed to produce thermal plasma hotter than kT=3 keV. Moreover, a kT=8 keV plasma is too hot to be bound to the Galactic center, and therefore would form a slow wind or fountain of plasma. The energy required to sustain such a freely-expanding plasma within the inner 20 pc of the Galaxy is ~10^40 erg/s, which corresponds to the entire kinetic energy of one supernova every 3000 y. This rate is unreasonably high. However, alternative explanations for the kT=8 keV diffuse emission are equally unsatisfying. We are left to conclude that either the diffuse emission is heated by an unanticipated source of energy, or that a population of faint (< 10^31 erg/s), hard X-ray sources that are a factor of 10 more numerous than CVs remains to be discovered. (Abridged)

The extremely sensitive Chandra Deep Field-North (CDF-N) pencil-beam X-ray survey is used to identify and characterize the X-ray emission from old high-latitude main sequence Galactic stars. Our principal goal is to investigate the expected long-term decay of magnetic activity of late-type stars due to the gradual spindown of stellar rotation from a magnetized stellar wind. Eleven X-ray sources constitute a well-defined sample of 2 G, 2 K0-K4, and 7 M2-M5 stars with median distance around 300 pc. X-ray luminosities are typically log Lx ~ 27 erg/s and is dominated by flares rather than quiescent coronal emission. Models of the population indicates that the CDF-N stars are the most magnetically active old disk stars. A substantial decline in X-ray luminosities over the 1<t<11 Gyr age interval is required. This is the first demonstration that the coronal and flaring components of stellar magnetic activity -- and presumably the interior magnetic dynamos responsible for the reconnecting fields at the stellar surface -- exhibit long-term decay over the age of the Galactic disk. The model that best fits the magnitudes, spectral types and X-ray luminosities of the sample has Lx ~ 1/t^2 erg/s which is faster than the 1/t decay rate predicted from widely accepted rotational spindown rates and X-ray-activity relations.

LkHA312 has been observed serendipitously with the ACIS-I detector on board Chandra with 26h continuous exposure. This H_alpha emission line star belongs to the star-forming region M78 (NGC2068). From the optical and NIR data, we show that it is a pre-main sequence (PMS) low-mass star with a weak NIR excess. This genuine T Tauri star displayed an X-ray flare with an unusual long rise phase (~8h). The X-ray emission was nearly constant during the first 18h of the observation, and then increased by a factor of 13 during a fast rise phase (~2h), and reached a factor of 16 above the quiescent X-ray level at the end of a gradual phase (~6h) showing a slower rise. To our knowledge this flare, with \~0.4-~0.5 cts/s, has the highest count rate observed so far with Chandra from a PMS low-mass star. By chance, the source position, 8.2' off-axis, protected this observation from pile-up. We make a spectral analysis of the X-ray emission versus time, showing that the plasma temperature of the quiescent phase and the flare peak reaches 29MK and 88MK, respectively. The quiescent and flare luminosities in the energy range 0.5--8keV corrected from absorption (N_H~1.7E21 cm^{-2}) are 6E30erg/s and ~1E32erg/s, respectively. The ratio of the quiescent X-ray luminosity on the LkHA312 bolometric luminosity is very high with log(L_X/L_bol)= -2.9, implying that the corona of LkHA312 reached the `saturation' level. The X-ray luminosity of the flare peak reaches ~2% of the stellar bolometric luminosity. The different phases of this flare are finally discussed in the framework of solar flares, which leads to the magnetic loop height from 3.1E10 to 1E11 cm (0.2-0.5 R*, i.e., 0.5-1.3 R_sun).

We further investigate subpixel event repositioning (SER) algorithms in application to Chandra X-ray Observatory (CXO) CCD imaging. SER algorithms have been applied to backside illuminated (BI) Advanced CCD Imaging Spectrometer (ACIS) devices, and demonstrate spatial resolution improvements in Chandra/ACIS observations. Here a new SER algorithm that is charge split dependent is added to the SER family. We describe the application of SER algorithms to frontside illuminated (FI) ACIS devices. The results of SER for FI CCDs are compared with those obtained from SER techniques applied to BI CCD event data. Both simulated data and Chandra/ACIS observations of the Orion Nebular Cluster were used to test and evaluate the achievement of the various SER techniques.

Despite centuries of close association, statistics and astronomy are surprisingly distant today. Most observational astronomical research relies on an inadequate toolbox of methodological tools. Yet the needs are substantial: astronomy encounters sophisticated problems involving sampling theory, survival analysis, multivariate classification and analysis, time series analysis, wavelet analysis, spatial point processes, nonlinear regression, bootstrap resampling and model selection. We review the recent resurgence of astrostatistical research, and outline new challenges raised by the emerging Virtual Observatory. Our essay ends with a list of research challenges and infrastructure for astrostatistics in the coming decade.

We report new results in our continuing study of the unique compact (1 pc extent), nearby (d = 97 pc), young (t = 9 Myr) stellar cluster dominated by the B9 star eta Chamaeleontis. An optical photometric survey spanning 1.3 x 1.3 pc adds two M5-M5.5 weak-lined T Tauri (WTT) stars to the cluster inventory which is likely to be significantly complete for primaries with masses > 0.15 M_sun. The cluster now consists of 17 primaries and approximately 9 secondaries lying within 100 AU of their primaries. The apparent distribution of 9:7:1 single:binary:triple systems shows 2-4 x higher multiplicity than in the field main sequence stars, and is comparable to that seen in other pre-main sequence (PMS) populations. The initial mass function (IMF) is consistent with that of rich young clusters and field stars. By extending the cluster IMF to lower masses, we predict 10-14 additional low mass stars with 0.08 < M < 0.15 M-sun and 10-15 brown dwarfs with 0.025 < M < 0.08 M_sun remain to be discovered. The eta Cha cluster extends the established stellar density and richness relationship for young open clusters. The radial distribution of stars is consistent with an isothermal sphere, but mass segregation is present with > 50 percent of the stellar mass residing in the inner 6 arcmin (0.17 pc). Considering that the eta Cha cluster is sparse, diffuse and young, the cluster may be an ideal laboratory for distinguishing between mass segregation that is primordial in nature, or arising from dynamical interaction processes.

We present the outcomes of a Chandra X-ray Observatory snapshot study of five nearby Herbig Ae/Be (HAeBe) stars which are kinematically linked with the Oph-Sco-Cen Association (OSCA). Optical photometric and spectroscopic followup was conducted for the HD 104237 field. The principal result is the discovery of a compact group of pre-main sequence (PMS) stars associated with HD 104237 and its codistant, comoving B9 neighbor epsilon Chamaeleontis AB. We name the group after the most massive member. The group has five confirmed stellar systems ranging from spectral type B9-M5, including a remarkably high degree of multiplicity for HD 104237 itself. The HD 104237 system is at least a quintet with four low mass PMS companions in nonhierarchical orbits within a projected separation of 1500 AU of the HAeBe primary. Two of the low-mass members of the group are actively accreting classical T Tauri stars. The Chandra observations also increase the census of companions for two of the other four HAeBe stars, HD 141569 and HD 150193, and identify several additional new members of the OSCA. We discuss this work in light of several theoretical issues: the origin of X-rays from HAeBe stars; the uneventful dynamical history of the high-multiplicity HD 104237 system; and the origin of the epsilon Cha group and other OSCA outlying groups in the context of turbulent giant molecular clouds. Together with the similar eta Cha cluster, we paint a portrait of sparse stellar clusters dominated by intermediate-mass stars 5-10 Myr after their formation.

(abridged) We present a catalog of 2357 point sources detected during 590 ks of Chandra observations of the 17-by-17 arcminute field around Sgr A*. This field encompasses a physical area of 40 by 40 pc at a distance of 8 kpc. The completeness limit of the sample at the Galactic center is 10^{31} erg s^{-1} (2.0--8.0 keV), while the detection limit is an order of magnitude lower. The 281 sources detected below 1.5 keV are mainly in the foreground of the Galactic center, while comparisons to the Chandra deep fields at high Galactic latitudes suggest that only about 100 of the observed sources are background AGN. The surface density of absorbed sources (not detected below 1.5 keV) falls off as 1/theta away from Sgr A*, in agreement with the distribution of stars in infrared surveys. Point sources brighter than our completeness limit produce 10% of the flux previously attributed to diffuse emission. The log(N)-log(S) distribution of the Galactic center sources is extremely steep (power-law slope alpha = 1.7). If this distribution extends down to a flux of 10^{-17} erg cm^{-1} s^{-1} (10^{29} erg s^{-1} at 8 kpc, 2.0--8.0 keV) with the same slope, then point sources would account for all of the previously reported diffuse emission. Therefore, the 2.0--8.0 keV luminosity distribution must flatten between 10^{29} - 10^{31} erg s^{-1}. Finally, the spectra of more than half of the Galactic center sources are very hard, and can be described by a power law ($E^{-Gamma}) with photon index Gamma < 1. Such hard spectra have been seen previously only from magnetically accreting white dwarfs and wind-accreting neutron stars, suggesting that there are large numbers of these systems in our field.

We present a comparison of the gas and dust properties of the dense interstellar matter in six nearby star-forming regions (d<500 pc): rho Oph, Cha I, R CrA, IC 348, NGC 1333, and Orion. We measure from Chandra and XMM-Newton observations the X-ray absorption toward pre-main sequence stars (PMS) without accretion disks (i.e., Class III sources) to obtain the total hydrogen column density N_{H,X}. For these sources we take from the literature the corresponding dust extinction in the near-infrared, A_J, or when unavailable we derive it from SED fitting using the available DENIS, 2MASS, ISOCAM and other data. We then compare N_{H,X} and A_J for each object, up to unprecedently high extinction. For the rho Oph dark cloud with a relatively large sample of 20 bona-fide Class III sources, we probe the extinction up to A_J <~ 14 (A_V <~ 45), and find a best-fit linear relation N_{H,X}/A_J = 5.6 (+/- 0.4)x10^{21} cm^{-2} mag^{-1}, adopting standard ISM abundances. The other regions reveal a large dispersion in the N_{H,X}/A_J ratio for each source but for lack of adequate IR data these studies remain limited to moderate extinctions (A_J <~ 1.5 or A_V <~5). For rho Oph, the N_{H,X}/A_J ratio is significantly lower (>~2 sigma) than the galactic value, derived using the standard extinction curve (R_V = 3.1). This result is consistent with the recent downwards revision of the metallicity of the Sun and stars in the solar vicinity. We find that the rho Oph dense cloud has the same metallicity than the local ISM when assuming that the galactic gas-to-dust ratio remains unchanged. The difference between galactic and local values of the gas-to-dust ratio can thus be attributed entirely to a difference in metallicity.

We present the first high-spatial-resolution X-ray images of two high-mass star forming regions, the Omega Nebula (M 17) and the Rosette Nebula (NGC 2237--2246), obtained with the Chandra X-ray Observatory Advanced CCD Imaging Spectrometer (ACIS) instrument. The massive clusters powering these H II regions are resolved at the arcsecond level into >900 (M 17) and >300 (Rosette) stellar sources similar to those seen in closer young stellar clusters. However, we also detect soft diffuse X-ray emission on parsec scales that is spatially and spectrally distinct from the point source population. The diffuse emission has luminosity L_x ~ 3.4e33 ergs/s in M~17 with plasma energy components at kT ~0.13 and ~0.6 keV (1.5 and 7 MK), while in Rosette it has L_x \~6e32 ergs/s with plasma energy components at kT ~0.06 and ~0.8 keV (0.7 and 9 MK). This extended emission most likely arises from the fast O-star winds thermalized either by wind-wind collisions or by a termination shock against the surrounding media. We establish that only a small portion of the wind energy and mass appears in the observed diffuse X-ray plasma; in these blister H II regions, we suspect that most of it flows without cooling into the low-density interstellar medium. These data provide compelling observational evidence that strong wind shocks are present in H II regions.

We report the X-ray detection of the brown dwarf (BD) companion TWA 5B in a $\simeq 12$ Myr old pre-main sequence binary system. We clearly resolve the faint companion (35 photons) separated from the X-ray luminous primary by 2 arcsec in a {\it Chandra} ACIS image. TWA 5B shows a soft X-ray spectrum with a low plasma temperature of only 0.3 keV and a constant flux during the 3 hour observation, of which the characteristics are commonly seen in the solar corona. The X-ray luminosity is 4$\times10^{27}$ erg s$^{-1}$ (0.1--10 keV band) or $\log L_X/L_{bol} = -3.4$. Comparing these properties to both younger and older BDs, we discuss the evolution of the X-ray emission in BDs. During their first few Myr, they exhibit high levels of X-ray activity as seen in higher mass pre-main sequence stars. The level in TWA 5B is still high at $t \simeq 12$ Myr in $\log L_X/L_{bol}$ while $kT$ has already substantially cooled.

A recent Chandra/ACIS observation of the Orion Nebula Cluster detected 1075 sources (Feigelson et al. 2002), providing a uniquely large and well-defined sample to study the dependence of magnetic activity on bulk properties for stars descending the Hayashi tracks. The following results are obtained: (1) X-ray luminosities L_t in the 0.5-8 keV band are strongly correlated with bolometric luminosity with <log L_t/L_bol> = -3.8 for stars with masses 0.7<M<2 Mo, an order of magnitude below the main sequence saturation level; (2) the X-ray emission drops rapidly below this level in some or all stars with 2<M<3 Mo; (3) the presence or absence of infrared circumstellar disks has no apparent relation to X-ray levels; and (4) X-ray luminosities exhibit a slight rise as rotational periods increase from 0.4 to 20 days. This last finding stands in dramatic contrast to the strong anticorrelation between X-rays and period seen in main sequence stars. The absence of a strong X-ray/rotation relationship in PMS stars, and particularly the high X-ray values seen in some very slowly rotating stars, is a clear indication that the mechanisms of magnetic field generation differ from those operating in main sequence stars. The most promising possibility is a turbulent dynamo distributed throughout the deep convection zone, but other models such as alpha-Omega dynamo with `supersaturation' or relic core fields are not immediately excluded. The drop in magnetic activity in intermediate-mass stars may reflect the presence of a significant radiative core. The evidence does not support X-ray production in large-scale star-disk magnetic fields.

We have analyzed JHKL observations of the stellar population of the ~9 Myr-old eta Chamaeleontis cluster. Using infrared (IR) colour-colour and colour-excess diagrams, we find the fraction of stellar systems with near-IR excess emission is 0.60 pm 0.13 (2_sigma). This results implies considerably longer disc lifetimes than found in some recent studies of other young stellar clusters. For the classical T Tauri (CTT) and weak-lined T Tauri (WTT) star population, we also find a strong correlation between the IR excess and H_alpha emission. The IR excesses of these stars indicate a wide range of star-disc activity; from a CTT star showing high levels of accretion, to CTT - WTT transition objects with evidence for some on-going accretion, and WTT stars with weak or absent IR excesses. Of the 15 known cluster members, 4 stars with IR excesses delta(K-L) > 0.4 mag are likely experiencing on-going accretion owing to strong or variable optical emission. The resulting accretion fraction (0.27 pm 0.13; 2_sigma) shows that the accretion phase, in addition to the discs themselves, can endure for at least ~10 Myr.

Star-forming regions have been the targets of X-ray observations since the dawn of satellite X-ray astronomy. The increase in sensitivity and/or spatial resolution offered by XMM-Newton and Chandra allows a dramatic improvement, both qualitative and quantitative, on our knowledge of high-energy phenomena in these regions and the underlying physical processes. We summarize here some recent developments: the Orion Nebula Cluster and its 1000+ stellar X-ray sources; Herbig-Haro objects and their high-speed shocks; protostars, brown dwarfs and their unusual magnetic activity; and the discovery of diffuse X-ray emission from HII regions, presumably related to strong winds from massive stars. The role that future X-ray missions may play in the field is already starting to be visible.

NGC 1333, a highly active star formation region within the Perseus molecular cloud complex, has been observed with the ACIS-I detector on board the Chandra X-ray Observatory. In our image with a sensitivity limit of ~ 1e28 erg/s, we detect 127 X-ray sources, of which, 95 are known cluster members. The X-ray luminosity function of the discovered YSO population spans a range of 1e28-3.2e31 erg/s in the 0.5-8 keV band. Comparison with K-band source counts indicates that we detect all of the known cluster members with K < 12, and about half of members with K > 12. We detect seven X-ray emitting YSOs which drive optically visible jets as well as one deeply embedded object that has not been detected in near-IR observations. The presence or absence of an outflows does not appear to produce any difference in X-ray properties of YSOs. We also find no systematic differences in X-ray luminosity distributions between two complete subsamples of CTTs and WTTs. That suggests that there is no difference in the astrophysical mechanism of WTT and CTT X-ray emission production. Additional results include: the X-ray emission from two late-B stars which illuminate the reflection nebula originates from unresolved late-type companions; two T Tauri stars are discovered as previously unknown components of visual binaries; and the X-ray counterpart of SVS 16 has the column density much lower than that expected from near-IR photometry and thus its X-ray luminosity is not anomalously high, as has been previously suggested.

The Orion Nebula Cluster and the molecular cloud in its vicinity have been observed with the ACIS-I detector on board the Chandra X-ray Observatory with 23 hours exposure. We detect 1075 X-ray sources: 91% are spatially associated with known stellar members of the cluster, and 7% are newly identified deeply embedded cloud members. This provides the largest X-ray study of a pre-main sequence stellar population. We examine here the X-ray properties of Orion young stars as a function of mass. Results include: (a) the discovery of rapid variability in the O9.5 31 M_o star \theta^2A Ori, and several early B stars, inconsistent with the standard model of X-ray production in small wind shocks; (b) support for the hypothesis that intermediate-mass mid-B through A type stars do not themselves produce significant X-ray emission; (c) confirmation that low-mass G- through M-type T Tauri stars exhibit powerful flaring but typically at luminosities considerably below the `saturation' level; (d) confirmation that the presence or absence of a circumstellar disk has no discernable effect on X-ray emission; (e) evidence that T Tauri plasma temperatures are often very high with T >= 100 MK, even when luminosities are modest and flaring is not evident; and (f) detection of the largest sample of pre-main sequence very low mass objects showing high flaring levels and a decline in magnetic activity as they evolve into L- and T-type brown dwarfs.

We report the discovery of soft (0.25--2 keV) x-ray emission from the Galilean satellites Io and Europa, probably Ganymede, and from the Io Plasma Torus (IPT). Bombardment by energetic (>10 keV) H, O, and S ions from the region of the IPT seems the likely source of the x-ray emission from the Galilean satellites. According to our estimates, fluorescent x-ray emission excited by solar x-rays, even during flares from the active Sun, charge-exchange processes, previously invoked to explain Jupiter's x-ray aurora and cometary x-ray emission, and ion stripping by dust grains fail to account for the observed emission. On the other hand, bremsstrahlung emission of soft X-rays from non-thermal electrons in the few hundred to few thousand eV range may account for a substantial fraction of the observed x-ray flux from the IPT.

To address the role of energetic processes in the solar nebula, we provide a detailed characterization of magnetic flaring in stellar analogs of the pre-main sequence Sun based on 23 hours observations of 43 analogs of the young Sun in the Orion Nebula Cluster obtained with the Chandra X-ray Observatory. We find the X-ray luminosities are strongly elevated over main sequence levels with average <log L_x> = 30.3 erg/s and <log L_x/L_*> = -3.9 (0.5-8 keV), and dozens of flares are present. Analogs of the <= 1 My old pre-main sequence Sun exhibit X-ray flares 10^{1.5} times more powerful and 10^{2.5} times more frequent than the most powerful flares seen on the contemporary Sun. Extrapolating the solar relationship between X-ray luminosity and proton fluence, we infer that the young Sun exhibited a 10^5-fold enhancement in energetic protons compared to contemporary levels. Unless the flare geometries are unfavorable, this inferred proton flux on the disk is sufficient to produce the observed meteoritic abundances of several important short-lived radioactive isotopes. Our study thus strengthens the astronomical foundation for local proton spallation models of isotopic anomalies in carbonaceous chondritic meteorites. The radiation, particles and shocks produced by the magnetic reconnection flares seen with Chandra may also have flash melted meteoritic chondrules and produced excess 21-Ne seen in meteoritic grains.

A limited-area survey of the eta Chamaeleontis cluster has identified 2 new late-type members. The more significant of these is ECHA J0843.3-7905 (= IRAS F08450-7854), a slowly-rotating (P = 12 d) M2 classical T Tauri (CTT) star with a spectrum dominated by Balmer emission. At a distance of 97 pc and cluster age of approx 9 Myr, the star is a nearby rare example of an 'old' CTT star and promises to be a rewarding laboratory for the study of disk structure and evolution in pre-main sequence (PMS) stars. The other new member is the M4 weak-lined T Tauri (WTT) star ECHA J0841.5-7853, which is the lowest mass (M = 0.2 M_sun) primary known in the cluster.

Most galactic nuclei are now believed to harbour supermassive black holes. Studies of stellar motions in the central few light-years of our Milky Way Galaxy indicate the presence of a dark object with a mass of about 2.6 million solar masses. This object is spatially coincident with Sagittarius A* (Sgr A*), the unique compact radio source located at the dynamical centre of our Galaxy. By analogy with distant quasars and nearby active galactic nuclei (AGN), Sgr A* is thought to be powered by the gravitational potential energy released by matter as it accretes onto a supermassive black hole. However, Sgr A* is much fainter than expected in all wavebands, especially in X-rays, casting some doubt on this model. Recently, we reported the first strong evidence of X-ray emission from Sgr A*. Here we report the discovery of rapid X-ray flaring from the direction of Sgr A*. These data provide compelling evidence that the X-ray emission is coming from accretion onto a supermassive black hole at the Galactic Centre, and the nature of the variations provides strong constraints on the astrophysical processes near the event horizon of the black hole.

An extremely deep X-ray survey (about 1 Ms) of the Hubble Deep Field North and its environs (about 450 arcmin^2) has been performed with the Advanced CCD Imaging Spectrometer on board the Chandra X-ray Observatory. This is one of the two deepest X-ray surveys ever performed; for point sources near the aim point it reaches 0.5-2.0 keV and 2-8 keV flux limits of 3 x 10^{-17} erg/cm^2/s and 2 x 10^{-16} erg/cm^2/s, respectively. Here we provide source catalogs along with details of the observations, data reduction, and technical analysis. Observing conditions, such as background, were excellent for almost all of the exposure. We have detected 370 distinct point sources: 360 in the 0.5-8.0 keV band, 325 in the 0.5-2.0 keV band, 265 in the 2-8 keV band, and 145 in the 4-8 keV band. Two new Chandra sources in the HDF-N itself are reported and discussed. Source positions are accurate to within 0.6-1.7 arcsec (at 90% confidence) depending mainly on the off-axis angle. We also detect two highly significant extended X-ray sources and several other likely extended X-ray sources. We present basic number count results for sources located near the center of the field. Source densities of 7100^{+1100}_{-940} deg^{-2} (at 4.2 x 10^{-17} erg/cm^2/s) and 4200^{+670}_{-580} deg^{-2} (at 3.8 x 10^{-16} erg/cm^2/s) are observed in the soft and hard bands, respectively.

We review topics related to the dispersal of young stars from their birth-sites, and focus in particular on the entourage of young stars related to the ongoing star-formation event in the Sco-Cen OB association. We conduct a follow-up kinematic study to that presented in Mamajek, Lawson, & Feigelson (2000; ApJ 544, 356) amongst nearby, isolated, young stars. In addition to the eta Cha and TW Hya groups, we find several more intriguing Sco-Cen outlier candidates: most notably beta Pic, PZ Tel, HD 199143, and HD 100546. We discuss the connection between Sco-Cen and the southern ``150 pc Conspiracy'' molecular clouds, and in particular, Corona Australis. The kinematic evidence suggests that many of the nearby, isolated ~10 Myr-old stars were born near Sco-Cen during the UCL and LCC starbursts 10-15 Myr ago. We hypothesize that these stars inherited 5-10 km/s velocities moving away from Sco-Cen, either through molecular cloud turbulence, or through formation in molecular clouds associated with the expanding Sco-Cen superbubbles (e.g. Loop I).

We present results from a 479.7 ks Chandra exposure of the Hubble Deep Field North (HDF-N) and its immediate vicinity. In our X-ray image, the deepest ever reported with a 0.5-2.0 keV flux limit of about 4.9 x 10^{-17} erg/cm^2/s, we detect four new HDF-N X-ray sources bringing the total number of such sources to 12. The new sources include two optically bright (R=18.3-18.8), low-redshift (z<0.15) galaxies, a Fanaroff-Riley I radio galaxy, and an edge-on spiral hosting either a powerful starburst or a low-luminosity active galactic nucleus (AGN). Notably, we have now detected X-ray emission from all luminous galaxies (M_V<-18) with z<0.15 known in the HDF-N. We have also detected the remarkable microJy radio source VLA J123642.09+621331.4, which is located just outside the HDF-N and has a likely redshift of z=4.424. We have detected X-ray variability from two of the previously known HDF-N X-ray sources, and spectral fitting shows clear evidence for X-ray absorption in the brightest X-ray source in the HDF-N, a z=0.960 broad-line AGN with associated Mg II absorption. Stacking analyses of optically bright HDF-N galaxies not individually detected in X-rays have provided estimates of their average X-ray fluxes, and we find that the X-ray luminosities of `normal' spirals at z~0.5 are not more than a factor of \~2 larger (per unit B-band luminosity) than those of spirals in the local Universe (z<0.01). This constrains models for the evolution of low-mass X-ray binary populations in galaxies in response to the declining cosmic star-formation rate. Monte-Carlo simulations support the validity of the stacking analyses and show that the Chandra Advanced CCD Imaging Spectrometer (ACIS) performs source detection well even with effective exposure times of ~8 Ms. (Abridged)

We report preliminary results from an arcsecond-resolution X-ray survey of nearby galaxies using ACIS on board Chandra. The total sample consists of 41 low-luminosity AGNs, including Seyferts, LINERs, and LINER/H II transition objects. In the initial subsample of 24 objects observed thus far, we detect in 62% of the objects a compact, pointlike source astrometrically coincident with either the optical or radio position of the nucleus. The high detection rate strongly suggests that the majority of the objects do contain weakly active, AGN-like cores, presumably powered by central massive black holes. The 2-10 keV luminosities of the nuclear sources range from 10^38 to 10^41 erg/s, with a median value of 2x10^38 erg/s. Our detection limit corresponds to 8x10^37 erg/s for the typical sample distance of 12 Mpc; this limit is two orders of magnitude fainter than the weakest sources of this kind previously studied using ASCA or BeppoSAX. The new data extend toward lower luminosities the known linear correlation between hard X-ray and H-alpha luminosity for broad-line AGNs. Many narrow-line objects do contain X-ray cores, consistent with either weak AGNs or X-ray binary systems, but they have X-ray luminosities a factor of 10 below the L_X-L(Halpha) relation of the broad-line sources. Their distributions of photon energies show no indication of exceptionally high absorption. The optical line emission in these nuclei is likely powered, at least in part, by stellar processes.

Observations of nearby star forming clouds with imaging X-ray telescopes have revealed that X-ray emission is elevated $10^1-10^4$ above main sequence levels in low-mass pre-main sequence (PMS) stars. The variability and spectral X-ray evidence, together with circularly polarized radio continuum flares seen in a few cases, strongly argues for an origin in magnetic reconnection flares. These high levels of magnetic activity are present from the earliest protostellar phase to the main sequence. After a brief review of past observations, three astrophysical issues are raised: the location of the flaring magnetic fields, the origin of these fields, and the effects of flare high-energy photons and particles on the environs. New results from Chandra observations of a well-defined samples of PMS solar analogues are presented, giving an improved view of magnetic flaring in the early Sun.

We report on the X-ray emission from the shell-like, non-thermal radio source Sgr A East located in the inner few parsecs of the Galaxy based on observations made with the ACIS detector on board Chandra. This is the first time Sgr A East has been clearly resolved from other complex structures in the region. The X-ray emitting region is concentrated within the central $\simeq 2$ pc of the larger radio shell. The spectrum shows strong K$\alpha$ lines from highly ionized ions of S, Ar, Ca, and Fe. A simple isothermal plasma model gives electron temperature $\sim 2$ keV, absorption column $\sim 1 \times 10^{23}$ H/cm^2, luminosity $\sim 8 \times 10^{34}$ ergs/s in the 2--10 keV band, and gas mass $\sim 2\eta^{1/2}$ M$_{\odot}$ with a filling factor $\eta$. The plasma appears to be rich in heavy elements, over-abundant by roughly a factor of four with respect to solar abundances, and shows a spatial gradient of elemental abundance: the spatial distribution of iron is more compact than that of the lighter elements. The gas mass and elemental abundance of the X-ray emission support the long-standing hypothesis that Sgr A East is a supernova remnant (SNR), maybe produced by the Type II supernova explosion of a massive star with a main-sequence mass of 13--20 M$_\odot$. The combination of the radio and X-ray morphologies classifies Sgr A East as a new metal-rich ``mixed morphology'' (MM) SNR. The size of the Sgr A East radio shell is the smallest of the known MM SNRs, which strongly suggests that the ejecta have expanded into a very dense interstellar medium. The ejecta-dominated chemical compositions of the plasma indicate that the ambient materials should be highly homogeneous. The relation between Sgr A East and the massive blackhole Sgr A* is discussed.

We present results of our Chandra observation with the ACIS-I instrument centered on the position of Sagittarius A* (Sgr A*), the compact nonthermal radio source associated with the massive black hole (MBH) at the dynamical center of the Milky Way Galaxy. We have obtained the first high-spatial-resolution (~1 arcsec), hard X-ray (0.5-7 keV) image of the central 40 pc (17 arcmin) of the Galaxy and have discovered an X-ray source, CXOGC J174540.0-290027, coincident with the radio position of Sgr A* to within 0.35 arcsec, corresponding to a maximum projected distance of 16 light-days for an assumed distance to the center of the Galaxy of 8.0 kpc. We received 222 +/-17 (1 sigma) net counts from the source in 40.3 ks. Due to the low number of counts, the spectrum is well fit either by an absorbed power-law model with photon index Gamma = 2.7 (1.8-4.0) and column density NH = (9.8 [6.8-14.2]) x 10^22 cm^-2 (90% confidence interval) or by an absorbed optically thin thermal plasma model with kT = 1.9 (1.4-2.8) keV and NH = (11.5 [8.4-15.9]) x 10^22 cm^-2. Using the power-law model, the measured (absorbed) flux in the 2-10 keV band is (1.3 [1.1-1.7]) x 10^-13 ergs cm^-2 s^-1, and the absorption-corrected luminosity is (2.4 [1.8-5.4]) x 10^33 ergs s^-1. We also briefly discuss the complex structure of the X-ray emission from the Sgr A radio complex and along the Galactic plane and present morphological evidence that Sgr A* and Sgr A West lie within the hot plasma in the central cavity of Sgr A East.

We present results from a 221.9 ks Chandra exposure of the HDF-N and its vicinity, concentrating on the 8.6' X 8.7' area covered by the Caltech Faint Field Galaxy Redshift Survey (the `Caltech area'). The minimum detectable fluxes in the 0.5-2 keV and 2-8 keV bands are 1.3e-16 cgs and 6.5e-16 cgs, respectively and a total of 82 sources are detected. More than 80% of the extragalactic X-ray background in the 2-8 keV band is resolved. Redshifts are available for 96% of the sources with R<23; the redshift range is 0.1-3.5 with most sources having z < 1.5. Eight of the X-ray sources are located in the HDF-N itself, including two not previously reported. A population of X-ray faint, optically bright, nearby galaxies emerges at soft-band fluxes of ~< 3e-16 cgs. We set the tightest constraints to date on the X-ray emission properties of microJy radio sources, mid-infrared sources detected by ISO, and very red (R-K_s > 5.0) objects. Where both the infrared and the X-ray coverage are deepest, 75% of the X-ray sources are detected by ISO; the high X-ray to infrared matching rate bodes well for future sensitive infrared observations of faint X-ray sources. Four of the 33 very red objects that have been identified in the Caltech area by Hogg et al. (2000) are detected in X-rays; these four are among our hardest Chandra sources, and we argue that they contain moderately luminous obscured AGN. Overall, however, the small Chandra detection fraction suggests a relatively small AGN content in the optically selected very red object population. (Abridged)

When observed with sufficiently high spatial resolution and sensitivity, star formation regions are unusually complex X-ray sources. Low-mass protostars and T Tauri stars, Herbig Ae/Be stars, OB and Wolf-Rayet stars are seen at levels $28 \leq \log L_x \leq 34$ erg/s. High-amplitude variability from magnetic reconnection flares are often present. From past star formation episodes, supernova remnants and X-ray binary systems may dominate the emission on large scales. Astrophysically, ionization of molecular material from yount stellar X-rays which may have important consequences for circumstellar disk evolution, bipolar flow ejection and star formation. We report here early results from Chandra X-ray Observatory studies of a wide range of star forming regions using the ACIS camera. They include portions of the nearby Ophiuchus and Perseus clouds, the Orion Nebula Cluster (ONC) and elsewhere in the Orion giant molecular clouds, molecular clouds and star clusters around the Galactic Center, the 30 Doradus region of the LMC, and the prototype starburst galaxy M 82.

General properties of the Orion star-forming region are discussed, with a focus on the dense Orion Nebula Cluster (ONC). This cluster contains between 2500 and 4500 objects located within a few parsecs of the eponymous Trapezium stars. Its members are aged <1 to a few Myr and encompass the full spectrum of stellar masses <50 Msun, as well as brown dwarfs detected with masses as low as <0.02 Msun (20 Mjup) thusfar. Recent results from optical, near-infrared, and x-ray studies of the stellar/sub-stellar population associated with this cluster are summarized.

About 1000 X-ray emitting young pre-main sequence (PMS) stars distributed in mass from 0.05 to 50 solar masses are detected in an image of the Orion Nebula obtained with the Advanced CCD Imaging Spectrometer on board the Chandra X-ray Observatory. This is the richest field of sources ever obtained in X-ray astronomy. ACIS sources include 85-90% of V<20 stars, plus a substantial fraction of deeply embedded stars with extinctions as high as A_V~60. The relationships between X-ray and other PMS stellar properties suggest that X-ray luminosity of lower-mass PMS stars depends more on mass, and possibly stellar rotation, than on bolometric luminosity as widely reported. A subsample of unabsorbed solar-mass stars show high X-ray luminosities around L_x ~ 2 x 10^30 erg/s for the first ~2 My but diverge during 2-10 My with X-ray emission plummeting in some stars but remaining high in others. This behavior is consistent with theories of stellar rotational history and magnetic dynamos. The Becklin-Neugebauer/Kleinman-Low (BN/KL) region of massive star formation shows both unabsorbed and embedded low-mass members of the Orion Nebula Cluster, the luminous infrared Source n, and a class of sources without photospheric counterparts that may be new magnetically active embedded PMS stars. Several X-ray sources are also variable radio emitters. Faint X-ray emission is seen close to, but apparently not coincident with the BN object. Its nature is not clear.

A young, nearby compact aggregate of X-ray emitting pre-main sequence stars was recently discovered in the vicinity of eta Cha (B8V). In this paper, we further investigate this cluster: its membership, its environs and origins. ROSAT HRI X-ray data for the cluster's T Tauri stars show high levels of magnetic activity and variability. The cluster has an anomalous X-ray luminosity function compared to other young clusters, deficient in stars with low, but detectable X-ray luminosities. This suggests that many low-mass members have escaped the surveyed core region. Photographic photometry from the USNO-A2.0 catalog indicates that additional, X-ray-quiet members exist in the cluster core region. The components of the eclipsing binary RS Cha, previously modeled in the literature as post-MS with discordant ages, are shown to be consistent with being coeval pre-MS stars. We compute the Galactic motion of the cluster from Hipparcos data, and compare it to other young stars and associations in the fourth Galactic quadrant. The kinematic study shows that the eta Cha cluster, the TW Hya association, and a new group near epsilon Cha, probably originated near the giant molecular cloud complex that formed the two oldest subgroups of the Sco-Cen OB association roughly 10-15 Myr ago. Their dispersal is consistent with the velocity dispersions seen in giant molecular clouds. A large H I filament and dust lane located near eta Cha has been identified as part of a superbubble formed by Sco-Cen OB winds and supernova remnants. The passage of the superbubble may have terminated star-formation in the eta Cha cluster and dispersed its natal molecular gas.

We have obtained two deep exposures of the rho Oph cloud core region with the ROSAT HRI. The improved position accuracy (1"-6") with respect to previous recent X-ray observations (ROSAT PSPC, and ASCA) allows us to remove positional ambiguities for the detected sources. We also cross-correlate the X-ray positions with IR sources found in the ISO-ISOCAM survey of the same region at 6.7 and 14.3 \mum, in addition to sources known from ground-based observations, which are young stars. We find that there is no statistically significant difference between the X-ray luminosity functions of HRI-detected Class II and Class III sources, i.e., T Tauri stars with and without disks, confirming that the contribution of these disks to X-ray emission or absorption, must be small. Most of the sources are variable, and their variability is consistent with a solar-like (hence magnetic) flare origin. We use the information given both by the ISOCAM survey and by our HRI deep exposure to study the T Tauri star population of the rho Oph dense cores. We confirm that essentially all Class II and Class III sources are X-ray emitters, and that a strong correlation exists between their X-ray luminosity, L_X, and their stellar luminosity, L_*, with L_X/L_* ~ 1E-4. Most of the new ISOCAM Class II sources are not detected, however, which we explain by the fact that their X-ray luminosities ``predicted'' on the basis of this correlation are too faint to be detected by the HRI. We predict that ~40 unknown faint or embedded Class III sources remain to be discovered in X-rays in the HRI/ISOCAM overlapping area, down to a limit of L_X ~ 3 \times 1E28 erg/s. We show that the bulk of these unknown Class III sources should be made of low- to very low-mass stars. Prospects for future detections with XMM-Newton and Chandra are discussed.

The quasar PKS 0637-753, the first celestial X-ray target of the Chandra X-ray Observatory, has revealed asymmetric X-ray structure extending from 3 to 12 arcsec west of the quasar, coincident with the inner portion of the jet previously detected in a 4.8 GHz radio image (Tingay et al. 1998). At a redshift of z=0.651, the jet is the largest (~100 kpc) and most luminous (~10^{44.6} ergs/s) of the few so far detected in X-rays. This letter presents a high resolution X-ray image of the jet, from 42 ks of data when PKS 0637-753 was on-axis and ACIS-S was near the optimum focus. For the inner portion of the radio jet, the X-ray morphology closely matches that of new ATCA radio images at 4.8 and 8.6 GHz. Observations of the parsec scale core using the VSOP space VLBI mission show structure aligned with the X-ray jet, placing important constraints on the X-ray source models. HST images show that there are three small knots coincident with the peak radio and X-ray emission. Two of these are resolved, which we use to estimate the sizes of the X-ray and radio knots. The outer portion of the radio jet, and a radio component to the east, show no X-ray emission to a limit of about 100 times lower flux. The X-ray emission is difficult to explain with models that successfully account for extra-nuclear X-ray/radio structures in other active galaxies. We think the most plausible is a synchrotron self-Compton (SSC) model, but this would imply extreme departures from the conventional minimum-energy and/or homogeneity assumptions. We also rule out synchrotron or thermal bremsstrahlung models for the jet X-rays, unless multicomponent or ad hoc geometries are invoked.

We present first results from an X-ray study of the Hubble Deep Field North (HDF-N) and its environs obtained using 166 ks of data collected by the Advanced CCD Imaging Spectrometer (ACIS) on board the Chandra X-ray Observatory. This is the deepest X-ray observation ever reported, and in the HDF-N itself we detect six X-ray sources down to a 0.5--8 keV flux limit of 4E-16 erg cm^-2 s^-1. Comparing these sources with objects seen in multiwavelength HDF-N studies shows positional coincidences with the extremely red object NICMOS J123651.74 +621221.4, an active galactic nucleus (AGN), three elliptical galaxies, and one nearby spiral galaxy. The X-ray emission from the ellipticals is consistent with that expected from a hot interstellar medium, and the spiral galaxy emission may arise from a `super-Eddington' X-ray binary or ultraluminous supernova remnant. Four of the X-ray sources have been detected at radio wavelengths. We also place X-ray upper limits on AGN candidates found in the HDF-N, and we present the tightest constraints yet on X-ray emission from the SCUBA submillimeter source population. None of the 10 high-significance submillimeter sources reported in the HDF-N and its vicinity is detected with Chandra ACIS. These sources appear to be dominated by star formation or have AGN with Compton-thick tori and little circumnuclear X-ray scattering.

Combining newly identified and previously known BL Lacs from the RASS-Green Bank (RGB) catalog, we present a sample of 127 BL Lacs, the largest ever derived from a single uniform survey. A Complete sample of 33 objects brighter than O=18.0 mag is also presented. These samples are compared to other known BL Lac samples and are generally found to exhibit properties intermediate between those of the previously disparate classes of High and Low energy peaked BL Lacs (HBLs and LBLs, respectively). This result is most dramatic in the distribution of the X-ray to radio logarithmic flux ratios, where the RGB BL Lacs are shown to peak precisely where the sharp dichotomy between the two subclasses was previously seen. The alpha_ro vs alpha_ox diagram also shows the RGB sample smoothly bridges the gap between the previously distinct subclasses of LBLs and HBLs. We also find a weak, but statistically significant correlation between the composite X-ray spectral index alpha_xox and alpha_ro. This implies that the more LBL-like RGB BL Lacs have a secondary source of X-ray emission, possibly from an inverse Compton component. We also present both the X-ray and radio logN-logS distributions for which the competing HBL/LBL unification scenarios have differing predictions. The unknown effects of the triple flux limit inherent in the RGB Complete sample makes quantitative analysis uncertain, but the characteristics of the RGB sample compare well both with results obtained from previous samples and with general theoretical predictions based on a simple Monte Carlo simulation. Our analysis indicates that the unimodal distribution of BL Lac properties found in the RGB sample likely reliably reflect the underlying population, while the bimodal distribution found in earlier studies arose primarily from observational selection effects.

The advent of large area deep radio and X-ray surveys is leading to the creation of many new BL Lac samples. In particular, the ROSAT All-Sky, Green Bank and FIRST surveys are proving to be rich sources of new BL Lacs. We will discuss the methods used in four independent BL Lac searches based on these surveys. Comparison of the broadband spectral energy distributions of these BL Lacs with those of previously known objects clearly points to the existence of a large previously unrecognized population of objects with characteristics intermediate between those exhibited by Low and High energy peaked BL Lacs.

A multivariate analysis of gamma-ray burst (GRB) bulk properties is presented to discriminate between distinct classes of GRBs. Several variables representing burst duration, fluence and spectral hardness are considered. Two multivariate clustering procedures are used on a sample of 797 bursts from the Third BATSE Catalog: a nonparametric average linkage hierarchical agglomerative clustering procedure validated with Wilks' $\Lambda^*$ and other MANOVA tests; and a parametric maximum likelihood model-based clustering procedure assuming multinormal populations calculated with the EM Algorithm and validated with the Bayesian Information Criterion. The two methods yield very similar results. The BATSE GRB population consists of three classes with the following Duration/Fluence/Spectrum bulk properties: Class I with long/bright/intermediate bursts, Class II with short/hard/faint bursts, and Class III with intermediate/intermediate/soft bursts. One outlier with poor data is also present. Classes I and II correspond to those reported by Kouveliotou et al. (1993), but Class III is clearly defined here for the first time.

We present new spectroscopic identifications for 169 objects in the RASS-Green Bank (RGB) catalog of radio- and X-ray-emitting AGN. These data significantly increase the fraction of bright RGB objects with classifications. Specifically, we report and discuss the classification of 66 radio-loud quasars, 53 BL Lacs, 33 Broad Line Radio Galaxies, 5 Narrow Line Radio Galaxies, 1 Seyfert I galaxy and 11 galaxies or galaxies in clusters. Over 78% of the identifications we present here are new. The observations we report were undertaken as part of our targeted search program to identify a new, large unbiased sample of BL Lac Objects and we therefore discuss the BL Lac sample extensively. Unlike many previous surveys, we impose no selection criteria based on optical morphology, color or broadband spectral energy distribution. Our classifications are based solely on a carefully defined set of self-consistent spectroscopic classification criteria. We show the 53 RGB presented here exhibit transitional properties between normal galaxies and BL Lacs discovered previously. We show there is no clear separation in CaII break strength between RGB BL Lacs and galaxies, with the distribution of break strengths varying smoothly between 0% and 50%. We also show that the newly discovered RGB BL Lacs reside in a "zone of avoidance" in the log(S_x/S_r) vs. log(S_o/S_r) diagram. This has important implications for BL Lac search strategies since it shows that RASS BL Lac samples will be severely incomplete if candidates are chosen only from among those objects with the highest S_x/S_r flux ratios.

The Rosat - Green Bank BL Lac sample consists of 119 objects and smoothly bridges the gap between the previously disparate subclasses of radio- and X-ray-selected objects. Further study of this sample will provide useful constraints to the unified scheme and help determine what, if any, modifications are necessary.

We present 5 GHz high resolution VLA observations of 2,127 radio- and X-ray-emitting sources found in both the Green Bank (GB) 5 GHz radio catalog and the ROSAT All-Sky Survey (RASS). We report core flux densities and positions accurate to +/-0.5". Combined with the GB measuremens of the total radio emission, we derive the core-to-lobe ratio of objects in our sample and discuss their core-dominance relative to samples of radio galaxies and BL Lacertae objects. Our results show the ROSAT/GreenBank (RGB) sample is approximately an order of magnitude more core-dominated than the radio galaxy sample, but is more than an order of magnitude less core-dominated than highly beamed BL Lacertae object samples. Using simple beaming models, this indicates the typical object in the RGB catalog exhibits beamed radio emission and is oriented at a modest angle to the line-of-sight (theta ~ 25-30 degrees). The case of the origin of the X-ray emission is not as clear; the data are consistent with either an anisotropic unbeamed or moderately beamed X-ray component.