Source author record

David H. Cohen

David H. Cohen 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

20works
5topics
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

20 published item(s)

preprint2022arXiv

Helium-like X-ray line complexes show that the hottest plasma on the O supergiant zeta Puppis is in its wind

We present an analysis of Chandra grating spectra of key helium-like line complexes to put constraints on the location with respect to the photosphere of the hottest ($T \gtrsim{6 \times 10^6}$ K) plasma in the wind of the O supergiant zeta Pup and to explore changes in the 18 years between two sets of observations of this star. We fit two models -- one empirical and one wind-shock-based -- to the S XV, Si XIII, and Mg XI line complexes and show that an origin in the wind flow, above $r \approx 1.5$ R$_{\ast}$, is strongly favored over an origin less than 0.3 R$_{\ast}$ above the photosphere ($r \lesssim 1.3$ R$_{\ast}$), especially in the more recent, very long-exposure data set. There is a modest increase in the line and continuum fluxes, line widths, wind absorption signatures, and of the hot plasma's distance from the photosphere in the 18 years since the first Chandra grating observation of zeta Pup. Both modes of modeling include the effects of dielectronic recombination satellite emission line blending on the helium-like complexes -- the first time this has been accounted for in the analysis of He-like line ratios in O stars.

preprint2016arXiv

A JVLA survey of the high frequency radio emission of the massive magnetic B- and O-type stars

We conducted a survey of seven magnetic O and eleven B-type stars with masses above $8M_{\odot}$ using the Very Large Array in the 1cm, 3cm and 13cm bands. The survey resulted in a detection of two O and two B-type stars. While the detected O-type stars - HD 37742 and HD 47129 - are in binary systems, the detected B-type stars, HD 156424 and ALS 9522, are not known to be in binaries. All four stars were detected at 3cm, whereas three were detected at 1cm and only one star was detected at 13cm. The detected B-type stars are significantly more radio luminous than the non-detected ones, which is not the case for O-type stars. The non-detections at 13cm are interpreted as due to thermal free-free absorption. Mass-loss rates were estimated using 3cm flux densities and were compared with theoretical mass-loss rates, which assume free-free emission. For HD 37742, the two values of the mass-loss rates were in good agreement, possibly suggesting that the radio emission for this star is mainly thermal. For the other three stars, the estimated mass-loss rates from radio observations were much higher than those expected from theory, suggesting either a possible contribution from non- thermal emission from the magnetic star or thermal or non-thermal emission due to interacting winds of the binary system, especially for HD 47129. All the detected stars are predicted to host centrifugal magnetospheres except HD 37742, which is likely to host a dynamical magnetosphere. This suggests that non-thermal radio emission is favoured in stars with centrifugal magnetospheres.

preprint2016arXiv

An `Analytic Dynamical Magnetosphere' formalism for X-ray and optical emission from slowly rotating magnetic massive stars

Slowly rotating magnetic massive stars develop "dynamical magnetospheres" (DM's), characterized by trapping of stellar wind outflow in closed magnetic loops, shock heating from collision of the upflow from opposite loop footpoints, and subsequent gravitational infall of radiatively cooled material. In 2D and 3D magnetohydrodynamic (MHD) simulations the interplay among these three components is spatially complex and temporally variable, making it difficult to derive observational signatures and discern their overall scaling trends.Within a simplified, steady-state analysis based on overall conservation principles, we present here an "analytic dynamical magnetosphere" (ADM) model that provides explicit formulae for density, temperature and flow speed in each of these three components -- wind outflow, hot post-shock gas, and cooled inflow -- as a function of colatitude and radius within the closed (presumed dipole) field lines of the magnetosphere. We compare these scalings with time-averaged results from MHD simulations, and provide initial examples of application of this ADM model for deriving two key observational diagnostics, namely hydrogen H-alpha emission line profiles from the cooled infall, and X-ray emission from the hot post-shock gas. We conclude with a discussion of key issues and advantages in applying this ADM formalism toward derivation of a broader set of observational diagnostics and scaling trends for massive stars with such dynamical magnetospheres.

preprint2016arXiv

KELT-17b: A hot-Jupiter transiting an A-star in a misaligned orbit detected with Doppler tomography

We present the discovery of a hot-Jupiter transiting the V=9.23 mag main-sequence A-star KELT-17 (BD+14 1881). KELT-17b is a 1.31 -0.29/+0.28 Mj, 1.525 -0.060/+0.065 Rj hot-Jupiter in a 3.08 day period orbit misaligned at -115.9 +/- 4.1 deg to the rotation axis of the star. The planet is confirmed via both the detection of the radial velocity orbit, and the Doppler tomographic detection of the shadow of the planet over two transits. The nature of the spin-orbit misaligned transit geometry allows us to place a constraint on the level of differential rotation in the host star; we find that KELT-17 is consistent with both rigid-body rotation and solar differential rotation rates (alpha < 0.30 at 2 sigma significance). KELT-17 is only the fourth A-star with a confirmed transiting planet, and with a mass of 1.635 -0.061/+0.066 Msun, effective temperature of 7454 +/- 49 K, and projected rotational velocity v sin I_* = 44.2 -1.3/+1.5 km/s; it is amongst the most massive, hottest, and most rapidly rotating of known planet hosts.

preprint2015arXiv

X-ray, UV and optical analysis of supergiants: $ε$ Ori

We present a multi-wavelength (X-ray to optical) analysis, based on non-local thermodynamic equilibrium photospheric+wind models, of the B0 Ia-supergiant: $ε$~Ori. The aim is to test the consistency of physical parameters, such as the mass-loss rate and CNO abundances, derived from different spectral bands. The derived mass-loss rate is $\dot{M}/\sqrt{f_\infty}\sim$1.6$\times$10$^{-6}$ M$_\odot$ yr$^{-1}$ where $f_\infty$ is the volume filling factor. However, the S IV $λλ$1062,1073 profiles are too strong in the models; to fit the observed profiles it is necessary to use $f_\infty<$0.01. This value is a factor of 5 to 10 lower than inferred from other diagnostics, and implies $\dot{M} \lesssim1 \times 10^{-7}$ M$_\odot$ yr$^{-1}$. The discrepancy could be related to porosity-vorosity effects or a problem with the ionization of sulfur in the wind. To fit the UV profiles of N V and O VI it was necessary to include emission from an interclump medium with a density contrast ($ρ_{cl}/ρ_{ICM}$) of $\sim$100. X-ray emission in H-He like and Fe L lines was modeled using four plasma components located within the wind. We derive plasma temperatures from $1 \times 10^{6}$ to $7\times 10^{6}$ K, with lower temperatures starting in the outer regions (R$_0\sim$3-6 R$_*$), and a hot component starting closer to the star (R$_0\lesssim$2.9 R$_*$). From X-ray line profiles we infer $\dot{M} <\, 4.9\times10^{-7}$ M$_\odot$ yr$^{-1}$. The X-ray spectrum ($\geq$0.1 kev) yields an X-ray luminosity $L_{\rm X}\sim 2.0\times10^{-7} L_{\rm bol}$, consistent with the superion line profiles. X-ray abundances are in agreement with those derived from the UV and optical analysis: $ε$ Ori is slightly enhanced in nitrogen and depleted in carbon and oxygen, evidence for CNO processed material.

preprint2014arXiv

Measuring mass-loss rates and constraining shock physics using X-ray line profiles of O stars from the Chandra archive

We quantitatively investigate the extent of wind absorption signatures in the X-ray grating spectra of all non-magnetic, effectively single O stars in the Chandra archive via line profile fitting. Under the usual assumption of a spherically symmetric wind with embedded shocks, we confirm previous claims that some objects show little or no wind absorption. However, many other objects do show asymmetric and blue shifted line profiles, indicative of wind absorption. For these stars, we are able to derive wind mass-loss rates from the ensemble of line profiles, and find values lower by an average factor of 3 than those predicted by current theoretical models, and consistent with H-alpha if clumping factors of f_cl ~ 20 are assumed. The same profile fitting indicates an onset radius of X-rays typically at r ~ 1.5 R_star, and terminal velocities for the X-ray emitting wind component that are consistent with that of the bulk wind. We explore the likelihood that the stars in the sample that do not show significant wind absorption signatures in their line profiles have at least some X-ray emission that arises from colliding wind shocks with a close binary companion. The one clear exception is zeta Oph, a weak-wind star that appears to simply have a very low mass-loss rate. We also reanalyse the results from the canonical O supergiant zeta Pup, using a solar-metallicity wind opacity model and find Mdot = 1.8 \times 10^{-6} M_sun/yr, consistent with recent multi-wavelength determinations.

preprint2014arXiv

Measuring the shock-heating rate in the winds of O stars using X-ray line spectra

We present a new method for using measured X-ray emission line fluxes from O stars to determine the shock-heating rate due to instabilities in their radiation-driven winds. The high densities of these winds means that their embedded shocks quickly cool by local radiative emission, while cooling by expansion should be negligible. Ignoring for simplicity any non-radiative mixing or conductive cooling, the method presented here exploits the idea that the cooling post-shock plasma systematically passes through the temperature characteristic of distinct emission lines in the X-ray spectrum. In this way, the observed flux distribution among these X-ray lines can be used to construct the cumulative probability distribution of shock strengths that a typical wind parcel encounters as it advects through the wind. We apply this new method (Gayley 2014) to Chandra grating spectra from five O stars with X-ray emission indicative of embedded wind shocks in effectively single massive stars. Correcting for wind absorption of the X-ray line emission is a crucial component of our analysis, and we use wind optical depth values derived from X-ray line-profile fitting (Cohen et al. 2014) in order to make that correction. The shock-heating rate results we derive for all the stars are quite similar: the average wind mass element passes through roughly one shock that heats it to at least $10^6$ K as it advects through the wind, and the cumulative distribution of shock strengths is a strongly decreasing function of temperature, consistent with a negative power-law of index $n \approx 3$, implying a marginal distribution of shock strengths that scales as $T^{-4}$, and with hints of an even steeper decline or cut-off above $10^7$ K.

preprint2013arXiv

Constraints on porosity and mass loss in O-star winds from modeling of X-ray emission line profile shapes

We fit X-ray emission line profiles in high resolution XMM-Newton and Chandra grating spectra of the early O supergiant Zeta Pup with models that include the effects of porosity in the stellar wind. We explore the effects of porosity due to both spherical and flattened clumps. We find that porosity models with flattened clumps oriented parallel to the photosphere provide poor fits to observed line shapes. However, porosity models with isotropic clumps can provide acceptable fits to observed line shapes, but only if the porosity effect is moderate. We quantify the degeneracy between porosity effects from isotropic clumps and the mass-loss rate inferred from the X-ray line shapes, and we show that only modest increases in the mass-loss rate (<~ 40%) are allowed if moderate porosity effects (h_infinity <~ R_*) are assumed to be important. Large porosity lengths, and thus strong porosity effects, are ruled out regardless of assumptions about clump shape. Thus, X-ray mass-loss rate estimates are relatively insensitive to both optically thin and optically thick clumping. This supports the use of X-ray spectroscopy as a mass-loss rate calibration for bright, nearby O stars.

preprint2012arXiv

An X-ray Survey of Colliding Wind Binaries

We have compiled a list of 36 O+O and 89 Wolf-Rayet binary candidates in the Milky Way and Magellanic clouds detected with the Chandra, XMM-Newton and ROSAT satellites to probe the connection between their X-ray properties and their system characteristics. Of the WR binaries with published parameters, all but two have kT > 0.9 keV. The most X-ray luminous WR binaries are typically very long period systems. The WR binaries show a nearly four-order of magnitude spread in X-ray luminosity, even among among systems with very similar WR primaries. Among the O+O binaries, short-period systems generally have soft X-ray spectra and longer period systems show harder X-ray spectra, again with a large spread in LX/Lbol.

preprint2012arXiv

Thin-shell mixing in radiative wind-shocks and the Lx-Lbol scaling of O-star X-rays

X-ray satellites since Einstein have empirically established that the X-ray luminosity from single O-stars scales linearly with bolometric luminosity, Lx ~ 10^{-7} Lbol. But straightforward forms of the most favored model, in which X-rays arise from instability-generated shocks embedded in the stellar wind, predict a steeper scaling, either with mass loss rate Lx ~ Mdot ~ Lbol^{1.7} if the shocks are radiative, or with Lx ~ Mdot^{2} ~ Lbol^{3.4} if they are adiabatic. This paper presents a generalized formalism that bridges these radiative vs. adiabatic limits in terms of the ratio of the shock cooling length to the local radius. Noting that the thin-shell instability of radiative shocks should lead to extensive mixing of hot and cool material, we propose that the associated softening and weakening of the X-ray emission can be parametrized as scaling with the cooling length ratio raised to a power m$, the "mixing exponent". For physically reasonable values m ~= 0.4, this leads to an X-ray luminosity Lx ~ Mdot^{0.6} ~ Lbol that matches the empirical scaling. To fit observed X-ray line profiles, we find such radiative-shock-mixing models require the number of shocks to drop sharply above the initial shock onset radius. This in turn implies that the X-ray luminosity should saturate and even decrease for optically thick winds with very high mass-loss rates. In the opposite limit of adiabatic shocks in low-density winds (e.g., from B-stars), the X-ray luminosity should drop steeply with Mdot^2. Future numerical simulation studies will be needed to test the general thin-shell mixing ansatz for X-ray emission.

preprint2011arXiv

A generalised porosity formalism for isotropic and anisotropic effective opacity and its effects on X-ray line attenuation in clumped O star winds

We present a generalised formalism for treating the porosity-associated reduction in continuum opacity that occurs when individual clumps in a stochastic medium become optically thick. We consider geometries resulting in either isotropic or anisotropic effective opacity, and, in addition to an idealised model in which all clumps have the same local overdensity and scale, we also treat an ensemble of clumps with optical depths set by Markovian statistics. This formalism is then applied to the specific case of bound-free absorption of X-rays in hot star winds, a process not directly affected by clumping in the optically thin limit. We find that the Markov model gives surprisingly similar results to those found previously for the single clump model, suggesting that porous opacity is not very sensitive to details of the assumed clump distribution function. Further, an anisotropic effective opacity favours escape of X-rays emitted in the tangential direction (the `venetian blind' effect), resulting in a 'bump' of higher flux close to line centre as compared to profiles computed from isotropic porosity models. We demonstrate how this characteristic line shape may be used to diagnose the clump geometry, and we confirm previous results that for optically thick clumping to significantly influence X-ray line profiles, very large porosity lengths, defined as the mean free path between clumps, are required. Moreover, we present the first X-ray line profiles computed directly from line-driven instability simulations using a 3-D patch method, and find that porosity effects from such models also are very small. This further supports the view that porosity has, at most, a marginal effect on X-ray line diagnostics in O stars, and therefore that these diagnostics do indeed provide a good `clumping insensitive' method for deriving O star mass-loss rates.

preprint2011arXiv

Chandra X-ray spectroscopy of the very early O supergiant HD 93129A: constraints on wind shocks and the mass-loss rate

We present analysis of both the resolved X-ray emission line profiles and the broadband X-ray spectrum of the O2 If* star HD 93129A, measured with the Chandra HETGS. This star is among the earliest and most massive stars in the Galaxy, and provides a test of the embedded wind shock scenario in a very dense and powerful wind. A major new result is that continuum absorption by the dense wind is the primary cause of the hardness of the observed X-ray spectrum, while intrinsically hard emission from colliding wind shocks contributes less than 10% of the X-ray flux. We find results consistent with the predictions of numerical simulations of the line-driving instability, including line broadening indicating an onset radius of X-ray emission of several tenths Rstar. Helium-like forbidden-to-intercombination line ratios are consistent with this onset radius, and inconsistent with being formed in a wind-collision interface with the star's closest visual companion at a distance of ~100 AU. The broadband X-ray spectrum is fit with a dominant emission temperature of just kT = 0.6 keV along with significant wind absorption. The broadband wind absorption and the line profiles provide two independent measurements of the wind mass-loss rate: Mdot = 5.2_{-1.5}^{+1.8} \times 10^{-6} Msun/yr and Mdot = 6.8_{-2.2}^{+2.8} \times 10^{-6} Msun/yr, respectively. This is the first consistent modeling of the X-ray line profile shapes and broadband X-ray spectral energy distribution in a massive star, and represents a reduction of a factor of 3 to 4 compared to the standard H-alpha mass-loss rate that assumes a smooth wind.

preprint2010arXiv

A mass-loss rate determination for zeta Puppis from the quantitative analysis of X-ray emission line profiles

We fit every emission line in the high-resolution Chandra grating spectrum of zeta Pup with an empirical line profile model that accounts for the effects of Doppler broadening and attenuation by the bulk wind. For each of sixteen lines or line complexes that can be reliably measured, we determine a best-fitting fiducial optical depth, tau_* = kappa*Mdot/4{pi}R_{\ast}v_{\infty}, and place confidence limits on this parameter. These sixteen lines include seven that have not previously been reported on in the literature. The extended wavelength range of these lines allows us to infer, for the first time, a clear increase in tau_* with line wavelength, as expected from the wavelength increase of bound-free absorption opacity. The small overall values of tau_*, reflected in the rather modest asymmetry in the line profiles, can moreover all be fit simultaneously by simply assuming a moderate mass-loss rate of 3.5 \pm 0.3 \times 10^{-6} Msun/yr, without any need to invoke porosity effects in the wind. The quoted uncertainty is statistical, but the largest source of uncertainty in the derived mass-loss rate is due to the uncertainty in the elemental abundances of zeta Pup, which affects the continuum opacity of the wind, and which we estimate to be a factor of two. Even so, the mass-loss rate we find is significantly below the most recent smooth-wind H-alpha mass-loss rate determinations for zeta Pup, but is in line with newer determinations that account for small-scale wind clumping. If zeta Pup is representative of other massive stars, these results will have important implications for stellar and galactic evolution.

preprint2010arXiv

X-ray spectral diagnostics of activity in massive stars

X-rays give direct evidence of instabilities, time-variable structure, and shock heating in the winds of O stars. The observed broad X-ray emission lines provide information about the kinematics of shock-heated wind plasma, enabling us to test wind-shock models. And their shapes provide information about wind absorption, and thus about the wind mass-loss rates. Mass-loss rates determined from X-ray line profiles are not sensitive to density-squared clumping effects, and indicate mass-loss rate reductions of factors of 3 to 6 over traditional diagnostics that suffer from density-squared effects. Broad-band X-ray spectral energy distributions also provide mass-loss rate information via soft X-ray absorption signatures. In some cases, the degree of wind absorption is so high that the hardening of the X-ray SED can be quite significant. We discuss these results as applied to the early O stars zeta Pup (O4 If), 9 Sgr (O4 V((f))), and HD 93129A (O2 If*).

preprint2009arXiv

X-ray Spectroscopy of the Radiation-Driven Winds of Massive Stars: Line Profile and Line Ratio Diagnostics

Massive stars drive powerful, supersonic winds via the radiative momentum associated with the thermal UV emission from their photospheres. Shock phenomena are ubiquitous in these winds, heating them to millions, and sometimes tens of millions, of degrees. The emission line spectra from the shock-heated plasma provide powerful diagnostics of the winds' physical conditions, which in turn provide constraints on models of wind shock heating. Here I show how x-ray line transfer is affected by photoelectric absorption in the partially ionized component of the wind and how it can be modeled to determine the astrophysically important mass-loss rates of these stellar winds. I also discuss how photoexcitation out of metastable excited levels of helium-like ions can provide critical information about the location of the hot plasma in magnetically channeled massive star winds.

preprint2006arXiv

Measurements and analysis of helium-like triplet ratios in the X-ray spectra of O-type stars

We discuss new methods of measuring and interpreting the forbidden-to-intercombination line ratios of helium-like triplets in the X-ray spectra of O-type stars, including accounting for the spatial distribution of the X-ray emitting plasma and using the detailed photospheric UV spectrum. Measurements are made for four O stars using archival Chandra HETGS data. We assume an X-ray emitting plasma spatially distributed in the wind above some minimum radius R_0. We find minimum radii of formation typically in the range of 1.25 < R_0 / R_* < 1.67, which is consistent with results obtained independently from line profile fits. We find no evidence for anomalously low f/i ratios and we do not require the existence of X-ray emitting plasmas at radii that are too small to generate sufficiently strong shocks.

preprint2005arXiv

Chandra HETGS Multi-Phase Spectroscopy of the Young Magnetic O Star theta^1 Orionis C

We report on four Chandra grating observations of the oblique magnetic rotator theta^1 Ori C (O5.5 V) covering a wide range of viewing angles with respect to the star's 1060 G dipole magnetic field. We employ line-width and centroid analyses to study the dynamics of the X-ray emitting plasma in the circumstellar environment, as well as line-ratio diagnostics to constrain the spatial location, and global spectral modeling to constrain the temperature distribution and abundances of the very hot plasma. We investigate these diagnostics as a function of viewing angle and analyze them in conjunction with new MHD simulations of the magnetically channeled wind shock mechanism on theta^1 Ori C. This model fits all the data surprisingly well, predicting the temperature, luminosity, and occultation of the X-ray emitting plasma with rotation phase.

preprint2003arXiv

High-Resolution {\it Chandra} Spectroscopy of tau Scorpii: A Narrow-Line X-ray Spectrum From a Hot Star

Long known to be an unusual early-type star by virtue of its hard and strong X-ray emission, tau Scorpii poses a severe challenge to the standard picture of O-star wind-shock X-ray emission. The Chandra HETGS spectrum now provides significant direct evidence that this B0.2 star does not fit this standard wind-shock framework. The many emission lines detected with the Chandra gratings are significantly narrower than what would be expected from a star with the known wind properties of tau Sco, although they are broader than the corresponding lines seen in late-type coronal sources. While line ratios are consistent with the hot plasma on this star being within a few stellar radii of the photosphere, from at least one He-like complex there is evidence that the X-ray-emitting plasma is located more than a stellar radius above the photosphere. The Chandra spectrum of tau Sco is harder and more variable than those of other hot stars, with the exception of the young magnetized O star theta Ori C. We discuss these new results in the context of wind, coronal, and hybrid wind-magnetic models of hot-star X-ray emission.

preprint2003arXiv

X-ray Emission Line Profile Modeling of O stars: Fitting a Spherically-Symmetric Analytic Wind-Shock Model to the Chandra Spectrum of zeta Puppis

X-ray emission line profiles provide the most direct insight into the dynamics and spatial distribution of the hot, X-ray-emitting plasma above the surfaces of OB stars. The O supergiant zeta Puppis shows broad, blueshifted, and asymmetric line profiles, generally consistent with the wind-shock picture. We model the profiles of eight lines in the Chandra HETGS spectrum of this prototypical hot star. The fitted lines indicate that the plasma is distributed throughout the wind starting close to the photosphere, that there is significantly less attenuation of the X-rays by the overlying wind than is generally supposed, and that there is not a strong trend in wind absorption with wavelength.

preprint2002arXiv

X ray emission line profile modeling of hot stars

The launch of high-spectral-resolution x-ray telescopes (Chandra, XMM) has provided a host of new spectral line diagnostics for the astrophysics community. In this paper we discuss Doppler-broadened emission line profiles from highly supersonic outflows of massive stars. These outflows, or winds, are driven by radiation pressure and carry a tremendous amount of kinetic energy, which can be converted to x rays by shock-heating even a small fraction of the wind plasma. The unshocked, cold wind is a source of continuum opacity to the x rays generated in the shock-heated portion of the wind. Thus the emergent line profiles are affected by transport through a two-component, moving, optically thick medium. While complicated, the interactions among these physical effects can provide quantitative information about the spatial distribution and velocity of the x-ray-emitting and absorbing plasma in stellar winds. We present quantitative models of both a spherically-symmetric wind and a wind with hot plasma confined in an equatorial disk by a dipole magnetic field.