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.

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.

The Tianlai Pathfinder is designed to demonstrate the feasibility of using a wide field of view radio interferometers to map the density of neutral hydrogen in the Universe after the Epoch of Reionizaton. This approach, called 21~cm intensity-mapping, promises an inexpensive means for surveying the large-scale structure of the cosmos. The Tianlai Pathfinder presently consists of an array of three, 15~m $\times$ 40~m cylinder telescopes and an array of sixteen, 6~m diameter dish antennas located in a radio-quiet part of western China. The two types of arrays were chosen to determine the advantages and disadvantages of each approach. The primary goal of the Pathfinder is to make 3D maps by surveying neutral hydrogen over large areas of the sky %$20,000 {\rm deg}^2$ in two different redshift ranges: first at $1.03 > z > 0.78$ ($700 - 800$~MHz) and later at $0.21 > z > 0.12$ ($1170 - 1270$~MHz). The most significant challenge to $21$~cm intensity-mapping is the removal of strong foreground radiation that dwarfs the cosmological signal. It requires exquisite knowledge of the instrumental response, i.e. calibration. In this paper, we provide an overview of the status of the Pathfinder and discuss the details of some of the analysis that we have carried out to measure the beam function of both arrays. We compare electromagnetic simulations of the arrays to measurements, discuss measurements of the gain and phase stability of the instrument, and provide a brief overview of the data processing pipeline.

The lunar technique is a method for maximising the collection area for ultra-high-energy (UHE) cosmic ray and neutrino searches. The method uses either ground-based radio telescopes or lunar orbiters to search for Askaryan emission from particles cascading near the lunar surface. While experiments using the technique have made important advances in the detection of nanosecond-scale pulses, only at the very highest energies has the lunar technique achieved competitive limits. This is expected to change with the advent of the Square Kilometre Array (SKA), the low-frequency component of which (SKA-low) is predicted to be able to detect an unprecedented number of UHE cosmic rays. In this contribution, the status of lunar particle detection is reviewed, with particular attention paid to outstanding theoretical questions, and the technical challenges of using a giant radio array to search for nanosecond pulses. The activities of SKA's High Energy Cosmic Particles Focus Group are described, as is a roadmap by which this group plans to incorporate this detection mode into SKA-low observations. Estimates for the sensitivity of SKA-low phases 1 and 2 to UHE particles are given, along with the achievable science goals with each stage. Prospects for near-future observations with other instruments are also described.

Observations of free-free continuum radio emission of four young main-sequence solar-type stars (EK Dra, Pi1 UMa, Chi1 Ori; and Kappa1 Cet) are studied to detect stellar winds or at least to place upper limits on their thermal radio emission, which is dominated by the ionized wind. These stars are excellent proxies for representing the young Sun. Upper limits on mass loss rates are calculated using their observational radio emission. Our aim is to re-examine the faint young Sun paradox by assuming that the young Sun was more massive in its past, and hence to find a possible solution for this famous problem. The observations of our sample are performed with the Karl G. Jansky VLA with excellent sensitivity, using the C-band and the Ku-band. ALMA observations are performed at 100 GHz. For the estimation of the mass loss limits, spherically symmetric winds and stationary, anisotropic, ionized winds are assumed. We compare our results to 1) mass loss rate estimates of theoretical rotational evolution models, and 2) to results of the indirect technique of determining mass loss rates: Lyman-alpha absorption. We are able to derive the most stringent direct upper limits on mass loss so far from radio observations. Two objects, EK Dra and Chi1 Ori, are detected at 6 and 14 GHz down to an excellent noise level. These stars are very active and additional radio emission identified as non-thermal emission was detected, but limits for the mass loss rates of these objects are still derived. The stars Pi1 UMa and Kappa1 Cet were not detected in either C-band or in Ku-band. For these objects we give upper limits to their radio free-free emission and calculate upper limits to their mass loss rates. Finally, we reproduce the evolution of the Sun and derive an estimate for the solar mass of the Sun at a younger age.

The lunar Askaryan technique is a method to study the highest-energy cosmic rays, and their predicted counterparts, the ultra-high-energy neutrinos. By observing the Moon with a radio telescope, and searching for the characteristic nanosecond-scale Askaryan pulses emitted when a high-energy particle interacts in the outer layers of the Moon, the visible lunar surface can be used as a detection area. Several previous experiments, at Parkes, Goldstone, Kalyazin, Westerbork, the ATCA, Lovell, LOFAR, and the VLA, have developed the necessary techniques to search for these pulses, but existing instruments have lacked the necessary sensitivity to detect the known flux of cosmic rays from such a distance. This will change with the advent of the SKA. The Square Kilometre Array (SKA) will be the world's most powerful radio telescope. To be built in southern Africa, Australia and New Zealand during the next decade, it will have an unsurpassed sensitivity over the key 100 MHz to few-GHZ band. We introduce a planned experiment to use the SKA to observe the highest-energy cosmic rays and, potentially, neutrinos. The estimated event rate will be presented, along with the predicted energy and directional resolution. Prospects for directional studies with phase 1 of the SKA will be discussed, as will the major technical challenges to be overcome to make full use of this powerful instrument. Finally, we show how phase 2 of the SKA could provide a vast increase in the number of detected cosmic rays at the highest energies, and thus to provide new insight into their spectrum and origin.

To investigate if magnetic fields are present in Wolf-Rayet stars, we selected a few stars in the Galaxy and one in the Large Magellanic Cloud (LMC). We acquired low-resolution spectropolarimetric observations with the ESO FORS2 instrument during two different observing runs. During the first run in visitor mode, we observed the LMC Wolf-Rayet star BAT99 7 and the stars WR6, WR7, WR18, and WR23 in our Galaxy. The second run in service mode was focused on monitoring the star WR6. Linear polarization was recorded immediately after the observations of circular polarization. During our visitor observing run, the magnetic field for the cyclically variable star WR6 was measured at a significance level of 3.3sigma (<B_z> = 258+-78G). Among the other targets, the highest value for the longitudinal magnetic field, <B_z> = 327+-141G, was measured in the LMC star BAT99 7. Spectropolarimetric monitoring of the star WR6 revealed a sinusoidal nature of the <B_z> variations with the known rotation period of 3.77d, significantly adding to the confidence in the detection. The presence of the rotation-modulated magnetic variability is also indicated in our frequency periodogram. The reported field magnitude suffers from significant systematic uncertainties at the factor 2 level, in addition to the quoted statistical uncertainties, owing to the theoretical approach used to characterize it. Linear polarization measurements showed no line effect in the stars, apart from WR6. BAT99 7, WR7, and WR,23 do not show variability of the linear polarization over two nights.

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.

The lunar Askaryan technique, which involves searching for Askaryan radio pulses from particle cascades in the outer layers of the Moon, is a method for using the lunar surface as an extremely large detector of ultra-high-energy particles. The high time resolution required to detect these pulses, which have a duration of around a nanosecond, puts this technique in a regime quite different from other forms of radio astronomy, with a unique set of associated technical challenges which have been addressed in a series of experiments by various groups. Implementing the methods and techniques developed by these groups for detecting lunar Askaryan pulses will be important for a future experiment with the Square Kilometre Array (SKA), which is expected to have sufficient sensitivity to allow the first positive detection using this technique. Key issues include correction for ionospheric dispersion, beamforming, efficient triggering, and the exclusion of spurious events from radio-frequency interference. We review the progress in each of these areas, and consider the further progress expected for future application with the SKA.

We present an overview of four phase-constrained Chandra HETGS X-ray observations of Delta Ori A. Delta Ori A is actually a triple system which includes the nearest massive eclipsing spectroscopic binary, Delta Ori Aa, the only such object which can be observed with little phase-smearing with the Chandra gratings. Since the fainter star, Delta Ori Aa2, has a much lower X-ray luminosity than the brighter primary, Delta Ori A provides a unique system with which to test the spatial distribution of the X-ray emitting gas around Delta Ori Aa1 via occultation by the photosphere of and wind cavity around the X-ray dark secondary. Here we discuss the X-ray spectrum and X-ray line profiles for the combined observation, having an exposure time of nearly 500 ksec and covering nearly the entire binary orbit. Companion papers discuss the X-ray variability seen in the Chandra spectra, present new space-based photometry and ground-based radial velocities simultaneous with the X-ray data to better constrain the system parameters, and model the effects of X-rays on the optical and UV spectrum. We find that the X-ray emission is dominated by embedded wind shock emission from star Aa1, with little contribution from the tertiary star Ab or the shocked gas produced by the collision of the wind of Aa1 against the surface of Aa2. We find a similar temperature distribution to previous X-ray spectrum analyses. We also show that the line half-widths are about $0.3-0.5\times$ the terminal velocity of the wind of star Aa1. We find a strong anti-correlation between line widths and the line excitation energy, which suggests that longer-wavelength, lower-temperature lines form farther out in the wind. Our analysis also indicates that the ratio of the intensities of the strong and weak lines of \ion{Fe}{17} and \ion{Ne}{10} are inconsistent with model predictions, which may be an effect of resonance scattering

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.

The origin of the most energetic particles in nature, the ultra-high-energy (UHE) cosmic rays, is still a mystery. Only the most energetic of these have sufficiently small angular deflections to be used for directional studies, and their flux is so low that even the 3,000 km^2 Pierre Auger detector registers only about 30 cosmic rays per year of these energies. A method to provide an even larger aperture is to use the lunar Askaryan technique, in which ground-based radio telescopes search for the nanosecond radio flashes produced when a cosmic ray interacts with the Moon's surface. The technique is also sensitive to UHE neutrinos, which may be produced in the decays of topological defects from the early universe. Observations with existing radio telescopes have shown that this technique is technically feasible, and established the required procedure: the radio signal should be searched for pulses in real time, compensating for ionospheric dispersion and filtering out local radio interference, and candidate events stored for later analysis. For the Square Kilometre Array (SKA), this requires the formation of multiple tied-array beams, with high time resolution, covering the Moon, with either SKA1-LOW or SKA1-MID. With its large collecting area and broad bandwidth, the SKA will be able to detect the known flux of UHE cosmic rays using the visible lunar surface - millions of square km - as the detector, providing sufficient detections of these extremely rare particles to address the mystery of their origin.

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.

In collisional ionization equilibrium (CIE), the X-ray spectrum from a plasma depends on the distribution of emission measure over temperature (DEM). Due to the well-known ill conditioning problem, no precisely resolved DEM can be inverted directly from the spectrum, so often only a gross parametrization of the DEM is used to approximate the data, in hopes that the parametrization can provide useful model-independent constraints on the heating process. However, ill conditioning also introduces ambiguity into the various different parametrizations that could approximate the data, which may spoil the perceived advantages of model independence. Thus, this paper instead suggests a single parametrization for both the heating mechanism and the X-ray sources, based on a model of impulsive heating followed by complete cooling. This approach is similar to a ``cooling flow'' approach, but allows injection at multiple initial temperatures, and applies even when the steady state is distribution of different shock strengths, as for a standing shock with a range of obliquities, or for embedded stochastic shocks that are only steady in a statistical sense. This produces an alternative parametrization for X-ray spectra that is especially streamlined for higher density plasmas with efficient radiative cooling, and provides internal consistency checks on the assumption of impulsive heating followed by complete cooling. The result is no longer model independent, but the results are more directly interpretable in terms of useful physical constraints on the impulsive heating distribution.

We obtained four pointings of over 100 ks each of the well-studied Wolf-Rayet star WR 6 with the XMM-Newton satellite. With a first paper emphasizing the results of spectral analysis, this follow-up highlights the X-ray variability clearly detected in all four pointings. However, phased light curves fail to confirm obvious cyclic behavior on the well-established 3.766 d period widely found at longer wavelengths. The data are of such quality that we were able to conduct a search for "event clustering" in the arrival times of X-ray photons. However, we fail to detect any such clustering. One possibility is that X-rays are generated in a stationary shock structure. In this context we favor a co-rotating interaction region (CIR) and present a phenomenological model for X-rays from a CIR structure. We show that a CIR has the potential to account simultaneously for the X-ray variability and constraints provided by the spectral analysis. Ultimately, the viability of the CIR model will require both intermittent long-term X-ray monitoring of WR 6 and better physical models of CIR X-ray production at large radii in stellar winds.

The active young protostar DG Tau has an extended jet that has been well studied at radio, optical, and X-ray wavelengths. We report sensitive new VLA full-polarization observations of the core and jet between 5 GHz and 8 GHz. Our high angular resolution observation at 8 GHz clearly shows an unpolarized inner jet with a size 42 AU (0.35") extending along a position angle similar to the optical-X ray outer jet. Using our nearly coeval 2012 VLA observations, we find a spectral-index=+0.46+/-0.05, which combined with the lack of polarization, is consistent with bremsstrahlung (free-free) emission, with no evidence for a non-thermal coronal component. By identifying the end of the radio jet as the optical depth unity surface, and calculating the resulting emission measure, we find our radio results are in agreement with previous optical line studies of electron density and consequent mass-loss rate. We also detect a weak radio knot at 5 GHz located 7" from the base of the jet, coincident with the inner radio knot detected by Rodriguez et al. (2012) in 2009 but at lower surface brightness. We interpret this as due to expansion of post-shock ionized gas in the three years between observations.

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.

We present the first high-resolution X-ray spectrum of a putatively single Wolf-Rayet star. 400 ks observations of WR 6 by the XMM-Newton-telescope resulted in a superb quality high-resolution X-ray spectrum. Spectral analysis reveals that the X-rays originate far out in the stellar wind, more than 30 stellar radii from the photosphere, and thus outside the wind acceleration zone where the line-driving instability could create shocks. The X-ray emitting plasma reaches temperatures up to 50\,MK, and is embedded within the un-shocked, "cool" stellar wind as revealed by characteristic spectral signatures. We detect a fluorescent Fe line at approx 6.4 keV. The presence of fluorescence is consistent with a two-component medium, where the cool wind is permeated with the hot X-ray emitting plasma. The wind must have a very porous structure to allow the observed amount of X-rays to escape. We find that neither the line-driving instability nor any alternative binary scenario can explain the data. We suggest a scenario where X-rays are produced when the fast wind rams into slow "sticky clumps" that resist acceleration. Our new data show that the X-rays in single WR-star are generated by some special mechanism different from the one operating in the O-star winds.

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 ~ Lx ~ Mdot^2 ~ Lbol^3.4 if they are adiabatic. We present here 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 then propose that the associated softening and weakening of the X-ray emission can be parametrized by 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. We conclude by noting that such thin-shell mixing may also be important for X-rays from colliding wind binaries, and that future numerical simulation studies will be needed to test this thin-shell mixing ansatz for X-ray emission.

During the past decade there have been several attempts to detect cosmogenic ultra high energy (UHE) neutrinos by searching for radio Cerenkov bursts resulting from charged impact showers in terrestrial ice or the lunar regolith. So far these radio searches have yielded no detections, but the inferred flux upper limits have started to constrain physical models for UHE neutrino generation. For searches which use the Moon as a target, we summarize the physics of the interaction, properties of the resulting Cerenkov radio pulse, detection statistics, effective aperture scaling laws, and derivation of upper limits for isotropic and point source models. We report on initial results from the RESUN search, which uses the Expanded Very Large Array configured in multiple sub-arrays of four antennas at 1.45 GHz pointing along the lunar limb. We detected no pulses of lunar origin during 45 observing hours. This implies upper limits to the differential neutrino flux E^2 dN/dE < 0.003 EeV km^{-2} s^{-1} sr^{-1} and < 0.0003 EeV km$^{-2} s^{-1} at 90% confidence level for isotropic and sampled point sources respectively, in the neutrino energy range 10^{21.6} < E(eV) < 10^{22.6}. The isotropic flux limit is comparable to the lowest published upper limits for lunar searches. The full RESUN search, with an additional 200 hours observing time and an improved data acquisition scheme, will be be an order of magnitude more sensitive in the energy range 10^{21} < E(eV) < 10^{22} than previous lunar-target searches, and will test Z burst models of neutrino generation.

By considering the advection and interaction of the vector momentum flux in highly supersonic spherically diverging winds, we derive a simple analytic description of the asymptotic opening angle of a wind-collision shock cone, in the approximation that the shocked gas is contained in a cone streaming out along a single characteristic opening angle. Both highly radiative and highly adiabatic limits are treated, and their comparison is the novel result. Analytic closed-form expressions are obtained for the inferred wind momentum ratios as a function of the observed shock opening angle, allowing the conspicuous shape of the asymptotic bow shock to be used as a preliminary constraint on more detailed modeling of the colliding winds. In the process, we explore from a general perspective the limitations in applying to the global shock geometry the so-called Dyson approximation, which asserts a local balance in the perpendicular ram pressure across the shock.

We derive analytic expressions, and approximate them in closed form, for the effective detection aperture for Cerenkov radio emission from ultra-high-energy neutrinos striking the Moon. The resulting apertures are in good agreement with recent Monte Carlo simulations and support the conclusion of James & Protheroe (2009)that neutrino flux upper limits derived from the GLUE search (Gorham et al.2004) were too low by an order of magnitude. We also use our analytic expressions to derive scaling laws for the aperture as a function of observational and lunar parameters. We find that at low frequencies downward-directed neutrinos always dominate, but at higher frequencies, the contribution from upward-directed neutrinos becomes increasingly important, especially at low neutrino energies. Detecting neutrinos from Earth near the GZK regime will likely require radio telescope arrays with extremely large collecting area and hundreds of hour of exposure time. Higher energy neutrinos are most easily detected using lower frequencies. Lunar surface roughness is a decisive factor for obtaining detections at higher frequencies and higher energies.

Results from the most extensive study of the time-evolving dust structure around the prototype "Pinwheel" nebula WR 104 are presented. Encompassing 11 epochs in three near-infrared filter bandpasses, a homogeneous imaging data set spanning more than 6 years (or 10 orbits) is presented. Data were obtained from the highly successful Keck Aperture Masking Experiment, which can recover high fidelity images at extremely high angular resolutions, revealing the geometry of the plume with unprecedented precision. Inferred properties for the (unresolved) underlying binary and wind system are orbital period 241.5 +/- 0.5 days and angular outflow velocity of 0.28 +/- 0.02 mas/day. An optically thin cavity of angular size 13.3 +/- 1.4 mas was found to lie between the central binary and the onset of the spiral dust plume. Rotational motion of the wind system induced by the binary orbit is found to have important ramifications: entanglement of the winds results in strong shock activity far downstream from the nose of the bowshock. The far greater fraction of the winds participating in the collision may play a key role in gas compression and the nucleation of dust at large radii from the central binary and shock stagnation point. Investigation of the effects of radiative braking pointed towards significant modifications of the simple hydrostatic colliding wind geometry, extending the relevance of this phenomena to wider binary systems than previously considered. Limits placed on the maximum allowed orbital eccentricity of e < 0.06 argue strongly for a prehistory of tidal circularization in this system. Finally we discuss the implications of Earth's polar (i < 16 deg) vantage point onto a system likely to host supernova explosions at future epochs.

Given that dynamically significant magnetic fields in at least some massive stars have now been measured, our contribution addresses the question, to what extent can fields be directly detected in circumstellar gas? The question speaks directly to the very interesting topic of line-driving physics coupled with magnetized plasmas, and how this coupling produces structure in the wind flow. We focus our attention on weak-field diagnostics. These come in two main types: the Hanle effect, which pertains to coherence effects for linear polarization from line scattering, and the weak longitudinal Zeeman effect, which pertains to circular polarization in lines.

To better understand Wolf-Rayet stars as progenitors of gamma-ray bursts, an understanding of the effect metallicity has on Wolf-Rayet mass loss is needed. Using simple analytic models, we study the Mdot - Z relation of a WN star and compare the results to similar models. We find that Mdot roughly follows a power law in Z with index 0.88 from -2.5 < log Z/Z_sun < -1 and appears to flatten by log Z/Z_sun ~ -0.5.

The hypothesis that CAK-type line driving is responsible for the large observed Wolf-Rayet (W-R) mass-loss rates has been called into question in recent theoretical studies. The purpose of this paper is to reconsider the plausibility of line driving of W-R winds within the standard approach using the Sobolev approximation while advancing the conceptual understanding of this topic. Due to the multiple scattering required in this context, of particular importance is the role of photon frequency redistribution into spectral gaps, which in the extreme limit yields the statistical Sobolev-Rosseland (SSR) mean approximation. Interesting limits to constrain are the extremes of no frequency redistribution, wherein the small radii and corresponding high W-R surface temperature induces up to twice the mass-loss rate relative to cooler stars, and the SSR limit, whereby the reduced efficiency of the driving drops the mass flux by as much as an order of magnitude whenever there exist significant gaps in the spectral line distribution. To see how this efficiency drop might be sufficiently avoided to permit high W-R mass loss, we explore the suggestion that ionization stratification may serve to fill the gaps globally over the wind. We find that global ionization changes can only fill the gaps sufficiently to cause about a 25% increase in the mass-loss rate over the local SSR limit. Higher temperatures and more ionization states (especially of iron) may be needed to achieve optically thick W-R winds, unless strong clumping corrections eliminate the need for such winds.

We examine radiatively driven mass loss from stars near and above the Eddington limit (Ledd). We begin by reviewing the instabilities that are expected to form extensive structure near Ledd. We investigate how this "porosity" can reduce the effective coupling between the matter and radiation. Introducing a new "porosity-length'' formalism, we derive a simple scaling for the reduced effective opacity, and use this to derive an associated scaling for the porosity-moderated, continuum-driven mass loss rate from stars that formally exceed Ledd. For a simple super-Eddington model with a single porosity length that is assumed to be on the order of the gravitational scale height, the overall mass loss is similar to that derived in previous porosity work. This is much higher than is typical of line-driven winds, but is still only a few percent of the photon tiring limit--for which the luminosity becomes insufficient to carry the flow out of the gravitational potential. To obtain still stronger mass loss that approaches observationally inferred values near this limit, we introduce a power-law-porosity model in which the associated structure has a broad range of scales. We show that the mass loss rate can be enhanced by a factor that increases with the Eddington parameter Gamma, such that for moderately large Gamma (> 3-4), mass loss rates could approach the photon tiring limit. Together with the ability to drive quite fast outflow speeds, the derived mass loss could explain the large inferred mass loss and flow speeds of giant outbursts in eta Carinae and other LBV stars.

We consider the consequences of appreciable line optical depth for the profile shape of X-ray emission lines formed in stellar winds. The hot gas is thought to arise in distributed wind shocks, and the line formation is predominantly via collisional excitation followed by radiative decay. Such lines are often modelled as optically thin, but the theory has difficulty matching resolved X-ray line profiles. We suggest that for strong lines of abundant metals, newly created photons may undergo resonance scattering, modifying the emergent profile. Using Sobolev theory in a spherically symmetric wind, we show that thick-line resonance scattering leads to emission profiles that still have blueshifted centroids like the thin lines, but which are considerably less asymmetric in appearance. We focus on winds in the constant-expansion domain, and derive an analytic form for the profile shape in the limit of large line and photoabsorptive optical depths. Our theory is applied to published {\it Chandra} observations of the O star $\zeta$ Pup.

We examine the effects of optically thick line forces on orbiting circumstellar disks, such as occur around Be stars. For radially streaming radiation, line forces are only effective if there is a strong radial velocity gradient, as occurs, for example, in a line-driven stellar wind. However, within an orbiting disk, the radial shear of the azimuthal velocity leads to strong line-of-sight velocity velocity gradients along nonradial directions. As such, in the proximity of a stellar surface extending over a substantial cone angle, the nonradial stellar radiation can impart a significant line force, even in the case of purely circular orbits. Given the highly supersonic nature of orbital velocity variations, we use the Sobolev approximation, thereby extending to the disk case the standard CAK formalism developed for line-driven winds. We delineate the parameter regimes for which radiative forces might alter disk properties; but even when radiative forces are small, we analytically quantify higher-order effects in the linear limit, including the precession of weakly elliptical orbits. We find that optically thick line forces can have observable implications for the dynamics of disks around Be stars, including the generation of either prograde or retrograde precession in slightly eccentric orbits. However, our analysis suggests a net retrograde effect, in apparent contradiction with observed long-term variations of violet/red line profile asymmetries from Be stars, which are generally thought to result from prograde propagation of a so-called ``one arm mode''. We also conclude that radiative forces may alter the dynamical properties at the surface of the disk where disk winds originate, and may even make low-density disks vulnerable to being blown away.

We analyse two ASCA observations of the highly eccentric O9III + B1III binary Iota Orionis obtained at periastron and apastron. Based on the assumption of a strong colliding winds shock between the stellar components, we expected to see significant variation in the X-ray emission between these phases. The observations proved otherwise: the X-ray luminosities and spectral distributions were remarkably similar. The only noteworthy feature was the hint of a proximity effect during periastron passage, supported also in the optical. We discuss the accuracy of our results, and also analyse archival ROSAT observations. We investigate why we do not see a clear colliding winds signature. A simple model shows that the wind attenuation to the expected position of the shock apex is negligible throughout the orbit, which poses the puzzling question of why the expected 1/D variation (ie. a factor of 7.5) in the intrinsic luminosity is not seen in the data. Two scenarios are proposed: either the colliding winds emission is unexpectedly weak such that intrinsic shocks in the winds dominate the emission, or, alternatively, that the emission observed is colliding winds emission but in a more complex form than we would naively expect. Complex hydrodynamical models are then analyzed. Despite strongly phase-variable emission from the models, both were consistent with the observations. We find that if the mass-loss rates of the stars are low then intrinsic wind shocks could dominate the emission. However, when we assume higher mass-loss rates of the stars, we find that the observed emission could also be consistent with a purely colliding winds origin. To distinguish between the different models X-ray observations with improved phase coverage will be necessary.

Redshifted L$\alpha$ absorption toward $\alpha$ Cen has been interpreted by Linsky & Wood (1996) and Frisch et al. (1996) as evidence for decelerated interstellar hydrogen piled up on the upstream side of the heliosphere. We utilize newly developed two-dimensional multi-fluid models of the solar wind interaction with the ISM to corroborate this interpretation by synthesizing the L$\alpha$ absorption profile predicted for this ``hydrogen wall''. Both subsonic and supersonic inflow into the heliosphere are considered, corresponding to one-shock and two-shock global morphologies, respectively. It is found that these two extremes give observably different redward absorption characteristics in the L$\alpha$ profiles, and our preliminary conclusion is that the L$\alpha$ profiles seen toward $\alpha$ Cen favor a barely subsonic model (Mach number 0.9). For such a model to hold, additional interstellar pressure terms, such as cosmic ray or magnetic pressures, must contribute. To make this conclusion more certain, an extended model-parameter survey is required, coupled with L$\alpha$ data along other lines of sight.