BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Jeremy J. Drake

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

ResearcherUnclaimed source record
astro-ph.SRastro-ph.HEastro-ph.GAastro-phastro-ph.EPastro-ph.IMApplicationsastro-ph.COphysics.data-an

Catalog footprint

What is connected

53works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

53 published item(s)

preprint2026arXiv

The Advanced X-ray Imaging Satellite (AXIS) Community Science Book

The AXIS Community Science Book represents the collective effort of 592 scientists worldwide to define the transformative science enabled by the Advanced X-ray Imaging Satellite (AXIS), a next-generation X-ray mission selected by NASA's Astrophysics Probe Program for Phase A study. AXIS will advance the legacy of high-angular-resolution X-ray astronomy with ~1.5'' imaging over a wide 24' field of view and an order of magnitude greater collecting area than Chandra in the 0.3-12 keV band. Combining sharp imaging, high throughput, and rapid response capabilities, AXIS will open new windows on virtually every aspect of modern astrophysics, exploring the birth and growth of supermassive black holes, the feedback processes that shape galaxies, the life cycles of stars and exoplanet environments, and the nature of compact stellar remnants, supernova remnants, and explosive transients. This book compiles 138 community-contributed science cases developed by five Science Working Groups focused on AGN and supermassive black holes, galaxy evolution and feedback, compact objects and supernova remnants, stellar physics and exoplanets, and time-domain and multi-messenger astrophysics. Together, these studies establish the scientific foundation for next-generation X-ray exploration in the 2030s and highlight strong synergies with facilities of the 2030s, such as JWST, Roman, Rubin/LSST, SKA, ALMA, ngVLA, and next-generation gravitational-wave and neutrino networks.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

New Clues to the Evolution of Dwarf Carbon Stars From Their Variability and X-ray Emission

As main-sequence stars with C$>$O, dwarf carbon (dC) stars are never born alone but inherit carbon-enriched material from a former asymptotic giant branch (AGB) companion. In contrast to M dwarfs in post-mass transfer binaries, C$_2$ and/or CN molecular bands allow dCs to be identified with modest-resolution optical spectroscopy, even after the AGB remnant has cooled beyond detectability. Accretion of substantial material from the AGB stars should spin up the dCs, potentially causing a rejuvenation of activity detectable in X-rays. Indeed, a few dozen dCs have recently been found to have photometric variability with periods under a day. However, most of those are likely post-common-envelope binaries (PCEBs), spin-orbit locked by tidal forces, rather than solely spun-up by accretion. Here, we study the X-ray properties of a sample of the five nearest known dCs with $Chandra$. Two are detected in X-rays, the only two for which we also detected short-period photometric variability. We suggest that the coronal activity detected so far in dCs is attributable to rapid rotation due to tidal locking in short binary orbits after a common-envelope phase, late in the thermally pulsing (TP) phase of the former C-AGB primary (TP-AGB).

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Simulating the Space Weather in the AU Mic System: Stellar Winds and Extreme Coronal Mass Ejections

Two close-in planets have been recently found around the M-dwarf flare star AU Microscopii (AU Mic). These Neptune-sized planets (AU Mic b and c) seem to be located very close to the so-called "evaporation valley" in the exoplanet population, making this system an important target for studying atmospheric loss on exoplanets. This process, while mainly driven by the high-energy stellar radiation, will be strongly mediated by the space environment surrounding the planets. Here we present an investigation on this last area, performing 3D numerical modeling of the quiescent stellar wind from AU Mic, as well as time-dependent simulations describing the evolution of a highly energetic Coronal Mass Ejection (CME) event in this system. Observational constraints on the stellar magnetic field and properties of the eruption are incorporated in our models. We carry out qualitative and quantitative characterizations of the stellar wind, the emerging CMEs, as well as the expected steady and transient conditions along the orbit of both exoplanets. Our results predict an extreme space weather for AU Mic and its planets. This includes sub-Alfvénic regions for the large majority of the exoplanet orbits, very high dynamic and magnetic pressure values in quiescence (varying within $10^{2} - 10^{5}$ times the dynamic pressure experienced by the Earth), and an even harsher environment during the passage of any escaping CME associated with the frequent flaring observed in AU Mic. These space weather conditions alone pose an immense challenge for the survival of the exoplanetary atmospheres (if any) in this system.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Factory and the Beehive. IV. A Comprehensive Study of the Rotation X-ray Activity Relation in Praesepe and the Hyades

X-ray observations of low-mass stars in open clusters are critical to understanding the dependence of magnetic activity on stellar properties and their evolution. Praesepe and the Hyades, two of the nearest, most-studied open clusters, are among the best available laboratories for examining the dependence of magnetic activity on rotation for stars with masses lower than $\approx 1\ M_{\odot}$. We present an updated study of the rotation X-ray activity relation in the two clusters. We updated membership catalogs that combine pre-Gaia catalogs with new catalogs based on Gaia Data Release 2. The resulting catalogs are the most inclusive ones for both clusters: 1739 Praesepe and 1315 Hyades stars. We collected X-ray detections for cluster members, for which we analyzed, re-analyzed, or collated data from ROSAT, the Chandra X-ray Observatory, the Neil Gehrels Swift Observatory, and XMM-Newton. We have detections for 326 Praesepe and 462 Hyades members, of which 273 and 164, respectively, have rotation periods, an increase of 6$\times$ relative to what was previously available. We find that at $\approx$700 Myr, only M dwarfs remain saturated in X-rays, with only tentative evidence for supersaturation. We also find a tight relation between the Rossby number and fractional X-ray luminosity $L_\mathrm{X}/L_\mathrm{bol}$ in unsaturated single members, suggesting a power-law index between $-3.2$ and $-3.9$. Lastly, we find no difference in the coronal parameters between binary and single members. These results provide essential insight into the relative efficiency of magnetic heating of the stars' atmospheres, thereby informing the development of robust age-rotation-activity relations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

X-rays in Cepheids: Identifying Low-Mass Companions of Intermediate-Mass Stars

X-ray observations have been made of a sample of 20 classical Cepheids, including two new observations (Polaris and {\it l} Car) reported here. The occurrence of X-ray flux around the pulsation cycle is discussed. Three Cepheids are detected ($δ$ Cep, $β$ Dor, and Polaris). X-rays have also been detected from the low--mass F, G, and K companions of 4 Cepheids (V473 Lyr, R Cru, V659 Cen, and W Sgr), and one hot companion (S Mus). Upper limits on the X-ray flux of the remaining Cepheids provide an estimate that 28\% have low mass companions. This fraction of low--mass companions in intermediate mass Cepheids is significantly lower than expected from random pairing with the field IMF. Combining the companion fraction from X-rays with that from ultraviolet observations results in a binary/multiple fraction of 57\% $\pm$12\% for Cepheids with the ratios q $>$ 0.1 and separations a $>$ 1 au. This is a lower limit since M stars are not included. X-ray observations detect less massive companions than other existing studies of intermediate mass stars. Our measured occurrence rate of unresolved, low-mass companions to Cepheids suggests that intermediate-period binaries derive from a combination of disk and core fragmentation and accretion. This yields a hybrid mass-ratio distribution that is skewed toward small values compared to a uniform distribution but is still top-heavy compared to random pairings drawn from the IMF.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Coronal Mass Ejections and Exoplanets: A Numerical Perspective

Coronal mass ejections (CMEs) are more energetic than any other class of solar phenomena. They arise from the rapid release of up to $10^{33}$ erg of magnetic energy mainly in the form of particle acceleration and bulk plasma motion. Their stellar counterparts, presumably involving much larger energies, are expected to play a fundamental role in shaping the environmental conditions around low-mass stars, in some cases perhaps with catastrophic consequences for planetary systems due to processes such as atmospheric erosion and depletion. Despite their importance, the direct observational evidence for stellar CMEs is almost non-existent. In this way, numerical simulations constitute extremely valuable tools to shed some light on eruptive behavior in the stellar regime. Here we review recent results obtained from realistic modeling of CMEs in active stars, highlighting their key role in the interpretation of currently available observational constraints. We include studies performed on M-dwarf stars, focusing on how emerging signatures in different wavelengths related to these events vary as a function of the magnetic properties of the star. Finally, the implications and relevance of these numerical results are discussed in the context of future characterization of host star-exoplanet systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Chandra High Energy Transmission Gratings Spectra of V3890 Sgr

The recurrent nova (RN) V3890 Sgr was observed during the 7th day after the onset of its most recent outburst, with the Chandra ACIS-S camera and High Energy Transmission Gratings (HETG). A rich emission line spectrum was detected, due to transitions of Fe-L and K-shell ions ranging from neon to iron. The measured absorbed flux is $\approx 10^{-10}$ erg cm$^{-2}$ s$^{-1}$ in the 1.4-15 Angstrom range (0.77-8.86 keV). The line profiles are asymmetric, blue-shifted and skewed towards the blue side, as if the ejecta moving towards us are less absorbed than the receding ones. The full width at half maximum of most emission lines is 1000-1200 km s$^{-1}$, with some extended blue wings. The spectrum is thermal and consistent with a plasma in collisional ionization equilibrium with column density 1.3 $\times 10^{22}$ cm$^{-2}$ and at least two components at temperatures of about 1 keV and 4 keV, possibly a forward and a reverse shock, or regions with differently mixed ejecta and red giant wind. The spectrum is remarkably similar to the symbiotic RNe V745 Sco and RS Oph, but we cannot distinguish whether the shocks occurred at a distance of few AU from the red giant, or near the giant's photosphere, in a high density medium containing only a small mass. The ratios of the flux in lines of aluminum, magnesium and neon relative to the flux in lines of silicon and iron probably indicate a carbon-oxygen white dwarf (CO WD).

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Investigating Opacity Modifications and Reaction Rate Uncertainties to Resolve the Cepheid Mass Discrepancy

Cepheid masses derived from pulsations or binary dynamics are generally lower than those derived from stellar evolution models. Recent efforts have been dedicated to investigating the effects of abundances, mass loss, rotation, convection and overshooting prescriptions for modifying the evolution tracks to reduce or remove this Cepheid mass discrepancy. While these approaches are promising, either alone or in combination, more work is required to distinguish between possible solutions. Here we investigate nuclear reaction rate and opacity modifications on Cepheid evolution using the MESA code. We discuss the effects of opacity increases at envelope temperatures of 200,000-400,000 K proposed to explain the pulsation properties of hybrid main-sequence beta Cep/Slowly Pulsating B (SPB) variables which will evolve into Cepheids. We make use of the RSP nonlinear radial pulsation modeling capability in MESA to calculate periods and radial velocity amplitudes of Galactic Cepheids V1334 Cyg, Polaris, and delta Cep.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Phase-modulated X-ray Emission from Cepheids due to Pulsation-Driven Shocks

Cepheids are pulsating variable stars with a periodic chromospheric response at UV wavelengths close to their minimum radius phase. Recently, an X-ray variable signature was captured in observations during the \emph{maximum} radius phase. This X-ray emission came as a surprise and is not understood. In this work, we use the modern astrophysical code, PLUTO, to investigate the effects of pulsations on Cepheid X-ray emission. We run a number of hydrodynamic numerical simulations with a variety of initial and boundary conditions in order to explore the capability of shocks to produce the observed phase-dependent X-ray behavior. Finally we use the Simulated Observations of X-ray Sources (SOXS) package to create synthetic spectra for each simulation case and link our simulations to observables. We show that, for certain conditions, we can reproduce observed X-ray fluxes at phases 0.4--0.8 when the Cepheid is at maximum radius. Our results span a wide range of mass-loss rates, $2\times10^{-13}$--$3\times10^{-8}$ $M_\odot$ yr$^{-1}$, and peak X-ray luminosities, $5\times10^{-17}$--$1.4\times10^{-12}$ erg cm$^{-2}$ s$^{-1}$. We conclude that Cepheids exhibit \textit{two component emission} with (a) shock waves being responsible for the phase dependent variable emission (phases 0.2 - 0.6), and (b) a separate quiescent mechanism being the dominant emission mechanism for the remaining phases.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Tuning the Exo-Space Weather Radio for Stellar Coronal Mass Ejections

Coronal mass ejections (CMEs) on stars other than the Sun have proven very difficult to detect. One promising pathway lies in the detection of type II radio bursts. Their appearance and distinctive properties are associated with the development of an outward propagating CME-driven shock. However, dedicated radio searches have not been able to identify these transient features in other stars. Large Alfvén speeds and the magnetic suppression of CMEs in active stars have been proposed to render stellar eruptions "radio-quiet". Employing 3D magnetohydrodynamic simulations, we study here the distribution of the coronal Alfvén speed, focusing on two cases representative of a young Sun-like star and a mid-activity M-dwarf (Proxima Centauri). These results are compared with a standard solar simulation and used to characterize the shock-prone regions in the stellar corona and wind. Furthermore, using a flux-rope eruption model, we drive realistic CME events within our M-dwarf simulation. We consider eruptions with different energies to probe the regimes of weak and partial CME magnetic confinement. While these CMEs are able to generate shocks in the corona, those are pushed much farther out compared to their solar counterparts. This drastically reduces the resulting type II radio burst frequencies down to the ionospheric cutoff, which impedes their detection with ground-based instrumentation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

X-ray Observations of the Peculiar Cepheid V473 Lyr Identify A Low-Mass Companion

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

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Deep Chandra survey of the Small Magellanic Cloud. III. Formation efficiency of High-Mass X-ray binaries

We have compiled the most complete census of High-Mass X-ray Binaries (HMXBs) in the Small Magellanic Cloud with the aim to investigate the formation efficiency of young accreting binaries in its low metallicity environment. In total, we use 127 X-ray sources with detections in our \chandra X-ray Visionary Program (XVP), supplemented by 14 additional (likely and confirmed) HMXBs identified by \cite{2016A&A...586A..81H} that fall within the XVP area, but are not either detected in our survey (9 sources) or matched with any XVP source that has at least one OB counterpart in the OGLE-III catalog (5 sources). Specifically, we examine the number ratio of the HMXBs [N(HMXBs)] to {\it (a)} the number of OB stars, {\it (b)} the local star-formation rate (SFR), and {\it (c)} the stellar mass produced during the specific star-formation burst, all as a function of the age of their parent stellar populations. Each of these indicators serves a different role, but in all cases we find that the HMXB formation efficiency increases as a function of time (following a burst of star formation) up to $\sim$40--60\,Myr, and then gradually decreases. The peak formation efficiency N(HMXB)/SFR is (49 $\pm$ 14) $[10^{-5}~{\rm M_{\odot}/yr}]^{-1}$, in good agreement with previous estimates of the average formation efficiency in the broad $\sim$20--60\,Myr age range. The frequency of HMXBs is a factor of 8$\times$ higher than at $\sim$10\,Myr, and 4$\times$ higher than at $\sim$260\,Myr, i.e. at earlier and later epochs, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Pointing Chandra Toward the Extreme Ultraviolet Fluxes of Very Low-Mass Stars

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

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

X-rays from Stars and Planetary Systems

The Chandra X-ray Observatory has completed a remarkable twenty years in orbit. A large part of the science program of Chandra during this time has involved the study of stars and their planetary systems. This primer aims to give the reader a taste of the enormous range of stellar and planetary astrophysics that Chandra has enabled. Beginning with a tour of the X-ray solar system zoo, including the stunning pulsating X-ray aurorae of Jupiter, we then move on to the hot million-degree outer atmospheres of stars like our own Sun, whose X-ray emission is driven by an internal magnetic dynamo. The same emission processes are also vigorously present in the youngest stars, and we highlight some Chandra observations and results on nascent stellar and planetary systems. Chandra surveys and high-resolution spectroscopy of massive stars have provided a new window on the means by which they scavenge X-ray emission from their radiatively-driven winds, sometimes modulating this output by strong underlying stellar magnetic fields. We touch upon the evanescent X-radiation from intermediate-mass stars before arriving at the inevitable evolutionary endpoints of all but massive stars, first in energized X-ray emitting planetary nebulae, then in the slowly cooling, soft-X-ray emitting photospheres of white dwarfs. We conclude with white dwarfs in close binary systems, rejuvenated by interaction with a companion and where accretion gives play to a new range of energetic behavior even more spectacular and cataclysmic than the coruscant astrophysical road down which they have travelled.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Collimation and asymmetry of the hot blast wave from the recurrent nova V745 Scorpii

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Solar-type dynamo behaviour in fully convective stars without a tachocline

In solar-type stars (with radiative cores and convective envelopes), the magnetic field powers star spots, flares and other solar phenomena, as well as chromospheric and coronal emission at ultraviolet to X-ray wavelengths. The dynamo responsible for generating the field depends on the shearing of internal magnetic fields by differential rotation. The shearing has long been thought to take place in a boundary layer known as the tachocline between the radiative core and the convective envelope. Fully convective stars do not have a tachocline and their dynamo mechanism is expected to be very different, although its exact form and physical dependencies are not known. Here we report observations of four fully convective stars whose X-ray emission correlates with their rotation periods in the same way as in Sun-like stars. As the X-ray activity - rotation relationship is a well-established proxy for the behaviour of the magnetic dynamo, these results imply that fully convective stars also operate a solar-type dynamo. The lack of a tachocline in fully convective stars therefore suggests that this is not a critical ingredient in the solar dynamo and supports models in which the dynamo originates throughout the convection zone.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Stellar flares and the dark energy of CMEs

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

SXP214, an X-ray Pulsar in the Small Magellanic Cloud, Crossing the Circumstellar Disk of the Companion

Located in the Small Magellanic Cloud (SMC), SXP214 is an X-ray pulsar in a high mass X-ray binary system with a Be-star companion. A recent survey of the SMC under a Chandra X-ray Visionary program found the source in a transition when the X-ray flux was on a steady rise. The Lomb-Scargle periodogram revealed a pulse period of 211.49 +/- 0.42 s, which is significantly (>5sigma) shorter than the previous measurements with XMM-Newton and RXTE. This implies that the system has gone through sudden spin-up episodes recently. The pulse profile shows a sharp eclipse-like feature with a modulation amplitude of >95%. The linear rise of the observed X-ray luminosity from <~2x to 7x10^35 erg s^-1 is correlated with steady softening of the X-ray spectrum, which can be described by the changes in the local absorption from N_H ~ 10^24 to <~10^20 cm^-2 for an absorbed power-law model. The soft X-ray emission below 2 keV was absent in the early part of the observation when only the pulsating hard X-ray component was observed, whereas at later times both soft and hard X-ray components were observed pulsating. A likely explanation is that the neutron star was initially hidden in the circumstellar disk of the companion, and later came out of the disk with the accreted material that continued fueling the observed pulsation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Missing Magnetic Morphology Term in Stellar Rotation Evolution

This study examines the relationship between magnetic field complexity and mass and angular momentum losses. Observations of open clusters have revealed a bimodal distribution of the rotation periods of solar-like stars that has proven difficult to explain under the existing rubric of magnetic braking. Recent studies suggest that magnetic complexity can play an important role in controlling stellar spin-down rates. However, magnetic morphology is still neglected in most rotation evolution models due to the difficulty of properly accounting for its effects on wind driving and angular momentum loss. Using state-of-the-art magnetohydrodynamical magnetized wind simulations we study the effect that different distributions of the magnetic flux at different levels of geometrical complexity have on mass and angular momentum loss rates. Angular momentum loss rates depend strongly on the level of complexity of the field but are independent of the way this complexity is distributed. We deduce the analytical terms representing the magnetic field morphology dependence of mass and angular momentum loss rates. We also define a parameter that best represents complexity for real stars. As a test, we use these analytical methods to estimate mass and angular momentum loss rates for $8$ stars with observed magnetograms and compare them to the simulated results. Magnetic field complexity provides a natural physical basis for stellar rotation evolution models requiring a rapid transition between weak and strong spin-down modes.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Space Weather of Proxima Centauri b

A planet orbiting in the "habitable zone" of our closest neighboring star, Proxima Centauri, has recently been discovered, and the next natural question is whether or not Proxima b is "habitable". Stellar winds are likely a source of atmospheric erosion that could be particularly severe in the case of M dwarf habitable zone planets that reside close to their parent star. Here we study the stellar wind conditions that Proxima b experiences over its orbit. We construct 3-D MHD models of the wind and magnetic field around Proxima Centauri using a surface magnetic field map for a star of the same spectral type and scaled to match the observed ~ 600 G surface magnetic field strength of Proxima. We examine the wind conditions and dynamic pressure over different plausible orbits that sample the constrained parameters of the orbit of Proxima b. For all the parameter space explored, the planet is subject to stellar wind pressures of more than 2000 times those experienced by Earth from the solar wind. During an orbit, Proxima b is also subject to pressure changes of 1 to 3 orders of magnitude on timescales of a day. Its magnetopause standoff distance consequently undergoes sudden and periodic changes by a factor of 2 to 5. Proxima b will traverse the interplanetary current sheet twice each orbit, and likely crosses into regions of subsonic wind quite frequently. These effects should be taken into account in any physically realistic assessment or prediction of its atmospheric reservoir, characteristics and loss.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Flows of X-ray gas reveal the disruption of a star by a massive black hole

Tidal forces close to massive black holes can violently disrupt stars that make a close approach. These extreme events are discovered via bright X-ray and optical/UV flares in galactic centers. Prior studies based on modeling decaying flux trends have been able to estimate broad properties, such as the mass accretion rate. Here we report the detection of flows of highly ionized X-ray gas in high-resolution X-ray spectra of a nearby tidal disruption event. Variability within the absorption-dominated spectra indicates that the gas is relatively close to the black hole. Narrow line widths indicate that the gas does not stretch over a large range of radii, giving a low volume filling factor. Modest outflow speeds of a few hundred kilometers per second are observed, significantly below the escape speed from the radius set by variability. The gas flow is consistent with a rotating wind from the inner, super-Eddington region of a nascent accretion disk, or with a filament of disrupted stellar gas near to the apocenter of an elliptical orbit. Flows of this sort are predicted by fundamental analytical theory and more recent numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Magnetic complexity as an explanation for bimodal rotation populations among young stars

Observations of young open clusters have revealed a bimodal distribution of fast and slower rotation rates that has proven difficult to explain with predictive models of spin down that depend on rotation rates alone. The Metastable Dynamo Model proposed recently by Brown, employing a stochastic transition probability from slow to more rapid spin down regimes, appears to be more successful but lacks a physical basis for such duality. Using detailed 3D MHD wind models computed for idealized multipole magnetic fields, we show that surface magnetic field complexity can provide this basis. Both mass and angular momentum losses decline sharply with increasing field complexity. Combined with observation evidence for complex field morphologies in magnetically active stars, our results support a picture in which young, rapid rotators lose angular momentum in an inefficient way because of field complexity. During this slow spin-down phase, magnetic complexity is eroded, precipitating a rapid transition from weak to strong wind coupling.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Simulating the sensitivity to stellar point sources of Chandra X-ray observations

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The dependence of stellar mass and angular momentum losses on latitude and on active region and dipolar magnetic fields

Rotation evolution of late-type stars is dominated by magnetic braking and the underlying factors that control this angular momentum loss are important for the study of stellar spin-down. In this work, we study angular momentum loss as a function of two different aspects of magnetic activity using a calibrated Alfvén wave-driven magnetohydrodynamic wind model: the strengths of magnetic spots and their distribution in latitude. By driving the model using solar and modified solar surface magnetograms, we show that the topology of the field arising from the net interaction of both small-scale and large-scale field is important for spin-down rates and that angular momentum loss is not a simple function of large scale magnetic field strength. We find that changing the latitude of magnetic spots can modify mass and angular momentum loss rates by a factor of two. The general effect that causes these differences is the closing down of large-scale open field at mid- and high-latitudes by the addition of the small-scale field. These effects might give rise to modulation of mass and angular momentum loss through stellar cycles, and present a problem for ab initio attempts to predict stellar spin-down based on wind models. For all the magnetogram cases considered here, from dipoles to various spotted distributions, we find that angular momentum loss is dominated by the mass loss at mid-latitudes. The spin-down torque applied by magnetized winds therefore acts at specific latitudes and is not evenly distributed over the stellar surface, though this aspect is unlikely to be important for understanding spin-down and surface flows on stars.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

X-ray Evidence for a Pole-Dominated Corona on AB Dor

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A 33 year constancy of the X-ray coronae of AR Lac and eclipse diagnosis of scale height

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Grid of MHD Models for Stellar Mass Loss and Spin-down Rates of Solar Analogs

Stellar winds are believed to be the dominant factor in spin down of stars over time. However, stellar winds of solar analogs are poorly constrained due to the challenges in observing them. A great improvement has been made in the last decade in our understanding of the mechanisms responsible for the acceleration of the solar wind and in the development of numerical models for solar and stellar winds. In this paper, we present a grid of Magnetohydrodynamic (MHD) models to study and quantify the values of stellar mass-loss and angular momentum loss rates as a function of the stellar rotation period, magnetic dipole component, and coronal base density. We derive simple scaling laws for the loss rates as a function of these parameters, and constrain the possible mass-loss rate of stars with thermally-driven winds. Despite the success of our scaling law in matching the results of the model, we find a deviation between the "solar dipole" case and a real case based on solar observations that overestimates the actual solar mass-loss rate by a factor of 3. This implies that the model for stellar fields might require a further investigation with higher complexity which might include the use of a filling factor for active regions, as well as the distribution of the strength of the small-scale fields. Mass loss rates in general are largely controlled by the magnetic field strength, with the wind density varying in proportion to the confining magnetic pressure $B^2$. We also find that the mass-loss rates obtained using our grid models drop much faster with the increase in rotation period than scaling laws derived using observed stellar activity. For main-sequence solar-like stars, our scaling law for angular momentum loss vs. poloidal magnetic field strength retrieves the well-known Skumanich decline of angular velocity with time.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A magnetic accretion switch in pre-cataclysmic binaries?

We have investigated the mass accretion rate implied by published surface abundances of Si and C in the white dwarf component of the 3.62 hr period pre-cataclysmic binary and planet host candidate QS Vir (DA+M2-4). Diffusion timescales for gravitational settling imply $\dot{M} \sim 10^{-16}M_\odot$ yr$^{-1}$ for the 1999 epoch of the observations, which is three orders of magnitude lower than measured from a 2006 {\it XMM-Newton} observation. This is the first time that large accretion rate variations have been seen in a detached pre-CV. A third body in a 14 yr eccentric orbit suggested in a recent eclipse timing study is too distant to perturb the central binary sufficiently to influence accretion. A hypothetical coronal mass ejection just prior to the {\it XMM-Newton} observation might explain the higher accretion rate, but the implied size and frequency of such events appear too great. We suggest accretion is most likely modulated by a magnetic cycle on the secondary acting as a wind "accretion switch", a mechanism that can be tested by X-ray and ultraviolet monitoring. If so, QS Vir and similar pre-CVs could provide powerful insights into hitherto inscrutable cataclysmic variable and M dwarf magnetospheres, and mass and angular momentum loss rates.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Burn Out or Fade Away? On the X-ray and Magnetic Death of Intermediate Mass Stars

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Magnetic Modulation of Stellar Angular Momentum Loss

Angular Momentum Loss is important for understanding astrophysical phenomena such as stellar rotation, magnetic activity, close binaries, and cataclysmic variables. Magnetic breaking is the dominant mechanism in the spin down of young late-type stars. We have studied angular momentum loss as a function of stellar magnetic activity. We argue that the complexity of the field and its latitudinal distribution are crucial for angular momentum loss rates. In this work we discuss how angular momentum is modulated by magnetic cycles, and how stellar spin down is not just a simple function of large scale magnetic field strength.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Pan-chromatic observations of the remarkable nova LMC 2012

We present the results of an intensive multiwavelength campaign on nova LMC 2012. This nova evolved very rapidly in all observed wavelengths. The time to fall two magnitudes in the V band was only 2 days. In X-rays the super soft phase began 13$\pm$5 days after discovery and ended around day 50 after discovery. During the super soft phase, the \Swift/XRT and \Chandra\ spectra were consistent with the underlying white dwarf being very hot, $\sim$ 1 MK, and luminous, $\sim$ 10$^{38}$ erg s$^{-1}$. The UV, optical, and near-IR photometry showed a periodic variation after the initial and rapid fading had ended. Timing analysis revealed a consistent 19.24$\pm$0.03 hr period in all UV, optical, and near-IR bands with amplitudes of $\sim$ 0.3 magnitudes which we associate with the orbital period of the central binary. No periods were detected in the corresponding X-ray data sets. A moderately high inclination system, $i$ = 60$\pm$10$^{\arcdeg}$, was inferred from the early optical emission lines. The {\it HST}/STIS UV spectra were highly unusual with only the \ion{N}{5} (1240Å) line present and superposed on a blue continuum. The lack of emission lines and the observed UV and optical continua from four epochs can be fit with a low mass ejection event, $\sim$ 10$^{-6}$ M$_{\odot}$, from a hot and massive white dwarf near the Chandrasekhar limit. The white dwarf, in turn, significantly illuminated its subgiant companion which provided the bulk of the observed UV/optical continuum emission at the later dates. The inferred extreme white dwarf characteristics and low mass ejection event favor nova LMC 2012 being a recurrent nova of the U Sco subclass.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Chandra Cygnus OB2 Legacy Survey: Design and X-ray Point Source Catalog

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Constraints on Massive Star Formation: Cygnus OB2 was always an Association

We examine substructure and mass segregation in the massive OB association Cygnus OB2 to better understand its initial conditions. Using a well understood Chandra X-ray selected sample of young stars we find that Cyg OB2 exhibits considerable physical substructure and has no evidence for mass segregation, both indications that the association is not dynamically evolved. Combined with previous kinematical studies we conclude that Cyg OB2 is dynamically very young, and what we observe now is very close to its initial conditions: Cyg OB2 formed as a highly substructured, unbound association with a low volume density (< 100 stars/pc^3). This is inconsistent with the idea that all stars form in dense, compact clusters. The massive stars in Cyg OB2 show no evidence for having formed particularly close to one another, nor in regions of higher than average density. Since Cyg OB2 contains stars as massive as ~100 Mo this result suggests that very massive stars can be born in relatively low-density environments. This would imply that the massive stars in Cyg OB2 did not form by competitive accretion, or by mergers.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Implications of mass and energy loss due to coronal mass ejections on magnetically-active stars

Analysis of a database of solar coronal mass ejections (CMEs) and associated flares over the period 1996-2007 finds well-behaved power law relationships between the 1-8 AA flare X-ray fluence and CME mass and kinetic energy. We extrapolate these relationships to lower and higher flare energies to estimate the mass and energy loss due to CMEs from stellar coronae, assuming that the observed X-ray emission of the latter is dominated by flares with a frequency as a function of energy dn/dE=kE^-alpha. For solar-like stars at saturated levels of X-ray activity, the implied losses depend fairly weakly on the assumed value of alpha and are very large: M_dot ~ 5x10^-10 M_sun/yr and E_dot ~ 0.1L_sun. In order to avoid such large energy requirements, either the relationships between CME mass and speed and flare energy must flatten for X-ray fluence >~ 10^31 erg, or the flare-CME association must drop significantly below 1 for more energetic events. If active coronae are dominated by flares, then the total coronal energy budget is likely to be up to an order of magnitude larger than the canonical 10^-3 L_bol X-ray saturation threshold. This raises the question of what is the maximum energy a magnetic dynamo can extract from a star? For an energy budget of 1% of L_bol, the CME mass loss rate is about 5x10^-11 M_sun/yr.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Ionized Nebula surrounding the Red Supergiant W26 in Westerlund 1

We present Hαimages of an ionized nebula surrounding the M2-5Ia red supergiant (RSG) W26 in the massive star cluster Westerlund 1. The nebula consists of a circumstellar shell or ring ~0.1pc in diameter and a triangular nebula ~0.2pc from the star that in high-resolution Hubble Space Telescope images shows a complex filamentary structure. The excitation mechanism of both regions is unclear since RSGs are too cool to produce ionizing photons and we consider various possibilities. The presence of the nebula, high stellar luminosity and spectral variability suggest that W26 is a highly evolved RSG experiencing extreme levels of mass-loss. As the only known example of an ionized nebula surrounding a RSG W26 deserves further attention to improve our understanding of the final evolutionary stages of massive stars.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Detection of accretion X-rays from QS Vir: cataclysmic or a lot of hot air?

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Multi-wavelength Observations of the Radio Magnetar PSR J1622-4950 and Discovery of its Possibly Associated Supernova Remnant

We present multi-wavelength observations of the radio magnetar PSR J1622-4950 and its environment. Observations of PSR J1622-4950 with Chandra (in 2007 and 2009) and XMM (in 2011) show that the X-ray flux of PSR J1622-4950 has decreased by a factor of ~50 over 3.7 years, decaying exponentially with a characteristic time of 360 +/- 11 days. This behavior identifies PSR J1622-4950 as a possible addition to the small class of transient magnetars. The X-ray decay likely indicates that PSR J1622-4950 is recovering from an X-ray outburst that occurred earlier in 2007, before the 2007 Chandra observations. Observations with the Australia Telescope Compact Array show strong radio variability, including a possible radio flaring event at least one and a half years after the 2007 X-ray outburst that may be a direct result of this X-ray event. Radio observations with the Molonglo Observatory Synthesis Telescope reveal that PSR J1622-4950 is 8' southeast of a diffuse radio arc, G333.9+0.0, which appears non-thermal in nature and which could possibly be a previously undiscovered supernova remnant. If G333.9+0.0 is a supernova remnant then the estimates of its size and age, combined with the close proximity and reasonable implied velocity of PSR J1622-4950, suggests that these two objects could be physically associated.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Photoevaporating Proplyd-like objects in Cygnus OB2

We report the discovery of ten proplyd-like objects in the vicinity of the massive OB association Cygnus OB2. They were discovered in IPHAS H-Alpha images and are clearly resolved in broad-band HST/ACS, near-IR and Spitzer mid-IR images. All exhibit the familiar tadpole shape seen in photoevaporating objects such as the Orion proplyds, with a bright ionization front at the head facing the central cluster of massive stars, and a tail stretching in the opposite direction. Many also show secondary ionization fronts, complex tail morphologies or multiple heads. We consider the evidence that these are either proplyds or `evaporating gaseous globules' (EGGs) left over from a fragmenting molecular cloud, but find that neither scenario fully explains the observations. Typical sizes are 50,000--100,000 AU, larger than the Orion proplyds, but in agreement with the theoretical scaling of proplyd size with distance from the ionizing source. These objects are located at projected separations of 6-14pc from the OB association, compared to 0.1pc for the Orion proplyds, but are clearly being photoionized by the 65 O-type stars in Cyg OB2. Central star candidates are identified in near- and mid-IR images, supporting the proplyd scenario, though their large sizes and notable asymmetries is more consistent with the EGG scenario. A third possibility is therefore considered, that these are a unique class of photoevaporating partially-embedded young stellar objects that have survived the destruction of their natal molecular cloud. This has implications for the properties of stars that form in the vicinity of massive stars.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Dynamics of an Expanding OB Association

We present 3-dimensional kinematical observations of the massive OB association Cygnus OB2 to identify the mechanisms responsible for disrupting young star clusters. The picture revealed by these observations is of a highly-substructured, dynamically unmixed OB association that does not exhibit the position-velocity correlations predicted by the theories of infant mortality or tidal stripping. These observations would appear to support a picture of hierarchical star formation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Stellar Activity - Rotation Relationship

Using a new catalog of 824 solar and late-type stars with X-ray luminosities and rotation periods we have studied the relationship between rotation and stellar activity. From an unbiased subset of this sample the power law slope of the unsaturated regime, $L_X/L_{bol}\propto Ro^β$, is fit as $β=-2.70\pm0.13$. This is inconsistent with the canonical $β=-2$ slope to a confidence of 5$σ$ and argues for an interface-type dynamo. Super-saturation is observed for the fastest rotators in our sample and its parametric dependencies are explored. Significant correlations are found with both the corotation radius and the excess polar updraft, the latter theory being supported by other observations. We also present a new X-ray population synthesis model of the mature stellar component of our Galaxy and use it to reproduce deep observations of a high Galactic latitude field. The model, XStar, can be used to test models of stellar spin-down and dynamo decay, as well as for estimating stellar X-ray contamination rates for non-stellar studies.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Accounting for Calibration Uncertainties in X-ray Analysis: Effective Areas in Spectral Fitting

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Chandra observations of Cygnus OB2

Cygnus OB2 is the nearest example of a massive star forming region, containing over 50 O-type stars and hundreds of B-type stars. We have analyzed two Chandra pointings in Cyg OB2, detecting ~1700 X-ray sources, of which ~1450 are thought to be members of the association. Optical and near-IR photometry has been obtained for ~90% of these sources from recent deep Galactic plane surveys. We have performed isochrone fits to the near-IR color-magnitude diagram, deriving ages of 3.5(+0.75,-1.0) and 5.25(+1.5,-1.0) Myrs for sources in the two fields, both with considerable spreads around the pre-MS isochrones. The presence of a second population in the region, somewhat older than the present-day O-type stars, has been suggested by other authors and fits with the ages derived here. The fraction of sources with inner circumstellar disks (as traced by the K-band excess) is found to be very low, but appropriate for a population of age ~5 Myrs. We measure the stellar mass functions and find a power-law slope of Gamma = -1.09 +/- 0.13, in good agreement with the global mean value estimated by Kroupa. A steepening of the mass function at high masses is observed and we suggest this is due to the presence of the previous generation of stars that have lost their most massive members. Finally, combining our mass function and an estimate of the radial density profile of the association suggests a total mass of Cyg OB2 of ~30,000 Msun, similar to that of many of our Galaxy's most massive star forming regions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Modeling the Young Sun's Solar Wind and its Interaction with Earth's Paleomagnetosphere

We present a focused parameter study of solar wind - magnetosphere interaction for the young Sun and Earth, $~3.5$ Ga ago, that relies on magnetohydrodynamic (MHD) simulations for both the solar wind and the magnetosphere. By simulating the quiescent young Sun and its wind we are able to propagate the MHD simulations up to Earth's magnetosphere and obtain a physically realistic solar forcing of it. We assess how sensitive the young solar wind is to changes in the coronal base density, sunspot placement and magnetic field strength, dipole magnetic field strength and the Sun's rotation period. From this analysis we obtain a range of plausible solar wind conditions the paleomagnetosphere may have been subject to. Scaling relationships from the literature suggest that a young Sun would have had a mass flux different from the present Sun. We evaluate how the mass flux changes with the aforementioned factors and determine the importance of this and several other key solar and magnetospheric variables with respect to their impact on the paleomagnetosphere. We vary the solar wind speed, density, interplanetary magnetic field strength and orientation as well as Earth's dipole magnetic field strength and tilt in a number of steady-state scenarios that are representative of young Sun-Earth interaction. This study is done as a first step of a more comprehensive effort towards understanding the implications of Sun-Earth interaction for planetary atmospheric evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Swift X-Ray Observations of Classical Novae. II. The Super Soft Source sample

The Swift GRB satellite is an excellent facility for studying novae. Its rapid response time and sensitive X-ray detector provides an unparalleled opportunity to investigate the previously poorly sampled evolution of novae in the X-ray regime. This paper presents Swift observations of 52 Galactic/Magellanic Cloud novae. We included the XRT (0.3-10 keV) X-ray instrument count rates and the UVOT (1700-8000 Angstroms) filter photometry. Also included in the analysis are the publicly available pointed observations of 10 additional novae the X-ray archives. This is the largest X-ray sample of Galactic/Magellanic Cloud novae yet assembled and consists of 26 novae with super soft X-ray emission, 19 from Swift observations. The data set shows that the faster novae have an early hard X-ray phase that is usually missing in slower novae. The Super Soft X-ray phase occurs earlier and does not last as long in fast novae compared to slower novae. All the Swift novae with sufficient observations show that novae are highly variable with rapid variability and different periodicities. In the majority of cases, nuclear burning ceases less than 3 years after the outburst begins. Previous relationships, such as the nuclear burning duration vs. t_2 or the expansion velocity of the eject and nuclear burning duration vs. the orbital period, are shown to be poorly correlated with the full sample indicating that additional factors beyond the white dwarf mass and binary separation play important roles in the evolution of a nova outburst. Finally, we confirm two optical phenomena that are correlated with strong, soft X-ray emission which can be used to further increase the efficiency of X-ray campaigns.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Testing EUV/X-ray Atomic Data for the Solar Dynamics Observatory

The Atmospheric Imaging Assembly (AIA) and the Exteme-ultraviolet Variability Experiment (EVE) onboard the Solar Dynamics Observatory include spectral windows in the X-ray/EUV band. Accuracy and completeness of the atomic data in this wavelength range is essential for interpretation of the spectrum and irradiance of the solar corona, and of SDO observations made with the AIA and EVE instruments. Here we test the X-ray/EUV data in the CHIANTI database to assess their completeness and accuracy in the SDO bands, with particular focus on the 94A and 131A AIA passbands. Given the paucity of solar observations adequate for this purpose, we use high-resolution X-ray spectra of the low-activity solar-like corona of Procyon obtained with the Chandra Low Energy Transmission Grating Spectrometer (LETGS). We find that while spectral models overall can reproduce quite well the observed spectra in the soft X-ray range ll < 50A, and at the EUV wavelengths ll >130A, they significantly underestimate the observed flux in the 50-130A wavelength range. The model underestimates the observed flux by a variable factor ranging from \approx 1.5, at short wavelengths below \sim50A, up to \approx5-7 in the \sim 70-125A range. In the AIA bands covered by LETGS, i.e. 94A and 131A, we find that the observed flux can be underestimated by large factors (\sim 3 and \sim 1.9 respectively, for the case of Procyon presented here). We discuss the consequences for analysis of AIA data and possible empirical corrections to the AIA responses to model more realistically the coronal emission in these passbands.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The stellar activity-rotation relationship and the evolution of stellar dynamos

We present a sample of 824 solar and late-type stars with X-ray luminosities and rotation periods. This is used to study the relationship between rotation and stellar activity and derive a new estimate of the convective turnover time. From an unbiased subset of this sample the power law slope of the unsaturated regime, L_X / L_bol = Ro^β, is fit as β= -2.70 +/- 0.13. This is inconsistent with the canonical β=-2 slope to a confidence of 5 sigma, and argues for an additional term in the dynamo number equation. From a simple scaling analysis this implies ΔΩ/ Ω= Ω^0.7, i.e. the differential rotation of solar-type stars gradually declines as they spin down. Super-saturation is observed for the fastest rotators in our sample and its parametric dependencies are explored. Significant correlations are found with both the corotation radius and the excess polar updraft, the latter theory providing a stronger dependence and being supported by other observations. We estimate mass-dependent empirical thresholds for saturation and super-saturation and map out three regimes of coronal emission. Late F-type stars are shown never to pass through the saturated regime, passing straight from super-saturated to unsaturated X-ray emission. The theoretical threshold for coronal stripping is shown to be significantly different from the empirical saturation threshold (Ro < 0.13), suggesting it is not responsible. Instead we suggest that a different dynamo configuration is at work in stars with saturated coronal emission. This is supported by a correlation between the empirical saturation threshold and the time when stars transition between convective and interface sequences in rotational spin-down models.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Stellar X-ray sources in the Chandra COSMOS survey

We present an analysis of the X-ray properties of a sample of solar- and late-type field stars identified in the Chandra Cosmic Evolution Survey (COSMOS), a deep (160ks) and wide (0.9 deg2) extragalactic survey. The sample of 60 sources was identified using both morphological and photometric star/galaxy separation methods. We determine X-ray count rates, extract spectra and light curves and perform spectral fits to determine fluxes and plasma temperatures. Complementary optical and near-IR photometry is also presented and combined with spectroscopy for 48 of the sources to determine spectral types and distances for the sample. We find distances ranging from 30pc to ~12kpc, including a number of the most distant and highly active stellar X-ray sources ever detected. This stellar sample extends the known coverage of the L_X-distance plane to greater distances and higher luminosities, but we do not detect as many intrinsically faint X-ray sources compared to previous surveys. Overall the sample is typically more luminous than the active Sun, representing the high-luminosity end of the disk and halo X-ray luminosity functions. The halo population appears to include both low-activity spectrally hard sources that may be emitting through thermal bremsstrahlung, as well as a number of highly active sources in close binaries.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Massive Star Forming Region Cygnus OB2. II. Integrated Stellar Properties and the Star Formation History

Cygnus OB2 is the nearest example of a massive star forming region, containing over 50 O-type stars and hundreds of B-type stars. We have analysed the properties of young stars in two fields in Cyg OB2 using the recently published deep catalogue of Chandra X-ray point sources with complementary optical and near-IR photometry. Our sample is complete to 1 Msun (excluding A and B-type stars that do not emit X-rays), making this the deepest study of the stellar properties and star formation history in Cyg OB2 to date. From Siess et al. (2000) isochrone fits to the near-IR color-magnitude diagram, we derive ages of 3.5 (+0.75/-1.0) and 5.25 (+1.5/-1.0) Myrs for sources in the two fields, both with considerable spreads around the pre-MS isochrones. The presence of a stellar population somewhat older than the present-day O-type stars, also fits in with the low fraction of sources with inner circumstellar disks (as traced by the K-band excess) that we find to be very low, but appropriate for a population of age ~5 Myrs. We also find that the region lacks a population of highly embedded sources that is often observed in young star forming regions, suggesting star formation in the vicinity has declined. We measure the stellar mass functions in this limit and find a power-law slope of Gamma = -1.09 +/- 0.13, in good agreement with the global mean value estimated by Kroupa (2002). A steepening of the slope at higher masses is observed and suggested as due to the presence of the previous generation of stars that have lost their most massive members. Finally, combining our mass function and an estimate of the radial density profile of the association suggests a total mass of Cyg OB2 of ~3 x 10^4 Msun, similar to that of many of our Galaxy's most massive star forming regions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

X-ray Flares of EV Lac: Statistics, Spectra, Diagnostics

We study the spectral and temporal behavior of X-ray flares from the active M-dwarf EV Lac in 200 ks of exposure with the Chandra/HETGS. We derive flare parameters by fitting an empirical function which characterizes the amplitude, shape, and scale. The flares range from very short (<1 ks) to long (10 ks) duration events with a range of shapes and amplitudes for all durations. We extract spectra for composite flares to study their mean evolution and to compare flares of different lengths. Evolution of spectral features in the density-temperature plane shows probable sustained heating. The short flares are significantly hotter than the longer flares. We determined an upper limit to the Fe K fluorescent flux, the best fit value being close to what is expected for compact loops.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

X-ray irradiated protoplanetary disk atmospheres II: Predictions from models in hydrostatic equilibrium

We present new models for the X-ray photoevaporation of circumstellar discs which suggest that the resulting mass loss (occurring mainly over the radial range 10-40 AU) may be the dominant dispersal mechanism for gas around low mass pre-main sequence stars, contrary to the conclusions of previous workers. Our models combine use of the MOCASSIN Monte Carlo radiative transfer code and a self-consistent solution of the hydrostatic structure of the irradiated disc. We estimate the resulting photoevaporation rates assuming sonic outflow at the surface where the gas temperature equals the local escape temperature and derive mass loss rates of ~10^{-9} M_sun/yr, typically a factor 2-10 times lower than the corresponding rates in our previous work (Ercolano et al., 2008) where we did not adjust the density structure of the irradiated disc. The somewhat lower rates, and the fact that mass loss is concentrated towards slightly smaller radii, result from the puffing up of the heated disc at a few AU which partially screens the disc at tens of AU. (.....) We highlight the fact that X-ray photoevaporation has two generic advantages for disc dispersal compared with photoevaporation by extreme ultraviolet (EUV) photons that are only modestly beyond the Lyman limit: the demonstrably large X-ray fluxes of young stars even after they have lost their discs and the fact that X-rays are effective at penetrating much larger columns of material close to the star (abridged).

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

X-ray Photoevaporation-starved T Tauri Accretion

X-ray luminosities of accreting T Tauri stars are observed to be systematically lower than those of non-accretors. There is as yet no widely accepted physical explanation for this effect, though it has been suggested that accretion somehow suppresses, disrupts or obscures coronal X-ray activity. Here, we suggest that the opposite might be the case: coronal X-rays modulate the accretion flow. We re-examine the X-ray luminosities of T Tauri stars in the Orion Nebula Cluster and find that not only are accreting stars systematically fainter, but that there is a correlation between mass accretion rate and stellar X-ray luminosity. We use the X-ray heated accretion disk models of Ercolano et al. to show that protoplanetary disk photoevaporative mass loss rates are strongly dependent on stellar X-ray luminosity and sufficiently high to be competitive with accretion rates. X-ray disk heating appears to offer a viable mechanism for modulating the gas accretion flow and could be at least partially responsible for the observed correlation between accretion rates and X-ray luminosities of T Tauri stars.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

X-Ray Spectroscopy of the Classical Nova V458 Vulpeculae with Suzaku

We conducted a target of opportunity X-ray observation of the classical nova V458 Vulpeculae 88 days after the explosion using the Suzaku satellite. With a 20 ks exposure, the X-ray Imaging Spectrometer detected X-ray emission significantly harder than typical super-soft source emission. The X-ray spectrum shows K lines from N, Ne, Mg, Si, and S, and L-series emission from Fe in highly ionized states. The spectrum can be described by a single temperature (0.64 keV) thin thermal plasma model in collisional equilibrium with a hydrogen-equivalent extinction column density of ~3e21/cm2, a flux of ~1e-12 erg/s/cm2, and a luminosity of ~6e34 erg/s in the 0.3-3.0 keV band at an assumed distance of 13 kpc. We found a hint of an enhancement of N and deficiencies of O and Fe relative to other metals. The observed X-ray properties can be interpreted as the emission arising from shocks of ejecta from an ONe-type nova.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

EUV Flare Activity in Late-Type Stars

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

0 citations0 reviewsNo code linkedOpen access