Search

We investigated the temporal and spectral features of $\gamma$ Cassiopeiae's X-ray emission within the context of the white dwarf accretion hypothesis. We find that the variabilities present in the X-ray data show two different signals, one primarily due to absorption and the other due to flickering like in non-magnetic cataclysmic variables. We then use this two-component insight to investigate previously un-reported simultaneous XMM and NuSTAR data. The model fitting results find white dwarf properties consistent with optical studies alongside a significant secondary, thermal source. We propose a secondary shock between the Be decretion disk and white dwarf accretion disk as the source. Finally, we analyzed a unique, low-count rate event of the XMM light curve as potential evidence for the white dwarf encountering Be decretion disk structures., Comment: Accepted for publication in ApJ

The Polstar small explorer concept is for an ultraviolet (UV) spectropolarimetry space telescope mission with a focus on massive star astrophysics. The instrument waveband will be from 115 nm - 286 nm for spectroscopy and 122 nm - 286 nm for polarimetry. All 4 Stokes parameters, IQUV, will be measured at a resolving power of R=20,000 (15 km/s velocity resolution). The telescope aperture will be 40 cm with an effective area of about 22 cm^2 at a reference wavelength of 150 nm. The thrust of the science goals will be to determine the astrophysics of angular momentum exchange and transport, and consequences for massive star properties and evolution. This includes the effects of rapid to critical rotation for individual stars (magnetic and non-magnetic), and the effects of mass transfer for massive binaries, including identification of stripped core stars. If selected by the NASA/SMEX program, Polstar would launch around 2031 and observe ~300 stars to achieve science goals. The mission will include a Guest Observer program to advance discovery in other areas of astrophysics., Comment: accepted to the journal, Bulletin de la Societe Royale des Sciences

High-resolution X-ray spectra of $\pi\,$Aqr, a $\gamma\,$Cas-type star, obtained with the Chandra/HETG grating spectrometer, revealed emission lines of H-like ions of Mg, Si, S, and Fe, a strong, hard continuum, and a lack of He-like ions, indicating the presence of very hot thermal plasma. The X-ray light curve showed significant fluctuations, with coherent variability at period of about 3400 seconds in one observation. The hardness ratio was relatively constant except for one observation in which the spectrum was much harder and more absorbed. We interpret the X-ray emission as arising from accretion onto the secondary, which is likely a magnetic white dwarf, an intermediate polar system., Comment: 7 pages, 3 figures, 3 tables

The bright WN4 star EZ CMa exhibits a 3.77 day periodicity in photometry, spectroscopy, and polarimetry but the variations in the measurements are not strictly phase-locked, exhibiting changes in reference times, amplitudes, and the shape of the variability happening over times as short as a few weeks. Recently, 137 days of contiguous, variable photometry from BRITE-Constellation was interpreted as caused either by large-scale dense wind structures modulated by rotation, or by a fast-precessing binary having a slightly shorter 3.626 day orbital period and a fast apsidal motion rate of $1315^\circ\,\text{yr}^{-1}$. We aim at testing the latter hypothesis through analysis of spectroscopy and focus on the N\,{\sc v} $\lambda\,4945$ line. We derive an orbital solution for the system and reject the 3.626 day period to represent the variations in the radial velocities of EZ CMa. An orbital solution with an orbital period of 3.77 days was obtained but at the cost of an extremely high and thus improbable apsidal motion rate. Our best orbital solution yields a period of $3.751\pm0.001$\,days with no apsidal motion. We place our results in the context of other variability studies and system properties. While we cannot fully reject the precessing binary model, we find that the corotating interaction region (CIR) hypothesis is better supported by these and other data through qualitative models of CIRs., Comment: accepted to MNRAS

Identifying trends between observational data and the range of physical parameters of massive stars is a critical step to the still-elusive full understanding of the source, structure, and evolution of X-ray emission from the stellar winds, requiring a substantial sample size and systematic analysis methods. The \emph{Chandra} data archive as of 2022 contains 37 high resolution spectra of O, B, and WR stars, observed with the \emph{Chandra}/HETGS and of sufficient quality to fit the continua and emission line profiles. Using a systematic approach to the data analysis, we explore morphological trends in the line profiles (i.e., O, Ne, Mg, Si) and find that the centroid offsets of resolved lines versus wavelength can be separated in three empirically-defined groups based on the amount of line broadening and centroid offset. Using \ion{Fe}{17} (15.01 \AA, 17.05 \AA) and \ion{Ne}{10} $\alpha$ (12.13 \AA) lines which are prevalent among the sample stars, we find a well-correlated linear trend of increasing Full Width Half Maximum (FWHM) with faster wind terminal velocity. The H-like/He-like total line flux ratio for strong lines displays different trends with spectral class depending on ion species. Some of the sources in our sample have peculiar properties (e.g., magnetic and $\gamma$ Cas-analogue stars) and we find that these sources stand out as outliers from more regular trends. Finally, our spectral analysis is presented summarily in terms of X-ray spectral energy distributions in specific luminosity for each source, plus tables of line identifications and fluxes., Comment: 78 pages, 46 figures, 41 tables, Accepted for publication in ApJ

The 2021 outburst of the symbiotic recurrent nova RS Oph was monitored with the Neutron Star Interior Composition Explorer Mission (NICER) in the 0.2-12 keV range from day one after the optical maximum, until day 88, producing an unprecedented, detailed view of the outburst development. The X-ray flux preceding the supersoft X-ray phase peaked almost 5 days after optical maximum and originated only in shocked ejecta for 21 to 25 days. The emission was thermal; in the first 5 days only a non-collisional-ionization equilibrium model fits the spectrum, and a transition to equilibrium occurred between days 6 and 12. The ratio of peak X-rays flux measured in the NICER range to that measured with Fermi in the 60 MeV-500 GeV range was about 0.1, and the ratio to the peak flux measured with H.E.S.S. in the 250 GeV-2.5 TeV range was about 100. The central supersoft X-ray source (SSS), namely the shell hydrogen burning white dwarf (WD), became visible in the fourth week, initially with short flares. A huge increase in flux occurred on day 41, but the SSS flux remained variable. A quasi-periodic oscillation every ~35 s was always observed during the SSS phase, with variations in amplitude and a period drift that appeared to decrease in the end. The SSS has characteristics of a WD of mass >1 M(solar). Thermonuclear burning switched off shortly after day 75, earlier than in 2006 outburst. We discuss implications for the nova physics., Comment: Accepted for publication in the Astrophysical Journal

The semi-regular variable star V Canum Venaticorum (V CVn) is well-known for its unusual linear polarization position angle (PA). Decades of observing V CVn reveal a nearly constant PA spanning hundreds of pulsation cycles. This phenomenon has persisted through variability that has ranged by 2 magnitudes in optical brightness and through variability in the polarization amplitude over 0.3% and 6.9%. Additionally, the polarization fraction of V CVn varies inversely with brightness. This paper presents polarization measurements obtained over three pulsation cycles. We find that the polarization maximum does not always occur precisely at the same time as the brightness minimum. Instead, we observe a small lead or lag in relation to the brightness minimum, spanning a period of a few days up to three weeks. Furthermore, the PA sometimes exhibits a non-negligible rotation, especially at lower polarization levels. To elucidate the unusual optical behavior of V CVn, we present a list of literature sources that also exhibit polarization variability with a roughly fixed PA. We find this correlation occurs in stars with high tangential space velocities, i.e., "runaway" stars, suggesting that the long-term constant PA is related to how the circumstellar gas is shaped by the star's high-speed motion through the interstellar medium., Comment: 9 pages + appendices, 7 figures, accepted for publication in A&A

We report precision ground-based broadband optical intensity and linear-polarisation light-curves for the sky's brightest WN8 star, WR40. WN8 stars are notorious for their high level of variability, stemming from stochastic clumps in their strong winds that are slower and less hot than the winds of most other Wolf-Rayet (WR) stars. We confirm previous results that many WR stars display an amplitude of variability that is an order-of-magnitude higher in photometry than in polarimetry. For the first time, the unique nature of near simultaneity of our photometric and polarimetric observations of WR40 allows us to check whether the two types of variability show correlated behaviour, of which we find none. Assuming simple temporal functions for the brightness and polarisation of individual clumps, a model for simulated light curves is found to reproduce the properties of the observations, specifically the absence of correlation between photometric and polarimetric variations, the ratio of standard deviations for photometric and polarimetric variability, and the ratio of the average intrinsic polarisation relative to its standard deviation. Mapping the solution for time variability to a spatial coordinate suggests that the wind clumps of WR40 are free-free emitting in addition to light scattering., Comment: to appear in MNRAS

The winds of massive stars are important for their direct impact on the interstellar medium, and for their influence on the final state of a star prior to it exploding as a supernova. However, the dynamics of these winds is understood primarily via their illumination from a single central source. The Doppler shift seen in resonance lines is a useful tool for inferring these dynamics, but the mapping from that Doppler shift to the radial distance from the source is ambiguous. Binary systems can reduce this ambiguity by providing a second light source at a known radius in the wind, seen from orbitally modulated directions. From the nature of the collision between the winds, a massive companion also provides unique additional information about wind momentum fluxes. Since massive stars are strong ultraviolet (UV) sources, and UV resonance line opacity in the wind is strong, UV instruments with a high resolution spectroscopic capability are essential for extracting this dynamical information. Polarimetric capability also helps to further resolve ambiguities in aspects of the wind geometry that are not axisymmetric about the line of sight, because of its unique access to scattering direction information. We review how the proposed MIDEX-scale mission Polstar can use UV spectropolarimetric observations to critically constrain the physics of colliding winds, and hence radiatively-driven winds in general. We propose a sample of 20 binary targets, capitalizing on this unique combination of illumination by companion starlight, and collision with a companion wind, to probe wind attributes over a range in wind strengths. Of particular interest is the hypothesis that the radial distribution of the wind acceleration is altered significantly, when the radiative transfer within the winds becomes optically thick to resonance scattering in multiple overlapping UV lines., Comment: 26 pages, 12 figures, Review in a topical collection series of Astrophysics and Space Sciences on the proposed Polstar satellite. arXiv admin note: substantial text overlap with arXiv:2111.11552

Massive star binaries are critical laboratories for measuring masses and stellar wind mass-loss rates. A major challenge is inferring viewing inclination and extracting information about the colliding wind interaction (CWI) region. Polarimetric variability from electron scattering in the highly ionized winds provides important diagnostic information about system geometry. We combine for the first time the well-known generalized treatment of \citet{brown_polarisation_1978} for variable polarization from binaries with the semi-analytic solution for the geometry and surface density CWI shock interface between the winds based on Canto et al 1996. Our calculations include some simplifications in the form of inverse square-law wind densities and the assumption of axisymmetry, but in so doing arrive at several robust conclusions. One is that when the winds are nearly equal (e.g., O\,+\,O binaries), the polarization has a relatively mild decline with binary separation. Another is that despite Thomson scattering being a gray opacity, the continuum polarization can show chromatic effects at ultraviolet wavelengths but will be mostly constant at longer wavelengths. Finally, when one wind dominates the other, as for example in WR+OB binaries, the polarization is expected to be larger at wavelengths where the OB component is more luminous, and generally smaller at wavelengths where the WR component is more luminous. This behavior arises because from the perspective of the WR star, the distortion of the scattering envelope from spherical is a minor perturbation situated far from the WR star. By contrast, the polarization contribution from the OB star is dominated by the geometry of the CWI shock., Comment: accepted for publication in the Astrophysical Journal. arXiv admin note: text overlap with arXiv:2111.11552

The mass-loss rates from single massive stars are high enough to form radio photospheres at large distances from the stellar surface where the wind is optically thick to (thermal) free-free opacity. Here we calculate the far-infrared, millimeter, and radio band spectral energy distributions (SEDs) that can result from the combination of free-free processes and synchrotron emission, to explore the conditions for non-thermal SEDs. Simplifying assumptions are adopted in terms of scaling relations for the magnetic field strength and the spatial distribution of relativistic electrons. The wind is assumed to be spherically symmetric, and we consider the effect of Razin suppression on the synchrotron emission. Under these conditions, long-wavelength SEDs with synchrotron emission can be either more steep or more shallow than the canonical asymptotic power-law SED from a non-magnetic wind. When non-thermal emission is present, the resultant SED shape is generally not a power-law; however, the variation in slope can change slowly with wavelength. Consequently, over a limited range of wavelengths, the SED can masquerade as approximately a power law. While most observed non-thermal long-wavelength spectra are associated with binarity, synchroton emission can have only a mild influence on single-star SEDs, requiring finer levels of wavelength sampling for detection of the effect., Comment: 15 pages, 3 figures, accepted for publication in The Astrophysical Journal

Understanding the evolution of massive binary stars requires accurate estimates of their masses. This understanding is critically important because massive star evolution can potentially lead to gravitational wave sources such as binary black holes or neutron stars. For Wolf-Rayet stars with optically thick stellar winds, their masses can only be determined with accurate inclination angle estimates from binary systems which have spectroscopic $M \sin i$ measurements. Orbitally-phased polarization signals can encode the inclination angle of binary systems, where the Wolf-Rayet winds act as scattering regions. We investigated four Wolf-Rayet + O star binary systems, WR 42, WR 79, WR 127, and WR 153, with publicly available phased polarization data to estimate their masses. To avoid the biases present in analytic models of polarization while retaining computational expediency, we used a Monte Carlo radiative transfer model accurately emulated by a neural network. We used the emulated model to investigate the posterior distribution of parameters of our four systems. Our mass estimates calculated from the estimated inclination angles put strong constraints on existing mass estimates for three of the systems, and disagrees with the existing mass estimates for WR 153. We recommend a concerted effort to obtain polarization observations that can be used to estimate the masses of Wolf-Rayet binary systems and increase our understanding of their evolutionary paths., Comment: 24 pages, 13 figures, accepted to ApJ

One objective of the Polstar spectropolarimetry mission is to characterize the degree of nonconservative mass transfer that occurs at various stages of binary evolution, from the initial mass reversal to the late Algol phase. The proposed instrument combines spectroscopic and polarimetric capabilities, where the spectroscopy can resolve Doppler shifts in UV resonance lines with 10 km/s precision, and polarimetry can resolve linear polarization with 1e-3 precision or better. The spectroscopy will identify absorption by mass streams seen in projection against the stellar disk as a function of orbital phase, hot accretion spots, as well as scattering from extended splash structures, circumbinary disks, and other flows in and above/below the orbital plane (e.g. jets) that fail to be transferred conservatively. The polarimetry affects more the light coming from material not seen against the stellar disk, allowing the geometry of the scattering to be tracked, resolving ambiguities left by the spectroscopy and light-curve information. For example, nonconservative mass streams ejected in the polar direction will produce polarization of the opposite sign from conservative transfer accreting in the orbital plane. Also, time domain coverage over a range of phases of the binary orbit are well supported by the Polstar observing strategy. Combining these elements will significantly improve our understanding of the mass transfer process and the amount of mass that can escape from the system, an important channel for changing the final mass, and ultimate supernova, of the large number of massive stars found in binaries at close enough separation to undergo interaction., Comment: Refereed paper in Ap&SS 367:Topical Collection, 23 pages, 10 figures

The most massive stars are thought to lose a significant fraction of their mass in a steady wind during the main-sequence and blue supergiant phases. This in turn sets the stage for their further evolution and eventual supernova, with consequences for ISM energization and chemical enrichment. Understanding these processes requires accurate observational constraints on the mass-loss rates of the most luminous stars, which can also be used to test theories of stellar wind generation. In the past, mass-loss rates have been characterized via collisional emission processes such as H$\alpha$ and free-free radio emission, but these so-called "density squared" diagnostics require correction in the presence of widespread clumping. Recent observational and theoretical evidence points to the likelihood of a ubiquitously high level of such clumping in hot-star winds, but quantifying its effects requires a deeper understanding of the complex dynamics of radiatively driven winds. Furthermore, large-scale structures arising from surface anisotropies and propagating throughout the wind can further complicate the picture by introducing further density enhancements, affecting mass-loss diagnostics. Time series spectroscopy of UV resonance lines with high resolution and high signal-to-noise are required to better understand this complex dynamics, and help correct "density squared" diagnostics of mass-loss rates. The proposed Polstar mission easily provides the necessary resolution at the sound-speed scale of 20 km s$^{-1}$, on three dozen bright targets with signal-to noise an order of magnitude above that of the celebrated IUE MEGA campaign, via continuous observations that track structures advecting through the wind in real time., Comment: white paper

As sources of chemical enrichment, ionizing radiation and energetic feedback, massive stars drive the ecology of their host galaxies despite their relative rarity, additionally to yielding compact remnants, which can generate gravitational waves. The evolution of massive stars is crucially informed by their detailed mass-loss history; however, wind structures on a variety of scales cause important uncertainties on their mass-loss rates. Binary systems can place further constraints on the mass-loss properties of massive stars, especially colliding-wind binaries. In this paper, we review how the proposed MIDEX-scale mission Polstar can critically constrain the physics of colliding winds (and hence radiatively-driven winds in general) with ultraviolet spectropolarimetric observations, providing an unprecedented improvement on the accuracy of the determination of both mass-loss rates and the velocity structure of the winds of massive stars. We propose a sample of 17 targets that will allow us to study a variety of wind-colliding systems spanning a large parameter space using the spatial information yielded by both spectroscopic and polarimetric data obtained with Polstar., Comment: 21 pages, 12 figures 1 table, Whitepaper

Polstar is a proposed NASA MIDEX mission that would feature a high resolution UV spectropolarimeter capable of measure all four Stokes parameters onboard a 60cm telescope. The mission would pioneer the field of time-domain UV spectropolarimetry. Time domain UV spectropolarimetry offers the best resource to determine the geometry and physical conditions of protoplanetary disks from the stellar surface to <5 AU. We detail two key objectives that a dedicated time domain UV spectropolarimetry survey, such as that enabled by Polstar, could achieve: 1) Test the hypothesis that magneto-accretion operating in young planet-forming disks around lower-mass stars transitions to boundary layer accretion in planet-forming disks around higher mass stars; and 2) Discriminate whether transient events in the innermost regions of planet-forming disks of intermediate mass stars are caused by inner disk mis-alignments or from stellar or disk emissions., Comment: 16 pages, 7 figures

The Polstar mission will provide for a space-borne 60cm telescope operating at UV wavelengths with spectropolarimetric capability capturing all four Stokes parameters (intensity, two linear polarization components, and circular polarization). Polstar's capabilities are designed to meet its goal of determining how circumstellar gas flows alter massive stars' evolution, and finding the consequences for the stellar remnant population and the stirring and enrichment of the interstellar medium, by addressing four key science objectives. In addition, Polstar will determine drivers for the alignment of the smallest interstellar grains, and probe the dust, magnetic fields, and environments in the hot diffuse interstellar medium, including for the first time a direct measurement of the polarized and energized properties of intergalactic dust. Polstar will also characterize processes that lead to the assembly of exoplanetary systems and that affect exoplanetary atmospheres and habitability. Science driven design requirements include: access to ultraviolet bands: where hot massive stars are brightest and circumstellar opacity is highest; high spectral resolution: accessing diagnostics of circumstellar gas flows and stellar composition in the far-UV at 122-200nm, including the NV, SiIV, and CIV resonance doublets and other transitions such as NIV, AlIII, HeII, and CIII; polarimetry: accessing diagnostics of circumstellar magnetic field shape and strength when combined with high FUV spectral resolution and diagnostics of stellar rotation and distribution of circumstellar gas when combined with low near-UV spectral resolution; sufficient signal-to-noise ratios: ~1000 for spectropolarimetric precisions of 0.1% per exposure; ~100 for detailed spectroscopic studies; ~10 for exploring dimmer sources; and cadence: ranging from 1-10 minutes for most wind variability studies., Comment: 13 pages, 6 figures, presented at SPIE Optics and Photonics 2021

One possibility for detecting low-amplitude pulsational variations is through gravitational microlensing. During a microlensing event, the temporary brightness increase leads to improvement in the signal-to-noise ratio, and thereby better detectability of pulsational signatures in light curves. We explore this possibility under two primary considerations. The first is when the standard point-source and point-lens approximation applies. In this scenario, dividing the observed light curve by the best-fitted microlensing model leads to residuals that result in pulsational features with improved uncertainties. The second is for transit events (single lens) or caustic crossing (binary lens). The point-source approximation breaks down, and residuals relative to a simple best-fitted microlensing model display more complex behavior. We employ a Monte-Carlo simulation of microlensing of pulsating variables toward the Galactic bulge for the surveys of OGLE and of KMTNet. We demonstrate that the efficiency for detecting pulsational signatures with intrinsic amplitudes of $<0.25$ mag during single and binary microlensing events, at differences in $\chi^{2}$ of $\Delta \chi^{2} >350$, is $\sim 50-60\%$. The maximum efficiency occurs for pulsational periods $P \simeq 0.1-0.3$ days. We also study the possibility that high-magnification microlensing events of non-radially pulsating stars (NRPs) could be misinterpreted as planetary or binary microlensing events. We conclude that small asymmetric features around lightcurve peaks due to stellar pulsations could be misdiagnosed with crossing (or passing close to) small caustic curves., Comment: 11 pages, 4 Figures, 4 Tables

The eclipsing binary epsilon Aur is unique in being a very long-period binary involving an evolved, variable F star and a suspected B main-sequence star enshrouded in an opaque circumstellar disk. The geometrical arrangement is that the disk is viewed almost perfectly edge-on, with alignment leading to a partial eclipse of the F star. Despite a global observing campaign for the 2009-11 eclipse, there remain outstanding questions about the nature of the binary, its components, the disk, and the evolutionary state of the system. We analyze optical-band polarimetry in conjunction with broad-band color variations to interpret brightness variations across the surface of the F star. We model this both during and after the 1982-84 eclipse for which an extensive and dense data set exists. We develop a model in terms of surface temperature variations characterized by a small global variation overlaid with a temperature variation described with low-order spherical harmonics. While not providing a detailed fit to the dataset, our modeling captures the overall characterization of the color and polarimetric variability. In particular, we are able to recover the gross behavior of the polarimetric excursion in the Q-U plane as observed during eclipse of the F star when compared to post-eclipse behavior., Comment: accepted version for Astronomical Journal

WR 25 is a colliding-wind binary star system comprised of a very massive O2.5If*/WN6 primary and an O-star secondary in a 208-day period eccentric orbit. These hot stars have strong, highly-supersonic winds which interact to form a bright X-ray source from wind-collision-shocks whose conditions change with stellar separation. Different views through the WR and O star winds are afforded with orbital phase as the stars move about their orbits, allowing for exploration of wind structure in ways not easy or even possible for single stars. We have analyzed an on-axis Chandra/HETGS spectrum of WR 25 obtained shortly before periastron when the X-rays emanating from the system are the brightest. From the on-axis observations, we constrain the line fluxes, centroids, and widths of various emission lines, including He-triplets of Si XIII and Mg XI. We have also been able to include several serendipitous off-axis HETG spectra from the archive and study their flux variation with phase. This is the first report on high-resolution spectral studies of WR 25 in X-rays., Comment: 27 pages, 10 figures, submitted to ApJ, comments welcome

Analysis of the recent long exposure Chandra X-ray observation of the early-type O star zeta Pup shows clear variability with a period previously reported in optical photometric studies. These 813 ks of HETG observations taken over a roughly one year time span have two signals of periodic variability: a high significance period of 1.7820 +/- 0.0008 day, and a marginal detection of periodic behavior close to either 5 day or 6 day period. A BRITE-Constellation nanosatellite optical photometric monitoring, using near-contemporaneous observations to the Chandra data, confirms a 1.78060 +/- 0.00088 day period for this star. The optical period coincides with the new Chandra period within their error ranges, demonstrating a link between these two wavebands and providing a powerful lever for probing the photosphere/wind connection in this star. The phase lag of the X-ray maximum relative to the optical maximum is approximately phi=0.45, but consideration of secondary maxima in both datasets indicates possibly two hot spots on the star with an X-ray phase lag of phi=0.1 each. The details of this periodic variation of the X-rays are probed by displaying a phased and trailed X-ray spectrum and by constructing phased light curves for wavelength bands within the HETG spectral coverage, ranging down to bands encompassing groups of emission lines. We propose that the 1.78 day period is the stellar rotation period and explore how stellar bright spots and associated co-rotating interacting regions or CIRs could explain the modulation of the optical and X-ray output for this star and their phase difference., Comment: 17 pages; 10 figures

We study the polarization produced by scattering from dust in a bow shock-shaped region of enhanced density surrounding a stellar source, using the Monte Carlo radiative transfer code SLIP. Bow shocks are structures formed by the interaction of the winds of fast-moving stars with the interstellar medium. Our previous study focused on the polarization produced in these structures by electron scattering; we showed that polarization is highly dependent on inclination angle and that multiple scattering changes the shape and degree of polarization. In contrast to electron scattering, dust scattering is wavelength-dependent, which changes the polarization behaviour. Here we explore different dust particle sizes and compositions and generate polarized spectral energy distributions for each case. We find that the polarization SED behaviour depends on the dust composition and grain size. Including dust emission leads to polarization changes with temperature at higher optical depth in ways that are sensitive to the orientation of the bow shock. In various scenarios and under certain assumptions, our simulations can constrain the optical depth and dust properties of resolved and unresolved bow shock-shaped scattering regions.Constraints on optical depth can provide estimates of local ISM density for observed bow shocks. We also study the impact of dust grains filling the region between the star and bow shock. We see that as the density of dust between the star and bow shock increases, the resulting polarization is suppressed for all the optical depth regimes., Comment: 21 pages, accepted for publication in MNRAS

We study the microlensing of Non-Radially Pulsating (NRP) stars. Pulsations are formulated for stellar radius and temperature using spherical harmonic functions with different values of l,m. The characteristics of the microlensing light curves from NRP stars are investigated in relation to different pulsation modes. For the microlensing of NRP stars, the light curve is not a simple multiplication of the magnification curve and the intrinsic luminosity curve of the source star, unless the effect of finite source size can be ignored. Three main conclusions can be drawn from the simulated light curves. First, for modes with $m\neq0$ and when the viewing inclination is more nearly pole-on, the stellar luminosity towards the observer changes little with pulsation phase. In this case, high-magnification microlensing events are chromatic and can reveal the variability of these source stars. Second, some combinations of pulsation modes produce nearly degenerate luminosity curves (e.g., (l,m)=(3,0), (5,0)). The resulting microlensing light curves are also degenerate, unless the lens crosses the projected source. Finally, for modes involving m=1, the stellar brightness centre does not coincide with the coordinate centre, and the projected source brightness centre moves in the sky with pulsation phase. As a result of this time-dependent displacement in the brightness centroid, the time of the magnification peak coincides with the closest approach of the lens to the brightness centre as opposed to the source coordinate centre. Binary microlensing of NRP stars and in caustic-crossing features are chromatic., Comment: 14 pages, 23 figures, accepted for publication by MNRAS

KQ Vel is a peculiar A0p star with a strong surface magnetic field of about 7.5 kG. It has a slow rotational period of nearly 8 years. Bailey et al. (2015) detected a binary companion of uncertain nature, and suggested it could be a neutron star or a black hole. In this Letter we analyze X-ray data obtained by the Chandra telescope to ascertain information about the stellar magnetic field and interaction between the star and its companion. We confirm previous X-ray detection of KQ Vel with a relatively large X-ray luminosity of 2\times 10^{30} erg/s. X-ray spectra suggest the presence of hot gas at > 20MK and, possibly, of a non-thermal component. X-ray light curves are variable, but better quality data are needed to determine periodicity if any. We interpret X-ray spectra as a combination of two components: the non-thermal emission arising from the aurora on the A0p star and the hot thermal plasma filling the extended shell surrounding the "propelling" neutron star. We explore various alternatives, but a hybrid model involving the stellar magnetosphere along with a hot shell around the propelling neutron star seems most plausible. We speculate that KQ Vel was originally a triple system, and the Ap star is a merger product. We conclude that KQ Vel is an intermediate-mass binary consisting of a strongly magnetic main sequence star and a neutron star., Comment: accepted to A&AL. 8 pages including the appendix where the formalism describing the X-ray emission from a hot shell around a propelling neutron star in settling accretion regime is presented

The structured winds of single massive stars can be classified into two broad groups: stochastic structure and organized structure. While the former is typically identified with clumping, the latter is typically associated with rotational modulations, particularly the paradigm of Co-rotating Interaction Regions (CIRs). While CIRs have been explored extensively in the UV band, and moderately in the X-ray and optical, here we evaluate radio variability from CIR structures assuming free-free opacity in a dense wind. Our goal is to conduct a broad parameter study to assess the observational feasibility, and to this end, we adopt a phenomenological model for a CIR that threads an otherwise spherical wind. We find that under reasonable assumptions, it is possible to obtain radio variability at the 10% level. The detailed structure of the folded light curve depends not only on the curvature of the CIR, the density contrast of the CIR relative to the wind, and viewing inclination, but also on wavelength. Comparing light curves at different wavelengths, we find that the amplitude can change, that there can be phase shifts in the waveform, and the the entire waveform itself can change. These characterstics could be exploited to detect the presence of CIRs in dense, hot winds., Comment: to appear inn MNRAS

Here, we study the microlensing of radially pulsating stars. Discerning and characterizing the properties of distant, faint pulsating stars is achievable through high-cadence microlensing observations. Combining stellar variability period with microlensing gives the source distance, type, and radius and helps better determine the lens parameters. Considering periodically variations in their radius and surface temperature, their microlensing light curves are resulted from multiplication of the magnification factor with variable finite size effect by the intrinsic brightness curves of pulsing source. The variable finite source size due to pulsation can be significant for transit and single microlensing and while caustic-crossing features. This kind of deviation in the magnification factor is considerable when the ratio of the source radius to the projected lens-source distance is in the range of $\rho_{\star}/u \in[0.4,10]$ and its duration is short and in the same order of the time of crossing the source radius. Other deviations due to variable source intensity and its area make colored and periodic deviations which are asymmetric with respect to the signs of pulsation phase. The positive phases makes deviations with larger amplitude that negative phase. These deviations dominate in filters with short wave lengths (e.g., $B-$band). The position of magnification peaks in microlensing of variable stars varies and this displacement differs in different filters., Comment: 10 pages, 6 figures; Accepted to be published by MNRAS

We have obtained a very deep exposure (813 ks) of $\zeta\,$Puppis (O4 supergiant) with the Chandra/HETG Spectrometer. Here we report on analysis of the 1-9 \r{A} region, especially well suited for Chandra, which has a significant contribution from continuum emission between well separated emission lines from high-ionization species. These data allow us to study the hottest plasma present through the continuum shape and emission line strengths. Assuming a powerlaw emission measure distribution which has a high-temperature cut-off, we find that the emission is consistent with a thermal spectrum having a maximum temperature of 12 MK. This implies an effective wind shock velocity of $900\,\mathrm{km\,s^{-1}}$, well below the wind terminal speed of $2250\,\mathrm{km\,s^{-1}}$. For X-ray emission which forms close to the star, the speed and X-ray flux are larger than can be easily reconciled with strictly self-excited line-deshadowing-instability models, suggesting a need for a fraction of the wind to be accelerated extremely rapidly right from the base. This is not so much a dynamical instability as a nonlinear response to changing boundary conditions., Comment: 18 pages, 4 figures; accepted for publication in the Astrophysical Journal

PolStar is an Explorer-class far ultraviolet (FUV) spectropolarimetry mission designed to target massive stars and their environments. PolStar will take advantage of resonance lines only available in the FUV to measure for the first time the magnetic and wind environment around massive stars to constrain models of rotation and mass loss.

Among Wolf-Rayet stars, those of subtype WN8 are the intrinsically most variable. We have explored the long-term photometric variability of the brightest known WN8 star, WR 40, through four contiguous months of time-resolved, single-passband optical photometry with the BRIght Target Explorer (BRITE) nanosatellite mission. The Fourier transform of the observed light-curve reveals that the strong light variability exhibited by WR 40 is dominated by many randomly-triggered, transient, low-frequency signals. We establish a model in which the whole wind consists of stochastic clumps following an outflow visibility promptly rising to peak brightness upon clump emergence from the optically thick pseudo-photosphere in the wind, followed by a gradual decay according to the right-half of a Gaussian. Free electrons in each clump scatter continuum light from the star. We explore a scenario where the clump size follows a power-law distribution, and another one with an ensemble of clumps of constant size. Both scenarios yield simulated light curves morphologically resembling the observed light curve remarkably well, indicating that one cannot uniquely constrain the details of clump size distribution with only a photometric light curve. Nevertheless, independent evidence favours a negative-index power law, as seen in many other astrophysical turbulent media., Comment: 11 pages, 8 figures, 3 tables; Monthly Notices of the Royal Astronomical Society (MNRAS), in press

We use a 3D Monte Carlo radiative transfer code to study the polarimetric and photometric variability from stationary corotating interaction regions (CIR) in the wind of massive stars. Our CIRs are approximated by Archimedean spirals of higher (or lower) density formed in a spherical wind originating from the star and we also made allowance for a bright gaussian spot at the base of the CIR. Comparing results from our code to previous analytical calculations in the optically thin case, we find differences which we attribute mainly to a better estimation of the total unpolarized flux reaching the observer. In the optically thick case, the differences with the analytical calculations are much larger, as multiple scattering introduces extra complexities including occultation effects. The addition of a gaussian spot does not alter the shape of the polarization curve significantly but does create a small excess in polarization. On the other hand, the effect can be larger on the light curve and can become dominant over the resulting CIR, depending on the spot parameters and density of the wind., Comment: 16 pages, 17 figures

We present new, contemporaneous HST STIS and XMM observations of the O7 III(n)((f)) star xi Per. We supplement the new data with archival IUE spectra, to analyze the variability of the wind lines and X-ray flux of xi Pper. The variable wind of this star is known to have a 2.086 day periodicity. We use a simple, heuristic spot model which fits the low velocity (near surface) IUE wind line variability very well, to demonstrate that the low velocity absorption in the new STIS spectra of N IV 1718 and Si IV 1402 vary with the same 2.086 day period. It is remarkable that the period and amplitude of the STIS data agree with those of the IUE spectra obtained 22 years earlier. We also show that the time variability of the new XMM fluxes are also consistent with the 2.086 day period. Thus, our new, multi-wavelength coordinated observations demonstrate that the mechanism which causes the UV wind line variability is also responsible for a significant fraction of the X-rays in single O stars. The sequence of events for the multi-wavelength light curve minima is: Si IV 1402, N IV 1718, and X-ray flux, each separated by a phase of about 0.06 relative to the 2.086 day period. Analysis of the X-ray fluxes shows that they become softer as they weaken. This is contrary to expectations if the variability is caused by periodic excess absorption. Further, the high resolution X-ray spectra suggest that the individual emission lines at maximum are more strongly blue shifted. If we interpret the low velocity wind line light curves in terms of our model, it implies that there are two bright regions, i.e., regions with less absorption, separated by 180 deg, on the surface of the star. We note that the presence and persistent of two spots separated by 180 deg suggests that a weak dipole magnetic field is responsible for the variability of the UV wind line absorption and X-ray flux in xi Per., Comment: 23 pages, 14 figures

We examine over 160 archival H$\alpha$ spectra from the Ritter Observatory for the interacting binary $\beta$~Lyr obtained between 1996 and 2000. The emission is characteristically double-peaked, but asymmetric, and with an absorption feature that is persistently blueshifted. Using a set of simplifying assumptions, phase varying emission line profiles are calculated for H$\alpha$ formed entirely in a Keplerian disk, and separately for the line formed entirely from an off-center bipolar flow. However, a dynamic spectrum of the data indicate the blueshifted feature is not always present, and the data are even suggestive of a drift of the feature in velocity shift. We explore whether a circumbinary envelope, hot spot on the accretion disk, or accretion stream could explain the observations. While none are satisfactory, an accretion stream explanation is somewhat promising., Comment: accepted to Astronomical Journal

Bow shocks and related density enhancements produced by the winds of massive stars moving through the interstellar medium provide important information regarding the motions of the stars, the properties of their stellar winds, and the characteristics of the local medium. Since bow shocks are aspherical structures, light scattering within them produces a net polarization signal even if the region is spatially unresolved. Scattering opacity arising from free electrons and dust leads to a distribution of polarized intensity across the bow shock structure. That polarization encodes information about the shape, composition, opacity, density, and ionisation state of the material within the structure. In this paper we use the Monte Carlo radiative transfer code SLIP to investigate the polarization created when photons scatter in a bow shock-shaped region of enhanced density surrounding a stellar source. We present results assuming electron scattering, and investigate the polarization behaviour as a function of optical depth, temperature, and source of photons for two different cases: pure scattering and scattering with absorption. In both regimes we consider resolved and unresolved cases. We discuss the implication of these results as well as their possible use along with observational data to constrain the properties of observed bow shock systems. In different situations and under certain assumptions, our simulations can constrain viewing angle, optical depth and temperature of the scattering region, and the relative luminosities of the star and shock., Comment: 20 pages, 13 figures, 1 table; Accepted for publication in MNRAS

The intriguing WN4b star WR6 has been known to display epoch-dependent spectroscopic, photometric and polarimetric variability for several decades. In this paper, we set out to verify if a simplified analytical model in which Corotating Interaction Regions (CIRs) threading an otherwise spherical wind is able to reproduce the many broadband continuum light curves from the literature with a reasonable set of parameters. We modified the optically thin model we developed in Ignace, St-Louis & Proulx-Giraldeau (2015) to approximately account for multiple scattering and used it to fit 13 separate datasets of this star. By including two CIRs in the wind, we obtained reasonable fits for all datasets with coherent values for the inclination of the rotation axis ($i_0=166^{\circ}$) and for its orientation in the plane of the sky, although in the latter case we obtained two equally acceptable values ($\psi=63^{\circ}$ and $\psi=152^{\circ}$) from the polarimetry. Additional line profile variation simulations using the Sobolev approximation for the line transfer allowed us to eliminate the $\psi=152^{\circ}$ solution. With the adopted configuration ($i_0=166^{\circ}$ and $\psi=63^{\circ}$), we were able to reproduce all datasets relatively well with two CIRs located near the stellar equator and always separated by $\sim 90^{\circ}$ in longitude. The epoch-dependency comes from the fact that these CIRs migrate along the surface of the star. Density contrasts smaller than a factor of two and large opening angles for the CIR ($\beta \geq 35^{\circ}$) were found to best reproduce the type of spectroscopic variability reported in the literature., Comment: 15 pages, 5 figures,5 tables

From $5.5$ months of dual-band optical photometric monitoring at the $1$ mmag level, BRITE-Constellation has revealed two simultaneous types of variability in the O4I(n)fp star $\zeta$ Puppis: one single periodic non-sinusoidal component superimposed on a stochastic component. The monoperiodic component is the $1.78$ d signal previously detected by Coriolis/SMEI, but this time along with a prominent first harmonic. The shape of this signal changes over time, a behaviour that is incompatible with stellar oscillations but consistent with rotational modulation arising from evolving bright surface inhomogeneities. By means of a constrained non-linear light curve inversion algorithm we mapped the locations of the bright surface spots and traced their evolution. Our simultaneous ground-based multi-site spectroscopic monitoring of the star unveiled cyclical modulation of its He II $\lambda4686$ wind emission line with the $1.78$-day rotation period, showing signatures of Corotating Interaction Regions (CIRs) that turn out to be driven by the bright photospheric spots observed by BRITE. Traces of wind clumps are also observed in the He II $\lambda4686$ line and are correlated with the amplitudes of the stochastic component of the light variations probed by BRITE at the photosphere, suggesting that the BRITE observations additionally unveiled the photospheric drivers of wind clumps in $\zeta$ Pup and that the clumping phenomenon starts at the very base of the wind. The origins of both the bright surface inhomogeneities and the stochastic light variations remain unknown, but a subsurface convective zone might play an important role in the generation of these two types of photospheric variability., Comment: 44 pages, 28 figures, 6 Tables; To appear in the Monthly Notices of the Royal Astronomical Society (MNRAS)

The transit method, employed by MOST, \emph{Kepler}, and various ground-based surveys has enabled the characterization of extrasolar planets to unprecedented precision. These results are precise enough to begin to measure planet atmosphere composition, planetary oblateness, star spots, and other phenomena at the level of a few hundred parts-per-million. However, these results depend on our understanding of stellar limb darkening, that is, the intensity distribution across the stellar disk that is sequentially blocked as the planet transits. Typically, stellar limb darkening is assumed to be a simple parameterization with two coefficients that are derived from stellar atmosphere models or fit directly. In this work, we revisit this assumption and compute synthetic planetary transit light curves directly from model stellar atmosphere center-to-limb intensity variations (CLIV) using the plane-parallel \textsc{Atlas} and spherically symmetric \textsc{SAtlas} codes. We compare these light curves to those constructed using best-fit limb-darkening parameterizations. We find that adopting parametric stellar limb-darkening laws lead to systematic differences from the more geometrically realistic model stellar atmosphere CLIV of about 50 -- 100 ppm at the transit center and up to 300 ppm at ingress/egress. While these errors are small they are systematic, and appear to limit the precision necessary to measure secondary effects. Our results may also have a significant impact on transit spectra., Comment: 12 pages, 14 figures, accepted for publication in ApJ after revisions

The influence of macroclumps for free-free spectral energy distributions (SEDs) of ionized winds is considered. The goal is to emphasize distinctions between microclumping and macroclumping effects. Microclumping can alter SED slopes and flux levels if the volume filling factor of the clumps varies with radius; however, the modifications are independent of the clump geometry. To what extent does macroclumping alter SED slopes and flux levels? In addressing the question, two specific types of macroclump geometries are explored: shell fragments ("pancake"-shaped) and spherical clumps. Analytic and semi-analytic results are derived in the limiting case that clumps never obscure one another. Numerical calculations based on a porosity formalism is used when clumps do overlap. Under the assumptions of a constant expansion, isothermal, and fixed ionization wind, the fragment model leads to results that are essentially identical to the microclumping result. Mass-loss rate determinations are not affected by porosity effects for shell fragments. By contrast, spherical clumps can lead to a reduction in long-wavelength fluxes, but the reductions are only significant for extreme volume filling factors., Comment: to appear in MNRAS

The Chandra and XMM-Newton X-ray telescopes have led to numerous advances in the study and understanding of astrophysical X-ray sources. Particularly important has been the much increased spectral resolution of modern X-ray instrumentation. Wind-broadened emission lines have been spectroscopically resolved for many massive stars. This contribution reviews approaches to the modeling of X-ray emission line profile shapes from single stars, including smooth winds, winds with clumping, optically thin versus thick lines, and the effect of a radius-dependent photoabsorption coefficient., Comment: to appear in Advances in Space Research

With a deep Chandra/HETGS exposure of WR 6, we have resolved emission lines whose profiles show that the X-rays originate from a uniformly expanding spherical wind of high X-ray-continuum optical depth. The presence of strong helium-like forbidden lines places the source of X-ray emission at tens to hundreds of stellar radii from the photosphere. Variability was present in X-rays and simultaneous optical photometry, but neither were correlated with the known period of the system or with each other. An enhanced abundance of sodium revealed nuclear processed material, a quantity related to the evolutionary state of the star. The characterization of the extent and nature of the hot plasma in WR 6 will help to pave the way to a more fundamental theoretical understanding of the winds and evolution of massive stars., Comment: Accepted by the Astrophysical Journal

The $\beta$ Cephei stars represent an important class of massive star pulsators probing the evolution of B-type stars and the transition from main sequence to hydrogen-shell burning evolution. By understanding $\beta$ Cep stars, we gain insights into the detailed physics of massive star evolution such as rotational mixing, convective core overshooting, magnetic fields and stellar winds, all of which play important roles. Similarly, modeling their pulsation provides additional information into their interior structures. Furthermore, measurements of the rate of change of pulsation period offer a direct measure of $\beta$ Cephei stellar evolution. In this work, we compute state-of-the-art stellar evolution models assuming different amounts of initial rotation and convective core overshoot and measure theoretical rates of period change for which we compare to rates previously measured for a sample of $\beta$ Cephei stars. The results of this comparison are mixed. For three stars, the rates are too small to infer any information from stellar evolution models, whereas for three other stars the rates are too large. We infer stellar parameters, such as mass and age, for two $\beta$ Cephei stars: $\xi^1$ CMa and $\delta$ Cet, that agree well with independent measurements. We explore ideas for why models may not predict the larger rates of period change. In particular, period drifts in $\beta$ Cep stars can artificially lead to overestimated rates of secular period change., Comment: 6 pages, 2 figures 2 tables, accepted to A&A

Classical Cepheids are key probes of both stellar astrophysics and cosmology as standard candles and pulsating variable stars. It is important to understand Cepheids in unprecedented detail in preparation for upcoming GAIA, JWST and extremely-large telescope observations. Cepheid eclipsing binary stars are ideal tools for achieving this goal, however there are currently only three known systems. One of those systems, OGLE-LMC-CEP1812, raises new questions about the evolution of classical Cepheids because of an apparent age discrepancy between the Cepheid and its red giant companion. We show that the Cepheid component is actually the product of a stellar merger of two main sequence stars that has since evolved across the Hertzsprung gap of the HR diagram. This post-merger product appears younger than the companion, hence the apparent age discrepancy is resolved. We discuss this idea and consequences for understanding Cepheid evolution., Comment: 5 pages, 3 figures, accepted to A&A

We report on both high-precision photometry from the MOST space telescope and ground-based spectroscopy of the triple system delta Ori A consisting of a binary O9.5II+early-B (Aa1 and Aa2) with P = 5.7d, and a more distant tertiary (O9 IV P > 400 yrs). This data was collected in concert with X-ray spectroscopy from the Chandra X-ray Observatory. Thanks to continuous coverage for 3 weeks, the MOST light curve reveals clear eclipses between Aa1 and Aa2 for the first time in non-phased data. From the spectroscopy we have a well constrained radial velocity curve of Aa1. While we are unable to recover radial velocity variations of the secondary star, we are able to constrain several fundamental parameters of this system and determine an approximate mass of the primary using apsidal motion. We also detected second order modulations at 12 separate frequencies with spacings indicative of tidally influenced oscillations. These spacings have never been seen in a massive binary, making this system one of only a handful of such binaries which show evidence for tidally induced pulsations., Comment: Accepted to Apj. Part of a 4 paper series. 11 figures, 4 tables

Polarization is a powerful tool for understanding stellar atmospheres and circumstellar environments. Mira and semi-regular variable stars have been observed for decades and some are known to be polarimetrically variable, however, the semi-regular variable V Canes Venatici displays an unusually large, unexplained amount of polarization. We present ten years of optical polarization observations obtained with the HPOL instrument, supplemented by published observations spanning a total interval of about forty years for V CVn. We find that V CVn shows large polarization variations ranging from 1 - 6%. We also find that for the past forty years the position angle measured for V CVn has been virtually constant suggesting a long-term, stable, asymmetric structure about the star. We suggest that this asymmetry is caused by the presence of a stellar wind bow shock and tail, consistent with the star's large space velocity., Comment: 8 pages, 5 figures, 1 table. Accepted for publication in A&A