We present a comprehensive analysis of Chandra X-ray observations of 15 young radio quasars at redshifts $4.5 < z < 5.0$. All sources are detected in the $0.5-7.0$ keV energy band. Emission spectra are extracted, and the average photon index for the sample is measured to be $1.5\pm0.1$. Unabsorbed rest-frame $2-10$ keV luminosities are found to range between $(0.5-23.2) \times 10^{45}$ erg s$^{-1}$. The optical-X-ray power-law spectral index $\alpha_{ox}$ is calculated for each source using optical/UV data available in the literature. The $\alpha_{ox}$-UV relationship is compared with other quasar surveys, and an anti-correlation is observed that agrees with independent estimates. Rest-frame radio and X-ray luminosities are established for the sample, and a correlation between the luminosities is detected. These multiwavelength results reinforce a lack of spectral evolution for quasars over a broad redshift range. We additionally identify three quasars from our multiwavelength analysis that are statistically significant outliers, with one source being a Compton-thick candidate in the early Universe, and discuss each in detail.

We perform detailed spectroscopy of the X-ray brightest supernova remnant (SNR) in the Large Magellanic Cloud (LMC), N132D, using Chandra archival observations. By analyzing the spectra of the entire well-defined rim, we determine the mean abundances for O, Ne, Mg, Si, S and Fe for the local LMC environment. We find evidence of enhanced O on the north-western and S on the north-eastern blast wave. By analyzing spectra interior to the remnant, we confirm the presence of a Si-rich relatively hot plasma (> 1.5 kev) that is also responsible for the Fe K emission. Chandra images show that the Fe K emission is distributed throughout the interior of the southern half of the remnant but does not extend out to the blast wave. We estimate the progenitor mass to be $15\pm5\,M_{\odot}$ using abundance ratios in different regions that collectively cover a large fraction of the remnant, as well as from the radius of the forward shock compared with models of an explosion in a cavity created by stellar winds. We fit ionizing and recombining plasma models to the Fe K emission and find that the current data cannot distinguish between the two, hence the origin of the high-temperature plasma remains uncertain. Our analysis is consistent with N132D being the result of a core-collapse supernova in a cavity created by its intermediate mass progenitor.

We summarize the 14th meeting of the International Astronomical Consortium for High Energy Calibration (IACHEC) held at \textit{Shonan Village} (Kanagawa, Japan) in May 2019. Sixty scientists directly involved in the calibration of operational and future high-energy missions gathered during 3.5 days to discuss the status of the cross-calibration between the current international complement of X-ray observatories, and the possibilities to improve it. This summary consists of reports from the various WGs with topics ranging from the identification and characterization of standard calibration sources, multi-observatory cross-calibration campaigns, appropriate and new statistical techniques, calibration of instruments and characterization of background, communication and preservation of knowledge, and results for the benefit of the astronomical community.

V473 Lyr is a classical Cepheid which is unique in having substantial amplitude variations with a period of approximately 3.3 years, thought to be similar to the Blazhko variations in RR Lyrae stars. We obtained an {\it XMM-Newton} observation of this star to followup a previous detection in X-rays. Rather than the X-ray burst and rapid decline near maximum radius seen in $\delta$ Cephei itself, the X-ray flux in V473 Lyr remained constant for a third of the pulsation cycle covered by the observation. Thus the X-rays are most probably not produced by the changes around the pulsation cycle. The X-ray spectrum is soft (kT = 0.6 keV), with X-ray properties which are consistent with a young low mass companion. Previously there was no evidence of a companion in radial velocities or in {\it Gaia} and {\it Hipparcos} proper motions. While this rules out companions which are very close or very distant, a binary companion at a separation between 30 and 300 AU is possible. This is an example of an X-ray observation revealing evidence of a low mass companion, which is important in completing the mass ratio statistics of binary Cepheids. Furthermore, the detection of a young X-ray bright companion is a further indication that the Cepheid (primary) is a Population I star, even though its pulsation behavior differs from other classical Cepheids.

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

Over the last two decades, the discovery of exoplanets has fundamentally changed our perception of the universe and humanity's place within it. Recent work indicates that a solar system's X-ray and high energy particle environment is of fundamental importance to the formation and development of the atmospheres of close-in planets such as hot Jupiters, and Earth-like planets around M stars. X-ray imaging and spectroscopy provide powerful and unique windows into the high energy flux that an exoplanet experiences, and X-ray photons also serve as proxies for potentially transfigurative coronal mass ejections. Finally, if the host star is a bright enough X-ray source, transit measurements akin to those in the optical and infrared are possible and allow for direct characterization of the upper atmospheres of exoplanets. In this brief white paper, we discuss contributions to the study of exoplanets and their environs which can be made by X-ray data of increasingly high quality that are achievable in the next 10--15 years.

It is now recognized that energetic stellar photon and particle radiation evaporates and erodes planetary atmospheres and controls upper atmospheric chemistry. Key exoplanet host stars will be too faint at X-ray wavelengths for accurate characterization using existing generation and future slated X-ray telescopes. Observation of stellar coronal mass ejections and winds are also beyond current instrumentation. In line with theCommittee on an Exoplanet Science Strategy recognition that holistic observational approaches are needed, we point out here that a full understanding of exoplanet atmospheres, their evolution and determination of habitability requires a powerful high-resolution X-ray imaging and spectroscopic observatory. This is the only capability that can: (1) characterize by proxy the crucial, difficult to observe, EUV stellar flux, its history and its variations for planet hosting stars; (2) observe the stellar wind; (3) detect the subtle Doppler signatures of coronal mass ejections.

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.

All cool stars show magnetic activity, and X-ray emission is the hallmark of this activity. Gaining an understanding of activity aids us in answering fundamental questions about stellar astrophysics and in determining the impact of activity on the exoplanets that orbit these stars. Stellar activity is driven by magnetic fields, which are ultimately powered by convection and stellar rotation. However, the resulting dynamo properties heavily depend on the stellar interior structure and are far from being understood. X-ray radiation can evaporate exoplanet atmospheres and damage organic materials on the planetary surface, reducing the probability that life can form or be sustained, but also provides an important source of energy for prebiotic chemical reactions. Over the next two years, the TESS mission will deliver a catalog of the closest exoplanets, along with rotation periods and activity diagnostics for millions of stars, whether or not they have a planet. We propose to include all cool stars that are TESS targets and bright enough for Chandra observations, as determined by their detection and flux in the ROSAT all-sky survey (RASS), to the list of Chandra Cool Attitude Targets (CATs). For each target, the signal will be sufficient to fit the coronal plasma with at least two temperature components, and compare abundances of groups of elements with low, medium, or high first ionization potential. Similar to the known relation between X-ray luminosity $L_X$ and rotation period, we can correlate stellar properties with coronal temperatures and abundances to constrain models for stellar activity, coronal heating, and stellar dynamos. Detailed X-ray characterization for even a subset of planet-hosting systems would dramatically advance our knowledge of what impact these emissions have on orbiting planets.

We propose a new soft clustering scheme for classifying galaxies in different activity classes using simultaneously 4 emission-line ratios; log([NII ]/Ha), log([SII]/Ha), log([OI]/Ha) and log([OIII]/Hb). We fit 20 multivariate Gaussian distributions to the 4-dimensional distribution of these lines obtained from the Sloan Digital Sky Survey (SDSS) in order to capture local structures and subsequently group the multivariate Gaussian distributions to represent the complex multi-dimensional structure of the joint distribution of galaxy spectra in the 4 dimensional line ratio space. The main advantages of this method are the use of all four optical-line ratios simultaneously and the adoption of a clustering scheme. This maximises the available information, avoids contradicting classifications, and treats each class as a distribution resulting in soft classification boundaries and providing the probability for an object to belong to each class. We also introduce linear multi-dimensional decision surfaces using support vector machines based on the classification of our soft clustering scheme. This linear multi-dimensional hard clustering technique shows high classification accuracy with respect to our soft-clustering scheme.

The X-ray and EUV emission of stars plays a key role in the loss and evolution of the atmospheres of their planets. The coronae of dwarf stars later than M6 appear to behave differently to those of earlier spectral types and are more X-ray dim and radio bright. Too faint to have been observed by the Extreme Ultraviolet Explorer, their EUV behavior is currently highly uncertain. We have devised a method to use the Chandra X-ray Observatory High Resolution Camera to provide a measure of EUV emission in the 50-170 \AA\ range and have applied it to the M6.5 dwarf LHS 248 in a pilot 10 ks exposure. Analysis with model spectra using simple, idealised coronal emission measure distributions inspired by an analysis of Chandra HETG spectra of the M5.5 dwarf Proxima Cen and results from the literature, finds greatest consistency with a very shallow emission measure distribution slope, $DEM \propto T^{3/2}$ or shallower, in the range $\log T=5.5$-$6.5$. Within $2\sigma$ confidence, a much wider range of slopes can be accommodated. Model spectra constrained by this method can provide accurate (within a factor of 2-4) synthesis and extrapolation of EUV spectra for wavelengths $<400$-500 \AA. At longer wavelengths models are uncertain by an order of magnitude or more, and depend on the details of the emission measure distribution at temperatures $\log T < 5.5$. The method is sensitive to possible incompleteness of plasma radiative loss models in the 30 170 \AA\ range for which re-examination would be warranted.

We summarize the outcome of the 13th meeting of the International Astronomical Consortium for High Energy Calibration (IACHEC), held at Tenuta dei Ciclamini (Avigliano Umbro, Italy) in April 2018. Fifty-one scientists directly involved in the calibration of operational and future high-energy missions gathered during 3.5 days to discuss the current status of the X-ray payload inter-calibration and possible approaches to improve it. This summary consists of reports from the various working groups with topics ranging from the identification and characterization of standard calibration sources, multi-observatory cross-calibration campaigns, appropriate and new statistical techniques, calibration of instruments and characterization of background, and communication and preservation of knowledge and results for the benefit of the astronomical community.

We analyze the X-ray spectra of the $\sim$8000 sources detected in the Cygnus OB2 Chandra Legacy Survey (Drake et al., this issue), with the goals of characterizing the coronal plasma of the young low-mass stars in the region and estimating their intrinsic X-ray luminosities. We adopt two different strategies for X-ray sources for which more or less than 20 photons were detected. For the brighter sample we fit the spectra with absorbed isothermal models. In order to limit uncertainties, for most of the fainter Cygnus OB2 members in this sample, we constrain the spectral parameters to characteristic ranges defined from the brightest stars. For X-ray sources with $<$20 net photons we adopt a conversion factor from detected photon flux to intrinsic flux. This was defined, building on the results for the previous sample, as a function of the 20% quantile of the detected photon energy distributions, which we prove to also correlate well with extinction. We then use the X-ray extinction from the spectral fits to constrain the ratio between optical and X-ray extinction toward Cygnus OB2, finding it consistent with standard "Galactic"' values, when properly accounting for systematics. Finally we exploit the large number of sources to constrain the average coronal abundances of several elements, through two different ensemble analyses of the X-ray spectra of low-mass Cygnus OB2 members. We find the pattern of abundances to be largely consistent with that derived for the young stellar coronae in the Orion Nebula Cluster.

Information about the physical properties of astrophysical objects cannot be measured directly but is inferred by interpreting spectroscopic observations in the context of atomic physics calculations. Ratios of emission lines, for example, can be used to infer the electron density of the emitting plasma. Similarly, the relative intensities of emission lines formed over a wide range of temperatures yield information on the temperature structure. A critical component of this analysis is understanding how uncertainties in the underlying atomic physics propagates to the uncertainties in the inferred plasma parameters. At present, however, atomic physics databases do not include uncertainties on the atomic parameters and there is no established methodology for using them even if they did. In this paper we develop simple models for the uncertainties in the collision strengths and decay rates for Fe XIII and apply them to the interpretation of density sensitive lines observed with the EUV Imagining spectrometer (EIS) on Hinode. We incorporate these uncertainties in a Bayesian framework. We consider both a pragmatic Bayesian method where the atomic physics information is unaffected by the observed data, and a fully Bayesian method where the data can be used to probe the physics. The former generally increases the uncertainty in the inferred density by about a factor of 5 compared with models that incorporate only statistical uncertainties. The latter reduces the uncertainties on the inferred densities, but identifies areas of possible systematic problems with either the atomic physics or the observed intensities.

Observations of young open clusters show a bimodal distribution of rotation periods that has been difficult to explain with existing stellar spin-down models. Detailed MHD stellar wind simulations have demonstrated that surface magnetic field morphology has a strong influence on wind-driven angular momentum loss. Observations suggest that faster rotating stars store a larger fraction of their magnetic flux in higher-order multipolar components of the magnetic field. In this work, we present a new model for stellar spin-down that, for the first time, accounts for the stellar surface magnetic field configuration. We show how a magnetic complexity that evolves from complex toward simple configurations as a star spins down can explain the salient features of stellar rotation evolution, including the bimodal distribution of both slow and fast rotators seen in young open clusters.

We present a large-scale study of diffuse X-ray emission in the nearby massive stellar association Cygnus OB2 as part of the Chandra Cygnus OB2 Legacy Program. We used 40 Chandra X-ray ACIS-I observations covering $\sim$1.0 deg$^2$. After removing 7924 point-like sources detected in our survey, background-corrected X-ray emission, the adaptive smoothing reveals large-scale diffuse X-ray emission. Diffuse emission was detected in the sub-bands Soft [0.5 : 1.2] and Medium [1.2 : 2.5], and marginally in the Hard [2.5 : 7.0] keV band. From X-ray spectral analysis of stacked spectra we compute a total [0.5 : 7.0 keV] diffuse X-ray luminosity of L$_{\rm x}^{\rm diff}\approx$4.2$\times$10$^{\rm 34}$ erg s$^{-1}$, characterized with plasma temperature components at kT$\approx$ 0.11, 0.40 and 1.18 keV, respectively. The HI absorption column density corresponding to these temperatures has a distribution consistent with N$_{\rm H}$ = 0.43, 0.80 and 1.39 $\times$10$^{22}$ cm$^{-2}$. The extended medium band energy emission likely arises from O-type stellar winds thermalized by wind-wind collisions in the most populated regions of the association, while the soft band emission probably arises from less energetic termination shocks against the surrounding Interstellar-Medium. Super-soft and Soft diffuse emission appears more widely dispersed and intense than the medium band emission. The diffuse X-ray emission is generally spatially coincident with low-extinction regions that we attribute to the ubiquitous influence of powerful stellar winds from massive stars and their interaction with the local Interstellar-Medium. Diffuse X-ray emission is volume-filling, rather than edge-brightened, oppositely to other star-forming regions. We reveal the first observational evidence of X-ray haloes around some evolved massive stars.

We study the thermal structure and energetics of the point-like EUV brightenings within a system of fan loops observed in the active region \textsl{AR~11520}. These brightenings were simultaneously observed on 2012 July 11 by the HIgh-resolution Coronal (Hi-C) imager and the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO). We identified 27 brightenings by automatically determining intensity enhancements in both Hi-C and AIA ~193~{\AA} light curves. The energetics of these brightenings were studied by using the Differential Emission Measure (DEM) diagnostics. The DEM weighted temperatures of these transients are in the range $\log{T} (K) = 6.2 - 6.6$ with radiative energies ${\approx}10^{24-25}$~ergs and densities ${\approx}$ a few times $10^{9}$~cm$^{-3}$. To the best of our knowledge, these are the smallest brightenings in EUV ever detected. We used these results to determine the mechanism of energy loss in these brightenings. Our analysis reveals that the dominant mechanism of energy loss for all the identified brightenings is conduction rather than radiation.

We report the simultaneous presence of chromospheric umbral flashes and associated umbral waves, and propagating coronal disturbances, in a sunspot and related active region. We have analyzed time-distance maps obtained using the observations from Atmospheric Imaging Assembly (AIA) on-board Solar Dynamics Observatory (SDO). These maps show the simultaneous occurrence of different sunspot oscillation and waves such as umbral flashes, umbral waves, and coronal waves. Analysis of the original light curves, i.e., without implementing any Fourier filtering on them, show that the amplitudes of different sunspot waves observed at different atmospheric layers change in synchronization with the light curves obtained from the umbral flash region, thus demonstrating that these oscillations are modulated by umbral flashes. This study provides the first observational evidence of the influence of sunspot oscillations within the umbra on other sunspot waves extending up to the corona. The properties of these waves and oscillations can be utilized to study the inherent magnetic coupling among different layers of the solar atmosphere above sunspots.

The Large Synoptic Survey Telescope is designed to provide an unprecedented optical imaging dataset that will support investigations of our Solar System, Galaxy and Universe, across half the sky and over ten years of repeated observation. However, exactly how the LSST observations will be taken (the observing strategy or "cadence") is not yet finalized. In this dynamically-evolving community white paper, we explore how the detailed performance of the anticipated science investigations is expected to depend on small changes to the LSST observing strategy. Using realistic simulations of the LSST schedule and observation properties, we design and compute diagnostic metrics and Figures of Merit that provide quantitative evaluations of different observing strategies, analyzing their impact on a wide range of proposed science projects. This is work in progress: we are using this white paper to communicate to each other the relative merits of the observing strategy choices that could be made, in an effort to maximize the scientific value of the survey. The investigation of some science cases leads to suggestions for new strategies that could be simulated and potentially adopted. Notably, we find motivation for exploring departures from a spatially uniform annual tiling of the sky: focusing instead on different parts of the survey area in different years in a "rolling cadence" is likely to have significant benefits for a number of time domain and moving object astronomy projects. The communal assembly of a suite of quantified and homogeneously coded metrics is the vital first step towards an automated, systematic, science-based assessment of any given cadence simulation, that will enable the scheduling of the LSST to be as well-informed as possible.

We summarize the outcome of the 12th meeting of the International Astronomical Consortium for High Energy Calibration (IACHEC), held at the UCLA conference center in Lake Arrowhead (California) in March 2017. 56 scientists directly involved in the calibration of operational and future high-energy missions gathered during 3.5 days to discuss the status of the X-ray payload inter-calibration, as well as possible ways to improve it. The "Thermal Supernovas Remnant" (SNR) Working Group presented a recently published paper on 1E0102.2-7219 as a calibration standard in the 0.5-1.0 keV band. A new method to measure the high-energy spectrum of the Crab Nebula and pulsar with NuSTAR without using its optics may yield a new absolute flux standard in the 3-7 keV band. A new ACIS contamination model - released with CALDB version 4.7.3 - leads to a significant improvement in modeling the spectral, spatial, and temporal properties of the contaminant. The first calibration results of the scientific payload on board Hitomi confirm the excellent performance of the instruments before the spacecraft operation problems leading to its loss. Finally, the meeting discussed extensively a novel statistic approach to formally identify in which direction the effective areas of different instruments would need to be changed to bring them into concordance. This method could inform future further calibration efforts.

The gravitational field of a galaxy can act as a lens and deflect the light emitted by a more distant object such as a quasar. Strong gravitational lensing causes multiple images of the same quasar to appear in the sky. Since the light in each gravitationally lensed image traverses a different path length from the quasar to the Earth, fluctuations in the source brightness are observed in the several images at different times. The time delay between these fluctuations can be used to constrain cosmological parameters and can be inferred from the time series of brightness data or light curves of each image. To estimate the time delay, we construct a model based on a state-space representation for irregularly observed time series generated by a latent continuous-time Ornstein-Uhlenbeck process. We account for microlensing, an additional source of independent long-term extrinsic variability, via a polynomial regression. Our Bayesian strategy adopts a Metropolis-Hastings within Gibbs sampler. We improve the sampler by using an ancillarity-sufficiency interweaving strategy and adaptive Markov chain Monte Carlo. We introduce a profile likelihood of the time delay as an approximation of its marginal posterior distribution. The Bayesian and profile likelihood approaches complement each other, producing almost identical results; the Bayesian method is more principled but the profile likelihood is simpler to implement. We demonstrate our estimation strategy using simulated data of doubly- and quadruply-lensed quasars, and observed data from quasars Q0957+561 and J1029+2623.

Observations indicate that magnetic fields in rapidly rotating stars are very strong, on both small and large scales. What is the nature of the resulting corona? Here we seek to shed some light on this question. We use the results of an anelastic dynamo simulation of a rapidly rotating fully-convective M-star to drive a physics-based model for the stellar corona. We find that due to the several kilo Gauss large-scale magnetic fields at high latitudes, the corona and its X-ray emission are dominated by star-size large hot loops, while the smaller, underlying colder loops are not visible much in the X-ray. Based on this result we propose that, in rapidly rotating stars, emission from such coronal structures dominates the quiescent, cooler but saturated X-ray emission.

Flares we observe on stars in white light, UV or soft X-rays are probably harbingers of coronal mass ejections (CMEs). If we use the Sun as a guide, large stellar flares will dissipate two orders of magnitude less X-ray radiative energy than the kinetic energy in the associated CME. Since coronal emission on active stars appears to be dominated by flare activity, CMEs pose a quandary for understanding the fraction of their energy budget stars can spend on magnetic activity. One answer is magnetic suppression of CMEs, in which the strong large-scale fields of active stars entrap and prevent CMEs unless their free energy exceeds a critical value. The CME-less flaring active region NOAA 2192 presents a possible solar analogue of this. Monster CMEs will still exist, and have the potential to ravage planetary atmospheres.

Stars of stellar type later than about M3.5 are believed to be fully convective and therefore unable to support magnetic dynamos like the one that produces the 11-year solar cycle. Because of their intrinsic faintness, very few late M stars have undergone long-term monitoring to test this prediction, which is critical to our understanding of magnetic field generation in such stars. Magnetic activity is also of interest as the driver of UV and X-ray radiation, as well as energetic particles and stellar winds, that affect the atmospheres of close-in planets that lie within habitable zones, such as the recently discovered Proxima b. We report here on several years of optical, UV, and X-ray observations of Proxima Centauri (GJ 551; dM5.5e): 15 years of ASAS photometry in the V band (1085 nights) and 3 years in the I band (196 nights), 4 years of Swift XRT and UVOT observations (more than 120 exposures), and 9 sets of X-ray observations from other X-ray missions (ASCA, XMM-Newton, and three Chandra instruments) spanning 22 years. We confirm previous reports of an 83-day rotational period and find strong evidence for a 7-year stellar cycle, along with indications of differential rotation at about the solar level. X-ray/UV intensity is anti-correlated with optical V-band brightness for both rotational and cyclical variations. From comparison with other stars observed to have X-ray cycles we deduce a simple empirical relationship between X-ray cyclic modulation and Rossby number, and we also present Swift UV grism spectra covering 2300-6000 {\AA}.

The Large Synoptic Survey Telescope (LSST) will be a discovery machine for the astronomy and physics communities, revealing astrophysical phenomena from the Solar System to the outer reaches of the observable Universe. While many discoveries will be made using LSST data alone, taking full scientific advantage of LSST will require ground-based optical-infrared (OIR) supporting capabilities, e.g., observing time on telescopes, instrumentation, computing resources, and other infrastructure. This community-based study identifies, from a science-driven perspective, capabilities that are needed to maximize LSST science. Expanding on the initial steps taken in the 2015 OIR System Report, the study takes a detailed, quantitative look at the capabilities needed to accomplish six representative LSST-enabled science programs that connect closely with scientific priorities from the 2010 decadal surveys. The study prioritizes the resources needed to accomplish the science programs and highlights ways that existing, planned, and future resources could be positioned to accomplish the science goals.

We analyze Chandra X-ray images of a sample of 11 quasars that are known to contain kiloparsec scale radio jets. The sample consists of five high-redshift (z >= 3.6) flat-spectrum radio quasars, and six intermediate redshift (2.1 < z < 2.9) quasars. The dataset includes four sources with integrated steep radio spectra and seven with flat radio spectra. A total of 25 radio jet features are present in this sample. We apply a Bayesian multi-scale image reconstruction method to detect and measure the X-ray emission from the jets. We compute deviations from a baseline model that does not include the jet, and compare observed X-ray images with those computed with simulated images where no jet features exist. This allows us to compute p-value upper bounds on the significance that an X- ray jet is detected in a pre-determined region of interest. We detected 12 of the features unambiguously, and an additional 6 marginally. We also find residual emission in the cores of 3 quasars and in the background of 1 quasar that suggest the existence of unresolved X-ray jets. The dependence of the X-ray to radio luminosity ratio on redshift is a potential diagnostic of the emission mechanism, since the inverse Compton scattering of cosmic microwave background photons (IC/CMB) is thought to be redshift dependent, whereas in synchrotron models no clear redshift dependence is expected. We find that the high-redshift jets have X-ray to radio flux ratios that are marginally inconsistent with those from lower redshifts, suggesting that either the X-ray emissions is due to the IC/CMB rather than the synchrotron process, or that high redshift jets are qualitatively different.

This is the stars chapter of the Kavli workshop report, which resulted from the community-based study of needed US OIR capabilities in the LSST era. The full report, which will include this chapter, is anticipated to be available in Fall 2016. See NOAO website (http://www.noao.edu/meetings/lsst-oir-study/) for more details.

In our Galaxy, star formation occurs in a variety of environments, with a large fraction of stars formed in clusters hosting massive stars. OB stars have an important feedback on the evolution of protoplanetary disks around nearby young stars and likely on the process of planet formation occurring in them. The nearby massive association Cygnus OB2 is an outstanding laboratory to study this feedback. It is the closest massive association to our Sun, and hosts hundreds of massive stars and thousands of low mass members. In this paper, we analyze the spatial variation of the disk fraction in Cygnus OB2 and we study its correlation with the local values of Far and Extreme ultraviolet radiation fields and the local stellar surface density. We present definitive evidence that disks are more rapidly dissipated in the regions of the association characterized by intense local UV field and large stellar density. In particular, the FUV radiation dominates disks dissipation timescales in the proximity (i.e. within 0.5 pc) of the O stars. In the rest of the association, EUV photons potentially induce a significant mass loss from the irradiated disks across the entire association, but the efficiency of this process is reduced at increasing distances from the massive stars due to absorption by the intervening intracluster material. We find that disk dissipation due to close stellar encounters is negligible in Cygnus OB2, and likely to have affected 1% or fewer of the stellar population. Disk dissipation is instead dominated by photoevaporation. We also compare our results to what has been found in other young clusters with different massive populations, concluding that massive associations like Cygnus OB2 are potentially hostile to protoplanetary disks, but that the environments where disks can safely evolve in planetary systems are likely quite common in our Galaxy.

The recurrent symbiotic nova V745 Sco exploded on 2014 February 6 and was observed on February 22 and 23 by the Chandra X-ray Observatory Transmission Grating Spectrometers. By that time the supersoft source phase had already ended and Chandra spectra are consistent with emission from a hot, shock-heated circumstellar medium with temperatures exceeding 10^7K. X-ray line profiles are more sharply peaked than expected for a spherically-symmetric blast wave, with a full width at zero intensity of approximately 2400 km/s, a full width at half maximum of 1200 +/- 30 km/s and an average net blueshift of 165 +/- 10 km/s. The red wings of lines are increasingly absorbed toward longer wavelengths by material within the remnant. We conclude that the blast wave was sculpted by an aspherical circumstellar medium in which an equatorial density enhancement plays a role, as in earlier symbiotic nova explosions. Expansion of the dominant X-ray emitting material is aligned close to the plane of the sky and most consistent with an orbit seen close to face-on. Comparison of an analytical blast wave model with the X-ray spectra, Swift observations and near-infrared line widths indicates the explosion energy was approximately 10^43 erg, and confirms an ejected mass of approximately 10^-7 Msun. The total mass lost is an order of magnitude lower than the accreted mass required to have initiated the explosion, indicating the white dwarf is gaining mass and is a supernova Type 1a progenitor candidate.

We investigate the uncertainties of fitted X-ray model parameters and fluxes for relatively faint Chandra ACIS-I source spectra. Monte-Carlo (MC) simulations are employed to construct a large set of 150,000 fake X-ray spectra in the low photon count statistics regime (from 20 to 350 net counts) using the XSPEC spectral model fitting package. The simulations employed both absorbed thermal (APEC) and non-thermal (power-law) models, in concert with the Chandra ACIS-I instrument response and interstellar absorption. Simulated X-ray spectra were fit assuming a wide set of different input parameters and C-statistic minimization criteria to avoid numerical artifacts in the accepted solutions. Results provide an error estimate for each parameter (absorption, NH, plasma temperature, kT, or power-law slope, Gamma, and flux, and for different background contamination levels. The distributions of these errors are studied as a function of the 1 sigma quantiles and we show how these correlate with different model parameters, net counts in the spectra and relative background level. Maps of uncertainty in terms of the 1 sigma quantiles for parameters and flux are computed as a function of spectrum net counts. We find very good agreement between our estimated X-ray spectral parameter and flux uncertainties and those recovered from spectral fitting for a subset of the X-ray sources detected in the Chandra Cygnus OB2 Legacy Survey diagnosed to be Association members and that have between 20 and 350 net counts. Our method can provide uncertainties for spectral parameters whenever formal X-ray spectral fits cannot be well-constrained, or are unavailable, and predictions useful for computing Chandra ACIS-I exposure times for observation planning.

COCOA-PUFS is an energy-diverse, time-domain study of the ultra-fast spinning, heavily spotted, yellow giant FK Com (HD117555; G4 III). This single star is thought to be a recent binary merger, and is exceptionally active by measure of its intense ultraviolet and X-ray emissions, and proclivity to flare. COCOA-PUFS was carried out with Hubble Space Telescope in the UV (120-300 nm), using mainly its high-performance Cosmic Origins Spectrograph, but also high-precision Space Telescope Imaging Spectrograph; Chandra X-ray Observatory in the soft X-rays (0.5-10 keV), utilizing its High-Energy Transmission Grating Spectrometer; together with supporting photometry and spectropolarimetry in the visible from the ground. This is an introductory report on the project. FK Com displayed variability on a wide range of time scales, over all wavelengths, during the week-long main campaign, including a large X-ray flare; "super-rotational broadening" of the far-ultraviolet "hot-lines" (e.g., Si IV 139 nm (T~80,000 K) together with chromospheric Mg II 280 nm and C II 133 nm (10,000-30,000 K); large Doppler swings suggestive of bright regions alternately on advancing and retreating limbs of the star; and substantial redshifts of the epoch-average emission profiles. These behaviors paint a picture of a highly extended, dynamic, hot (10 MK) coronal magnetosphere around the star, threaded by cooler structures perhaps analogous to solar prominences, and replenished continually by surface activity and flares. Suppression of angular momentum loss by the confining magnetosphere could temporarily postpone the inevitable stellar spindown, thereby lengthening this highly volatile stage of coronal evolution.

Telescopes such as the Atmospheric Imaging Assembly aboard the Solar Dynamics Observatory, a NASA satellite, collect massive streams of high resolution images of the Sun through multiple wavelength filters. Reconstructing pixel-by-pixel thermal properties based on these images can be framed as an ill-posed inverse problem with Poisson noise, but this reconstruction is computationally expensive and there is disagreement among researchers about what regularization or prior assumptions are most appropriate. This article presents an image segmentation framework for preprocessing such images in order to reduce the data volume while preserving as much thermal information as possible for later downstream analyses. The resulting segmented images reflect thermal properties but do not depend on solving the ill-posed inverse problem. This allows users to avoid the Poisson inverse problem altogether or to tackle it on each of $\sim$10 segments rather than on each of $\sim$10$^7$ pixels, reducing computing time by a factor of $\sim$10$^6$. We employ a parametric class of dissimilarities that can be expressed as cosine dissimilarity functions or Hellinger distances between nonlinearly transformed vectors of multi-passband observations in each pixel. We develop a decision theoretic framework for choosing the dissimilarity that minimizes the expected loss that arises when estimating identifiable thermal properties based on segmented images rather than on a pixel-by-pixel basis. We also examine the efficacy of different dissimilarities for recovering clusters in the underlying thermal properties. The expected losses are computed under scientifically motivated prior distributions. Two simulation studies guide our choices of dissimilarity function. We illustrate our method by segmenting images of a coronal hole observed on 26 February 2015.

Variable-intensity astronomical sources are the result of complex and often extreme physical processes. Abrupt changes in source intensity are typically accompanied by equally sudden spectral shifts, i.e., sudden changes in the wavelength distribution of the emission. This article develops a method for modeling photon counts collected from observation of such sources. We embed change points into a marked Poisson process, where photon wavelengths are regarded as marks and both the Poisson intensity parameter and the distribution of the marks are allowed to change. To the best of our knowledge this is the first effort to embed change points into a marked Poisson process. Between the change points, the spectrum is modeled non-parametrically using a mixture of a smooth radial basis expansion and a number of local deviations from the smooth term representing spectral emission lines. Because the model is over parameterized we employ an $\ell_1$ penalty. The tuning parameter in the penalty and the number of change points are determined via the minimum description length principle. Our method is validated via a series of simulation studies and its practical utility is illustrated in the analysis of the ultra-fast rotating yellow giant star known as FK Com.

The Chandra Cygnus OB2 Legacy Survey is a wide and deep X-ray survey of the nearby and massive Cygnus OB2 association. The survey has detected ~8,000 X-ray sources, the majority of which are pre-main sequence X-ray emitting young stars in the association itself. To facilitate quantitative scientific studies of these sources as well as the underlying OB association it is important to understand the sensitivity of the observations and the level of completeness the observations have obtained. Here we describe the use of a hierarchical Monte Carlo simulation to achieve this goal by combining the empirical properties of the observations, analytic estimates of the source verification process, and an extensive set of source detection simulations. We find that our survey reaches a 90% completeness level for a pre-main-sequence population at the distance of Cyg OB2 at an X-ray luminosity of 4 x 10^30 ergs/s and a stellar mass of 1.3 Msun for a randomly distributed population. For a spatially clustered population such as Cyg~OB2 the 90% completeness level is reached at 1.1 Msun instead, as the sources are more concentrated in areas of our survey with a high exposure. These simulations can easily be adapted for use with other X-ray observations and surveys, and we provide X-ray detection efficiency curves for a very wide array of source and background properties to allow these simulations to be easily exploited by other users.

Unexpected structure in images of astronomical sources often presents itself upon visual inspection of the image, but such apparent structure may either correspond to true features in the source or be due to noise in the data. This paper presents a method for testing whether inferred structure in an image with Poisson noise represents a significant departure from a baseline (null) model of the image. To infer image structure, we conduct a Bayesian analysis of a full model that uses a multiscale component to allow flexible departures from the posited null model. As a test statistic, we use a tail probability of the posterior distribution under the full model. This choice of test statistic allows us to estimate a computationally efficient upper bound on a p-value that enables us to draw strong conclusions even when there are limited computational resources that can be devoted to simulations under the null model. We demonstrate the statistical performance of our method on simulated images. Applying our method to an X-ray image of the quasar 0730+257, we find significant evidence against the null model of a single point source and uniform background, lending support to the claim of an X-ray jet.

We present an analysis of the high-resolution Chandra observation of the multiple system, HD 155555 (an RS CVn type binary system, HD 155555 AB, and its spatially resolved low-mass companion HD 155555 C). This is an intriguing system which shows properties of both an active pre-main sequence star and a synchronised (main sequence) binary. We obtain the emission measure distribution, temperature structures, plasma densities, and abundances of this system and compare them with the coronal properties of other young/active stars. HD 155555 AB and HD 155555 C produce copious X-ray emission with log Lx of 30.54 and 29.30, respectively, in the 0.3-6.0 keV energy band. The light curves of individual stars show variability on timescales of few minutes to hours. We analyse the dispersed spectra and reconstruct the emission measure distribution using spectral line analysis. The resulting elemental abundances exhibit inverse first ionisation potential effect in both cases. An analysis of He-like triplets yields a range of coronal electron densities ~10^10-10^13 cm-3. Since HD 155555 AB is classified both as an RS CVn and a PMS star, we compare our results with those of other slightly older active main-sequence stars and T Tauri stars, which indicates that the coronal properties of HD 155555 AB closely resemble that of an older RS CVn binary rather than a younger PMS star. Our results also suggests that the properties of HD 155555 C is very similar to those of other active M dwarfs.

We present a powerful new algorithm that combines both spatial information (event locations and the point spread function) and spectral information (photon energies) to separate photons from overlapping sources. We use Bayesian statistical methods to simultaneously infer the number of overlapping sources, to probabilistically separate the photons among the sources, and to fit the parameters describing the individual sources. Using the Bayesian joint posterior distribution, we are able to coherently quantify the uncertainties associated with all these parameters. The advantages of combining spatial and spectral information are demonstrated through a simulation study. The utility of the approach is then illustrated by analysis of observations of FK Aqr and FL Aqr with the XMM-Newton Observatory and the central region of the Orion Nebula Cluster with the Chandra X-ray Observatory.

A fine analysis of spectral line widths and Doppler shifts employing Fourier transform and cross-correlation techniques has been applied to Chandra HETG spectra obtained in 1999 October of the rapidly rotating young star AB Doradus in order to investigate its coronal topology. The observation lasted 52.3ks, covering 1.2 rotation periods. The X-ray light curve obtained from integrating the dispersed signal revealed a moderate intensity flare midway through the exposure in which the count rate increased sharply by about 50% and subsequently decayed over the next 10ks. We find no significant Doppler shifts in the spectra or modulation of the light curve that could be attributed to rotation of dominant coronal structures at this epoch. Individual spectral line widths are statistically consistent with thermal broadening and formally require no rotational broadening, while the $1\sigma$ limit to rotational broadening corresponds to a compact corona restricted to latitudes $>57\deg$. Fourier analysis suggests a small amount of rotational broadening is present consistent with a corona restricted largely to the poles, and excludes models with surface rotational broadening or greater. These results present direct spectroscopic evidence that the dominant coronal activity on rapidly-rotating active stars is associated with the dark polar spots commonly seen in photospheric Doppler images, and support models in which these spots are of mixed magnetic polarity that forms closed loops.

The young massive OB association Cygnus OB2, in the Cygnus X complex, is the closest (1400 pc) star forming region to the Sun hosting thousands of young low mass stars and up to 1000 OB stars, among which are some of the most massive stars known in our Galaxy. This region holds great importance for several fields of modern astrophysics, such as the study of the physical properties of massive and young low-mass stars and the feedback provided by massive stars on star and planet formation process. Cygnus OB2 has been recently observed with Chandra/ACIS-I as part of the 1.08Msec Chandra Cygnus OB2 Legacy Project. This survey detected 7924 X-ray sources in a square degree area centered on Cyg OB2. Since a proper classification and study of the observed X-ray sources also requires the analysis of their optical and infrared counterparts, we combined a large and deep set of optical and infrared catalogs available for this region with our new X-ray catalog. In this paper we describe the matching procedure and present the combined catalog containing 5703 sources. We also briefly discuss the nature of the X-ray sources with optical and infrared counterparts using their position in the color-magnitude and color-color diagrams.

We present the results of the first strong lens time delay challenge. The motivation, experimental design, and entry level challenge are described in a companion paper. This paper presents the main challenge, TDC1, which consisted of analyzing thousands of simulated light curves blindly. The observational properties of the light curves cover the range in quality obtained for current targeted efforts (e.g.,~COSMOGRAIL) and expected from future synoptic surveys (e.g.,~LSST), and include simulated systematic errors. \nteamsA\ teams participated in TDC1, submitting results from \nmethods\ different method variants. After a describing each method, we compute and analyze basic statistics measuring accuracy (or bias) $A$, goodness of fit $\chi^2$, precision $P$, and success rate $f$. For some methods we identify outliers as an important issue. Other methods show that outliers can be controlled via visual inspection or conservative quality control. Several methods are competitive, i.e., give $|A|<0.03$, $P<0.03$, and $\chi^2<1.5$, with some of the methods already reaching sub-percent accuracy. The fraction of light curves yielding a time delay measurement is typically in the range $f = $20--40\%. It depends strongly on the quality of the data: COSMOGRAIL-quality cadence and light curve lengths yield significantly higher $f$ than does sparser sampling. Taking the results of TDC1 at face value, we estimate that LSST should provide around 400 robust time-delay measurements, each with $P<0.03$ and $|A|<0.01$, comparable to current lens modeling uncertainties. In terms of observing strategies, we find that $A$ and $f$ depend mostly on season length, while P depends mostly on cadence and campaign duration.

We present a rigorous description of the general problem of aperture photometry in high energy astrophysics photon-count images, in which the statistical noise model is Poisson, not Gaussian. We compute the full posterior probability density function for the expected source intensity for various cases of interest, including the important cases in which both source and background apertures contain contributions from the source, and when multiple source apertures partially overlap. A Bayesian approach offers the advantages that it allows one to (a) include explicit prior information on source intensities, (b) propagate posterior distributions as priors for future observations, and (c) use Poisson likelihoods, making the treatment valid in the low counts regime. Elements of this approach have been implemented in the Chandra Source Catalog.

The Cygnus OB2 association is the largest concentration of young and massive stars within 2 kpc of the Sun, including an estimated 65 O-type stars and hundreds of OB stars. The Chandra Cygnus OB2 Legacy Survey is a large imaging program undertaken with the Advanced CCD Imaging Spectrometer onboard the Chandra X-ray Observatory. The survey has imaged the central 0.5 deg^2 of the Cyg OB2 association with an effective exposure of 120ks and an outer 0.35 deg^2 area with an exposure of 60ks. Here we describe the survey design and observations, the data reduction and source detection, and present a catalog of 8,000 X-ray point sources. The survey design employs a grid of 36 heavily (~50%) overlapping pointings, a method that overcomes Chandra's low off-axis sensitivity and produces a highly uniform exposure over the inner 0.5 deg^2. The full X-ray catalog is described here and is made available online.

The nature of the mechanisms apparently driving X-rays from intermediate mass stars lacking strong convection zones or massive winds remains poorly understood, and the possible role of hidden, lower mass close companions is still unclear. A 20ks Chandra HRC-I observation of HR 4796A, an 8 Myr old main sequence A0 star devoid of close stellar companions, has been used to search for a signature or remnant of magnetic activity from the Herbig Ae phase. X-rays were not detected and the X-ray luminosity upper limit was L_X =< 1.3x10^27 erg/s. The result is discussed in the context of various scenarios for generating magnetic activity, including rotational shear and subsurface convection. A dynamo driven by natal differential rotation is unlikely to produce observable X rays, chiefly because of the difficulty in getting the dissipated energy up to the surface of the star. A subsurface convection layer produced by the ionisation of helium could host a dynamo that should be effective throughout the main-sequence but can only produce X-ray luminosities of order 10^25 erg/s. This luminosity lies only moderately below the current detection limit for Vega. Our study supports the idea that X-ray production in Herbig Ae/Be stars is linked largely to the accretion process rather than the properties of the underlying star, and that early A stars generally decline in X-ray luminosity at least 100,000 fold in only a few million years.

Abridged. Here we report on the X-ray activity of the primary star, HD189733 A, using a new XMM-Newton observation and a comparison with the previous X-ray observations. The spectrum in the quiescent intervals is described by two temperatures at 0.2 keV and 0.7 keV, while during the flares a third component at 0.9 keV is detected. We obtain estimates of the electron density in the range $n_e = 1.6 - 13 \times 10^{10}$ cm$^{-3}$ and thus the corona of HD189733 A appears denser than the solar one. {For the third time, we observe a large flare that occurred just after the eclipse of the planet. Together with the flares observed in 2009 and 2011, the events are restricted to a small planetary phase range of $\phi = 0.55-0.65$. Although we do not find conclusive evidence of a significant excess of flares after the secondary transits, we suggest that the planet might trigger such flares when it passes close to locally high magnetic field of the underlying star at particular combinations of stellar rotational phases and orbital planetary phases. For the most recent flares, a wavelet analysis of the light curve suggests a loop of length of four stellar radii at the location of the bright flare, and a local magnetic field of order of 40-100 G, in agreement with the global field measured in other studies. The loop size suggests an interaction of magnetic nature between planet and star, separated by only $\sim8 R_*$. We also detect the stellar companion (HD 189733 B, $\sim12"$ from the primary star) in this XMM observation. Its very low X-ray luminosity ($L_X = 3.4\times 10^{26}$ erg s$^{-1}$) confirms the old age of this star and of the binary system. The high activity of the primary star is best explained by a transfer of angular momentum from the planet to the star.

Extensive X-ray and extreme ultraviolet (EUV) photometric observations of the eclipsing RS CVn system AR Lac were obtained over the years 1997 to 2013 with the Chandra X-ray Observatory Extreme Ultraviolet Explorer. During primary eclipse, HRC count rates decrease by ~40%. A similar minimum is seen during one primary eclipse observed by EUVE but not in others owing to intrinsic source variability. Little evidence for secondary eclipses is present in either the X-ray or EUV data, reminiscent of earlier X-ray and EUV observations. Primary eclipses allow us to estimate the extent of a spherically symmetric corona on the primary G star of about 1.3Rsun, or 0.86Rstar, and indicate the G star is likely brighter than the K component by a factor of 2-5. Brightness changes not attributable to eclipses appear to be dominated by stochastic variability and are generally non-repeating. X-ray and EUV light curves cannot therefore be reliably used to reconstruct the spatial distribution of emission assuming only eclipses and rotational modulation are at work. Moderate flaring is observed, where count rates increase by up to a factor of three above quiescence. Combined with older ASCA, Einstein, EXOSAT, ROSAT and Beppo-SAX observations, the data show that the level of quiescent coronal emission at X-ray wavelengths has remained remarkably constant over 33 years, with no sign of variation due to magnetic cycles. Variations in base level X-ray emission seen by Chandra over 13 years are only ~10%, while variations back to pioneering Einstein observations in 1980 amount to a maximum of 45% and more typically about 15%.

This article offers a short guide to the steps scientists can take to ensure that their data and associated analyses continue to be of value and to be recognized. In just the past few years, hundreds of scholarly papers and reports have been written on questions of data sharing, data provenance, research reproducibility, licensing, attribution, privacy, and more, but our goal here is not to review that literature. Instead, we present a short guide intended for researchers who want to know why it is important to "care for and feed" data, with some practical advice on how to do that.

The formation of stars in massive clusters is one of the main modes of the star formation process. However, the study of massive star forming regions is hampered by their typically large distances to the Sun. One exception to this is the massive star forming region Cygnus OB2 in the Cygnus X region, at the distance of about 1400 pc. Cygnus OB2 hosts very rich populations of massive and low-mass stars, being the best target in our Galaxy to study the formation of stars, circumstellar disks, and planets in presence of massive stars. In this paper we combine a wide and deep set of photometric data, from the r band to 24 micron, in order to select the disk bearing population of stars in Cygnus OB2 and identify the class I, class II, and stars with transition and pre-transition disks. We selected 1843 sources with infrared excesses in an area of 1 degree x 1 degree centered on Cyg OB2 in several evolutionary stages: 8.4% class I, 13.1% flat-spectrum sources, 72.9% class II, 2.3% pre-transition disks, and 3.3% transition disks. The spatial distribution of these sources shows a central cluster surrounded by a annular overdensity and some clumps of recent star formation in the outer region. Several candidate subclusters are identified, both along the overdensity and in the rest of the association.

In order to fully understand the gravitational collapse of molecular clouds, the star formation process and the evolution of circumstellar disks, these phenomena must be studied in different Galactic environments with a range of stellar contents and positions in the Galaxy. The young massive association Cygnus OB2, in the Cygnus-X region, is an unique target to study how star formation and the evolution of circumstellar disks proceed in the presence of a large number of massive stars. We present a catalog obtained with recent optical observations in r,i,z filters with OSIRIS, mounted on the $10.4\,m$ GTC telescope, which is the deepest optical catalog of Cyg OB2 to date. The catalog consist of 64157 sources down to M=0.15 solar masses at the adopted distance and age of Cyg OB2. A total of 38300 sources have good photometry in all three bands. We combined the optical catalog with existing X-ray data of this region, in order to define the cluster locus in the optical diagrams. The cluster locus in the r-i vs. i-z diagram is compatible with an extinction of the optically selected cluster members in the 2.64<AV<5.57 range. We derive an extinction map of the region, finding a median value of AV=4.33 in the center of the association, decreasing toward the north-west. In the color-magnitude diagrams, the shape of the distribution of main sequence stars is compatible with the presence of an obscuring cloud in the foreground at about 850+/-25 pc from the Sun.

Sunspot numbers form a comprehensive, long-duration proxy of solar activity and have been used numerous times to empirically investigate the properties of the solar cycle. A number of correlations have been discovered over the 24 cycles for which observational records are available. Here we carry out a sophisticated statistical analysis of the sunspot record that reaffirms these correlations, and sets up an empirical predictive framework for future cycles. An advantage of our approach is that it allows for rigorous assessment of both the statistical significance of various cycle features and the uncertainty associated with predictions. We summarize the data into three sequential relations that estimate the amplitude, duration, and time of rise to maximum for any cycle, given the values from the previous cycle. We find that there is no indication of a persistence in predictive power beyond one cycle, and conclude that the dynamo does not retain memory beyond one cycle. Based on sunspot records up to October 2011, we obtain, for Cycle 24, an estimated maximum smoothed monthly sunspot number of 97 +- 15, to occur in January--February 2014 +- 6 months.

High-energy irradiation of exoplanets has been identified to be a key influence on the stability of these planets' atmospheres. So far, irradiation-driven mass-loss has been observed only in two Hot Jupiters, and the observational data remain even more sparse in the super-earth regime. We present an investigation of the high-energy emission in the CoRoT-7 system, which hosts the first known transiting super-earth. To characterize the high-energy XUV radiation field into which the rocky planets CoRoT-7b and CoRoT-7c are immersed, we analyzed a 25 ks XMM-Newton observation of the host star. Our analysis yields the first clear (3.5 sigma) X-ray detection of CoRoT-7. We determine a coronal temperature of ca. 3 MK and an X-ray luminosity of 3*10^28 erg/s. The level of XUV irradiation on CoRoT-7b amounts to ca. 37000 erg/cm^2/s. Current theories for planetary evaporation can only provide an order-of-magnitude estimate for the planetary mass loss; assuming that CoRoT-7b has formed as a rocky planet, we estimate that CoRoT-7b evaporates at a rate of about 1.3*10^11 g/s and has lost ca. 4-10 earth masses in total.

An XMM-Newton observation of the nearby "pre-cataclysmic" short-period (P_orb = 3.62 hr) binary QS Vir (EC 13471-1258) revealed regular narrow X-ray eclipses when the white dwarf passed behind its M2-4 dwarf companion. The X-ray emission provides a clear signature of mass transfer and accretion onto the white dwarf. The low-resolution XMM-Newton EPIC spectra are consistent with a cooling flow model and indicate an accretion rate of Mdot= 1.7\times10^-13M\odot/yr. At 48 pc distant, QS Vir is then the second nearest accreting cataclysmic variable known, with one of the lowest accretion rates found to date for a non-magnetic system. To feed this accretion through a wind would require a wind mass loss rate of Mdot ~ 2 \times 10^-12M\odot/yr if the accretion efficiency is of the order of 10%. Consideration of likely mass loss rates for M dwarfs suggests this is improbably high and pure wind accretion unlikely. A lack of accretion disk signatures also presents some difficulties for direct Roche lobe overflow. We speculate that QS Vir is on the verge of Roche lobe overflow, and that the observed mass transfer could be supplemented by upward chromospheric flows on the M dwarf, analogous to spicules and mottles on the Sun, that escape the Roche surface to be subsequently swept up into the white dwarf Roche lobe. If so, QS Vir would be in a rare evolutionary phase lasting only a million years. The X-ray luminosity of the M dwarf estimated during primary eclipse is L_X = 3 \times 10^28 erg/s, which is consistent with that of rapidly rotating "saturated" K and M dwarfs.

We study the effect of a white dwarf on the spin-down of a cataclysmic variable system using a three-dimensional magnetohydrodynamic numerical model. The model includes the stellar corona, the stellar wind, and the WD mass and magnetic field. The existence of the WD modifies the system spin-down by physically blocking the stellar wind, restructuring the wind, channeling the wind towards the WD surface, and by modifying the shape and size of the Alfv\'en surface. The combination of these processes differs among a set of simple test cases, and the resulting angular momentum loss rates vary by factors of 2-3, and by factors of two relative to a test model with a single M dwarf. While the model employs some simplifications, the results suggest angular momentum loss schemes currently employed in cataclysmic variable studies do not require drastic revision. Insights are also gained on wind accretion. We find that efficient accretion switches on quite rapidly with decreasing orbital separation. Accretion rates depend on magnetic field alignment and should be modulated by magnetic cycles on the M dwarf. For particular values of white dwarf magnetic field strength, an efficient syphoning of coronal plasma from the inward facing M dwarf hemisphere occurs. Wind accretion rates are expected to vary by factors of 10 or more between fairly similar close binaries, depending on magnetic field strengths and orbital separation.

We report on the follow-up XMM-Newton observation of the planet-hosting star HD 189733 we obtained in April 2011. We observe a flare just after the secondary transit of the hot Jupiter. This event shares the same phase and many of the characteristics of the flare we observed in 2009. We suggest that a systematic interaction between planet and stellar magnetic fields when the planet passes close to active regions on the star can lead to periodic variability phased with planetary motion. By mean of high resolution X-ray spectroscopy with RGS we determine that the corona of this star is unusually dense.

Aims: We carried out high-resolution spectroscopy and BV(I)_C photometric monitoring of the two fastest late-type rotators in the nearby Beta Pictoris moving group, HD199143 (F7V) and CD-641208 (K7V). The motivation for this work is to investigate the rotation periods and photospheric spot patterns of these very young stars, with a longer term view to probing the evolution of rotation and magnetic activity during the early phases of main-sequence evolution. We also aim to derive information on key physical parameters, such as rotational velocity and rotation period. Methods: We applied maximum entropy (ME) and Tikhonov regularizing (TR) criteria to derive the surface spot map distributions of the optical modulation observed in HD199143 (F7 V) and CD-641208 (K7V). We also used cross-correlation techniques to determine stellar parameters such as radial velocities and rotational velocities. Lomb-Scargle periodograms were used to obtain the rotational periods from differential magnitude time series. Results: We find periods and inclinations of 0.356 days and 21.5deg for HD199143, and 0.355 days and 50.1deg for CD-641208. The spot maps of HD199143 obtained from the ME and TR methods are very similar, although the latter gives a smoother distribution of the filling factor. Maps obtained at two different epochs three weeks apart show a remarkable increase in spot coverage amounting to ~7% of the surface of the photosphere over a time period of only ~20 days. The spot maps of CD-641208 from the two methods show good longitudinal agreement, whereas the latitude range of the spots is extended to cover the whole visible hemisphere in the TR map. The distributions obtained from the first light curve of HD199143 show the presence of an extended and asymmetric active longitude with the maximum filling factor at longitude ~325degree.

V1647 Ori is a young eruptive variable star, illuminating a reflection nebula (McNeil's Nebula). It underwent an outburst in 2003 before fading back to its pre-outburst brightness in 2006. In 2008, V1647 Ori underwent a new outburst. The observed properties of the 2003-2006 event are different in several respects from both the EXor and FUor type outbursts, and suggest that this star might represent a new class of eruptive young stars, younger and more deeply embedded than EXors, and exhibiting variations on shorter time scales than FUors. In outburst, the star lights up the otherwise invisible McNeil's nebular - a conical cloud likely accumulated from previous outbursts. We present follow-up photometric as well as optical and near-IR spectroscopy of the nebula obtainted during the 2008-2009 outburst. We will also present results from contemporaneous X-ray observations. These multi-wavelength observations of V1647 Ori, obtained at this key early stage of the outburst, provide a snapshot of the "lighting up" of the nebula, probe its evolution through the event, and enable comparison with the 2003-2006 outburst.

We carry out a numerical simulation depicting the effects of a Coronal Mass Ejection (CME) event on a close-in giant planet in an extrasolar system. We drive the CME in a similar manner as in simulations of space weather events on Earth. The simulation includes the planetary orbital motion, which leads to the forming of a comet-like planetary magnetotail which is oriented almost perpendicular to the direction of propagation of the CME. The combination of this feature and the fact that the CME does not expand much by the time it reaches the planet leads to a unique CME-magnetosphere interaction, where the CME itself is highly affected by the presence of the planetary magnetosphere. We find that the planet is well-shielded from CME penetration, even for a relatively weak internal magnetic field. The planetary angular momentum loss associated with such an event is negligible compared to the total planetary angular momentum. We also find that the energy which is deposited in the magnetosphere is much higher than in the case of the Earth, and our simulation suggests there is a large-scale change in the orientation of the magnetosphere-ionosphere current system during the CME event.

The frequency of heating events in the corona is an important constraint on the coronal heating mechanisms. Observations indicate that the intensities and velocities measured in active region cores are effectively steady, suggesting that heating events occur rapidly enough to keep high temperature active region loops close to equilibrium. In this paper, we couple observations of Active Region 10955 made with XRT and EIS on \textit{Hinode} to test a simple steady heating model. First we calculate the differential emission measure of the apex region of the loops in the active region core. We find the DEM to be broad and peaked around 3\,MK. We then determine the densities in the corresponding footpoint regions. Using potential field extrapolations to approximate the loop lengths and the density-sensitive line ratios to infer the magnitude of the heating, we build a steady heating model for the active region core and find that we can match the general properties of the observed DEM for the temperature range of 6.3 $<$ Log T $<$ 6.7. This model, for the first time, accounts for the base pressure, loop length, and distribution of apex temperatures of the core loops. We find that the density-sensitive spectral line intensities and the bulk of the hot emission in the active region core are consistent with steady heating. We also find, however, that the steady heating model cannot address the emission observed at lower temperatures. This emission may be due to foreground or background structures, or may indicate that the heating in the core is more complicated. Different heating scenarios must be tested to determine if they have the same level of agreement.

The first release of the Chandra Source Catalog (CSC) contains ~95,000 X-ray sources in a total area of ~0.75% of the entire sky, using data from ~3,900 separate ACIS observations of a multitude of different types of X-ray sources. In order to maximize the scientific benefit of such a large, heterogeneous data-set, careful characterization of the statistical properties of the catalog, i.e., completeness, sensitivity, false source rate, and accuracy of source properties, is required. Characterization efforts of other, large Chandra catalogs, such as the ChaMP Point Source Catalog (Kim et al. 2007) or the 2 Mega-second Deep Field Surveys (Alexander et al. 2003), while informative, cannot serve this purpose, since the CSC analysis procedures are significantly different and the range of allowable data is much less restrictive. We describe here the characterization process for the CSC. This process includes both a comparison of real CSC results with those of other, deeper Chandra catalogs of the same targets and extensive simulations of blank-sky and point source populations.

While considerable advance has been made to account for statistical uncertainties in astronomical analyses, systematic instrumental uncertainties have been generally ignored. This can be crucial to a proper interpretation of analysis results because instrumental calibration uncertainty is a form of systematic uncertainty. Ignoring it can underestimate error bars and introduce bias into the fitted values of model parameters. Accounting for such uncertainties currently requires extensive case-specific simulations if using existing analysis packages. Here we present general statistical methods that incorporate calibration uncertainties into spectral analysis of high-energy data. We first present a method based on multiple imputation that can be applied with any fitting method, but is necessarily approximate. We then describe a more exact Bayesian approach that works in conjunction with a Markov chain Monte Carlo based fitting. We explore methods for improving computational efficiency, and in particular detail a method of summarizing calibration uncertainties with a principal component analysis of samples of plausible calibration files. This method is implemented using recently codified Chandra effective area uncertainties for low-resolution spectral analysis and is verified using both simulated and actual Chandra data. Our procedure for incorporating effective area uncertainty is easily generalized to other types of calibration uncertainties.

We carry out the first time-dependent numerical MagnetoHydroDynamic modeling of an extrasolar planetary system to study the interaction of the stellar magnetic field and wind with the planetary magnetosphere and outflow. We base our model on the parameters of the HD 189733 system, which harbors a close-in giant planet. Our simulation reveals a highly structured stellar corona characterized by sectors with different plasma properties. The star-planet interaction varies in magnitude and complexity, depending on the planetary phase, planetary magnetic field strength, and the relative orientation of the stellar and planetary fields. It also reveals a long, comet-like tail which is a result of the wrapping of the planetary magnetospheric tail by its fast orbital motion. A reconnection event occurs at a specific orbital phase, causing mass loss from the planetary magnetosphere that can generate a hot spot on the stellar surface. The simulation also shows that the system has sufficient energy to produce hot-spots observed in Ca II lines in giant planet hosting stars. However, the short duration of the reconnection event suggests that such SPI cannot be observed persistently.

We analyze coordinated Hinode XRT and EIS observations of a non-flaring active region to investigate the thermal properties of coronal plasma taking advantage of the complementary diagnostics provided by the two instruments. In particular we want to explore the presence of hot plasma in non-flaring regions. Independent temperature analyses from the XRT multi-filter dataset, and the EIS spectra, including the instrument entire wavelength range, provide a cross-check of the different temperature diagnostics techniques applicable to broad-band and spectral data respectively, and insights into cross-calibration of the two instruments. The emission measure distribution, EM(T), we derive from the two datasets have similar width and peak temperature, but show a systematic shift of the absolute values, the EIS EM(T) being smaller than XRT EM(T) by approximately a factor 2. We explore possible causes of this discrepancy, and we discuss the influence of the assumptions for the plasma element abundances. Specifically, we find that the disagreement between the results from the two instruments is significantly mitigated by assuming chemical composition closer to the solar photospheric composition rather than the often adopted "coronal" composition (Feldman 1992). We find that the data do not provide conclusive evidence on the high temperature (log T[K] >~ 6.5) tail of the plasma temperature distribution, however, suggesting its presence to a level in agreement with recent findings for other non-flaring regions.

We present the science case for a broadband X-ray imager with high-resolution spectroscopy, including simulations of X-ray spectral diagnostics of both active regions and solar flares. This is part of a trilogy of white papers discussing science, instrument (Bandler et al. 2010), and missions (Bookbinder et al. 2010) to exploit major advances recently made in transition-edge sensor (TES) detector technology that enable resolution better than 2 eV in an array that can handle high count rates. Combined with a modest X-ray mirror, this instrument would combine arcsecondscale imaging with high-resolution spectra over a field of view sufficiently large for the study of active regions and flares, enabling a wide range of studies such as the detection of microheating in active regions, ion-resolved velocity flows, and the presence of non-thermal electrons in hot plasmas. It would also enable more direct comparisons between solar and stellar soft X-ray spectra, a waveband in which (unusually) we currently have much better stellar data than we do of the Sun.

We present a numerical Magnetohydrodynamic (MHD) study of the dependence of stellar mass and angular momentum- loss rates on the orbital distance to close-in giant planets. We find that the mass loss rate drops by a factor of $\approx$1.5-2, while the angular momentum loss rate drops by a factor of $\approx$4 as the distance decreases past the Alfv\'en surface. This reduction in angular momentum loss is due to the interaction between the stellar and planetary Alfv\'en surfaces, which modifies the global structure of the stellar corona and stellar wind on the hemisphere facing the planet, as well as the opposite hemisphere. The simulation also shows that the magnitude of change in angular momentum loss rate depends mostly on the strength of the planetary magnetic field and not on its polarity. The interaction however, begins at greater separation if the overall field topology of the star and the planet are of anti-aligned. Our results are consistent with evidence for excess angular momentum in stars harboring close-in giant planets, and show that the reduction in wind-driven angular momentum loss can compete with, and perhaps dominate, spin-up due to tidal interaction.

We report on two XMM-Newton observations of the planetary host star HD189733. The system has a close in planet and it can potentially affect the coronal structure via interactions with the magnetosphere. We have obtained X-ray spectra and light curves from EPIC and RGS on board XMM-Newton which we have analyzed and interpreted. We reduced X-ray data from primary transit and secondary eclipse occurred in April 17th 2007 and May 18th 2009, respectively. In the April 2007 observation only variability due to weak flares is recognized. In 2009 HD189733 exhibited a X-ray flux always larger than in the 2007 observation. The average flux in 2009 was higher than in 2007 observation by a factor of 45%. During the 2009 secondary eclipse we observed a softening of the X-ray spectrum significant at level of ~3 sigma. Further, we observed the most intense flare recorded at either epochs. This flare occurred 3 ks after the end of the eclipse.The flare decay shows several minor ignitions perhaps linked to the main event and hinting for secondary loops that emit triggered by the main loop. Magneto-Hydro-Dynamical (MHD) simulations show that the magnetic interaction between planet and star enhances the density and the magnetic field in a region comprised between the planet and the star because of their relative orbital/rotation motion. X-ray observations and model predictions are globally found in agreement, despite the quite simple MHD model and the lack of precise estimate of parameters including the alignment and the intensity of stellar and planetary magnetic fields. Future observations should confirm or disprove this hypothesis, by determining whether flares are systematically recurring in the light curve at the same planetary phase.

We present the results of a Spitzer photometric investigation into the IR excesses of close binary systems. In a sample of 10 objects, excesses in IRAC and MIPS24 bands implying the presence of warm dust are found for 3. For 2 objects we do not find excesses reported in earlier IRAS studies. We discuss the results in the context of the scenario suggested by Rhee and co-workers, in which warm dust is continuously created by destructive collisions between planetary bodies. A simple numerical model for the steady-state distribution of dust in one IR excess system shows a central clearing of radius 0.22 AU caused by dynamical perturbations from the binary star. This is consistent with the size of the central clearing derived from the Spitzer spectral energy distribution. We conclude that close binaries could be efficient "destroyers of worlds", lead to destabilize the orbits of their planetary progeny by magnetically-driven angular momentum loss and secular shrinkage of the binary separation.

We perform a numerical simulation of the corona of the young, rapidly rotating K0 dwarf AB Doradus using a global MHD model. The model is driven by a surface map of the radial magnetic field constructed using Zeeman-Doppler Imaging. We find that the global structure of the stellar corona is dominated by strong azimuthal tangling of the magnetic field due to the rapid rotation. The MHD solution enables us to calculate realistic Alfv\'en surfaces and we can therefore estimate the stellar mass loss rate and angular momentum loss rate without making undue theoretical simplifications. We consider three cases, parametrized by the base density of the corona, that span the range of possible solutions for the system. We find that overall, the mass and angular-momentum loss rates are higher than in the solar case; the mass loss rates are 10 to 500 times higher, and the angular momentum loss rate can be up to $3\times{10}^4$ higher than present day solar values. Our simulations show that this model can be use to constrain the wide parameter space of stellar systems. It also shows that an MHD approach can provide more information about the physical system over the commonly used potential field extrapolation.

A common problem in astrophysics is determining how bright a source could be and still not be detected. Despite the simplicity with which the problem can be stated, the solution involves complex statistical issues that require careful analysis. In contrast to the confidence bound, this concept has never been formally analyzed, leading to a great variety of often ad hoc solutions. Here we formulate and describe the problem in a self-consistent manner. Detection significance is usually defined by the acceptable proportion of false positives (the TypeI error), and we invoke the complementary concept of false negatives (the TypeII error), based on the statistical power of a test, to compute an upper limit to the detectable source intensity. To determine the minimum intensity that a source must have for it to be detected, we first define a detection threshold, and then compute the probabilities of detecting sources of various intensities at the given threshold. The intensity that corresponds to the specified TypeII error probability defines that minimum intensity, and is identified as the upper limit. Thus, an upper limit is a characteristic of the detection procedure rather than the strength of any particular source and should not be confused with confidence intervals or other estimates of source intensity. This is particularly important given the large number of catalogs that are being generated from increasingly sensitive surveys. We discuss the differences between these upper limits and confidence bounds. Both measures are useful quantities that should be reported in order to extract the most science from catalogs, though they answer different statistical questions: an upper bound describes an inference range on the source intensity, while an upper limit calibrates the detection process. We provide a recipe for computing upper limits that applies to all detection algorithms.

We present a three-dimensional simulation of the corona of an FK Com-type rapidly rotating G giant using a magnetohydrodynamic model that was originally developed for the solar corona in order to capture the more realistic, non-potential coronal structure. We drive the simulation with surface maps for the radial magnetic field obtained from a stellar dynamo model of the FK Com system. This enables us to obtain the coronal structure for different field topologies representing different periods of time. We find that the corona of such an FK Com-like star, including the large scale coronal loops, is dominated by a strong toroidal component of the magnetic field. This is a result of part of the field being dragged by the radial outflow, while the other part remains attached to the rapidly rotating stellar surface. This tangling of the magnetic field,in addition to a reduction in the radial flow component, leads to a flattening of the gas density profile with distance in the inner part of the corona. The three-dimensional simulation provides a global view of the coronal structure. Some aspects of the results, such as the toroidal wrapping of the magnetic field, should also be applicable to coronae on fast rotators in general, which our study shows can be considerably different from the well-studied and well-observed solar corona. Studying the global structure of such coronae should also lead to a better understanding of their related stellar processes, such as flares and coronal mass ejections, and in particular, should lead to an improved understanding of mass and angular momentum loss from such systems.

The Chandra Source Catalog (CSC) is a general purpose virtual X-ray astrophysics facility that provides access to a carefully selected set of generally useful quantities for individual X-ray sources, and is designed to satisfy the needs of a broad-based group of scientists, including those who may be less familiar with astronomical data analysis in the X-ray regime. The first release of the CSC includes information about 94,676 distinct X-ray sources detected in a subset of public ACIS imaging observations from roughly the first eight years of the Chandra mission. This release of the catalog includes point and compact sources with observed spatial extents <~ 30''. The catalog (1) provides access to the best estimates of the X-ray source properties for detected sources, with good scientific fidelity, and directly supports scientific analysis using the individual source data; (2) facilitates analysis of a wide range of statistical properties for classes of X-ray sources; and (3) provides efficient access to calibrated observational data and ancillary data products for individual X-ray sources, so that users can perform detailed further analysis using existing tools. The catalog includes real X-ray sources detected with flux estimates that are at least 3 times their estimated 1 sigma uncertainties in at least one energy band, while maintaining the number of spurious sources at a level of <~ 1 false source per field for a 100 ks observation. For each detected source, the CSC provides commonly tabulated quantities, including source position, extent, multi-band fluxes, hardness ratios, and variability statistics, derived from the observations in which the source is detected. In addition to these traditional catalog elements, for each X-ray source the CSC includes an extensive set of file-based data products that can be manipulated interactively.

Since the first discovery of an extrasolar planetary system more than a decade ago, hundreds more have been discovered. Surprisingly, many of these systems harbor Jupiter-class gas giants located close to the central star, at distances of 0.1 AU or less. Observations of chromospheric 'hot spots' that rotate in phase with the planetary orbit, and elevated stellar X-ray luminosities,suggest that these close-in planets significantly affect the structure of the outer atmosphere of the star through interactions between the stellar magnetic field and the planetary magnetosphere. Here we carry out the first detailed three-dimensional MagnetoHydroHynamics (MHD) simulation containing the two magnetic bodies and explore the consequences of such interactions on the steady-state coronal structure. The simulations reproduce the observable features of 1) increase in the total X-ray luminosity, 2) appearance of coronal hot spots, and 3) phase shift of these spots with respect to the direction of the planet. The proximate cause of these is an increase in the density of coronal plasma in the direction of the planet, which prevents the corona from expanding and leaking away this plasma via a stellar wind. The simulations produce significant low temperature heating. By including dynamical effects, such as the planetary orbital motion, the simulation should better reproduce the observed coronal heating.

We simulate the effect of latitudinal variations in the location of star spots, as well as their magnetic field strength, on stellar angular momentum loss to the stellar wind. We use the Michigan solar corona global MagnetoHydroDynamic model, which incorporates realistic relation between the magnetic field topology and the wind distribution. We find that the spots location significantly affects the stellar wind structure, and as a result, the total mass loss rate and angular momentum loss rate. In particular, we find that the angular momentum loss rate is controlled by the mass flux when spots are located at low latitudes but is controlled by an increased plasma density between the stellar surface and the Alfven surface when spots are located at high latitudes. Our results suggest that there might be a feedback mechanism between the magnetic field distribution, wind distribution, angular momentum loss through the wind, and the motions at the convection zone that generate the magnetic field. This feedback might explain the role of coronal magnetic fields in stellar dynamos.

The X-Ray Telescope (XRT) on the Japanese/USA/UK {\it Hinode (Solar-B)} spacecraft has detected emission from a quiescent active region core that is consistent with nanoflare heating. The fluxes from 10 broadband X-ray filters and filter combinations were used to constructed Differential Emission Measure (DEM) curves. In addition to the expected active region peak at Log T = 6.3-6.5, we find a high-temperature component with significant emission measure at Log T $>$ 7.0. This emission measure is weak compared to the main peak -- the DEM is down by almost three orders of magnitude -- which accounts of the fact that it has not been observed with earlier instruments. It is also consistent with spectra of quiescent active regions: no Fe XIX lines are observed in a CHIANTI synthetic spectrum generated using the XRT DEM distribution. The DEM result is successfully reproduced with a simple two-component nanoflare model.

We study the spectral energy distributions and evolution of a large sample of optically selected quasars from the Sloan Digital Sky Survey (SDSS) that were observed in 323 Chandra images analyzed by the Chandra Multiwavelength Project (ChaMP). Our highest-confidence matched sample includes 1135 X-ray detected quasars in the redshift range 0.2<z<5.4, representing some 36Msec of effective exposure. Spectroscopic redshifts are available for about 1/3 of the detected sample; elsewhere, redshifts are estimated photometrically. With 56 z>3 QSOs detected, we find no evidence for evolution out to z~5 for either the X-ray photon index Gamma or for the ratio of optical/UV to X-ray flux alpha_ox. About 10% of detected QSOs are obscured (Nh>1E22), but the fraction might reach ~1/3 if most non-detections are absorbed. We confirm a significant correlation between alpha_ox and optical luminosity, but it flattens or disappears for fainter AGN alone. Gamma hardens significantly both towards higher X-ray luminosity, and for relatively X-ray loud quasars. These trends may represent a relative increase in non-thermal X-ray emission, and our findings thereby strengthen analogies between Galactic black hole binaries and AGN.

We have carried out a survey of X-ray emission from stars with giant planets, combining both archival and targeted surveys. Over 230 stars have been currently identified as possessing planets, and roughly a third of these have been detected in X-rays. We carry out detailed statistical analysis on a volume limited sample of main sequence star systems with detected planets, comparing subsamples of stars that have close-in planets with stars that have more distant planets. This analysis reveals strong evidence that stars with close-in giant planets are on average more X-ray active by a factor ~4 than those with planets that are more distant. This result persists for various sample selections. We find that even after accounting for observational sample bias, a significant residual difference still remains. This observational result is consistent with the hypothesis that giant planets in close proximity to the primary stars influences the stellar magnetic activity.

We present 348 X-ray emitting stars identified from correlating the Extended Chandra Multiwavelength Project (ChaMP), a serendipitous wide-area X-ray survey, with the Sloan Digital Sky Survey (SDSS). We use morphological star/galaxy separation, an SDSS quasar catalog, an optical color-magnitude cut, and X-ray data quality tests to create our catalog, the ChaMP Extended Stellar Survey (ChESS), from a sample of 2121 matched ChaMP/SDSS sources. Our cuts retain 92% of the spectroscopically confirmed stars while excluding 99.6% of the 684 spectoscopically confirmed extragalactic sources. Fewer than 3% of the sources in our final catalog are previously identified stellar X-ray emitters; we expect ~10% of the catalog is composed by giants, and identify seven giant stars and three cataclysmic variables. We derive distances, X-ray and bolometric luminosities for these stars, revealing that this catalog fills the gap between the nearby stars identified by the ROSAT All-Sky Survey and the more distant stars detected in deep Chandra and XMM-Newton surveys. For 36 newly identified X-ray emitting M stars we calculate L_(Halpha)/L_(bol). L_(H\alpha)/L_(bol) and L_(X)/L_(bol) are linearly related below L_(X)/L_(bol) ~ 3 x 10^(-4), while L_(Halpha)/L_(bol) appears to turn over at larger L_(X)/L_(bol) values. Stars with reliable SDSS photometry have an ~0.1 mag blue excess in (u-g), likely due to increased chromospheric continuum emission. Photometric metallicity estimates suggest the sample is split between the young and old disk populations of the Galaxy; the lowest activity sources belong to the old disk population, a clear signature of the decay of magnetic activity with age. Future papers will present analyses of source variability and comparisons of this catalog to models of stellar activity in the Galactic disk.

We present the first observation of planet-induced stellar X-ray activity, identified for the HD 179949 system, using Chandra / ACIS-S. The HD 179949 system consists of a close-in giant planet orbiting an F9V star. Previous ground-based observations already showed enhancements in Ca II K in phase with the planetary orbit. We find an ~30% increase in the X-ray flux over quiescent levels coincident with the phase of the Ca II enhancements. There is also a trend for the emission to be hotter at increased fluxes, confirmed by modeling, showing the enhancement at ~1 keV compared to ~0.4 keV for the background star.

AB Dor, Speedy Mic and Rst137B are in their early post-T Tauri evolutionary phase (<100Myr), at the age of fastest rotation in the life of late-type stars. They straddle the coronal saturation-supersaturation boundary first defined by young stars in open clusters. High resolution Chandra X-ray spectra have been analysed to study their coronal properties as a function of coronal activity parameters Rossby number, $L_X/L_{bol}$ and a coronal temperature index. Plasma emission measure distributions as a function of temperature show broad peaks at T~10e7K. Differences between stars suggest that as supersaturation is reached the DEM slope below the temperature of peak DEM becomes shallower, while the DEM drop-off above this temperature becomes more pronounced. A larger sample comprising our three targets and 22 active stars studied in the recent literature reveals a general increase of plasma at T>10e7 toward the saturated-supersaturated boundary but a decline beyond this among supersaturated stars. All three of the stars studied in detail here show lower coronal abundances of the low FIP elements Mg, Si and Fe, relative to the high FIP elements S, O and Ne, as compared to the solar mixture. The coronal Fe abundances of the stellar sample are inversely correlated with Lx/Lbol, declining slowly with rising Lx/Lbol, but with a much more sharp decline at Lx/Lbol>3x10e-4. For dwarfs the Fe abundance is also well-correlated with Rossby number. The coronal O/Fe ratios for dwarfs show a clear increase with decreasing Rossby number, apparently reaching saturation at [O/Fe]=0.5 at the coronal supersaturation boundary. Similar increases in O/Fe with increasing coronal temperature and $L_X/L_{bol}$ are seen.

We present results from two Chandra imaging observations of Ross 154, a nearby flaring M dwarf star. During a 61-ks ACIS-S exposure, a very large flare occurred (the equivalent of a solar X3400 event, with L_X = 1.8x10^30 ergs/s) in which the count rate increased by a factor of over 100. The early phase of the flare shows evidence for the Neupert effect, followed by a further rise and then a two-component exponential decay. A large flare was also observed at the end of a later 48-ks HRC-I observation. Emission from the non-flaring phases of both observations was analyzed for evidence of low level flaring. From these temporal studies we find that microflaring probably accounts for most of the `quiescent' emission, and that, unlike for the Sun and the handful of other stars that have been studied, the distribution of flare intensities does not appear to follow a power-law with a single index. Analysis of the ACIS spectra, which was complicated by exclusion of the heavily piled-up source core, suggests that the quiescent Ne/O abundance ratio is enhanced by a factor of ~2.5 compared to the commonly adopted solar abundance ratio, and that the Ne/O ratio and overall coronal metallicity during the flare appear to be enhanced relative to quiescent abundances. Based on the temperatures and emission measures derived from the spectral fits, we estimate the length scales and plasma densities in the flaring volume and also track the evolution of the flare in color-intensity space. Lastly, we searched for a stellar-wind charge-exchange X-ray halo around the star but without success; because of the relationship between mass-loss rate and the halo surface brightness, not even an upper limit on the stellar mass-loss rate can be determined.

We have analyzed Chandra calibration observations of Betelgeuse (alpha Ori, M2Iab, m_V=0.58, 131 pc) obtained at the aimpoint locations of the HRC-I (8 ks), HRC-S (8 ks), and ACIS-I (5 ks). Betelgeuse is undetected in all the individual observations as well as cumulatively. We derive upper limits to the X-ray count rates and compute the corresponding X-ray flux and luminosity upper limits for coronal plasma that may potentially exist in the atmosphere of Betelgeuse over a range of temperatures, T=0.3-10 MK. We place a flux limit at the telescope of fx ~ 4x10^(-15) ergs s^(-1) cm^(-2) at T=1 MK. The upper limit is lowered by a factor of ~3 at higher temperatures, roughly an order of magnitude lower than that obtained previously. Assuming that the entire stellar surface is active, these fluxes correspond to a surface flux limit that ranges from 30-7000 ergs s^(-1) cm^(-2) at T=1 MK, to ~1 ergs s^(-1) cm^(-2) at higher temperatures, five orders of magnitude below the quiet Sun X-ray surface flux. We discuss the implications of our analysis in the context of models of a buried corona and a pervasive magnetic carpet. We rule out the existence of a solar-like corona on Betelgeuse, but cannot rule out the presence of low-level emission on the scale of coronal holes.

The high-resolution X-ray spectroscopy made possible by the 1999 deployment of the Chandra X-ray Observatory has revolutionized our understanding of stellar X-ray emission. Many puzzles remain, though, particularly regarding the mechanisms of X-ray emission from OB stars. Although numerous individual stars have been observed in high-resolution, realizing the full scientific potential of these observations will necessitate studying the high-resolution Chandra dataset as a whole. To facilitate the rapid comparison and characterization of stellar spectra, we have compiled a uniformly processed database of all stars observed with the Chandra High Energy Transmission Grating (HETG). This database, known as X-Atlas, is accessible through a web interface with searching, data retrieval, and interactive plotting capabilities. For each target, X-Atlas also features predictions of the low-resolution ACIS spectra convolved from the HETG data for comparison with stellar sources in archival ACIS images. Preliminary analyses of the hardness ratios, quantiles, and spectral fits derived from the predicted ACIS spectra reveal systematic differences between the high-mass and low-mass stars in the atlas and offer evidence for at least two distinct classes of high-mass stars. A high degree of X-ray variability is also seen in both high and low-mass stars, including Capella, long thought to exhibit minimal variability. X-Atlas contains over 130 observations of approximately 25 high-mass stars and 40 low-mass stars and will be updated as additional stellar HETG observations become public. The atlas has recently expanded to non-stellar point sources, and Low Energy Transmission Grating (LETG) observations are currently being added as well.

We analyze 205 ks of imaging data of the active binary, Capella, obtained with the Chandra High Resolution Camera Imager (HRC-I) to determine whether Capella shows any variability at timescales < 50 ks. We find that a clear signal for variability is present for timescales < 20 ks, and that the light curves show evidence for excess fluctuation over that expected from a purely Poisson process. This overdispersion is consistent with variability at the 2-7% level, and suggests that the coronae on the binary components of Capella are composed of low-density plasma and low-lying loops.

We present the Chandra Multiwavelength Project (ChaMP) X-ray point source catalog with ~6,800 X-ray sources detected in 149 Chandra observations covering \~10 deg^2. The full ChaMP catalog sample is seven times larger than the initial published ChaMP catalog. The exposure time of the fields in our sample ranges from 0.9 to 124 ksec, corresponding to a deepest X-ray flux limit of f_{0.5-8.0} = 9 x 10^{-16} erg/cm2/sec. The ChaMP X-ray data have been uniformly reduced and analyzed with ChaMP-specific pipelines, and then carefully validated by visual inspection. The ChaMP catalog includes X-ray photometric data in 8 different energy bands as well as X-ray spectral hardness ratios and colors. To best utilize the ChaMP catalog, we also present the source reliability, detection probability and positional uncertainty. To quantitatively assess those parameters, we performed extensive simulations. In particular, we present a set of empirical equations: the flux limit as a function of effective exposure time, and the positional uncertainty as a function of source counts and off axis angle. The false source detection rate is ~1% of all detected ChaMP sources, while the detection probability is better than ~95% for sources with counts >30 and off axis angle <5 arcmin. The typical positional offset between ChaMP X-ray source and their SDSS optical counterparts is 0.7+-0.4 arcsec, derived from ~900 matched sources.

A commonly used measure to summarize the nature of a photon spectrum is the so-called Hardness Ratio, which compares the number of counts observed in different passbands. The hardness ratio is especially useful to distinguish between and categorize weak sources as a proxy for detailed spectral fitting. However, in this regime classical methods of error propagation fail, and the estimates of spectral hardness become unreliable. Here we develop a rigorous statistical treatment of hardness ratios that properly deals with detected photons as independent Poisson random variables and correctly deals with the non-Gaussian nature of the error propagation. The method is Bayesian in nature, and thus can be generalized to carry out a multitude of source-population--based analyses. We verify our method with simulation studies, and compare it with the classical method. We apply this method to real world examples, such as the identification of candidate quiescent Low-mass X-ray binaries in globular clusters, and tracking the time evolution of a flare on a low-mass star.

(abridged) We have investigated 136 Chandra extragalactic sources without broad optical emission lines, including 93 galaxies with narrow emission lines (NELG) and 43 with only absorption lines (ALG). Based on fx/fo, Lx, X-ray spectral hardness and optical emission line diagnostics, we have conservatively classified 36 normal galaxies (20 spirals and 16 ellipticals) and 71 AGNs. We found no statistically significant evolution in Lx/LB, within the limited z range. We have built log(N)-log(S), after correcting for completeness based on a series of simulations. The best-fit slope is -1.5 for both S and B energy bands, which is considerably steeper than that of the AGN-dominated cosmic background sources, but slightly flatter than the previous estimate, indicating normal galaxies will not exceed the AGN population until fx ~ 2 x 10-18 erg s-1 cm-2 (a factor of ~5 lower than the previous estimate). A group of NELGs appear to be heavily obscured in X-rays, i.e., a typical type 2 AGN. After correcting for intrinsic absorption, their X-ray luminosities could be Lx > 10^44 erg s-1, making them type 2 quasar candidates. While most X-ray luminous ALGs (XBONG - X-ray bright, optically normal galaxy candidates) do not appear to be significantly absorbed, we found two heavily obscured objects, which could be as luminous as an unobscured broad-line quasar. Among 43 ALGs, we found two E+A galaxy candidates with strong Balmer absorption lines, but no [OII] line. The X-ray spectra of both galaxies are soft and one of them has a nearby close companion galaxy, supporting the merger/interaction scenario rather than the dusty starburst hypothesis.

$\beta$ Cet, 31 Com and $\mu$ Vel represent the main stages through which late-type giants evolve during their lifetime (the Hertzsprung gap (31 Com), the rapid braking zone ($\mu$ Vel) and the core helium burning ``clump'' phase ($\beta$ Cet)). An analysis of their high resolution {\it Chandra} X-ray spectra reveals similar coronal characteristics in terms of both temperature structure and element abundances for the more evolved stars ($\mu$ Vel and $\beta$ Cet) with slight differences for the `younger' giant (31 Com). The coronal temperature structure of 31 Com is significantly hotter showing a clear peak while $\beta$ Cet and $\mu$ Vel show a plateau. $\beta$ Cet and $\mu$ Vel show evidence for a FIP effect in which coronae are depleted in high FIP elements relative to their photospheres by a factor of $\sim 2$. In contrast, 31 Com is characterized by a lack of FIP effect. In other words, neither depletion nor enhancement relative to stellar photospheric values is found. We conclude that the structural changes during the evolution of late-type giants could be responsible for the observed differences in coronal abundances and temperature structure. In particular, the size of the convection zone coupled with the rotation rate seem obvious choices for playing a key role in determining coronal characteristics.

Observations of Fe XVIII and Fe XIX X-ray, EUV, and FUV line emission, formed at the peak of Capella's (alpha Aurigae) emission measure distribution and ubiquitous in spectra of many cool stars and galaxies, provide a unique opportunity to test the robustness of Fe XVIII and Fe XIX spectral models. The Astrophysical Plasma Emission Code (APEC) is used to identify over 35 lines from these two ions alone, and to compare model predictions with spectra obtained with the Chandra Low Energy Transmission Grating and High Energy Transmission Grating Spectrometers, the Far Ultraviolet Spectroscopic Explorer, and the Extreme Ultraviolet Explorer. Some flux discrepancies larger than factors of two are found between observations of Fe XVIII and Fe XIX lines and predictions by APEC and other models in common usage. In particular the X-ray resonance lines for both ions are stronger than predicted by all models relative to the EUV resonance lines. The multiwavelength observations demonstrate the importance of including dielectronic recombination and proton impact excitation, and of using accurate wavelengths in spectral codes. These ions provide important diagnostic tools for 10^7 K plasmas currently observed with Chandra, XMM-Newton, and FUSE.

The zero-age main-sequence star AB Dor and the K dwarf component of the V471 Tau close binary have essentially identical rotation rates and spectral types. An analysis of their high resolution {\it Chandra} X-ray spectra reveals remarkably similar coronal characteristics in terms of both temperature structure and element abundances. Both stars show depletions of low FIP elements by factors of $\sim 3$, with higher FIP elements showing more mild depletions. No evidence for enhancements of very low FIP ($< 7$ eV) elements, such as Na, Al and Ca, as compared to other low FIP elements was found. The abundance anomaly pattern for AB Dor and V471 Tau is similar to, though less extreme than, the abundance anomalies exhibited by active RS-CVn-type binaries. While we find statistically significant structure in the underlying differential emission measure distributions of these stars over narrow temperature intervals, this structure is strongly dependent on the lines used in the analysis and is probably spurious. Based on their X-ray similarities, we conclude that the exact evolutionary state of a star has little effect on coronal characteristics, and that the parameters that dominate coronal structure and composition are simply the rotation rate and spectral type.

We present the analysis and interpretation of an observation of the flare star AD Leo (dM3e) with Chandra/LETG. The high resolution X-ray spectrum allowed us to infer the plasma emission measure distribution (EMD) vs. temperature, and the abundances of individual elements in the corona of this magnetically-active star, during a typical state characterized by significant variability but no large flaring event. We have also measured plasma densities at various temperatures using spectroscopic diagnostics provided by He-like triplets and Fe XXI lines. We show that the present EMD is similar to those previously obtained from EUVE spectra during quiescent and flaring states, confirming the long-term stability of the corona of AD Leo. At variance with the case of other active stars, the EMD of AD Leo is characterized by a significantly shallower slope, compatible with that predicted by static models of isobaric loops with constant cross-section and uniform heating. We discuss such coronal modeling to infer the average properties of the corona in terms of loop populations, including estimates of the surface filling factor. We also show that the EMD is compatible with the model of a corona continuously heated by flares. The coronal composition is such that the element abundances, relative to solar values, tend to increase with the First Ionization Potential, with few exceptions. The line to continuum ratios suggest a nearly solar coronal metallicity.

In a study of Chandra High Energy Transmission Grating spectra of Algol, we clearly detect Doppler shifts caused by the orbital motion of Algol B. These data provide the first definitive proof that the X-ray emission of Algol is dominated by the secondary, in concordance with expectations that Algol A (B8) is X-ray dark. The measured Doppler shifts are slightly smaller than expected, implying an effective orbital radius of about 10 Rsolar, instead of 11.5 Rsolar for the Algol B center of mass. This could be caused by a small contribution of X-ray flux from Algol A (10-15%), possibly through accretion. The more likely explanation is an asymmetric corona biased toward the system center of mass by the tidal distortion of the surface of Algol B. Analysis of the strongest lines indicates excess line broadening of ~150 km/s above that expected from thermal motion and surface rotation. Possible explanations include turbulence, flows or explosive events, or rotational broadening from a radially extended corona. We favor the latter scenario and infer that a significant component of the corona at temperatures <10^7 K has a scale height of order the stellar radius. This is supported by the shape of the X-ray lightcurve and the shallow dip at secondary eclipse. We also examine the O VII intercombination and forbidden lines in a Low Energy Transmission Grating Spectrograph observation and find no change in their relative line fluxes as the system goes from quadrature to primary eclipse. Since these lines are strongly affected by UV irradiation from Algol A, this supports the conjecture that the corona of Algol B at temperatures of several million K must be significantly extended and/or located toward the poles to avoid being shadowed from Algol A during primary eclipse.

The Chandra Multi-wavelength Project (ChaMP) is a wide-area (~14 deg^2) survey of serendipitous Chandra X-ray sources, aiming to establish fair statistical samples covering a wide range of characteristics (such as absorbed AGNs, high z clusters of galaxies) at flux levels (fX ~ 10^-15 - 10^-14 erg sec-1 cm-2) intermediate between the Chandra Deep surveys and previous missions. We present the first ChaMP catalog, which consists of 991 near on-axis, bright X-ray sources obtained from the initial sample of 62 observations. The data have been uniformly reduced and analyzed with techniques specifically developed for the ChaMP and then validated by visual examination. To assess source reliability and positional uncertainty, we perform a series of simulations and also use Chandra data to complement the simulation study. The false source detection rate is found to be as good as or better than expected for a given limiting threshold. On the other hand, the chance of missing a real source is rather complex, depending on the source counts, off-axis distance (or PSF), and background rate. The positional error (95% confidence level) is usually < 1" for a bright source, regardless of its off-axis distance while it can be as large as 4" for a weak source (~20 counts) at a large off-axis distance (Doff-axis > 8'). We have also developed new methods to find spatially extended or temporary variable sources and those sources are listed in the catalog.

We present the first results of ChaMP X-ray source properties obtained from the initial sample of 62 observations. The data have been uniformly reduced and analyzed with techniques specifically developed for the ChaMP and then validated by visual examination. Utilizing only near on-axis, bright X-ray sources (to avoid problems caused by incompleteness and the Eddington bias), we derive the Log(N)-Log(S) relation in soft (0.5-2 keV) and hard (2-8 keV) energy bands. The ChaMP data are consistent with previous results of ROSAT, ASCA and Chandra deep surveys. In particular, our data nicely fill in the flux gap in the hard band between the Chandra Deep Field data and the previous ASCA data. We check whether there is any systematic difference in the source density between cluster and non-cluster fields and also search for field-to-field variations, both of which have been previously reported. We found no significant field-to-field cosmic variation in either test within the statistics (~1 sigma) across the flux levels included in our sample. In the X-ray color-color plot, most sources fall in the location characterized by photon index = 1.5-2 and NH = a few x 10^20 cm^2, suggesting that they are typical broad-line AGNs. There also exist a considerable number of sources with peculiar X-ray colors (e.g., highly absorbed, very hard, very soft). We confirm a trend that on average the X-ray color hardens as the count rate decreases. Since the hardening is confined to the softest energy band (0.3-0.9 keV), we conclude it is most likely due to absorption. We cross-correlate the X-ray sources with other catalogs and describe their properties in terms of optical color, X-ray-to-optical luminosity ratio and X-ray colors.

(Abridged) We investigate the EUV and X-ray flare rate distribution in radiated energy of the late-type active star AD Leo using long EUVE and BeppoSAX observations. We compare the observed light curves with light curves synthesized from model flares that are distributed in energy according to a power law with selectable index alpha (dN/dE ~ E^{-alpha}). Two methods are applied, the first comparing flux distributions of the binned data, and the second using the distributions of photon arrival time differences in the unbinned data. We find acceptable alpha values between 2.0-2.5 for the EUVE DS and the BeppoSAX LECS data. The BeppoSAX MECS data indicate a somewhat shallower energy distribution than the LECS data, which is attributed to the harder range of sensitivity of the MECS detector and the increasing peak temperatures of flares with increasing total (radiative) energy. The results suggest that flares can play an important role in the energy release of this active corona. We discuss caveats related to time variability, total energy, and multiple power-law distributions. Studying the limiting case of a corona that is entirely heated by a population of flares, we derive an expression for the time-averaged coronal differential emission measure distribution (DEM) that can be used as a diagnostic for the flare energy distribution. The shape of the analytical DEM agrees with previously published DEMs from observations of active stars.

The Chandra Multiwavength Plane (ChaMPlane) Survey of the galactic plane incorporates serendipitous sources from selected Chandra pointings in or near the galactic plane (b < 12deg; >20 ksec; lack of bright diffuse or point sources) to measure or constrain the luminosity function of low-luminosity accretion sources in the Galaxy. The primary goal is to detect and identify accreting white dwarfs (cataclysmic variables, with space density still uncertain by a factor of >10-100), neutron stars and black holes (quiescent low mass X-ray binaries) to constrain their space densities and thus origin and evolution. Secondary objectives are to identify Be stars in high mass X-ray binaries and constrain their space densities, and to survey the H-R diagram for stellar coronal sources. A parallel optical imaging under the NOAO Long Term Survey program provides deep optical images using the Mosaic imager on the CTIO and KPNO 4-m telescopes. The 36arcmin X 36arcmin optical images (Halpha, R, V and I) cover ~5X the area of each enclosed Chandra ACIS FOV, providing an extended survey of emission line objects for comparison with Chandra. Spectroscopic followup of optical counterparts is then conducted, thus far with WIYN and Magellan. The X-ray preliminary results from both the Chandra and optical surveys will be presented, including logN-logS vs. galactic position (l,b) and optical idenifications.

We investigate the contribution of very weak flares to the coronal luminosity of low-mass active stars. We analyze EUVE/DS events data from FK Aqr, V1054 Oph, and AD Leo and conclude that in all these cases the coronal emission is dominated by flares to such an extent that in some cases the entire emission may be ascribed to flare heating. We have developed a new method to directly model for the first time stochastically produced flare emission, including undetectable flares, and their effects on the observed photon arrival times. We find that the index of the power-law distribution of flare energies (dN/dE ~ E^{-alpha}) is 2.6+-0.34, 2.74+-0.35, and 2.03-2.32 for FK Aqr, V1054 Oph, and AD Leo respectively. We also find that the flare component accounts for a large fraction (generally >50 percent) of the total flux.

Deep Chandra LETG+HRC-S observations of the isolated neutron star candidate RX J1856.5-3754 have been analysed to search for metallic and resonance cyclotron spectral features and for pulsation behaviour. As found from earlier observations, the X-ray spectrum is well-represented by a ~ 60 eV (7e5 K) blackbody. No unequivocal evidence of spectral line or edge features has been found, arguing against metal-dominated models. The data contain no evidence for pulsation and we place a 99% confidence upper limit of 2.7% on the unaccelerated pulse fraction over a wide frequency range from 1e-4 to 100 Hz. We argue that the derived interstellar medium neutral hydrogen column density of 8e19 <= N_H <= 1.1e20 per sq. cm favours the larger distance from two recent HST parallax analyses, placing RX J1856.5-3754 at ~ 140 pc instead of ~ 60 pc, and in the outskirts of the R CrA dark molecular cloud. That such a comparatively rare region of high ISM density is precisely where an isolated neutron star re-heated by accretion of interstellar matter would be expected is either entirely coincidental, or current theoretical arguments excluding this scenario for RX J1856.5-3754 are premature. Taken at face value, the combined observational evidence -- a lack of spectral and temporal features and an implied radius at infinity of 3.8-8.2 km that is too small for current neutron star models -- points to a more compact object, such as allowed for quark matter equations of state.

Metallicity can be measured by analyzing the spectra in the X-ray region and comparing the flux in spectral lines to the flux in the underlying Bremsstrahlung continuum. In this paper we propose new Bayesian methods which directly model the Poisson nature of the data and thus are expected to exhibit improved sampling properties. Our model also accounts for the Poisson nature of background contamination of the observations, image blurring due to instrument response, and the absorption of photons in space. The resulting highly structured hierarchical model is fit using the Gibbs sampler, data augmentation and Metropolis-Hasting. We demonstrate our methods with the X-ray spectral analysis of several "Metal Abundance Deficient" stars. The model is designed to summarize the relative frequency of the energy of photons (X-ray or gamma-ray) arriving at a detector. Independent Poisson distributions are more appropriate to model the counts than the commonly used normal approximation. We model the high energy tail of the ASCA spectrum of each of the stars as a combination of a Bremsstrahlung continuum and ten narrow emission lines, included at positions of known strong lines. Statistical analysis is based on two source observations and one background observation. We use sequential Bayesian analysis for the two source observations; the posterior distribution from the first analysis is used to construct a prior for the second. Sensitivity of the final results to the choise of prior is investigated by altering the prior.

The likelihood ratio test (LRT) and the related $F$ test, do not (even asymptotically) adhere to their nominal $\chi^2$ and $F$ distributions in many statistical tests common in astrophysics, thereby casting many marginal line or source detections and non-detections into doubt. Although there are many legitimate uses of these statistics, in some important cases it can be impossible to compute the correct false positive rate. For example, it has become common practice to use the LRT or the $F$ test for detecting a line in a spectral model or a source above background despite the lack of certain required regularity conditions. In these and other settings that involve testing a hypothesis that is on the boundary of the parameter space, {\it contrary to common practice, the nominal $\chi^2$ distribution for the LRT or the $F$ distribution for the $F$ test should not be used}. In this paper, we characterize an important class of problems where the LRT and the $F$ test fail and illustrate this non-standard behavior. We briefly sketch several possible acceptable alternatives, focusing on Bayesian posterior predictive probability-values. We present this method in some detail, as it is a simple, robust, and intuitive approach. This alternative method is illustrated using the gamma-ray burst of May 8, 1997 (GRB 970508) to investigate the presence of an Fe K emission line during the initial phase of the observation.

Wavelets are scaleable, oscillatory functions that deviate from zero only within a limited spatial regime and have average value zero. In addition to their use as source characterizers, wavelet functions are rapidly gaining currency within the source detection field. Wavelet-based source detection involves the correlation of scaled wavelet functions with binned, two-dimensional image data. If the chosen wavelet function exhibits the property of vanishing moments, significantly non-zero correlation coefficients will be observed only where there are high-order variations in the data; e.g., they will be observed in the vicinity of sources. In this paper, we describe the mission-independent, wavelet-based source detection algorithm WAVDETECT, part of the CIAO software package. Aspects of our algorithm include: (1) the computation of local, exposure-corrected normalized (i.e. flat-fielded) background maps; (2) the correction for exposure variations within the field-of-view; (3) its applicability within the low-counts regime, as it does not require a minimum number of background counts per pixel for the accurate computation of source detection thresholds; (4) the generation of a source list in a manner that does not depend upon a detailed knowledge of the point spread function (PSF) shape; and (5) error analysis. These features make our algorithm considerably more general than previous methods developed for the analysis of X-ray image data, especially in the low count regime. We demonstrate the algorithm's robustness by applying it to various images.

In spring 1999, a long coordinated observing campaign was performed on the flare star AD Leo, including {\it EUVE}, {\it BeppoSAX}, the {\it VLA}, and optical telescopes. The campaign covered a total of 44 days. We obtained high-quality light curves displaying ongoing variability on various timescales, raising interesting questions on the role of flare-like events for coronal heating. We performed Kolmogorov-Smirnov tests to compare the observations with a large set of simulated light curves composed of statistical flares that are distributed in energy as a power law of the form ${\rm d}N/{\rm d}E \propto E^{-\alpha}$ with selectable index $\alpha$. We find best-fit $\alpha$ values slightly above a value of 2, indicating that the extension of the flare population toward small energies could be important for the generation of the overall X-ray emission.

The launch of the Chandra X-ray Observatory in July 1999 opened a new era in X-ray astronomy. Its unprecedented, <0.5" spatial resolution and low background are providing views of the X-ray sky 10-100 times fainter than previously possible. We have initiated a serendipitous survey (ChaMP) using Chandra archival data to flux limits covering the range between those reached by current satellites and those of the small area Chandra deep surveys. We estimate the survey will cover ~5 sq.deg./year to X-ray fluxes (2-10 keV) in the range 1E(-13)-6E(-16) erg/cm^2/s discovering ~2000 new X-ray sources, ~80% of which are expected to be AGN. The ChaMP has two parts, the extragalactic survey (ChaMP) and the galactic plane survey (ChaMPlane). ChaMP promises profoundly new science return on a number of key questions at the current frontier of many areas of astronomy including (1) locating and studying high redshift clusters and so constraining cosmological parameters (2) defining the true population of AGN, including those that are absorbed, and so constraining the accretion history of the universe, (3) filling in the gap in the luminosity/redshift plane between Chandra deep and previous surveys in studying the CXRB, (4) studying coronal emission from late-type stars and (5) search for CVs and quiescent Low-Mass X-ray Binaries (qLXMBs) to measure their luminosity functions. In this paper we summarize the status, predictions and initial results from the X-ray analysis and optical imaging.

Over the past 10 years Bayesian methods have rapidly grown more popular as several computationally intensive statistical algorithms have become feasible with increased computer power. In this paper, we begin with a general description of the Bayesian paradigm for statistical inference and the various state-of-the-art model fitting techniques that we employ (e.g., Gibbs sampler and Metropolis- Hastings). These algorithms are very flexible and can be used to fit models that account for the highly hierarchical structure inherent in the collection of high-quality spectra and thus can keep pace with the accelerating progress of new space telescope designs. The methods we develop, which will soon be available in the CIAO software package, explicitly model photon arrivals as a Poisson process and, thus, have no difficulty with high resolution low count X-ray and gamma-ray data. We expect these methods to be useful not only for the recently launched Chandra X-ray observatory and XMM but also new generation telescopes such as Constellation X, GLAST, etc. In the context of two examples (Quasar S5 0014+813 and Hybrid-Chromosphere Supergiant Star alpha TrA) we illustrate a new highly structured model and how Bayesian posterior sampling can be used to compute estimates, error bars, and credible intervals for the various model parameters.

\textit{Extreme Ultraviolet Explorer} Deep Survey observations of cool stars (spectral type F to M) have been used to investigate the distribution of coronal flare rates in energy and its relation to activity indicators and rotation parameters. Cumulative and differential flare rate distributions were constructed and fitted with different methods. Power laws are found to approximately describe the distributions. A trend toward flatter distributions for later-type stars is suggested in our sample. Assuming that the power laws continue below the detection limit, we have estimated that the superposition of flares with radiated energies of about $10^{29}-10^{31}$ergs could explain the observed radiative power loss of these coronae, while the detected flares are contributing only $\approx 10$%. While the power-law index is not correlated with rotation parameters (rotation period, projected rotational velocity, Rossby number) and only marginally with the X-ray luminosity, the flare occurrence rate is correlated with all of them. The occurrence rate of flares with energies larger than $10^{32}$ergs is found to be proportional to the average total stellar X-ray luminosity. Thus, energetic flares occur more often in X-ray bright stars than in X-ray faint stars. The normalized occurrence rate of flares with energies larger than $10^{32}$ergs increases with increasing \lxlbol and stays constant for saturated stars. A similar saturation is found below a critical Rossby number. The findings are discussed in terms of simple statistical flare models in an attempt to explain the previously observed trend for higher average coronal temperatures in more active stars. It is concluded that flares can contribute a significant amount of energy to coronal heating in active stars.

We present a deep ASCA observation of a Broad Absorption Line Quasar (BALQSO) PG0946+301. The source was clearly detected in one of the gas imaging spectrometers, but not in any other detector. If BALQSOs have intrinsic X-ray spectra similar to normal radio-quiet quasars, our observations imply that there is Thomson thick X-ray absorption (N_H >~10^{24} cm^{-2}) toward PG0946+301. This is the largest column density estimated so far toward a BALQSO. The absorber must be at least partially ionized and may be responsible for attenuation in the optical and UV. If the Thomson optical depth toward BALQSOs is close to one, as inferred here, then spectroscopy in hard X-rays with large telescopes like XMM would be feasible.

We present the first X-ray spectrum obtained by the Low Energy Transmission Grating Spectrometer (LETGS) aboard the Chandra X-ray Observatory. The spectrum is of Capella and covers a wavelength range of 5-175 A (2.5-0.07 keV). The measured wavelength resolution, which is in good agreement with ground calibration, is $\Delta \lambda \simeq$ 0.06 A (FWHM). Although in-flight calibration of the LETGS is in progress, the high spectral resolution and unique wavelength coverage of the LETGS are well demonstrated by the results from Capella, a coronal source rich in spectral emission lines. While the primary purpose of this letter is to demonstrate the spectroscopic potential of the LETGS, we also briefly present some preliminary astrophysical results. We discuss plasma parameters derived from line ratios in narrow spectral bands, such as the electron density diagnostics of the He-like triplets of carbon, nitrogen, and oxygen, as well as resonance scattering of the strong Fe XVII line at 15.014 A.

We compare the count-rates of active binaries observed by Einstein and ROSAT. We find significant evidence for short-term variability (approx 0.32) and marginal evidence for a larger variability at longer timescales (approx 0.38). Modeling this excess as cyclic variability, we place an upper limit on the relative amplitude of the cyclic component (I_cyc/I_min < 4) for active binaries, significantly less than the variations seen in the Solar X-ray output over the Solar cycle.