Engineering Papers⌕ Search

Engineering topics

Macfarlane, J. J.

Publications and source records attributed to Macfarlane, J. J..

ASCA solid state imaging spectrometer observations of O stars

We report ASCA Solid State Imaging Spectrometer (SIS) x-ray observations of the O stars delta Ori and lambda Ori. The energy resolution of the SIS allows us to resolve features in the O star x-ray spectra which are not apparent in spectra obtained by x-ray spectrometers with lower energy resolution. SIS spectra from both stars show evidence of line emission, suggesting the thermal nature of the x-ray source. However, the observed line strengths are different for the two stars. The observed stellar x-ray spectra are not well described by isothermal models although absorbed thermal emission models with two or more temperatures can provide an adequate fit to the data. For both stars we present evidence of absorbing columns significantly larger than the known ISM columns, indicative of absorption by a circumstellar medium, presumably the stellar winds. In addition, the lambda Ori spectrum shows the presence of emission at energies greater than 3 keV which is not seen in the delta Ori spectrum.

Corcoran, M. F.↗

EUVE spectroscopy of epsilon Canis Majoris (B2 II) from 70 to 730 A

We present spectra of the brightest stellar source of extreme ultraviolet (EUV) radiation longward of 400 A, the B2 II star, epsilon CMa. These data were taken with the three spectrometers aboard the NASA Extreme Ultraviolet Explorer satellite (EUVE) during the first cycle of pointed observations. We report on our initial studies of the continuum and line spectrum of the stellar photosphere in the 320 to 730 A region, and on the wind emission lines observed in the 170-375 A region. This is the first EUV spectrum of an early-type star, and thus makes epsilon CMa the most comprehensively observed B star from the X-ray to infrared regimes. The radiation in both the H Lyman continuum and He I continuum (shortward of 504 A) are found to be significantly greater than predicted by both Local Thermodynamic Equilibrium (LTE) and non-LTE model atmospheres. Since epsilon CMa also exhibits a mid-infrared excess, this points to the outer layers being warmer than the models indicate. The anomalously large Lyman continuum flux, combined with the very low column density measured in the direction toward this star implies that it is the dominant source of hydrogen ionization of the local interstellar medium in the immediate vicinity of the sun. All of the lines predicted to be strong from model atmospheres are present and several wind absorption features are also identified. We have detected emission lines from highly ionized iron that are consistent with the ROSAT Position Sensitive Proportional Counter (PSPC) observations if a multi-temperature emission model is used, and the assumption is made that there is significant absorption beyond that of the neutral phase of the ISM. The spectrum shows strong O III 374 A line emission produced by the Bowen flourescence mechanism, which has not previously been observed in the spectra of hot stars.

Cassinelli, J. P.↗

X-ray induced ionization in the winds of near-main-sequence O and B stars

We investigate the effects of X-rays on the ionization state of stellar winds for O and early-B stars along the main sequence. In our investigation, detailed statistical equilibrium, radiative transfer, and atomic physics models are used to compute ionization distributions for H, He, C, N, O, and Si. X-rays are modeled as a spatially distributed source within the wind, with parameters constrained by ROSAT and Einstein observations. Our results indicate that the ionization balance in the winds of stars with spectral type B0 V and later is significantly altered by the X-ray radiation field. Unlike the case of denser O star winds, where the X-rays tend to perturb the level of ionization, the ionization state of the bulk wind of early-B stars can be significantly increased by soft X-rays. We examine in detail the case of tau Sco (B0 V), which has been well-studied at UV and X-ray wavelengths. Comparisons are made between calculated ionization fractions and those deduced from UV observations. In addition, we address the sensitivity of our results to the X-ray source characteristics, the wind temperature, and the photospheric extreme ultraviolet (EUV) flux. Our results suggest the possibility that for early-B stars X-rays play a critical role in both influencing the radiation line driving force, as well as ionizing and heating the wind all the way down to the top of the photosphere.

Macfarlane, J. J.↗

ASCA solid state imaging spectrometer observations of O stars

We report ASCA Solid State Imaging Spectrometer (SIS) X-ray observations of the O stars delta Ori and lambda Ori. The energy resolution of the SIS allows us to resolve features in the O star X-ray spectra which are not apparent in spectra obtained by X-ray spectrometers with lower energy resolution. SIS spectra from both stars show evidence of line emission, suggesting the thermal nature of the X-ray source. However, the observed line strengths are different for the two stars. The observed stellar X-ray spectra are not well described by isothermal models although absorbed thermal emission models with two or more temperatures can provide an adequate fit to the data. For both stars we present evidence of absorbing columns significantly larger than the known Interstellar Medium (ISM) columns, indicative of absorption by a circumstellar medium, presumably the stellar winds. In addition, the lambda Ori spectrum shows the presence of emission at energies greater than 3 keV which is not seen in the delta Ori spectrum.

Corcoran, M. F.↗

Constraints on hot star X-ray source characteristics from combinded analysis of X-ray and UV observations

Results from wind ionization calculations are presented which show how the P-Cygni profiles of 'superionized' species such as O VI can provide information about the X-ray source characteristics of early-type stars. Using detailed radiative and atomic physics models, we find that a significant source of X-ray emission from zeta Pup comes from a region in the wind located within rougly 1 to 2 stellar radii of the photosphere. Our results suggest that X-rays sources in which emission occurs exclusively at large radii (r greater than or approximately equal to a few R(sub *)) are inconsistent with UV P-Cygni profiles for O VI. Instead, we find that X-ray emission from shocks distributed throughout the lower regions of the wind (r approximately equal to 1-2 R(sub *)) is consistent with both X-ray and UV data, as well as mass loss rates deduced from radio and H-alpha observations.

Macfarlane, J. J.↗

X-ray emission from near-main-sequence B stars

Results of a Roengten Satellite (ROSAT) X-ray survey of 12 nearby near-main-sequence B stars are presented. Objects with very low interstellar hydrogen column density were chosen to study the soft X-ray emission properties of the stars. All of the stars were detected at the 3 sigma level. Spectral fits to seven of the stars are presented, and temperatures and source emission measures are derived. The spectra and characterized by emission from gas at a temperature of about 2 x 10(exp 6) K, which is lower than that typically observed for O stars. The ratio L(sub x)/L(sub Bol) decreases sharply with spectral types later than B1, reaching values of about 10(exp -9) at B3 V; much less than the value of 10(exp -7) that holds for O stars and early-B stars. The survey includes four beta Cephei stars and two Be stars. Our observations are consistent with the finding that beta Cephei variables have softer spectra than other stars of the same spectral type. The X-ray emission measures for the stars are compared with upper and lower limits derived using mass loss rates and terminal velocities from line-driven wind theory. It is concluded that if shocks embedded in the winds are the source of the X-ray emission then a significant fraction of the winds in B stars must be hot. One to a few shocks can account for the observed emission measures. As a function of spectral type, there appears to be a transition in X-ray emission properties at about B1 to B1.5, with later type stars having a smaller L(sub x)/L(sub Bol) ratio. Coronal and shock models are discussed in regards to the BO V star tau Sco which has an anomalously high source temperature.

Cassinelli, J. P.↗

Non-LTE radiative transfer with lambda-acceleration - Convergence properties using exact full and diagonal lambda-operators

We investigate the convergence properties of Lambda-acceleration methods for non-LTE radiative transfer problems in planar and spherical geometry. Matrix elements of the 'exact' A-operator are used to accelerate convergence to a solution in which both the radiative transfer and atomic rate equations are simultaneously satisfied. Convergence properties of two-level and multilevel atomic systems are investigated for methods using: (1) the complete Lambda-operator, and (2) the diagonal of the Lambda-operator. We find that the convergence properties for the method utilizing the complete Lambda-operator are significantly better than those of the diagonal Lambda-operator method, often reducing the number of iterations needed for convergence by a factor of between two and seven. However, the overall computational time required for large scale calculations - that is, those with many atomic levels and spatial zones - is typically a factor of a few larger for the complete Lambda-operator method, suggesting that the approach should be best applied to problems in which convergence is especially difficult.

Macfarlane, J. J.↗

Statistical mechanics of light elements at high pressure. VIII - Thomas-Fermi-Dirac theory for binary mixtures of H with He, C, and O

We present three-dimensional Thomas-Fermi-Dirac calculations of lattice mixing energies of hydrogen with carbon and oxygen atoms, respectively. The results are used to derive effective interatomic potentials for use in liquid-state mixture calculations. We then use the potentials to derive analytic expressions for binary mixture-free energies and to map out the phase diagrams of mixtures of hydrogen with, respectively, helium, carbon, and oxygen, over a pressure range of about 5 to about 10 to the 3rd Mbar. Within this pressure range, all three of the latter elements are found to have unlimited solubility in metallic hydrogen over a temperature range which lies above their pure-element melting temperatures, and which includes likely interior temperatures in the Jovian planets.

Hubbard, W. B.↗

Statistical mechanics of light elements at high pressure. VI - Liquid-state calculations with Thomas-Fermi-Dirac theory

A model free energy is developed for hydrogen-helium mixtures based on solid-state Thomas-Fermi-Dirac calculations at pressures relevant to the interiors of giant planets. Using a model potential similar to that for a two-component plasma, effective charges for the nuclei (which are in general smaller than the actual charges because of screening effects) are parameterized, being constrained by calculations at a number of densities, compositions, and lattice structures. These model potentials are then used to compute the equilibrium properties of H-He fluids using a charged hard-sphere model. The results find critical temperatures of about 0 K, 500 K, and 1500 K, for pressures of 10, 100, and 1000 Mbar, respectively. These phase separation temperatures are considerably lower (approximately 6,000-10,000 K) than those found from calculations using free electron perturbation theory, and suggest that H-He solutions should be stable against phase separation in the metallic zones of Jupiter and Saturn.

Macfarlane, J. J.↗

Statistical mechanics of light elements at high pressure. V Three-dimensional Thomas-Fermi-Dirac theory

A numerical technique for solving the Thomas-Fermi-Dirac (TED) equation in three dimensions, for an array of ions obeying periodic boundary conditions, is presented. The technique is then used to calculate deviations from ideal mixing for an alloy of hydrogen and helium at zero temperature and high presures. Results are compared with alternative models which apply perturbation theory to calculation of the electron distribution, based upon the assumption of weak response of the electron gas to the ions. The TFD theory, which permits strong electron response, always predicts smaller deviations from ideal mixing than would be predicted by perturbation theory. The results indicate that predicted phase separation curves for hydrogen-helium alloys under conditions prevailing in the metallic zones of Jupiter and Saturn are very model dependent.

Macfarlane, J. J.↗

Internal structure of Uranus

An updated study is presented of Uranus interior models using current information about the planet's gravity field and rotation rate. The most plausible model, both from the point of view of recent data and cosmogony, has a central core of iron and magnesium silicates, an outer envelope of liquid water, methane, and ammonia, and a deep 'atmosphere' of almost four earth masses of hydrogen, helium, and methane. The 'atmosphere' contains a gravitationally nonnegligible amount of methane - about 40% by mass. All plausible models are most consistent with a rotation period of 15 to 16 hours.

Macfarlane, J. J.↗

Theoretical predictions of deuterium abundances in the Jovian planets

Current concepts for the origin of the Jovian planets and current constraints on their interior structure are used to support the argument that the presence of large amounts of 'ice' (H2O, CH4, and NH3) in Uranus and Neptune indicates temperature low enough to condense these species at the time Uranus and Neptune formed. Such low temperatures, however, imply orders-of-magnitude fractionation effects for deuterium into the 'ice' component if isotopic equilibration can occur. The present models thus imply that Uranus and Neptune should have D/H ratio at least four times primordial, contrary to observation for Uranus. It is found that the Jovian and Saturnian D/H should be close to primordial regardless of formation scenario.

Hubbard, W. B.↗

Interior structure of Saturn inferred from Pioneer 11 gravity data

The structure of Saturn is studied via a fourth-order theory for rotating planets and equations of state for the envelope which depend parametrically on the helium abundance, on the starting temperature for the adiabat, and on adopted forms of the pressure-density curve in the region of transition from molecular to metallic hydrogen. Models are constrained by the values of J2 and J4 obtained from the Pioneer-Saturn celestial mechanics experiment. Equations of state are tested by computing Jupiter models, which can now be subjected to a more stringent comparison with observed zonal harmonics. It is found that Saturn has a low-density hydrogen-helium envelope with no evidence for enhancement of H2O, CH4, or other abundant compounds. Such compounds are presumably located near the core. The helium mass abundance for Saturn's envelope appears to be in the range of approximately 0.12 to 0.19, but this result is very model-dependent. The helium abundance in the envelope of Jupiter is apparently very similar to that of Saturn.

Hubbard, W. B.↗

Pioneer Saturn celestial mechanics experiment

The paper investigates a continuous round-trip radio link at S band (2.2 GHz) that was maintained during the Pioneer Saturn encounter between stations of the Deep Space Network and the spacecraft. From an analysis of the Doppler shift in the radio carrier frequency, a number of gravitational effects on the trajectory were determined. Gravitational moments for Saturn were found from a preliminary analysis, as well as mass values for the Saturn satellites Rhea, Iapetus, and Titan. It was determined that the densities of all three satellites are low and consistent with the compositions of ices. Theoretical calculations for the Saturn interior are described which use the latest observational data, including Pioneer Saturn and state-of-the-art physics for the internal composition.

Anderson, J. D.↗

Structure and evolution of Uranus and Neptune

Three-layer interior models of Uranus and Neptune with central rocky cores, mantles of water, methane, and ammonia (the 'ices'), and outer envelopes primarily composed of hydrogen and helium are presented. The models incorporate a new H2O equation of state based on experimental data which is considerably 'softer' than previous H2O equations of state. Corrections for interior temperatures approximately 5000 K are included in the models, and the thermal evolution of both planets is investigated using recent heat flow measurements. It is found that the evolutionary considerations are consistent with gravitational field data in supporting models with approximately solar abundances of 'ice' and 'rock'. Evolutionary considerations indicate that initial temperatures and luminosities for Uranus and Neptune were not substantially higher than the present value. Both planets apparently have relatively small approximately 1-2 earth masses) hydrogen-helium envelopes, with Neptune's envelope smaller than Uranus'. A monotonic trend is evident among the Jovian planets: all have central rock-ice cores of approximately 15 earth masses, but with hydrogen-helium envelopes which decrease in mass from Jupiter to Saturn to Uranus to Neptune.

Hubbard, W. B.↗