Source author record

J. J. Drake

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

75works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

75 published item(s)

preprint2022arXiv

Space Weather-driven Variations in Ly$α$ Absorption Signatures of Exoplanet Atmospheric Escape: MHD Simulations and the Case of AU Mic b

We simulate the space environment around AU Microscopii b and the interaction between the magnetized stellar wind with a planetary atmospheric outflow for ambient stellar wind conditions and Coronal Mass Ejection (CME) conditions. We also calculate synthetic Ly$α$ absorption due to neutral hydrogen in the ambient and the escaping planetary atmosphere affected by this interaction. We find that the Ly$α$ absorption is highly variable due to the highly-varying stellar wind conditions. A strong Doppler blue-shift component is observed in the Ly$α$ profile, in contradiction to the actual escape velocity observed in the simulations themselves. This result suggest that the strong Doppler blue-shift is likely attributed to the stellar wind, not the escaping neutral atmosphere, either through its advection of neutral planetary gas, or through the creation of a fast neutral flow via charge exchange between the stellar wind ions and the planetary neutrals. Indeed, our CME simulations indicate a strong stripping of magnetospheric material from the planet, including some of the neutral escaping atmosphere. Our simulations show that the pressure around close-in exoplanets is not much lower, and may be even higher, than the pressure at the top of the planetary atmosphere. Thus, the neutral atmosphere is hydrodynamically escaping with a very small velocity ($<15~km~s^{-1}$). Moreover, our simulations show that an MHD treatment is essential in order to properly capture the coupled magnetized stellar wind and the escaping atmosphere, despite of the atmosphere being neutral. This coupling should be considered when interpreting Ly$α$observations in the context of exoplanets atmospheric escape.

preprint2022arXiv

The 2021 outburst of the recurrent nova RS Ophiuchi observed in X-rays by the Neil Gehrels Swift Observatory: a comparative study

On 2021 August 8, the recurrent nova RS Ophiuchi erupted again, after an interval of 15.5 yr. Regular monitoring by the Neil Gehrels Swift Observatory began promptly, on August 9.9 (0.37 day after the optical peak), and continued until the source passed behind the Sun at the start of November, 86 days later. Observations then restarted on day 197, once RS Oph emerged from the Sun constraint. This makes RS Oph the first Galactic recurrent nova to have been monitored by Swift throughout two eruptions. Here we investigate the extensive X-ray datasets from 2006 and 2021, as well as the more limited data collected by EXOSAT in 1985. The hard X-rays arising from shock interactions between the nova ejecta and red giant wind are similar following the last two eruptions. In contrast, the early super-soft source (SSS) in 2021 was both less variable and significantly fainter than in 2006. However, 0.3-1 keV light-curves from 2021 reveal a 35 s quasi-periodic oscillation consistent in frequency with the 2006 data. The Swift X-ray spectra from 2021 are featureless, with the soft emission typically being well parametrized by a simple blackbody, while the 2006 spectra showed much stronger evidence for superimposed ionized absorption edges. Considering the data after day 60 following each eruption, during the supersoft phase the 2021 spectra are hotter, with smaller effective radii and lower wind absorption, leading to an apparently reduced bolometric luminosity. We explore possible explanations for the gross differences in observed SSS behaviour between the 2006 and 2021 outbursts.

preprint2021arXiv

The Super-Soft Source Phase of the recurrent nova V3890 Sgr

The 30-year recurrent symbiotic nova V3890 Sgr exploded 2019 August 28 and was observed with multiple X-ray telescopes. An XMM-Newton observation during the SSS phase captured a high degree of X-ray variability including a deep dip in the middle of the observation, an initial rise of similar depth and shape and, after the deep dip, smaller dips of 10% amplitude, which might be periodic over 18.1-minutes. An eclipse model of the dips yields clump sizes and orbital radii of 0.5-8 and 5-150 white dwarf radii, respectively. The simultaneous UV light curve shows no significant variations beyond slow fading. The RGS spectrum contains both residual shock emission at short wavelengths and the SSS emission at longer wavelengths. The shock temperature has clearly decreased compared to an earlier Chandra observation (day 6). The dip spectrum is dominated by emission lines like in U Sco. The intensity of underlying blackbody-like emission is much lower with the blackbody normalisation yielding a similar radius as during the brighter phases, while the lower bolometric luminosity is ascribed to lower T_eff. This would be inconsistent with clump occultations unless Compton scattering of the continuum emission reduces the photon energies to mimic a lower effective temperature. However, systematic uncertainties are high. The absorption lines in the bright SSS spectrum are blue-shifted by 870+/-10 km/s before the dip and 900+/-10 km/s, after the dip. The reproduction of the observed spectrum is astonishing, especially that only a single absorbing layer is necessary while three such layers are needed to reproduce the RGS spectrum of V2491 Cyg. The ejecta of V3890 Sgr are thus more homogeneous than many other SSS spectra indicate. Abundance determination is in principle possible but highly uncertain. Generally, solar abundances are found except for N and possibly O higher by an order of magnitude.

preprint2020arXiv

COBRaS: The e-MERLIN 21 cm Legacy survey of Cygnus OB2

The role of massive stars is central to an understanding of galactic ecology. It is important to establish the details of how massive stars provide radiative, chemical, and mechanical feedback in galaxies. Central to these issues is an understanding of the evolution of massive stars, and the critical role of mass loss via strongly structured winds and stellar binarity. Ultimately, massive stellar clusters shape the structure and energetics of galaxies. We aim to conduct high-resolution, deep field mapping at 21cm of the core of the massive Cygnus OB2 association and to characterise the properties of the massive stars and colliding winds at this waveband. We used seven stations of the e-MERLIN radio facility, with its upgraded bandwidth and enhanced sensitivity to conduct a 21cm census of Cygnus OB2. Based on 42 hours of observations, seven overlapping pointings were employed over multiple epochs during 2014 resulting in 1 sigma sensitivities down to ~21microJy and a resolution of ~180mas. A total of 61 sources are detected at 21cm over a ~0.48deg x 0.48deg region centred on the heart of the Cyg OB2 association. Of these 61 sources, 33 are detected for the first time. We detect a number of previously identified sources including four massive stellar binary systems, two YSOs, and several known X-ray and radio sources. We also detect the LBV candidate (possible binary system) and blue hypergiant (BHG) star of Cyg OB2 #12. The 21cm observations secured in the COBRaS Legacy project provide data to constrain conditions in the outer wind regions of massive stars; determine the non-thermal properties of massive interacting binaries; examine evidence for transient sources, including those associated with young stellar objects; and provide unidentified sources that merit follow-up observations. The 21cm data are of lasting value and will serve in combination with other key surveys of Cyg OB2.

preprint2016arXiv

FK Comae Berenices, King of Spin: The COCOA-PUFS Project

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.

preprint2016arXiv

Giant coronal loops dominate the quiescent X-ray emission in rapidly rotating M stars

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.

preprint2016arXiv

Optical, UV, and X-Ray Evidence for a 7-Year Stellar Cycle in Proxima Centauri

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

preprint2016arXiv

Photoevaporation and close encounters: how the environment around Cygnus OB2 affects the evolution of protoplanetary disks

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.

preprint2016arXiv

Physical parameters of IPHAS-selected classical Be stars. (I. Determination procedure and evaluation of the results.)

We present a semi-automatic procedure to obtain fundamental physical parameters and distances of classical Be (CBe) stars, based on the Barbier-Chalonge-Divan (BCD) spectrophotometric system. Our aim is to apply this procedure to a large sample of CBe stars detected by the IPHAS photometric survey, to determine their fundamental physical parameters and to explore their suitability as galactic structure tracers. In this paper we describe the methodology used and the validation of the procedure by comparing our results with those obtained from different independent astrophysical techniques for subsamples of stars in common with other studies. We also present a test case study of the galactic structure in the direction of the Perseus Galactic Arm, in order to compare our results with others recently obtained with different techniques and the same sample of stars. We did not find any significant clustering of stars at the expected positions of the Perseus and Outer Galactic Arms, in agreement with previous studies in the same area that we used for verification.

preprint2016arXiv

Simulating the Environment Around Planet-Hosting Stars - I. Coronal Structure

We present the results of a detailed numerical simulation of the circumstellar environment around three exoplanet-hosting stars. A state-of-the-art global magnetohydrodynamic (MHD) model is considered, including Alfvén wave dissipation as a self-consistent coronal heating mechanism. This paper contains the description of the numerical set-up, evaluation procedure, and the simulated coronal structure of each system (HD 1237, HD 22049 and HD 147513). The simulations are driven by surface magnetic field maps, recovered with the observational technique of Zeeman Doppler Imaging (ZDI). A detailed comparison of the simulations is performed, where two different implementations of this mapping routine are used to generate the surface field distributions. Quantitative and qualitative descriptions of the coronae of these systems are presented, including synthetic high-energy emission maps in the Extreme Ultra-Violet (EUV) and Soft X-rays (SXR) ranges. Using the simulation results, we are able to recover similar trends as in previous observational studies, including the relation between the magnetic flux and the coronal X-ray emission. Furthermore, for HD 1237 we estimate the rotational modulation of the high-energy emission due to the various coronal features developed in the simulation. We obtain variations, during a single stellar rotation cycle, up to 15\% for the EUV and SXR ranges. The results presented here will be used, in a follow-up paper, to self-consistently simulate the stellar winds and inner astrospheres of these systems.

preprint2016arXiv

Simulating the environment around planet-hosting stars - II. Stellar winds and inner astrospheres

We present the results of a comprehensive numerical simulation of the environment around three exoplanet-host stars (HD 1237, HD 22049, and HD 147513). Our simulations consider one of the latest models currently used for space weather studies in the Heliosphere. Large-scale magnetic field maps, recovered with two implementations of the tomographic technique of Zeeman-Doppler imaging, serve to drive steady-state solutions in each system. This paper contains the description of the stellar wind and inner astrosphere, while the coronal structure was previously discussed in Alvarado-Gómez et al. (2016). The analysis includes the magneto-hydrodynamical properties of the stellar wind, the associated mass and angular momentum loss rates, as well as the topology of the astrospheric current sheet in each system. A systematic comparison among the considered cases is performed, including two reference solar simulations covering activity minimum and maximum. For HD 1237, we investigate the interactions between the structure of the developed stellar wind, and a possible magnetosphere around the Jupiter-mass planet in this system. We find that the process of particle injection into the planetary atmosphere is dominated by the density distribution rather than velocity profile of the stellar wind. In this context, we predict a maximum exoplanetary radio emission of 12 mJy at 40 MHz in this system, assuming the crossing of a high-density streamer during periastron passage. Furthermore, in combination with the analysis performed in Alvarado-Gómez et al. (2016), we obtain for the first time a fully simulated mass loss-activity relation, which is compared and discussed in the context of the relation based on astrospheric detections proposed by Wood et al. (2005a). Finally, we provide a characterisation of the 3D properties of the stellar wind of these systems, at the inner edges of their habitable zones.

preprint2016arXiv

The statistical uncertainties on X-ray flux and spectral parameters from Chandra ACIS-I observations of faint sources: Application to the Cygnus OB2 Association

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.

preprint2015arXiv

Activity and Magnetic Field Structure of the Sun-Like Planet Hosting Star HD 1237

We analyse the magnetic activity characteristics of the planet hosting Sun-like star, HD 1237, using HARPS spectro-polarimetric time-series data. We find evidence of rotational modulation of the magnetic longitudinal field measurements consistent with our ZDI analysis, with a period of 7 days. We investigate the effect of customising the LSD mask to the line depths of the observed spectrum and find that it has a minimal effect on shape of the extracted Stokes V profile but does result in a small increase in the S/N ($\sim$ 7%). We find that using a Milne-Eddington solution to describe the local line profile provides a better fit to the LSD profiles in this slowly rotating star, which also impacts the recovered ZDI field distribution. We also introduce a fit-stopping criterion based on the information content (entropy) of the ZDI maps solution set. The recovered magnetic field maps show a strong (+90 G) ring-like azimuthal field distribution and a complex radial field dominating at mid latitudes ($\sim$45 degrees). Similar magnetic field maps are recovered from data acquired five months apart. Future work will investigate how this surface magnetic field distribution impacts the coronal magnetic field and extended environment around this planet-hosting star.

preprint2015arXiv

Classical T Tauri stars with VPHAS$+$: I: H$α$ and $u$-band accretion rates in the Lagoon Nebula M8

We estimate the accretion rates of 235 Classical T Tauri star (CTTS) candidates in the Lagoon Nebula using $ugri$H$α$ photometry from the VPHAS+ survey. Our sample consists of stars displaying H$α$-excess, the intensity of which is used to derive accretion rates. For a subset of 87 stars, the intensity of the $u$-band excess is also used to estimate accretion rates. We find the mean variation in accretion rates measured using H$α$ and $u$-band intensities to be $\sim$ 0.17 dex, agreeing with previous estimates (0.04-0.4 dex) but for a much larger sample. The spatial distribution of CTTS align with the location of protostars and molecular gas suggesting that they retain an imprint of the natal gas fragmentation process. Strong accretors are concentrated spatially, while weak accretors are more distributed. Our results do not support the sequential star forming processes suggested in the literature.

preprint2015arXiv

Optical and infrared counterparts of the X-ray sources detected in the Chandra Cygnus OB2 Legacy Survey

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.

preprint2015arXiv

Short-period X-ray oscillations in super-soft novae and persistent SSS

Transient short-period <100s oscillations have been found in the X-ray light curves of three novae during their SSS phase and in one persistent SSS. We pursue an observational approach to determine possible driving mechanisms and relations to fundamental system parameters such as the white dwarf mass. We performed a systematic search for short-period oscillations in all available XMM-Newton and Chandra X-ray light curves of persistent SSS and novae during their SSS phase. To study time evolution, we divided each light curve into short time segments and computed power spectra. We then constructed dynamic power spectra from which we identified transient periodic signals even when only present for a short time. From all time segments of each system, we computed fractions of time when periodic signals were detected. In addition to the previously known systems with short-period oscillations, RS Oph (35s), KT Eri (35s), V339 Del (54s), and Cal 83 (67s), we found one additional system, LMC 2009a (33s), and also confirm the 35s period from Chandra data of KT Eri. The amplitudes of oscillations are of order <15% of the respective count rates and vary without any clear dependence on the X-ray count rate. The fractions of the time when the respective periods were detected at 2-sigma significance (duty cycle) are 11.3%, 38.8%, 16.9%, 49.2%, and 18.7% for LMC 2009a, RS Oph, KT Eri, V339 Del, and Cal 83, respectively. The respective highest duty cycles found in a single observation are 38.1%, 74.5%, 61.4%, 67.8%, and 61.8%.

preprint2015arXiv

Swift detection of the super-swift switch-on of the super-soft phase in nova V745 Sco (2014)

V745 Sco is a recurrent nova, with the most recent eruption occurring in February 2014. V745 Sco was first observed by Swift a mere 3.7 hr after the announcement of the optical discovery, with the super-soft X-ray emission being detected around four days later and lasting for only ~two days, making it both the fastest follow-up of a nova by Swift and the earliest switch-on of super-soft emission yet detected. Such an early switch-on time suggests a combination of a very high velocity outflow and low ejected mass and, together with the high effective temperature reached by the super-soft emission, a high mass white dwarf (>1.3 M_sun). The X-ray spectral evolution was followed from an early epoch where shocked emission was evident, through the entirety of the super-soft phase, showing evolving column density, emission lines, absorption edges and thermal continuum temperature. UV grism data were also obtained throughout the super-soft interval, with the spectra showing mainly emission lines from lower ionization transitions and the Balmer continuum in emission. V745 Sco is compared with both V2491 Cyg (another nova with a very short super-soft phase) and M31N 2008-12a (the most rapidly recurring nova yet discovered). The longer recurrence time compared to M31N 2008-12a could be due to a lower mass accretion rate, although inclination of the system may also play a part. Nova V745 Sco (2014) revealed the fastest evolving super-soft source phase yet discovered, providing a detailed and informative dataset for study.

preprint2015arXiv

The Interaction of Venus-like, M-dwarf Planets with the Stellar Wind of Their Host Star

We study the interaction between the atmospheres of Venus-like, non-magnetized exoplanets orbiting an M-dwarf star, and the stellar wind using a multi-species Magnetohydrodynaic (MHD) model. We focus our investigation on the effect of enhanced stellar wind and enhanced EUV flux as the planetary distance from the star decreases. Our simulations reveal different topologies of the planetary space environment for sub- and super-Alfvenic stellar wind conditions, which could lead to dynamic energy deposition in to the atmosphere during the transition along the planetary orbit. We find that the stellar wind penetration for non-magnetized planets is very deep, up to a few hundreds of kilometers. We estimate a lower limit for the atmospheric mass-loss rate and find that it is insignificant over the lifetime of the planet. However, we predict that when accounting for atmospheric ion acceleration, a significant amount of the planetary atmosphere could be eroded over the course of a billion years.

preprint2015arXiv

X-ray Fading and Expansion in the "Miniature Supernova Remnant" of GK Persei

We report on a second epoch of Chandra X-ray imaging spectroscopy of the spatially-resolved old nova remnant GK Persei. An ACIS-S3 observation of 97.4 ks was conducted in November 2013 after a lapse of 13.8 years from the last visit in 2000. The X-ray emitting nebula appeared more faint and patchy compared with the first epoch. The flux decline was particularly evident in fainter regions and the mean decline was 30-40 % in the 0.5-1.2 keV energy band. A typical expansion of the brightest part of the remnant was 1.9 arcsec, which corresponds to an expansion rate of 0.14 arcsec yr^{-1}. The soft X-ray spectra extracted from both the 2000 and 2013 data can be explained by a non-equilibrium ionization collisional plasma model convolved with interstellar absorption, though do not allow us to constrain the origin of the flux evolution. The plasma temperature has not significantly evolved since the 2000 epoch and we conclude that the fading of the X-ray emission is due largely to expansion. This implies that recent expansion has been into a lower density medium, a scenario that is qualitatively consistent with the structure of the circumstellar environment photographed soon after the initial explosion more than a century ago. Fainter areas are fading more quickly than brighter areas, indicating that they are fainter because of a lower ambient medium density and consequently more rapid expansion.

preprint2014arXiv

A deep catalogue of classical Be stars in the direction of the Perseus Arm: spectral types and interstellar reddenings

We present a catalogue of 247 photometrically and spectroscopically confirmed fainter classical Be stars (13 < r < 16) in the direction of the Perseus Arm of the Milky Way (-1 < b < +4, 120 < l < 140). The catalogue consists of 181 IPHAS-selected new classical Be stars, in addition to 66 objects that were studied by Raddi et al. (2013) more closely, and 3 stars identified as classical Be stars in earlier work. This study more than doubles the number known in the region. Photometry spanning 0.6 to 5 micron, spectral types, and interstellar reddenings are given for each object. The spectral types were determined from low-resolution spectra (lambda / Delta-lambda ~ 800-2000), to a precision of 1-3 subtypes. The interstellar reddenings are derived from the (r - i) colour, using a method that corrects for circumstellar disc emission. The colour excesses obtained range from E(B-V) = 0.3 up to 1.6 - a distribution that modestly extends the range reported in the literature for Perseus-Arm open clusters. For around half the sample, the reddenings obtained are compatible with measures of the total sightline Galactic extinction. Many of these are likely to lie well beyond the Perseus Arm.

preprint2014arXiv

Accretion and outflow in the proplyd-like objects near Cygnus OB2

Cygnus OB2 is the most massive association within 2kpc from the Sun, hosting hundreds of massive stars, thousands of young low mass members, and some sights of active star formation in the surrounding cloud. Recently, 10 photoevaporating proplyd-like objects with tadpole-shaped morphology were discovered in the outskirts of the OB association, approximately 6-14pc away from its center. The classification of these objects is ambiguous, being either evaporating residuals of the parental cloud which are hosting a protostar inside, or disk-bearing stars with an evaporating disk, such as the evaporating proplyds observed in the Trapezium Cluster in Orion. In this paper we present a study based on low resolution optical spectroscopic observations made with the Optical System for Imaging and low Resolution Integrated Spectroscopy (OSIRIS), mounted on the 10.4m Gran Telescopio CANARIAS (GTC), of two of these protostars. The spectrum of one of the objects shows evidence of accretion but not of outflows. In the latter object, the spectra show several emission lines indicating the presence of an actively accreting disk with outflow. We present estimates of the mass loss rate and the accretion rate from the disk, showing that the former exceeds the latter as observed in other known objects with evaporating disks. We also show evidence of a strong variability in the integrated flux observed in these objects, as well as in the accretion and outflow diagnostics.

preprint2014arXiv

Chasing the Identification of ASCA Galactic Objects (ChIcAGO) - An X-ray Survey of Unidentified Sources in the Galactic Plane I: Source Sample and Initial Results

We present the Chasing the Identification of ASCA Galactic Objects (ChIcAGO) survey, which is designed to identify the unknown X-ray sources discovered during the ASCA Galactic Plane Survey (AGPS). Little is known about most of the AGPS sources, especially those that emit primarily in hard X-rays (2-10 keV) within the F_x ~ 10^-13 to 10^-11 erg cm^-2 s^-1 X-ray flux range. In ChIcAGO, the subarcsecond localization capabilities of Chandra have been combined with a detailed multi-wavelength follow-up program, with the ultimate goal of classifying the >100 unidentified sources in the AGPS. Overall to date, 93 unidentified AGPS sources have been observed with Chandra as part of the ChIcAGO survey. A total of 253 X-ray point sources have been detected in these Chandra observations within 3' of the original ASCA positions. We have identified infrared and optical counterparts to the majority of these sources, using both new observations and catalogs from existing Galactic plane surveys. X-ray and infrared population statistics for the X-ray point sources detected in the Chandra observations reveal that the primary populations of Galactic plane X-ray sources that emit in the F_x ~ 10^-13 to 10^-11 erg cm^-2 s^-1 flux range are active stellar coronae, massive stars with strong stellar winds that are possibly in colliding-wind binaries, X-ray binaries, and magnetars. There is also a fifth population that is still unidentified but, based on its X-ray and infrared properties, likely comprise partly of Galactic sources and partly of active galactic nuclei.

preprint2014arXiv

Magnetospheric Structure and Atmospheric Joule Heating of Habitable Planets Orbiting M-dwarf Stars

We study the magnetospheric structure and the ionospheric Joule Heating of planets orbiting M-dwarf stars in the habitable zone using a set of magnetohydrodynamic (MHD) models. The stellar wind solution is used to drive a model for the planetary magnetosphere, which is coupled with a model for the planetary ionosphere. Our simulations reveal that the space environment around close-in habitable planets is extreme, and the stellar wind plasma conditions change from sub- to super-Alfvenic along the planetary orbit. As a result, the magnetospheric structure changes dramatically with a bow shock forming in the super-Alfvenic sectors, while no bow shock forms in the sub-Alfvenic sectors. The planets reside most of the time in the sub-Alfvenic sectors with poor atmospheric protection. A significant amount of Joule Heating is provided at the top of the atmosphere as a result of the planetary interaction with the stellar wind. For the steady-state solution, the heating is about 0.1-3\% of the total incoming stellar irradiation, and it is enhanced by 50\% for the time-dependent case. The significant Joule Heating obtained here should be considered in models for the atmospheres of habitable planets in terms of the thickness of the atmosphere, the top-side temperature and density, the boundary conditions for the atmospheric pressure, and particle radiation and transport.

preprint2014arXiv

The corona of HD 189733 and its X-ray activity

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 $ϕ= 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.

preprint2014arXiv

The Second Data Release of the INT Photometric H-Alpha Survey of the Northern Galactic Plane (IPHAS DR2)

The INT/WFC Photometric H-Alpha Survey of the Northern Galactic Plane (IPHAS) is a 1800 square degrees imaging survey covering Galactic latitudes |b| < 5 deg and longitudes l = 30 to 215 deg in the r, i and H-alpha filters using the Wide Field Camera (WFC) on the 2.5-metre Isaac Newton Telescope (INT) in La Palma. We present the first quality-controlled and globally-calibrated source catalogue derived from the survey, providing single-epoch photometry for 219 million unique sources across 92% of the footprint. The observations were carried out between 2003 and 2012 at a median seeing of 1.1 arcsec (sampled at 0.33 arcsec/pixel) and to a mean 5σ-depth of 21.2 (r), 20.0 (i) and 20.3 (H-alpha) in the Vega magnitude system. We explain the data reduction and quality control procedures, describe and test the global re-calibration, and detail the construction of the new catalogue. We show that the new calibration is accurate to 0.03 mag (rms) and recommend a series of quality criteria to select the most reliable data from the catalogue. Finally, we demonstrate the ability of the catalogue's unique (r-Halpha, r-i) diagram to (1) characterise stellar populations and extinction regimes towards different Galactic sightlines and (2) select H-alpha emission-line objects. IPHAS is the first survey to offer comprehensive CCD photometry of point sources across the Galactic Plane at visible wavelengths, providing the much-needed counterpart to recent infrared surveys.

preprint2014arXiv

The VST Photometric Halpha Survey of the Southern Galactic Plane and Bulge (VPHAS+)

The VST Photometric Halpha Survey of the Southern Galactic Plane and Bulge (VPHAS+) is surveying the southern Milky Way in u, g, r, i and Halpha at 1 arcsec angular resolution. Its footprint spans the Galactic latitude range -5 < b < +5 at all longitudes south of the celestial equator. Extensions around the Galactic Centre to Galactic latitudes +/-10 bring in much of the Galactic Bulge. This ESO public survey, begun on 28th December 2011, reaches down to 20th magnitude (10-sigma) and will provide single-epoch digital optical photometry for around 300 million stars. The observing strategy and data pipelining is described, and an appraisal of the segmented narrowband Halpha filter in use is presented. Using model atmospheres and library spectra, we compute main-sequence (u - g), (g - r), (r - i) and (r - Halpha) stellar colours in the Vega system. We report on a preliminary validation of the photometry using test data obtained from two pointings overlapping the Sloan Digital Sky Survey. An example of the (u - g, g - r) and (r - Halpha, r - i) diagrams for a full VPHAS+ survey field is given. Attention is drawn to the opportunities for studies of compact nebulae and nebular morphologies that arise from the image quality being achieved. The value of the u band as the means to identify planetary-nebula central stars is demonstrated by the discovery of the central star of NGC 2899 in survey data. Thanks to its excellent imaging performance, the VST/OmegaCam combination used by this survey is a perfect vehicle for automated searches for reddened early-type stars, and will allow the discovery and analysis of compact binaries, white dwarfs and transient sources.

preprint2013arXiv

First results of an H-alpha based search of classical Be stars in the Perseus Arm and beyond

We investigate a region of the Galactic plane, between 120 <= l <= 140 and -1 <= b <= +4, and uncover a population of moderately reddened (E(B-V) \sim 1) classical Be stars within and beyond the Perseus and Outer Arms. 370 candidate emission line stars (13 <= r <= 16) selected from the INT Photometric H-alpha Survey of the Northern Galactic plane (IPHAS) have been followed up spectroscopically. A subset of these, 67 stars with properties consistent with those of classical Be stars, have been observed at sufficient spectral resolution (Delta_lambda \sim 2 - 4 Angstrom) at blue wavelengths to narrow down their spectral types. We determine these to a precision estimated to be +/- 1 sub-type and then we measure reddenings via SED fitting with reference to appropriate model atmospheres. Corrections for contribution to colour excess from circumstellar discs are made using an established scaling to H-alpha emission equivalent width. Spectroscopic parallaxes are obtained after luminosity class has been constrained via estimates of distances to neighbouring A/F stars with similar reddenings. Overwhelmingly, the stars in the sample are confirmed as luminous classical Be stars at heliocentric distances ranging from 2 kpc up to \sim 12 kpc. However, the errors are presently too large to enable the cumulative distribution function with respect to distance to distinguish between models placing the stars exclusively in spiral arms, or in a smooth exponentially-declining distribution.

preprint2013arXiv

Obscuration effects in Super-Soft-Source X-ray spectra

Super-Soft-Source (SSS) X-ray spectra are blackbody-like spectra with effective temperatures ~3-7x10^5 K and luminosities of 10^{35-38} erg/s. SSS grating spectra display atmospheric absorption lines. Radiation transport atmosphere models can be used to derive physical parameters, but more sophisticated models are required. We bypass the complications of spectral models and concentrate on the data in a comparative, qualitative study. We inspect all available X-ray grating SSS spectra to determine systematic, model-independent trends. We use comparative plots of spectra of different systems to find common and different features. The results are interpreted in the context of system parameters obtained from the literature. We find two distinct types of SSS spectra which we name SSa and SSe. Their main observational characteristics are either clearly visible absorption lines or emission lines, respectively, while both types contain atmospheric continuum emission. SSe may be obscured SSa systems, which is supported by similarities between SSe and SSa with obscured and unobscured AGN, respectively. Further, we find all known or suspected high-inclination systems to emit permanently in an SSe state. Some sources are found to transition between SSa and SSe states, becoming SSe when fainter. SSS spectra are subject to various occultation processes. In Cal 87, the accretion disc blocks the central hot source when viewed edge on. In novae, the accretion disc may have been destroyed during the initial explosion but could have reformed by the time of the SSS phase. In addition, clumpy ejecta may lead to temporary obscuration events. The emission lines originate from reprocessed emission in the accretion disc, its wind or further out in clumpy ejecta while Thomson scattering allows continuum emission to be visible also during total obscuration of the central hot source.

preprint2013arXiv

The protoplanetary disks in the nearby massive star forming region Cygnus OB2

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.

preprint2013arXiv

X-ray Development of the Classical Nova V2672 Ophiuchi with Suzaku

We report the Suzaku detection of a rapid flare-like X-ray flux amplification early in the development of the classical nova V2672 Ophiuchi. Two target-of-opportunity ~25 ks X-ray observations were made 12 and 22 days after the outburst. The flux amplification was found in the latter half of day 12. Time-sliced spectra are characterized by a growing supersoft excess with edge-like structures and a relatively stable optically-thin thermal component with Ka emission lines from highly ionized Si. The observed spectral evolution is consistent with a model that has a time development of circumstellar absorption, for which we obtain the decline rate of ~10-40 % in a time scale of 0.2 d on day 12. Such a rapid drop of absorption and short-term flux variability on day 12 suggest inhomogeneous ejecta with dense blobs/holes in the line of sight. Then on day 22 the fluxes of both supersoft and thin-thermal plasma components become significantly fainter. Based on the serendipitous results we discuss the nature of this source in the context of both short- and long-term X-ray behavior.

preprint2013arXiv

X-ray Eclipse Diagnosis of the Evolving Mass Loss in the Recurrent Nova U Scorpii 2010

We report the Suzaku detection of the earliest X-ray eclipse seen in the recurrent nova U Scorpii 2010. A target-of-opportunity observation 15 days after the outburst found a 27+/-5% dimming in the 0.2-1.0 keV energy band at the predicted center of an eclipse. In comparison with the X-ray eclipse depths seen at two later epochs by XMM-Newton, the source region shrank by about 10-20% between days 15 and 35 after the outburst. The X-ray eclipses appear to be deeper than or similar to contemporaneous optical eclipses, suggesting the X-ray and optical source region extents are comparable on day 15. We raise the possibility of the energy dependency in the photon escape regions, and that this would be a result of the supersoft X-ray opacity being higher than the Thomson scattering optical opacity at the photosphere due to bound-free transitions in abundant metals that are not fully ionized. Assuming a spherically symmetric explosion model, we constrain the mass-loss rate as a function of time. For a ratio of actual to Thomson opacity of 10-100 in supersoft X-rays, we find a total ejecta mass of about 10^{-7}-10^{-6} M_{\odot}.

preprint2012arXiv

Chandra grating spectroscopy of three hot white dwarfs

High-resolution soft X-ray spectroscopic observations of single hot white dwarfs are scarce. With the Chandra Low-Energy Transmission Grating, we have observed two white dwarfs, one is of spectral type DA (LB 1919) and the other is a non-DA of spectral type PG1159 (PG 1520+525). The spectra of both stars are analyzed, together with an archival Chandra spectrum of another DA white dwarf (GD 246). The soft X-ray spectra of the two DA white dwarfs are investigated in order to study the effect of gravitational settling and radiative levitation of metals in their photospheres. LB 1919 is of interest because it has a significantly lower metallicity than DAs with otherwise similar atmospheric parameters. GD 246 is the only white dwarf known that shows identifiable individual iron lines in the soft X-ray range. For the PG1159 star, a precise effective temperature determination is performed in order to confine the position of the blue edge of the GW Vir instability region in the HRD. (abridged)

preprint2012arXiv

Chandra/ACIS-I study of the X-ray properties of the NGC 6611 and M16 stellar population

Mechanisms regulating the origin of X-rays in YSOs and the correlation with their evolutionary stage are under debate. Studies of the X-ray properties in young clusters allow to understand these mechanisms. One ideal target for this analysis is the Eagle Nebula (M16), with its central cluster NGC6611. At 1750 pc from the Sun, it harbors 93 OB stars, together with a population of low-mass stars from embedded protostars to disk-less Class III objects, with age <= 3Myrs. We study an archival 78 ksec Chandra/ACIS-I observation of NGC6611, and two new 80ksec observations of the outer region of M16, one centered on the Column V, and one on a region of the molecular cloud with ongoing star-formation. We detect 1755 point sources, with 1183 candidate cluster members (219 disk-bearing and 964 disk-less). We study the global X-ray properties of M16 and compare them with those of the Orion Nebula Cluster. We also compare the level of X-ray emission of Class II and Class III stars, and analyze the X-ray spectral properties of OB stars. Our study supports the lower level of X-ray activity for the disk-bearing stars with respect to the disk-less members. The X-ray Luminosity Function (XLF) of M16 is similar to that of Orion, supporting the universality of the XLF in young clusters. 85% of the O stars of NGC6611 have been detected in X-rays. With only one possible exception, they show soft spectra with no hard component, indicating that mechanisms for the production of hard X-ray emission in O stars are not operating in NGC 6611.

preprint2012arXiv

Discovery of X-ray Emission in the Old Classical Nova DK Lacertae

We report the discovery of X-ray emission at the position of the old classical nova DK Lacertae using the Swift satellite. Three observations were conducted using the X-ray telescope 62 years after the discovery of the nova, yielding 46 source signals in an exposure time of 4.8 ks. A background-subtracted count rate was 9+/-2x10^{-3} counts s^{-1}, corresponding to a detection significance level of 5-sigma. The X-ray spectrum was characterized by a continuum extending up to about 7 keV, which can be modeled by a power-law component with a photon index of 1.4--5.6, or by a thermal bremsstrahlung component with a temperature of 0.7--13.3 keV, convolved with interstellar absorption with an equivalent hydrogen column density of 0.3--2.4x10^{22} cm^{-2}. Assuming a distance of 3900 pc to the source, the luminosity was 10^{32}--10^{34} ergs s^{-1} in the 0.3--10 keV energy band. The origin of X-rays is considered to be either mass accretion on the white dwarf or adiabatic shocks in nova ejecta, with the former appearing much more likely. In either case, DK Lacertae represents a rare addition to the exclusive club of X-ray emitting old novae.

preprint2012arXiv

Optical photometric GTC/OSIRIS observations of the young massive association Cygnus OB2

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.

preprint2012arXiv

Soft coronal X-rays from β Pictoris

A type stars are expected to be X-ray dark, yet weak emission has been detected from several objects in this class. We present new Chandra/HRC-I observations of the A5 V star β Pictoris. It is clearly detected with a flux of 9+-2 10^{-4} counts/s. In comparison with previous data this constrains the emission mechanism and we find that the most likely explanation is an optically thin, collisionally dominated, thermal emission component with a temperature around 1.1 MK. We interpret this component as a very cool and dim corona, with \log L_X/L_{bol}=-8.2 (0.2-2.0 keV). Thus, it seems that β Pictoris shares more characteristics with cool stars than previously thought.

preprint2012arXiv

The Cosmic Ray Intensity Near the Archean Earth

We employ three-dimensional state of the art magnetohydrodynamic models of the early solar wind and heliosphere and a two-dimensional model for cosmic ray transport to investigate the cosmic ray spectrum and flux near the Archean Earth. We assess how sensitive the cosmic ray spectrum is to changes in the sunspot placement and magnetic field strength, the large scale dipole magnetic field strength, the wind ram pressure, and the Sun's rotation period. Overall, our results confirm earlier work that suggested the Archean Earth would have experienced a greatly reduced cosmic ray flux than is the case today. The cosmic ray reduction for the early Sun is mainly due to the shorter solar rotation period and tighter winding of the Parker spiral, and to the different surface distribution of the more active solar magnetic field. These effects lead to a global reduction of the cosmic ray flux at 1AU by up to two orders of magnitude or more. Variations in the sunspot magnetic field have more effect on the flux than variations in the dipole field component. The wind ram pressure affects the cosmic ray flux through its influence on the size of the heliosphere via the pressure balance with the ambient interstellar medium. Variations in the interstellar medium pressure experienced by the solar system in orbit through Galaxy could lead to order of magnitude changes in the cosmic ray flux at Earth on timescales of a few million years.

preprint2012arXiv

The effect of limited spatial resolution of stellar surface magnetic field maps on MHD wind and coronal X-ray emission models

We study the influence of the spatial resolution on scales of $5°$ and smaller of solar surface magnetic field maps on global magnetohydrodynamic solar wind models, and on a model of coronal heating and X-ray emission. We compare the solutions driven by a low-resolution Wilcox Solar Observatory magnetic map, the same map with spatial resolution artificially increased by a refinement algorithm, and a high-resolution Solar and Heliospheric Observatory Michelson Doppler Imager map. We find that both the wind structure and the X-ray morphology are affected by the fine-scale surface magnetic structure. Moreover, the X-ray morphology is dominated by the closed loop structure between mixed polarities on smaller scales and shows significant changes between high and low resolution maps. We conclude that three-dimensional modeling of coronal X-ray emission has greater surface magnetic field spatial resolution requirements than wind modeling, and can be unreliable unless the dominant mixed polarity magnetic flux is properly resolved.

preprint2012arXiv

The Role of White Dwarfs in Cataclysmic Variable Spin-down

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én 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.

preprint2011arXiv

Close to the Dredge: Precise X-ray C and N Abundances in lambda Andromeda and its Precocious RGB Mixing Problem

Chandra LETG+HRC-S and XMM-Newton RGS spectra of H-like C and N formed in the corona of the lambda And primary star, a mildly metal-poor G8 III-IV first ascent giant that completed dredge-up ~50 Myr ago, have been used to make a precise measurement of its surface C/N ratio. We obtain the formal result [C/N]=0.03+/-0.07, which is typical of old disk giants and in agreement with standard dredge-up theory for stars of about 1 M_sun or lower. In contrast, these stars as a group, including lambda And, have 12C/13C ~< 20, which is much lower than standard model predictions. We show that the abundances of the old disk giants are consistent with models including thermohaline mixing that begins at the red giant branch luminosity function "bump". Instead, lambda And indicates that the 12C/13C anomaly can be present immediately following dredge-up, contrary to current models of extra mixing on the red giant branch.

preprint2011arXiv

From X-ray dips to eclipse: Witnessing disk reformation in the recurrent nova USco

The 10th recorded outburst of the recurrent eclipsing nova USco was observed simultaneously in X-ray, UV, and optical by XMM-Newton on days 22.9 and 34.9 after outburst. Two full passages of the companion in front of the nova ejecta were observed, witnessing the reformation of the accretion disk. On day 22.9, we observed smooth eclipses in UV and optical but deep dips in the X-ray light curve which disappeared by day 34.9, then yielding clean eclipses in all bands. X-ray dips can be caused by clumpy absorbing material that intersects the line of sight while moving along highly elliptical trajectories. Cold material from the companion could explain the absence of dips in UV and optical light. The disappearance of X-ray dips before day 34.9 implies significant progress in the formation of the disk. The X-ray spectra contain photospheric continuum emission plus strong emission lines, but no clear absorption lines. Both continuum and emission lines in the X-ray spectra indicate a temperature increase from day 22.9 to day 34.9. We find clear evidence in the spectra and light curves for Thompson scattering of the photospheric emission from the white dwarf. Photospheric absorption lines can be smeared out during scattering in a plasma of fast electrons. We also find spectral signatures of resonant line scattering that lead to the observation of the strong emission lines. Their dominance could be a general phenomenon in high-inclination systems such as Cal87.

preprint2011arXiv

Modeling the 2010 blast wave of the symbiotic-like nova V407 Cygni

(Abridged) The symbiotic-like binary Mira and nova V407 Cyg was observed in outburst on March 2010 and monitored in several wavelength bands. Here we report on multi-dimensional hydrodynamic simulations describing the 2010 outburst of V407 Cyg, exploring the first 60 days of evolution. The model takes into account thermal conduction and radiative cooling; the pre-explosion system conditions included the companion star and a circumbinary density enhancement. The simulations showed that the blast and the ejecta distribution are both aspherical due to the inhomogeneous circumstellar medium in which they expand; in particular they are significantly collimated in polar directions (producing a bipolar shock morphology) if the circumstellar envelope is characterized by an equatorial density enhancement. The blast is partially shielded by the Mira companion, producing a wake with dense and hot post-shock plasma on the rear side of the companion star; most of the X-ray emission produced during the evolution of the blast arises from this plasma structure. The observed X-ray lightcurve can be reproduced, assuming values of outburst energy and ejected mass similar to those of RS Oph and U Sco, if a circumbinary gas density enhancement is included in the model. In this case, the 2010 blast propagated through a circumbinary gas density enhancement with radius of the order of 40 AU and gas density \approx 10^6 cm^{-3} and the mass of ejecta in the outburst was M_{ej} \approx 2\times 10^{-7} M_{\odot} with an explosion energy E_{0} \approx 2\times 10^{44} erg. Alternatively, the model can produce a similar X-ray lightcurve without the need of a circumbinary gas density enhancement only if the outburst energy and ejected mass were similar to those at the upper end of ranges for classical novae, namely M_{ej} \approx 5\times 10^{-5} M_{\odot} and E_{0} \approx 5\times 10^{46} erg.

preprint2011arXiv

Multi-Wavelength Study of the 2008-2009 Outburst of V1647 Ori

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.

preprint2011arXiv

Neon Insights from Old Solar X-rays: a Plasma Temperature Dependence of the Coronal Neon Content

An analysis using modern atomic data of fluxes culled from the literature for O VIII and Ne IX lines observed in solar active regions by the P78 and Solar Maximum Mission satellites confirms that the coronal Ne/O abundance ratio varies by a factor of two or more, and finds an increase in Ne/O with increasing active region plasma temperature. The latter is reminiscent of evidence for increasing Ne/O with stellar activity in low-activity coronae that reaches a "neon saturation" in moderately active stars at approximately twice the historically accepted solar value of about 0.15 by number. We argue that neon saturation represents the underlying stellar photospheric compositions, and that low activity coronae, including that of the Sun, are generally depleted in neon. The implication would be that the solar Ne/O abundance ratio should be revised upward by a factor of about two to n(Ne)/n(O)~0.3. Diverse observations of neon in the local cosmos provide some support for such a revision. Neon would still be of some relevance for reconciling helioseismology with solar models computed using recently advocated chemical mixtures with lower metal content.

preprint2011arXiv

Starspots on the fastest rotators in the Beta Pic moving group

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.

preprint2011arXiv

Swift observations of the March 2011 outburst of the cataclysmic variable NSV 1436: a probable dwarf nova (Research Note)

The March 2011 outburst of the poorly-studied cataclysmic variable NSV 1436 offered an opportunity to decide between dwarf nova and recurrent nova classifications. We use seven daily observations in the X-ray and UV by the Swift satellite, together with AAVSO V photometry, to characterise the outburst and decline behaviour. The short optical outburst coincided with a faint and relatively soft X-ray state, whereas in decline to fainter optical magnitudes the X-ray source was harder and brighter. These attributes, and the modest optical outburst amplitude, indicate that this was a dwarf nova outburst and not a recurrent nova. The rapid optical fading suggests an orbital period below 2 hours.

preprint2011arXiv

T Tauri candidates and accretion rates using IPHAS: method and application to IC1396

The INT Photometric H-Alpha Survey (IPHAS) is a 1800 deg^2 survey of the Northern Galactic Plane, reaching down to r'~21. We demonstrate how the survey can be used to (1) reliably select classical T Tauri star candidates and (2) constrain the mass accretion rates with an estimated relative uncertainty of 0.6 dex. IPHAS is a necessary addition to spectroscopic surveys because it allows large and uniform samples of accretion rates to be obtained with a precise handle on the selection effects. We apply the method on a region of 7 deg^2 towards the HII region IC1396 in Cepheus OB2 and identify 158 pre-main sequence candidates with masses between 0.2 and 2.0 Msun and accretion rates between 10^-9.2 and 10^-7.0 Msun yr^-1. We find a power-law dependency between the stellar mass and the accretion rates with a slope of α= 1.1 +/- 0.2, which is less steep than indicated by previous studies. We discuss the influence of method-dependent systematic effects on the determination of this relationship. The majority of our sample consists of faint, previously unknown, low-mass T Tauri candidates (56 per cent between 0.2 and 0.5 Msun). Many candidates are clustered in front of three bright-rimmed molecular clouds, which are being ionized by the massive star HD206267 (O6.5V). We discover a spatio-temporal gradient of increasing accretion rates, increasing Spitzer infrared excess, and younger ages away from the ionizing star, providing a strong indication that the formation of these clusters has been sequentially triggered by HD206267 during the last ~1 Myr.

preprint2011arXiv

The Dynamics of Stellar Coronae Harboring Hot-jupiters I. A Time-dependent MHD Simulation of the Interplanetary Environment in the HD 189733 Planetary System

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.

preprint2011arXiv

The Dynamics of Stellar Coronae Harboring Hot-jupiters II. A Space Weather Event on A Hot-jupiter

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.

preprint2011arXiv

X-ray Study of Rekindled Accretion in the Classical Nova V2491 Cygni

We conducted an X-ray spectroscopic study of the classical nova V2491 Cygni using our target-of-opportunity observation data with the Suzaku and XMM-Newton satellites as well as archived data with the Swift satellite. Medium-resolution (R~10-50) spectra were obtained using the X-ray CCD spectrometers at several post-nova epochs on days 9, 29, 40, 50, and 60-150 in addition to a pre-nova interval between days -322 and -100 all relative to the time when the classical nova was spotted. We found remarkable changes in the time series of the spectra: (a) In the pre-nova phase and on day 9, the 6.7 keV emission line from Fe XXV was significantly detected. (b) On day 29, no such emission line was found. (c) On day 40, the 6.7 keV emission line emerged again. (d) On days 50 and 60-150, three emission lines at 6.4, 6.7, and 7.0 keV respectively from quasi-neutral Fe, Fe XXV, and Fe XXVI were found. Statistically significant changes of the Fe K line intensities were confirmed between day 29 and 50. Based on these phenomena, we conclude that (1) the post-nova evolution can be divided into two different phases, (2) ejecta is responsible for the X-ray emission in the earlier phase, while rekindled accretion is for the later phase, and (3) the accretion process is considered to be reestablished as early as day 50 when the quasi-neutral Fe emission line emerged, which is a common signature of accretion from magnetic cataclysmic variables.

preprint2011arXiv

XMM-Newton X-ray and UV observations of the fast nova V2491 Cyg during the supersoft source phase

Two XMM observations of the fast classical nova V2491Cyg were carried out on days 39.93 and 49.62 after discovery, during the supersoft source (SSS) phase, yielding simultaneous X-ray and UV light curves and high-resolution X-ray spectra. The first X-ray light curve is highly variable with periodic oscillations (37.2 min) after an extended dip of factor of three lasting ~3 hours. The cause of the dip is currently unexplained and could have the same origin as similar events in V4743Sgr and RSOph, as it occurred on the same time scale. The 37-min period is not present during the dip and also not in the second observation. The UV light curves are variable but contain no dips and no period. High-resolution X-ray spectra are presented for 4 intervals of different intensity. All spectra are atmospheric continua with absorption lines and absorption edges. Interstellar lines of OI and NI are seen at their rest wavelengths, and a large number of high-ionization absorption lines are found at blue shifts indicating an expansion velocity of 3000-3400 km/s, which does not change significantly during the epochs of observation. Comparisons with the slower nova V4743Sgr and the symbiotic recurrent nova RSOph are presented. The SSS spectrum of V4743Sgr is much softer with broader and more complex photospheric absorption lines. Meanwhile, the absorption lines in RSOph are as narrow as in V2491Cyg, but they are less blue shifted. A remarkable similarity in the continua of V2491Cyg and RSOph is found. The only differences are smaller line shifts and additional emission lines in RSOph that are related to a dense stellar wind from the evolved companion. Three unidentified absorption lines are present in the X-ray spectra of all three novae, with rest wavelengths 26.05AA, 29.45AA, and 30.0AA. No satisfactory spectral model is currently available for the soft X-ray spectra of novae in outburst.

preprint2010arXiv

Close Binaries with Infrared Excess: Destroyers of Worlds?

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.

preprint2010arXiv

Identification of a Population of X-ray Emitting Massive Stars in the Galactic Plane

We present X-ray, infrared, optical and radio observations of four previously unidentified Galactic plane X-ray sources, AX J163252-4746, AX J184738-0156, AX J144701-5919 and AX J144547-5931. Detection of each source with the Chandra X-ray Observatory has provided sub-arcsecond localizations, which we use to identify bright infrared counterparts to all four objects. Infrared and optical spectroscopy of these counterparts demonstrate that all four X-ray sources are extremely massive stars, with spectral classifications Ofpe/WN9 (AX J163252-4746), WN7 (AX J184738-0156 = WR121a), WN7-8h (AX J144701-5919) and OIf+ (AX J144547-5931). AX J163252-4746 and AX J184738-0156 are both luminous, hard, X-ray emitters with strong Fe XXV emission lines in their X-ray spectra at ~6.7 keV. The multi-wavelength properties of AX J163252-4746 and AX J184738-0156 are not consistent with isolated massive stars or accretion onto a compact companion; we conclude that their X-ray emission is most likely generated in a colliding-wind binary system. For both AX J144701-5919 and AX J144547-5931, the X-ray emission is an order of magnitude less luminous and with a softer spectrum. These properties are consistent with a colliding-wind binary interpretation for these two sources also, but other mechanisms for the generation of X-rays cannot be excluded. There are many other as yet unidentified X-ray sources in the Galactic plane, with X-ray properties similar to those seen for AX J163252-4746, AX J184738-0156, AX J144701-5919 and AX J144547-5931. This may indicate a substantial population of X-ray-emitting massive stars and colliding-wind binaries in the Milky Way.

preprint2010arXiv

K-Shell Photoabsorption Studies of the Carbon Isonuclear Sequence

K-shell photoabsorption cross sections for the isonuclear C I - C IV ions have been computed using the R-matrix method. Above the K-shell threshold, the present results are in good agreement with the independent-particle results of Reilman & Manson (1979). Below threshold, we also compute the strong 1s -> np absorption resonances with the inclusion of important spectator Auger broadening effects. For the lowest 1s -> 2p, 3p resonances, comparisons to available C II, C III, and C IV experimental results show good agreement in general for the resonance strengths and positions, but unexplained discrepancies exist. Our results also provide detailed information on the C I K-shell photoabsorption cross section including the strong resonance features, since very limited laboratory experimental data exist. The resultant R-matrix cross sections are then used to model the Chandra X-ray absorption spectrum of the blazar Mkn 421.

preprint2010arXiv

Magnetic Structure of Rapidly Rotating FK Comae-Type Coronae

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.

preprint2010arXiv

The Coronal Structure of AB Doradus

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én 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.

preprint2010arXiv

The interactions of winds from massive young stellar objects

The supersonic stellar and disk winds possessed by massive young stellar objects will produce shocks when they collide against the interior of a pre-existing bipolar cavity (resulting from an earlier phase of jet activity). The shock heated gas emits thermal X-rays which may be observable by spaceborne observa- tories such as the Chandra X-ray Observatory. Hydrodynamical models are used to explore the wind-cavity interaction. Radiative transfer calculations are performed on the simulation output to produce synthetic X-ray observations, allowing constraints to be placed on model parameters through comparisons with observations. The model reveals an intricate interplay between the inflowing and outflowing material and is successful in reproducing the observed X-ray count rates from massive young stellar objects.

preprint2010arXiv

The orbital period of V458 Vulpeculae, a post double common-envelope nova

We present time-resolved optical spectroscopy of V458 Vulpeculae (Nova Vul 2007 No. 1) spread over a period of 15 months starting 301 days after its discovery. Our data reveal radial velocity variations in the HeII λ5412 and HeII λ4686 emission lines. A period analysis of the radial velocity curves resulted in a period of 98.09647 \pm 0.00025 min (0.06812255 \pm 0.00000017 d) which we identify with the orbital period of the binary system. V458 Vul is therefore the planetary nebula central binary star with the shortest period known. We explore the possibility of the system being composed of a relatively massive white dwarf (M1 \gsim 1.0 Msun) accreting matter from a post-asymptotic giant branch star which produced the planetary nebula observed. In this scenario, the central binary system therefore underwent two common-envelope episodes. A combination of previous photoionisation modelling of the nebular spectra, post-asymptotic giant branch evolutionary tracks and the orbital period favour a mass of M2 \sim 0.6 Msun for the donor star. Therefore, the total mass of the system may exceed the Chandrasekhar mass, which makes V458 Vul a Type Ia supernova progenitor candidate.

preprint2010arXiv

The Super-Soft X-ray Phase of Nova RS Ophiuchi 2006

Swift X-ray observations of the ~60 day super-soft phase of the recurrent nova RS Ophiuchi 2006 show the progress of nuclear burning on the white dwarf in exquisite detail. First seen 26 days after the optical outburst, this phase started with extreme variability likely due to variable absorption, although intrinsic white dwarf variations are not excluded. About 32 days later, a steady decline in count-rate set in. NLTE model atmosphere spectral fits during the super-soft phase show that the effective temperature of the white dwarf increases from ~65 eV to ~90 eV during the extreme variability phase, falling slowly after about day 60 and more rapidly after day 80. The bolometric luminosity is seen to be approximately constant and close to Eddington from day 45 up to day 60, the subsequent decline possibly signalling the end of extensive nuclear burning. Before the decline, a multiply-periodic, ~35 s modulation of the soft X-rays was present and may be the signature of a nuclear fusion driven instability. Our measurements are consistent with a white dwarf mass near the Chandrasekhar limit; combined with a deduced accumulation of mass transferred from its binary companion, this leads us to suggest RS Oph is a strong candidate for a future supernova explosion. The main uncertainty now is whether the WD is the CO type necessary for a SN Ia. This may be confirmed by detailed abundance analyses of spectroscopic data from the outbursts.

preprint2009arXiv

Hard X-ray flux from low-mass stars in the Cygnus OB2 Association

Context. The Cygnus OB2 association, the central engine of the Cygnus X star-forming region, is the subject of an extensive INTEGRAL Key Project that will accumulate 6Ms of observations. Analysis of 2Ms of observations by De Becker and co-workers provides the most sensitive limit yet obtained on hard X-ray emission from the cluster. Aims. We investigate the X-ray emission in the 20-40 keV band expected from the flaring low-mass stellar population in Cygnus OB2. We discuss whether such emission needs to be considered in the interpretation of existing and future X-ray observations of the region, and whether such observations might provide insight into the high-energy processes on low-mass pre-main sequence stars. Methods. The total hard X-ray flux from low-mass stars is estimated by assuming the observed soft X-ray emission stems from a superposition of flares. We further assume the ratio of hard X-ray to soft X-ray emission is described by a scaling found for solar flares by Isola and co-workers. Results. We estimate the low-mass stellar hard X-ray flux in the 20-40 keV band to lie in the range ~2x10^31-6x10^32 erg/s and discuss some potential biases that might affect this result. Conclusions. Hard X-ray emission could lie at a level not much below the current observed flux upper limits for Cygnus OB2. If this emission could be detected, it would provide insight into the hard X-ray production of large flares on pre-main sequence stars. We highlight the penetrating power of hard X-rays from low-mass stellar populations as a possible pointer to our Galaxy's hidden star-forming clusters and super-clusters using more sensitive observations from future missions.

preprint2009arXiv

Interactions of the magnetospheres of stars and close-in giant planets

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.

preprint2009arXiv

Suzaku Detection of Super-hard X-ray Emission from the Classical Nova V2491 Cygni

We report the detection of super-hard (>10 keV) X-ray emission extending up to 70 keV from the classical nova V2491 Cygni using the Suzaku observatory. We conducted two ~20 ks target-of-opportunity observations 9 and 29 days after the outburst on 2008 April 11, yielding wide energy range spectra by combining the X-ray Imaging Spectrometer and the Hard X-ray Detector. On day 9, a spectrum was obtained at 1.0-70 keV with the Fe XXV Kαline feature and a very flat continuum, which is explained by thermal plasma with a 3 keV temperature and power-law emission with a photon index of 0.1 attenuated by a heavy extinction of 1.4x10^{23} cm^{-2}. The 15-70 keV luminosity at 10.5 kpc is 6x10^{35} ergs s^{-1}. The super-hard emission was not present on day 29. This is the highest energy at which X-rays have been detected from a classical nova. We argue a non-thermal origin for the emission, which suggests the presence of accelerated charged particles in the nova explosion.

preprint2009arXiv

The interactions of winds from massive young stellar objects: X-ray emission, dynamics, and cavity evolution

2D axis-symmetric hydrodynamical simulations are presented which explore the interaction of stellar and disk winds with surrounding infalling cloud material. The star, and its accompanying disk, blow winds inside a cavity cleared out by an earlier jet. The collision of the winds with their surroundings generates shock heated plasma which reaches temperatures up to ~10^8 K. Attenuated X-ray spectra are calculated from solving the equation of radiative transfer along lines-of-sight. This process is repeated at various epochs throughout the simulations to examine the evolution of the intrinsic and attenuated flux. We find that the dynamic nature of the wind-cavity interaction fuels intrinsic variability in the observed emission on timescales of several hundred years. This is principally due to variations in the position of the reverse shock which is influenced by changes in the shape of the cavity wall. The collision of the winds with the cavity wall can cause clumps of cloud material to be stripped away. Mixing of these clumps into the winds mass-loads the flow and enhances the X-ray emission measure. The position and shape of the reverse shock plays a key role in determining the strength and hardness of the X-ray emission. In some models the reverse shock is oblique to much of the stellar and disk outflows, whereas in others it is closely normal over a wide range of polar angles. For reasonable stellar and disk wind parameters the integrated count rate and spatial extent of the intensity peak for X-ray emission agree with \textit{Chandra} observations of the deeply embedded MYSOs S106 IRS4, Mon R2 IRS3 A, and AFGL 2591.(abridged)

preprint2009arXiv

The Massive Star Forming Region Cygnus OB2. I. Chandra catalogue of association members

We present a catalogue of 1696 X-ray sources detected in the massive star forming region (SFR) Cygnus OB2 and extracted from two archival Chandra observations of the center of the region. A deep source extraction routine, exploiting the low background rates of Chandra observations was employed to maximize the number of sources extracted. Observations at other wavelengths were used to identify low count-rate sources and remove likely spurious sources. Monte Carlo simulations were also used to assess the authenticity of these sources. We used a Bayesian technique to identify optical or near-IR counterparts for 1501 (89%) of our sources, using deep observations from the INT Photometric H-alpha Survey, the Two Micron All Sky Survey, and the UKIRT Infrared Deep Sky Survey-Galactic plane Survey. 755 (45%) of these objects have six-band r', H-alpha, i', J, H, and K optical and near-IR photometry. From an analysis of the Poisson false-source probabilities for each source we estimate that our X-ray catalogue includes < 1% of false sources, and an even lower fraction when only sources with optical or near-IR associations are considered. A Monte Carlo simulation of the Bayesian matching scheme allows this method to be compared to more simplified matching techniques and enables the various sources of error to be quantified. The catalogue of 1696 objects presented here includes X-ray broad band fluxes, model fits, and optical and near-IR photometry in what is one of the largest X-ray catalogue of a single SFR to date. The high number of stellar X-ray sources detected from relatively shallow observations confirms the status and importance of Cygnus OB2 as one of our Galaxy's most massive SFRs.

preprint2009arXiv

The Massive Star Forming Region, Cygnus OB2

We present results from a catalogue of 1696 X-ray point sources detected in the massive star forming region (SFR) Cygnus OB2, the majority of which have optical or near-IR associations. We derive ages of 3.5 and 5.25 Myrs for the stellar populations in our two fields, in agreement with recent studies that suggest that the central 1-3 Myr OB association is surrounded and contaminated by an older population with an age of 5-10 Myrs. The fraction of sources with proto-planetary disks, as traced by K-band excesses, are unusually low. Though this has previously been interpreted as due to the influence of the large number of OB stars in Cyg OB2, contamination from an older population of stars in the region could also be responsible. An initial mass function is derived and found to have a slope of Gamma = -1.27, in agreement with the canonical value. Finally we introduce the recently approved Chandra Cygnus OB2 Legacy Survey that will image a 1 square degree area of the Cygnus OB2 association to a depth of 120 ks, likely detecting ~10,000 stellar X-ray sources.

preprint2009arXiv

The Structure of the Outer Galactic Disc as revealed by IPHAS early A Stars

This study is an investigation of the stellar density profile of the Galactic disc in the Anticentre direction. We select over 40,000 early A stars from IPHAS photometry in the Galactic longitude range 160 < l < 200 close to the equatorial plane (-1 < b < +1). We then compare their observed reddening-corrected apparent magnitude distribution with simulated photometry obtained from parameterised models in order to set constraints on the Anticentre stellar density profile. By selecting A stars, we are appraising the properties of a population only ~100 Myrs old. We find the stellar density profile of young stars is well fit to an exponential with length scale of (3020 \pm 120_{statistical} \pm 180_{systematic}) pc, which is comparable to that obtained in earlier studies, out to a Galactocentric radius of R_T = (13.0 \pm 0.5_{statistical} \pm 0.6_{systematic}) kpc. At larger radii the rate of decline appears to increase with the scale length dropping to (1200 \pm 300_{statistical} \pm 70_{systematic}) pc. This result amounts to a refinement of the conclusions reached in previous studies that the stellar density profile is abruptly truncated. The IPHAS A star data are not compatible with models that propose a sudden change in metallicity at R_G = 10 kpc.

preprint2009arXiv

X-ray optical depth diagnostics of T Tauri accretion shocks

In classical T Tauri stars, X-rays are produced by two plasma components: a hot low-density plasma, with frequent flaring activity, and a high-density lower temperature plasma. The former is coronal plasma related to the stellar magnetic activity. The latter component, never observed in non-accreting stars, could be plasma heated by the shock formed by the accretion process. However its nature is still being debated. Our aim is to probe the soft X-ray emission from the high-density plasma component in classical T Tauri stars to check whether this is plasma heated in the accretion shock or whether it is coronal plasma. High-resolution X-ray spectroscopy allows us to measure individual line fluxes. We analyze X-ray spectra of the classical T Tauri stars MP Muscae and TW Hydrae. Our aim is to evaluate line ratios to search for optical depth effects, which are expected in the accretion-driven scenario. We also derive the plasma emission measure distributions EMD, to investigate whether and how the EMD of accreting and non accreting young stars differ. The results are compared to those obtained for the non-accreting weak-line T Tauri star TWA 5. We find evidence of resonance scattering in the strongest lines of MP Mus, supporting the idea that soft X-rays are produced by plasma heated in the accretion shock. We also find that the EMD of MP Mus has two peaks: a cool peak at temperatures expected for plasma heated in the accretion shock, and a hot peak typical of coronal plasma. The shape of the EMD of MP Mus appears to be the superposition of the EMD of a pure coronal source, like TWA 5, and an EMD alike that of TW Hydrae, which is instead dominated by shock-heated plasma.

preprint2008arXiv

Coronal Structure and Abundances in Young Fast Rotators

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.

preprint2008arXiv

Gas Phase diagnostics of Protoplanetary disk extension

We investigate the potential of using ratios of fine structure and near-infrared forbidden line transitions of atomic carbon to diagnose protoplanetary disk extension. Using results from 2D photoionisation and radiative transfer modeling of a realistic protoplanetary disk structure irradiated by X-rays from a T Tauri star, we obtain theoretical emission maps from which we construct radial distributions of the strongest emission lines produced in the disk. We show that ratios of fine structure to near-infrared forbidden line emission of atomic carbon are especially promising to constrain the minimum size of gaseous protoplanetary disks. While theoretically viable, the method presents a number of observational difficulties that are also discussed here.

preprint2006arXiv

On the Relevance of Compton Scattering for the Soft X-ray Spectra of Hot DA White Dwarfs

We re-examine the effects of Compton scattering on the emergent spectra of hot DA white dwarfs in the soft X-ray range. Earlier studies have implied that sensitive X-ray observations at wavelengths $λ< 50$ Åmight be capable of probing the flux deficits predicted by the redistribution of electron-scattered X-ray photons toward longer wavelengths. We adopt two independent numerical approaches to the inclusion of Compton scattering in the computation of pure hydrogen atmospheres in hydrostatic equilibrium. One employs the Kompaneets diffusion approximation formalism, while the other uses the cross-sections and redistribution functions of Guilbert. Models and emergent spectra are computed for stellar parameters representative of HZ 43 and Sirius B, and for models with an effective temperature $T_{\rm eff} = 100 000$ K. The differences between emergent spectra computed for Compton and Thomson scattering cases are completely negligible in the case of both HZ 43 and Sirius B models, and are also negligible for all practical purposes for models with temperatures as high as $T_{\rm eff} = 100 000$ K. Models of the soft X-ray flux from these stars are instead dominated by uncertainties in their fundamental parameters.

preprint2006arXiv

Swift observations of the 2006 outburst of the recurrent nova RS Ophiuchi: I. Early X-ray emission from the shocked ejecta and red giant wind

RS Ophiuchi began its latest outburst on 2006 February 12. Previous outbursts have indicated that high velocity ejecta interact with a pre-existing red giant wind, setting up shock systems analogous to those seen in Supernova Remnants. However, in the previous outburst in 1985, X-ray observations did not commence until 55 days after the initial explosion. Here we report on Swift observations covering the first month of the 2006 outburst with the Burst Alert (BAT) and X-ray Telescope (XRT) instruments. RS Oph was clearly detected in the BAT 14-25 keV band from t=0 to $t\sim6$ days. XRT observationsfrom 0.3-10 keV, started at 3.17 days after outburst. The rapidly evolving XRT spectra clearly show the presence of both line and continuum emission which can be fitted by thermal emission from hot gas whose characteristic temperature, overlying absorbing column, $[N_H]_W$, and resulting unabsorbed total flux decline monotonically after the first few days. Derived shock velocities are in good agreement with those found from observations at other wavelengths. Similarly, $[N_H]_W$ is in accord with that expected from the red giant wind ahead of the forward shock. We confirm the basic models of the 1985 outburst and conclude that standard Phase I remnant evolution terminated by $t\sim10$ days and the remnant then rapidly evolved to display behaviour characteristic of Phase III. Around t=26 days however, a new, luminous and highly variable soft X-ray source began to appear whose origin will be explored in a subsequent paper.

preprint2004arXiv

X-ray spectroscopy of the unsteady quiescent corona of AD Leo with Chandra

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.

preprint2003arXiv

Chandra Multi-wavelength Project (ChaMP). I. First X-ray Source Catalog

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.

preprint2003arXiv

Chandra Multi-wavelength Project (ChaMP). II. First Results of X-ray Source Properties

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.

preprint2000arXiv

First Light Measurements of Capella with the Low Energy Transmission Grating Spectrometer aboard the Chandra X-ray Observatory

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 $Δλ\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.