An official website of the United States government
Abstract 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 (WR) 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 WR 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 the 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 disagree 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 WR binary systems and increase our understanding of their evolutionary paths.
more »
« less
Modelling variable linear polarization produced by Co-rotating Interaction Regions (CIRs) across optical recombination lines of Wolf–Rayet stars
ABSTRACT Massive star winds are structured both stochastically (‘clumps’) and often coherently (Co-rotation Interaction Regions, or CIRs). Evidence for CIRs threading the winds of Wolf–Rayet (WR) stars arises from multiple diagnostics including linear polarimetry. Some observations indicate changes in polarization position angle across optical recombination emission lines from a WR star wind but limited to blueshifted Doppler velocities. We explore a model involving a spherical wind with a single conical CIR stemming from a rotating star as qualitative proof-of-concept. To obtain a realistic distribution of limb polarization and limb darkening across the pseudo-photosphere formed in the optically thick wind of a WR star, we used Monte Carlo radiative transfer (MCRT). Results are shown for a parameter study. For line properties similar to WR 6 (EZ CMa; HD 50896), combining the MCRT results, a simple model for the CIR, and the Sobolev approximation for the line formation, we were able to reproduce variations in both polarization amplitude and position angle commensurate with observations. Characterizing CIRs in WR winds has added importance for providing stellar rotation periods since the vsin i values are unobtainable because the pseudo-photosphere forms in the wind itself.
more »
« less
- Award ID(s):
- 2009412
- PAR ID:
- 10465694
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 526
- Issue:
- 1
- ISSN:
- 0035-8711
- Format(s):
- Medium: X Size: p. 1298-1307
- Size(s):
- p. 1298-1307
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract 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 Brown et al. for variable polarization from binaries with the semianalytic solution for the geometry and surface density CWI shock interface between the winds based on Cantó et al. Our calculations include some simplifications in the form of inverse-square law wind densities and the assumption of axisymmetry, but in so doing they 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.more » « less
-
The effect of the turbulence that is associated with solar wind corotating interaction regions (CIRs) on transport of galactic cosmic rays remains an outstanding problem in space science. Observations show that the intensities of the plasma and magnetic fluctuations are enhanced within a CIR. The velocity shear layer between the slow and fast wind embedded in a CIR is thought to be responsible for this enhancement in turbulent energy. We perform physics-based magnetohydrodynamic simulations of the plasma background and turbulent fluctuations in the solar wind dominated by CIRs for radial distances between 0.3 and 5 au. A simple but effective approach is used to incorporate the inner boundary conditions for the solar wind and magnetic field for the periods 2007–2008 and 2017–2018. Legendre coefficients at the source surface obtained from the Wilcox Solar Observatory library are utilized for dynamic reconstructions of the current sheet and the fast and slow streams at the inner boundary. The dynamic inner boundary enables our simulations to generate CIRs that are reasonably comparable with observations near Earth. While the magnetic field structure is reasonably well reproduced, the enhancements in the turbulent energy at the stream interfaces are smaller than observed. A superposed epoch analysis is performed over several CIRs from the simulation and compared to the superposed epoch analysis of the observed CIRs. The results for the turbulent energy and correlation length are used to estimate the diffusion tensor of galactic cosmic rays. The derived diffusion coefficients could be used for more realistic modeling of cosmic rays in a dynamically evolving inner heliosphere.more » « less
-
null (Ed.)Massive Wolf-Rayet (WR) stars in binary systems may produce supernovae capable of emitting long duration gamma ray bursts. Characterizing the structure of the colliding winds in these systems may help constrain the mass loss and transfer properties and help predict their future evolution. I will present new spectropolarimetric data for the possible WR+O binary system WR 71, collected using RSS at the Southern African Large Telescope. WR 71 is a WN6 whose binary status is unknown, but it displays similar spectropolarimetric variations to the known WR+O binary system V444 Cygni. I investigate the orbital and rotational velocity of WR 71's winds by analyzing its polarized emission line profiles as a function of phase, the first analysis of its kind. I compare the line polarization behavior with predictive models of both colliding wind binaries and single stars with co-rotating interaction regions. Describing the wind structure of WR 71 will help determine the rate of mass loss from the system, an important indicator for LGRB progenitors, and shed light on its binary status.more » « less
-
ABSTRACT The bright WN4 star EZ CMa exhibits a 3.77 d 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 d 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 d 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 v λ 4945 line. We derive an orbital solution for the system and reject the 3.626 d period to represent the variations in the radial velocities of EZ CMa. An orbital solution with an orbital period of 3.77 d 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 ± 0.001 d 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.more » « less
