-
Deriving X-ray Line Profiles for Massive-Star Winds from Momentum-Conserving Dynamical Working Surface Solutions
Authors:
Sean J. Gunderson,
Kenneth G. Gayley
Abstract:
We present a general procedure for deriving a line profile model for massive star X-ray spectra that captures the dynamics of the wind more directly. The basis of the model is the analytic solution to the problem of variable jets in Herbig-Haro objects given by \citet{Canto2000}. In deriving our model, we generalize this jet solution to include flows with a prescribed nonzero acceleration for the…
▽ More
We present a general procedure for deriving a line profile model for massive star X-ray spectra that captures the dynamics of the wind more directly. The basis of the model is the analytic solution to the problem of variable jets in Herbig-Haro objects given by \citet{Canto2000}. In deriving our model, we generalize this jet solution to include flows with a prescribed nonzero acceleration for the context of radiatively driven winds. We provide example line profiles generated from our model for the case of sinusoidal velocity and mass ejection variations. The example profiles show the expected shape of massive star X-ray emission lines, as well as interesting but complicated trends with the model parameters. This establishes the possibility that observed X-rays could be a result of temporal variations seeded at the wind base, rather than purely generated intrinsically within the wind volume, and can be described via a quantitative language that connects with the physical attributes of those variations, consistently with the downstream momentum-conserving nature of radiatively cooled shocked radial flows.
△ Less
Submitted 9 April, 2024;
originally announced April 2024.
-
Observed epochal variations in X-ray lines from the O Supergiant $ζ$ Puppis do not require substantial changes in the wind mass flux
Authors:
Sean J. Gunderson,
Kenneth G. Gayley,
David P. Huenemoerder,
Pragati Pradhan,
Nathan A. Miller
Abstract:
We fit the high resolution \textit{Chandra} X-ray spectra of the O supergiant $ζ$ Puppis using the variable boundary condition (VBC) line model to test the stability of its mass-loss rate between two epochs of observation: 2000 March and 2018 July -- 2019 August. At issue is whether the observed variations are induced by global changes in the cool (unshocked) wind itself or are isolated to the loc…
▽ More
We fit the high resolution \textit{Chandra} X-ray spectra of the O supergiant $ζ$ Puppis using the variable boundary condition (VBC) line model to test the stability of its mass-loss rate between two epochs of observation: 2000 March and 2018 July -- 2019 August. At issue is whether the observed variations are induced by global changes in the cool (unshocked) wind itself or are isolated to the local pockets of hot gas (i.e., changes in the frequency and location of the shocks). Evidence in the literature favored the possibility of a 40 per cent increase in the mass flux of the entire stellar wind, based on X-ray reabsorption from a line-deshadowing-instability-inspired parameterization, whereas our fit parameters are consistent with a constant mass flux with a change in the velocity variations that determine the locations where shocks form. Our results suggest the shocks in the more recent data are formed at somewhat larger radii, mimicking the enhanced blueshifts and increased line fluxes interpreted in the previous analysis as being due to increases in both the X-ray generation and reabsorption from an overall stronger wind.
△ Less
Submitted 20 February, 2024;
originally announced February 2024.
-
Testing the Reliability of X-rays as a Tool for Constraining Mass-loss Rates of Hot Stars
Authors:
Sean J. Gunderson,
Kenneth G. Gayley,
Pragati Pradhan,
David P. Huenemoerder,
Nathan A. Miller
Abstract:
We fit a new line shape model to \textit{Chandra} X-ray spectra of the O supergiant $ζ$ Puppis to test the robustness of mass-loss rates derived from X-ray wind line profiles against different assumed heating models. Our goal is to track the hot gas by replacing the common assumption that it is proportional to the cool gas emission measure. Instead of assuming a turn-on radius for the hot gas (as…
▽ More
We fit a new line shape model to \textit{Chandra} X-ray spectra of the O supergiant $ζ$ Puppis to test the robustness of mass-loss rates derived from X-ray wind line profiles against different assumed heating models. Our goal is to track the hot gas by replacing the common assumption that it is proportional to the cool gas emission measure. Instead of assuming a turn-on radius for the hot gas (as appropriate for the line-deshadowing instability internal to the wind), we parametrize the hot gas in terms of a mean-free path for accelerated low-density gas to encounter slower high-density material. This alternative model is equally successful as previous approaches at fitting X-ray spectral lines in the 5 -- 17 Å wavelength range. We find that the characteristic radii where the hottest gas appears is inversely proportional to line formation temperature, suggesting that stronger shocks appear generally closer to the surface. This picture is more consistent with pockets of low-density, rapid acceleration at the lower boundary than with an internally generated wind instability. We also infer an overall wind mass-loss rate from the profile shapes with a technique used previously in the literature. In doing so, we find evidence that the mass-loss rate derived from X-ray wind line profiles is not robust with respect to changes in the specific heating picture used.
△ Less
Submitted 10 June, 2022;
originally announced June 2022.
-
Method and new tabulations for flux-weighted line-opacity and radiation line-force in supersonic media
Authors:
L. G. Poniatowski,
N. D. Kee,
J. O. Sundqvist,
F. A. Driessen,
N. Moens,
S. P. Owocki,
K. G. Gayley,
L. Decin,
A. de Koter,
H. Sana
Abstract:
In accelerating and supersonic media, the interaction of photons with spectral lines can be of ultimate importance. However, fully accounting for such line forces currently can only be done by specialised codes in 1-D steady-state flows. More general cases and higher dimensions require alternative approaches. We presented a comprehensive and fast method for computing the radiation line-force using…
▽ More
In accelerating and supersonic media, the interaction of photons with spectral lines can be of ultimate importance. However, fully accounting for such line forces currently can only be done by specialised codes in 1-D steady-state flows. More general cases and higher dimensions require alternative approaches. We presented a comprehensive and fast method for computing the radiation line-force using tables of spectral line-strength distribution parameters, which can be applied in arbitrary (multi-D, time-dependent) simulations, including those accounting for the line-deshadowing instability, to compute the appropriate opacities. We assumed local thermodynamic equilibrium to compute a flux-weighted line opacity from $>4$ million spectral lines. We derived the spectral line strength and tabulated the corresponding line-distribution parameters for a range of input densities $ρ\in[10^{-20},10^{-10}]gcm^{-3}$ and temperatures $T\in[10^4,10^{4.7}]K$. We found that the variation of the line distribution parameters plays an essential role in setting the wind dynamics in our models. In our benchmark study, we also found a good overall agreement between the O-star mass-loss rates of our models and those derived from steady-state studies using more detailed radiative transfer. Our models reinforce that self-consistent variation of the line-distribution parameters is important for the dynamics of line-driven flows. Within a well-calibrated O-star regime, our results support the proposed methodology. In practice, utilising the provided tables, yielded a factor $>100$ speed-up in computational time compared to specialised 1-D model-atmosphere codes of line-driven winds, which constitutes an important step towards efficient multi-D simulations. We conclude that our method and tables are ready to be exploited in various radiation-hydrodynamic simulations where the line force is important.
△ Less
Submitted 21 April, 2022;
originally announced April 2022.
-
Ultraviolet Spectropolarimetry: on the origin of rapidly rotating B stars
Authors:
C. E. Jones,
J. Labadie-Bartz,
D. V. Cotton,
Y. Nazé,
G. J. Peters,
D. J. Hillier,
C. Neiner,
N. D. Richardson,
J. L. Hoffman,
A. C. Carciofi,
J. P. Wisniewski,
K. G. Gayley,
M. W. Suffak,
R. Ignace,
P. A. Scowen
Abstract:
UV spectroscopy and spectropolarimetry hold the key to understanding certain aspects of massive stars that are largely inaccessible with optical or longer wavelength observations. This is especially true for the rapidly-rotating Be and Bn stars, owing to their high temperatures, geometric asymmetries, binary properties, and evolutionary history. UV spectropolarimetric observations are extremely se…
▽ More
UV spectroscopy and spectropolarimetry hold the key to understanding certain aspects of massive stars that are largely inaccessible with optical or longer wavelength observations. This is especially true for the rapidly-rotating Be and Bn stars, owing to their high temperatures, geometric asymmetries, binary properties, and evolutionary history. UV spectropolarimetric observations are extremely sensitive to the photospheric consequences of rapid rotation (i.e. oblateness, temperature, and surface gravity gradients). Our polarized radiative-transfer modelling predicts that with low-resolution UV spectropolarimetry covering 120 -- 300 nm the inclination angle of a rapid rotator can be determined to within 5 degrees, and the rotation rate to within 1%. The origin of rapid rotation in Be/n stars can be explained by either single-star or binary evolution, but their relative importance is largely unknown. Some Be stars have hot sub-luminous (sdO) companions, which at an earlier phase transferred their envelope (and with it mass and angular momentum) to the present-day rapid rotator. Through spectral modelling of a wide range of simulated Be/n+sdO configurations, we demonstrate that high-resolution high-SNR UV spectroscopy can detect an sdO star even when $\sim$1,000 times fainter in the UV than its Be/n star companion. This degree of sensitivity is needed to more fully explore the parameter space of Be/n+sdO binaries, which so far has been limited to about a dozen systems with relatively luminous sdO stars. We suggest that a UV spectropolarimetric survey of Be/n stars is the next step forward in understanding this population. Such a dataset would, when combined with population synthesis models, allow for the determination of the relative importance of the possible evolutionary pathways traversed by these stars, which is also crucial for understanding their future evolution and fate.
△ Less
Submitted 27 December, 2022; v1 submitted 15 November, 2021;
originally announced November 2021.
-
The Hydrogen Intensity and Real-time Analysis eXperiment: 256-Element Array Status and Overview
Authors:
Devin Crichton,
Moumita Aich,
Adam Amara,
Kevin Bandura,
Bruce A. Bassett,
Carlos Bengaly,
Pascale Berner,
Shruti Bhatporia,
Martin Bucher,
Tzu-Ching Chang,
H. Cynthia Chiang,
Jean-Francois Cliche,
Carolyn Crichton,
Romeel Dave,
Dirk I. L. de Villiers,
Matt A. Dobbs,
Aaron M. Ewall-Wice,
Scott Eyono,
Christopher Finlay,
Sindhu Gaddam,
Ken Ganga,
Kevin G. Gayley,
Kit Gerodias,
Tim Gibbon,
Austin Gumba
, et al. (75 additional authors not shown)
Abstract:
The Hydrogen Intensity and Real-time Analysis eXperiment (HIRAX) is a radio interferometer array currently in development, with an initial 256-element array to be deployed at the South African Radio Astronomy Observatory (SARAO) Square Kilometer Array (SKA) site in South Africa. Each of the 6m, $f/0.23$ dishes will be instrumented with dual-polarisation feeds operating over a frequency range of 40…
▽ More
The Hydrogen Intensity and Real-time Analysis eXperiment (HIRAX) is a radio interferometer array currently in development, with an initial 256-element array to be deployed at the South African Radio Astronomy Observatory (SARAO) Square Kilometer Array (SKA) site in South Africa. Each of the 6m, $f/0.23$ dishes will be instrumented with dual-polarisation feeds operating over a frequency range of 400-800 MHz. Through intensity mapping of the 21 cm emission line of neutral hydrogen, HIRAX will provide a cosmological survey of the distribution of large-scale structure over the redshift range of $0.775 < z < 2.55$ over $\sim$15,000 square degrees of the southern sky. The statistical power of such a survey is sufficient to produce $\sim$7 percent constraints on the dark energy equation of state parameter when combined with measurements from the Planck satellite. Additionally, HIRAX will provide a highly competitive platform for radio transient and HI absorber science while enabling a multitude of cross-correlation studies. In this paper, we describe the science goals of the experiment, overview of the design and status of the sub-components of the telescope system, and describe the expected performance of the initial 256-element array as well as the planned future expansion to the final, 1024-element array.
△ Less
Submitted 17 January, 2022; v1 submitted 28 September, 2021;
originally announced September 2021.
-
Clumping in the Winds of Wolf-Rayet Stars
Authors:
André-Nicolas Chené,
Nicole St-Louis,
Anthony F. J. Moffat,
Kenneth G. Gayley
Abstract:
We attempt to determine the driver for clumping in hot-star winds by extending the measure of the spectral variability level of Galactic Wolf-Rayet stars to by far the hottest known among them, the WN2 star WR 2 and the WO2 stars WR 102 and WR 142. These three stars have T* = 140 kK and 200 kK, the last two being well above the bulk of WR stars with T* ~ 40-120 kK. This full temperature range for…
▽ More
We attempt to determine the driver for clumping in hot-star winds by extending the measure of the spectral variability level of Galactic Wolf-Rayet stars to by far the hottest known among them, the WN2 star WR 2 and the WO2 stars WR 102 and WR 142. These three stars have T* = 140 kK and 200 kK, the last two being well above the bulk of WR stars with T* ~ 40-120 kK. This full temperature range for WR stars is much broader than that of their O-star progenitors (~30-50 kK), so is better suited to look for any temperature dependence of wind clumping. We have obtained multiple observations with high signal-to-noise, moderate-resolution spectroscopy in search of smallscale variability in the strong emission lines from the dense winds of these three extreme stars, and find a very low-level of variability in both stars. Temperature and terminal velocity are correlated, so faster winds show a lower variability, though this trend goes against any predictions made involving Line Deshadowing Instability (LDI) only, implying that instabilities intrinsic to LDI are not the main source of wind clumping. Instead, it could be taken as support for the suggestion that clumps are caused by a sub-surface convection zone (SSCZ) at T ~ 170 kK, since such an SSCZ would have little opportunity to operate under the hydrostatic surface of these hottest WR stars. It is still possible, however, that an SSCZ-related driver could interact with nonlinear line instability effects to enhance or possibly even produce clumps.
△ Less
Submitted 22 September, 2020;
originally announced September 2020.
-
A Deep Exposure in High Resolution X-Rays Reveals the Hottest Plasma in the $ζ\,$Puppis Wind
Authors:
David P. Huenemoerder,
Richard Ignace,
Nathan A. Miller,
Kenneth G. Gayley,
Wolf-Rainer Hamann,
Jennifer Lauer,
Anthony F. J. Moffat,
Yaël Nazé,
Joy S. Nichols,
Lidia Oskinova,
Noel D. Richardson,
Wayne Waldron
Abstract:
We have obtained a very deep exposure (813 ks) of $ζ\,$Puppis (O4 supergiant) with the Chandra/HETG Spectrometer. Here we report on analysis of the 1-9 Å 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 conti…
▽ More
We have obtained a very deep exposure (813 ks) of $ζ\,$Puppis (O4 supergiant) with the Chandra/HETG Spectrometer. Here we report on analysis of the 1-9 Å 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.
△ Less
Submitted 15 March, 2020;
originally announced March 2020.
-
Progress in the Construction and Testing of the Tianlai Radio Interferometers
Authors:
Santanu Das,
Christopher J. Anderson,
Reza Ansari,
Jean-Eric Campagne,
Daniel Charlet,
Xuelei Chen,
Zhiping Chen,
Aleksander J. Cianciara,
Pierre Colom,
Yanping Cong,
Kevin G. Gayley,
Jingchao Geng,
Jie Hao,
Qizhi Huang,
Celeste S. Keith,
Chao Li,
Jixia Li,
Yichao Li,
Chao Liu,
Tao Liu,
Christophe Magneville,
John P. Marriner,
Jean-Michel Martin,
Marc Moniez,
Trevor M. Oxholm
, et al. (22 additional authors not shown)
Abstract:
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,…
▽ More
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.
△ Less
Submitted 26 June, 2018; v1 submitted 12 June, 2018;
originally announced June 2018.
-
Radio emission and mass loss rate limits of four young solar-type stars
Authors:
Bibiana Fichtinger,
Manuel Güdel,
Robert L. Mutel,
Gregg Hallinan,
Eric Gaidos,
Stephen L. Skinner,
Christene Lynch,
Kenneth G. Gayley
Abstract:
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 usin…
▽ More
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.
△ Less
Submitted 27 February, 2017;
originally announced February 2017.
-
The lunar Askaryan technique with the Square Kilometre Array
Authors:
Clancy W. James,
Jaime Alvarez-Muniz,
Justin D. Bray,
Stijn Buitink,
Rustam D. Dagkesamanskii,
Ronald D. Ekers,
Heino Falcke,
Ken G. Gayley,
Tim Huege,
Maaijke Mevius,
Robert L. Mutel,
Raymond J. Protheroe,
Olaf Scholten,
Ralph E. Spencer,
Sander ter Veen
Abstract:
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 are…
▽ More
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.
△ Less
Submitted 8 August, 2016;
originally announced August 2016.
-
Searching for a magnetic field in Wolf-Rayet stars using FORS2 spectropolarimetry
Authors:
S. Hubrig,
K. Scholz,
W. -R. Hamann,
M. Schoeller,
R. Ignace,
I. Ilyin,
K. G. Gayley,
L. M. Oskinova
Abstract:
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…
▽ More
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.
△ Less
Submitted 4 March, 2016;
originally announced March 2016.
-
Probing Wolf-Rayet Winds: Chandra/HETG X-Ray Spectra of WR 6
Authors:
David P. Huenemoerder,
Kenneth G. Gayley,
Wolf-Rainer Hamann,
Richard Ignace,
Joy S. Nichols,
Lidia Oskinova,
Andrew M. T. Pollock,
Norbert S. Schulz,
Tomer Shenar
Abstract:
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…
▽ More
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.
△ Less
Submitted 5 November, 2015;
originally announced November 2015.
-
The lunar Askaryan technique: a technical roadmap
Authors:
J. D. Bray,
J. Alvarez-Muniz,
S. Buitink,
R. D. Dagkesamanskii,
R. D. Ekers,
H. Falcke,
K. G. Gayley,
T. Huege,
C. W. James,
M. Mevius,
R. L. Mutel,
R. J. Protheroe,
O. Scholten,
R. E. Spencer,
S. ter Veen
Abstract:
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 oth…
▽ More
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.
△ Less
Submitted 17 September, 2015;
originally announced September 2015.
-
Measuring the shock-heating rate in the winds of O stars using X-ray line spectra
Authors:
David H. Cohen,
Zequn Li,
Kenneth G. Gayley,
Stanley P. Owocki,
Jon O. Sundqvist,
Veronique Petit,
Maurice A. Leutenegger
Abstract:
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 coolin…
▽ More
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.
△ Less
Submitted 2 September, 2014;
originally announced September 2014.
-
Lunar detection of ultra-high-energy cosmic rays and neutrinos with the Square Kilometre Array
Authors:
J. D. Bray,
J. Alvarez-Muñiz,
S. Buitink,
R. D. Dagkesamanskii,
R. D. Ekers,
H. Falcke,
K. G. Gayley,
T. Huege,
C. W. James,
M. Mevius,
R. L. Mutel,
R. J. Protheroe,
O. Scholten,
R. E. Spencer,
S. ter Veen
Abstract:
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 a…
▽ More
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.
△ Less
Submitted 19 December, 2014; v1 submitted 26 August, 2014;
originally announced August 2014.
-
Parametrizing Impulsive X-ray Heating with a Cumulative Initial-Temperature Distribution
Authors:
K. G. Gayley
Abstract:
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 usef…
▽ More
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.
△ Less
Submitted 26 January, 2014;
originally announced January 2014.
-
The XMM-Newton EPIC X-ray Light Curve Analysis of WR 6
Authors:
R. Ignace,
K. G. Gayley,
W. R. Hamann,
D. P. Huenemoerder,
L. M. Oskinova,
A. M. T. Pollock,
M. McFall
Abstract:
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 long…
▽ More
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.
△ Less
Submitted 26 July, 2013;
originally announced July 2013.
-
VLA Observations of DG Tau's Radio Jet: A highly collimated thermal outflow
Authors:
C. Lynch,
R. L. Mutel,
M. Güdel,
T. Ray,
S. L. Skinner,
P. C. Schneider,
K. G. Gayley
Abstract:
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 r…
▽ More
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.
△ Less
Submitted 4 February, 2013;
originally announced February 2013.
-
Thin-shell mixing in radiative wind-shocks and the Lx-Lbol scaling of O-star X-rays
Authors:
Stanley P. Owocki,
Jon O. Sundqvist,
David H. Cohen,
Kenneth G. Gayley
Abstract:
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…
▽ More
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.
△ Less
Submitted 18 December, 2012;
originally announced December 2012.
-
High-resolution X-ray spectroscopy reveals the special nature of Wolf-Rayet star winds
Authors:
L. M. Oskinova,
K. G. Gayley,
W. -R. Hamann,
D. P. Huenemoerder,
R. Ignace,
A. M. T. Pollock
Abstract:
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-dr…
▽ More
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.
△ Less
Submitted 18 May, 2012; v1 submitted 7 February, 2012;
originally announced February 2012.
-
UHE neutrino searches using a Lunar target: First Results from the RESUN search
Authors:
T. R. Jaeger,
R. L. Mutel,
K. G. Gayley
Abstract:
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 se…
▽ More
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.
△ Less
Submitted 30 October, 2009;
originally announced October 2009.
-
Asymptotic Opening Angles for Colliding-Wind Bow Shocks: the Characteristic-Angle Approximation
Authors:
K. G. Gayley
Abstract:
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 limi…
▽ More
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.
△ Less
Submitted 9 May, 2009;
originally announced May 2009.
-
Analytic Aperture Calculation and Scaling Laws for Radio Detection of Lunar-Target UHE Neutrinos
Authors:
K. G. Gayley,
R. L. Mutel,
T. R. Jaeger
Abstract:
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.2…
▽ More
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.
△ Less
Submitted 1 May, 2009; v1 submitted 22 April, 2009;
originally announced April 2009.
-
Circumstellar Magnetic Field Diagnostics from Line Polarization
Authors:
Richard Ignace,
Kenneth G. Gayley
Abstract:
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 focu…
▽ More
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.
△ Less
Submitted 20 August, 2007; v1 submitted 14 August, 2007;
originally announced August 2007.
-
A Conceptual Analysis of Mass Loss and Frequency Redistribution in Wolf-Rayet Winds
Authors:
A. J. Onifer,
K. G. Gayley
Abstract:
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…
▽ More
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.
△ Less
Submitted 12 July, 2007;
originally announced July 2007.
-
Wolf-Rayet Mass-Loss Limits Due to Frequency Redistribution
Authors:
Andrew J. Onifer,
Kenneth G. Gayley
Abstract:
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 mul…
▽ More
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.
△ Less
Submitted 26 October, 2005; v1 submitted 13 October, 2005;
originally announced October 2005.
-
A Porosity-Length Formalism for Photon-Tiring-Limited Mass Loss from Stars Above the Eddington Limit
Authors:
Stanley P. Owocki,
Kenneth G. Gayley,
Nir J. Shaviv
Abstract:
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 opa…
▽ More
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.
△ Less
Submitted 23 September, 2004;
originally announced September 2004.
-
Profile Shapes for Optically Thick X-ray Emission Lines from Stellar Winds
Authors:
R. Ignace,
K. G. Gayley
Abstract:
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…
▽ More
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 $ζ$ Pup.
△ Less
Submitted 28 November, 2001;
originally announced November 2001.
-
Line Forces in Keplerian Circumstellar Disks and Precession of Nearly Circular Orbits
Authors:
K. G. Gayley,
R. Ignace,
S. P. Owocki
Abstract:
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 veloci…
▽ More
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.
△ Less
Submitted 24 March, 2001;
originally announced March 2001.
-
One- vs. Two-Shock Heliosphere: Constraining Models with GHRS LyA Spectra toward alpha Cen
Authors:
K. G. Gayley,
G. P. Zank,
H. L. Pauls,
P. C. Frisch,
D. E. Welty
Abstract:
Redshifted L$α$ absorption toward $α$ 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$α$ absorption profile pr…
▽ More
Redshifted L$α$ absorption toward $α$ 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$α$ 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$α$ profiles, and our preliminary conclusion is that the L$α$ profiles seen toward $α$ 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$α$ data along other lines of sight.
△ Less
Submitted 4 June, 1997;
originally announced June 1997.