The population of massive stars in R136 from HST/FOC UV observations

New ultraviolet (1300 A, 3400 A),HST FOC observations have been used to derive the UV color-magnitude diagram (CMD) of R136, with the main scientific goal of studying the upper end of the stellar mass function at ultraviolet wavelengths where the color degeneracy encountered in visual CMDs is less severe. The CMD has been compared to a set of theoretical isochrones, which have been computed using the latest generation of evolutionary models and model atmospheres for early type stars. Wolf-Rayet stars are included. Comparison of theTheoretical andobserved CMD suggests that there are no stars brighter than M₁₃₀–11. We use the observed main sequence turn-off and the known spectroscopic properties of the stellar population to derive constraints on the most probable age of R136. The presence of WNL stars and the lack of red supergiants suggests a most likely age of 3±1 Myr. A theoretical isochrone of 3±1 Myr is consistent with the observed stellar content of R136 if the most massive stars have initial masses around 50 M.Bases on Observations with the NASA/ESA Hubble Space Telescope, obtained at the STScI, which is operated by AURA, Inc., under NASA contract NAS5-26555.Astrophysics Division, Space Science Department, ESA     

Spectral analyses of Wolf-Rayet stars: Theory,results, conclusions

Stratified Non-LTE models for expanding atmospheres became available in the recent years. They are based on the idealizing assumptions of spherical symmetry, stationarity and radiative equilibrium. From a critical discussion we conclude that this standard model is basically adequate for describing real Wolf-Rayet atmospheres and hence can be applied for quantitative spectral analyses of their spectra.By means of these models, the fundamental parameters have been determined meanwhile for the majority of the known Galactic WR stars. Most of them populate a vertical strip in the Herzsprung-Russell diagram at effective temperatures of 35 kK, the luminosities ranging from 10^4.5 to 10^5.9 L . Only early-type WN stars with strong lines and WC stars are hotter. The chemical composition of WR atmospheres corresponds to nuclear-processed material (WN: hydrogen burning in the CNO cycle; WC: helium burning). Hydrogen is depleted but still detectable in the cooler part of the WN subclass.Different scenarios for the evolutionary formation of the Wolf-Rayet stars are discussed in the light of the empirical data provided from the spectral analyses. Post-red-supergiant evolution can principally explain the basic observational properties, except the rather low luminosities of a considerable fraction of WN stars. Among the alternative scenarios, close-binary evolution can theoretically produce the least-luminous WN stars. However, final conclusions about the evolutionary formation of the WR stars are not yet possible.     

The primordial Helium-4 abundance determination: systematic effects

By extrapolating to O/H = N/H = 0 the empirical correlations Y–O/H and Y–N/H defined by a relatively large sample of 45 Blue Compact Dwarfs (BCDs), we have obtained a primordial ⁴Helium mass fraction Y _p=0.2443±0.0015 with dY/dZ=2.4±1.0. This result is in excellent agreement with the average Y _p=0.2452±0.0015 determined in the two most metal-deficient BCDs known, I Zw 18 (Z /50) and SBS 0335–052 (Z /41), where the correction for He production is smallest. The quoted error (1) of 1% is statistical and does not include systematic effects. We examine various systematic effects including collisional excitation of hydrogen lines, ionization structure and temperature fluctuation effects, and underlying stellar Hei absorption, and conclude that combining all systematic effects, our Y _p may be underestimated by 2–4%. Taken at face value, our Y _p implies a baryon-to-photon number ratio =(4.7^+1.0 _–0.8)×10^–10 and a baryon mass fraction _b h ² ₁₀₀=0.017±0.005 (2), consistent with the values obtained from deuterium and Cosmic Microwave Background measurements. Correcting Y _p upward by 2–4% would make the agreement even better.     

Characteristics of nearby interstellar matter

There is a warm tenuous partially ionized cloud (T10⁴ K,n(HI)0.1 cm^–3,n(Hii 0.22–0.44 cm^–3) surrounding the solar system which regulates the environment of the solar system, determines the structure of the heliopause region, and feeds neutral interstellar gas into the inner solar system. The velocity (V–20 km s^–1 froml335°,b0° in the local standard of rest) and enhanced Caii and Feii abundances of this cloud suggest an origin as evaporated gas from cloud surfaces in the Scorpius-Centaurus Association. Although the soft X-ray emission attributed to the Local Bubble is enigmatic, optical and ultraviolet data are consistent with bubble formation caused by star formation epochs in the Scorpius-Centaurus Association as regulated by the nearby spiral arm configuration. The cloud surrounding the solar system (the local fluff) appears to be the leading region of an expanding interstellar structure (the squall line) which contains a magnetic field causing polarization of the light of nearby stars, and also absorption features in nearby upwind stars. The velocity vectors of the solar system and local fluff are perpendicular in the local standard of rest. Combining this information with the low column densities seen towards Sirius in the anti-apex direction, and the assumption that the cloud velocity vector is parallel to the surface normal, suggests that the Sun entered the local fluff within the historical past (less than 10 000 years ago) and is skimming the surface of the cloud. Comparison of magnesium absorption lines towards Sirius and anomalous cosmic-ray data suggest the local fluff is in ionization equilibrium.Reason has moons, but moons not hers, Lie mirror'd on her sea, Confounding her astronomers, But, O! delighting me.Ralph Hodgson     

The atmospheric composition,extinction and luminosity of the LBV R71

ESO 3.6m Caspec spectra of the LMC luminous blue variable (LBV) taken at minimum have been analysed using NLTE model atmospheres and line formation calculations to derive atmospheric parameters and chemical composition. Using the silicon ionization balance and the hydrogen Balmer lines we deriveT _eff =17250, log g=1.80 and a microturbulent velocity of 15–20 km/s. The analysis yields abundance ratios by number of approximately 0.43 for He/H, 0.03 for C/N and 0.14 for O/N, implying that enrichment of the atmosphere by processed material has taken place. We have re-evaluated the reddening of R71 using IUE low resolution data and published UBVRIJHKL photometry and derive a value for A _V of 0.63. We also construct an extinction curve using archive IUE data for mid-B LMC supergiants and show that the extinction is anomalous; the 2175A bump being almost absent and the far UV rise very pronounced. A comparison of our model flux in theV-band with the observed (dereddened)V magnitude and the D.M. of the LMC (18.45), implies that the bolometric magnitude or R71 is –9.9. This is significantly higher than the value of –9.0 usually adopted for R71 and suggests that this object may not in fact be a subluminous LBV.     

A New Mechanism of Coronal Heating

The interaction between network magnetic fields and emerging intranetwork fields may lead to magnetic reconnection and microflares, which generate fast shocks with an Alfvén Mach number M _A<2. Protons and less abundant ions in the solar corona are then heated and accelerated by fast shocks. Our study of shock heating shows that (a) the nearly nondeflection of ion motion across the shock ramp leads to a large perpendicular thermal velocity (v _th), which is an increasing function of the mass/charge ratio; (b) the heating by subcritical shocks with 1.1 M_A 1.5 leads to a large temperature anisotropy with T/T 50 for O⁵⁺ ions and a mild anisotropy with T_/T 1.2 for protons; (c) the large perpendicular thermal velocity of He⁺⁺ and O⁵⁺ ions can be converted to the radial outflow velocity (u) in the divergent coronal field lines; and (d) the heating and acceleration by shocks with 1.1 M_A 1.5 can lead to u(O⁵⁺) v _th(O⁵⁺) 460 km s^–1 for O⁵⁺ ions, u(He⁺⁺) v _th(He⁺⁺) 360 km s^–1 for He⁺⁺ ions, and u(H⁺) v _th(H⁺) 240 km s^–1 for protons at r=3–4 R . Our results can explain recent SOHO observations of the heating and acceleration of protons and heavier ions in the solar corona.     

Cosmic-ray Composition as Observed by Ulysses

The cosmic ray isotopic composition measurements from the High Energy Telescope (HET) on the Ulysses spacecraft are reviewed. The source isotopic composition of key elements is found to be surprisingly like the Solar system abundances with the notable exception of ²²Ne. The average density of interstellar material cosmic rays traverse is found to be 0.25 atom cm^–3, corresponding to a confinement time of 20 Myr. Vanadium isotopic abundances are shown to be consistent with weak cosmic-ray reacceleration. The implications of these measurements are discussed.     

The evolution of massive stars to explosion

We review the possible evolutionary paths from massive stars to explosive endpoints as various types of supernovae associated with Population I and hence with massive stars: Type II-P, Type II-L, Type Ib, Type Ic, and the hybrid events SN 1987K and SN 1993J. We identify SN 1954A as another hybrid event from the evidence for both H and He in its spectrum with velocities nearly the same as SN 1983J. Evidence for ejected⁵⁶Ni mass of 0.07 M suggests that SN II-P underwent standard iron core collapse, not collapse of an O–Ne–Mg core nor thermonuclear explosion of a C–O core. Most SN II-P presumably arise in single stars or wide binaries of 10–20 M. There may be indirect evidence for duplicity in some cases in the form of strong Ba II lines, such as characterized SN 1987A. SN II-L are recognizably distinct from typical SN II-P and must undergo a significantly different evolution. Despite indications that SN II-L have small envelopes that may be helium enriched, they are also distinct from events like SN 1993J that must have yet again a different evolution. The SN II-L that share a common Luminosity seem to have ejected a small nickel mass and hence may come from stars with O–Ne–Mg cores. The amount of nickel ejected by the exceptionally bright events, SN 1980K and SN 1979C, remains controversial. SN Ib require the complete loss of the H envelope, either to a binary companion or to a wind. The few identified have relatively large ejecta masses. It is not clear what evolutionary processes distinguish SN Ib's evolving in binary systems from hybrid events that retain some H in the envelope. SN Ic events are both H and He deficient. Binary models that can account for transfer of an extended helium envelope from low mass helium cores, 2 to 4 M, imply C–O core masses that are roughly consistent with that deduced from the ejecta mass plus a neutron star, 2 to 3 M. It is possible that the hybrid events are the result of Roche lobe overflow and that the pure events, SN Ib or SN Ic, result from common envelope evolution.     

Inside and outside of supernovae

A newly formed neutron star in a supernova finds itself in a dense environment, in which the gravitational energy of accreting matter can be lost to neutrinos. For the conditions in SN 1987A, 0.1M may have fallen back onto the central neutron star on a timescale of hours after the explosion, after which the accretion rate is expected to drop sharply. Radiation is trapped in the flow until the mass accretion rate drops to 2×10^–4 M yr^–1 at which point radiation can begin to escape from the shocked envelope at an Eddington limit luminosity. Between this neutrino limit and the Eddington limit, 3×10^–8 M yr^–1, there are no steady, spherical solutions for neutron star accretion. SN 1987A should have reached the neutrino limit within a year of the explosion; the current lack of an Eddington luminosity can be attributed to black hole formation or to a clearing of the neutron star envelope. There is no evidence for newly formed neutron stars in supernovae. Radio supernovae, which were initially interpreted as pulsar activity, probably involve circumstellar interaction; SN 1993J shows especially good evidence for outer shock phenomena.     

The Origin of Primary Cosmic Rays: Constraints from ACE Elemental and Isotopic Composition Observations

Cosmic-ray isotope observations from NASAs Advanced Composition Explorer (ACE) mission have been used to investigate the composition of cosmic-ray source material. Source abundances relative to ⁵⁶Fe are reported for eleven isotopes of Ca, Fe, Co, and Ni, including the very rare isotopes ⁴⁸Ca and ⁶⁴Ni. Although the source abundances range over a factor 10⁴, most of the ratios to ⁵⁶Fe are consistent with solar-system values to within 20%. However, there are some notable differences, the most significant being an excess of (70±30)% relative to the solar system for the cosmic-ray source ratio ⁵⁸Fe/⁵⁶Fe. The possible association of such an excess with a contribution to the cosmic-ray source from Wolf–Rayet star ejecta is discussed.     

Electric gaps in pulsar magnetospheres

Two contrasting models of the magnetosphere are studied: a classical, quasi-mhd model, in which it is the componentE that accelerates electrons to relativistic energies; and a model with e⁺-e^– production, in which the componentE is locally crucial, both in an acceleration domain near the star and in a dissipation domain beyond the light-cylinder.     

β Cephei stars from a photometric point of view

This is an observational review, with an emphasis on photometric data and their interpretation. Two lists are presented, one containing Cephei stars, and the other, Cephei suspects. These lists then serve as a basis for discussing such topics as the location of Cephei stars in the observational and theoretical H-R diagrams, the evolutionary state of these stars, the period-luminosity and period-luminosity-color relations, and observational identification of pulsation modes. The paper also includes references to recent work connected with the theoretical discovery that an opacity mechanism is responsible for the excitation of Cephei-star pulsations. Finally, observational programs for verifying the consequences of this discovery are suggested.Belgian Fund for Scientific Research (NFWO).     

The x-ray cluster baryon crisis

Nucleosynthesis in the standard hot big bang cosmology offers a successful account of the production of the light nuclides during the early evolution of the Universe. Consistency among the predicted and observed abundances of D,³He,⁴He and⁷Li leads to restrictive lower and upper bounds to the present density of nucleons. In particular, the upper bound ensures that nucleons cannot account for more than a small fraction (<0.06h ₅₀ ^–2 ) of the mass in a critical density (Einstein-de Sitter) Universe. In contrast, x-ray observations of rich clusters of galaxies suggest strongly that baryons (in hot gas) contribute a significant fraction of the total cluster mass (0.2h ₅₀ ^–3/2 ). If, indeed, clusters do provide a fair sample of the mass in the Universe, this crisis forces us to consider other ways of mitigating it, including the politically incorrect possibility that <1. The options, including magnetic or turbulent pressure, clumping and non-zero space curvature and/or cosmological constant, are discussed.     

Starburst galaxies

Bursts of massive star formation can temporarily dominate the luminosity of galaxies spanning a wide range of morphological types. This review is concerned primarily with such events in the central 1 kpc region of spiral galaxies which result from bar driven inflows of gas triggered by interactions or mergers. Most of the stellar radiant luminosity of such bursts is absorbed by dust and re-emitted in the far-infrared and is accompanied by radio and X-ray emission from supernova remnants which can also act collectively to drive galaxy scale outflows. Both evolutionary stellar models and estimates of the gas depletion times are consistent with typical burst durations of 10^7–8 yr. Spatially-resolved studies of nearby starburst galaxies reveal that the activity is distributed over many individual star forming complexes within rings and other structures organized by interactions between bars and the disc over a range of scales. More distant and extreme examples associated with mergers of massive spirals have luminosities > 10¹³ L and molecular gas masses > 10¹⁰ M implying star formation rates > 1000 M yr^–1 which can only be sustained for 10⁷ yr. In the most luminous merging systems, however, the relative importance of starburst and AGN activity and their possible evolutionary connection is still a hotly debated issue. Also controversial are suggestions that starbursts in addition to a black hole are required to account for the properties of AGNs or that starbursts alone may be sufficient under certain conditions. In a wider context, starbursts must clearly have played an important role in galaxy formation and evolution at earlier times. Recent detections of high redshift galaxies show that star formation was underway at z 4 but do not support a continuing increase of the strong evolution in the co-moving star formation density seen out to z l. Primeval starburst pre-cursors of spheroidal systems also remain elusive. The most distant candidates are radio galaxies and quasars at z = 4–5 and a possible population of objects responsible for an isotropic sub-mm wave background tentatively claimed to have been detected by the COBE satellite.     

The fundamental parameters of the central stars of eight WR ring nebulae

We describe work that has recently been completed on deriving the fundamental parameters of eight WR stars through the photoionization modelling of their surrounding nebulae using non-LTE WR flux distributions. The resulting effective temperatures range from 57 000–71 000 K for the WN4-5 stars and <30 000–42 000 K for the WN6-8 stars. The derived stellar parameters are compared with those obtained from stellar emission line modelling. We find good agreement for the hot early WN stars, indicating that the non-LTE WR flux distributions have essentially the correct shape in the crucial far-UV region. We find lower temperatures for the four cooler late WN stars, particularly for the two WN6 stars. For the nebulae surrounding these stars, we find that the model flux distributions produce too much nebular ionization. We suggest that these discrepancies arise because of the lack of line-blanketing in the WR atmospheres. For the WO1 central star of G2.4+1.4, with strong nebular He II 4686 A emission, we derive a temperature of 105 000 K, somewhat less than previous estimates. The positions of our eight WR stars on the H-R diagram are compared with the evolutionary tracks of Maeder (1990) for solar metallicity. In common with previous workers, we find that our derived luminosities are too low, giving an initial mass range of 25–40 M, below that expected for the majority of WR stars.     

Large scale structure of the universe

The current situation with the cosmological model and fundamental constants is briefly reviewed. Here, we concentrate on evolutionary effects of large-scale structure formation, in particular, the relationship with the quasar distribution and dynamics is discussed. We argue that groups of bright quasars with few or more than dozen of members within regions l _LS(100–150)h ^–1 Mpc found atz<2 may belong to concentrations of young rich clusters of galaxies, and thus be distant Great Attractors like the local GA or the Shapley concentration. These early large-scale galactic structures (i) provide a natural way to bias the distribution of Abell clusters, and (ii) suggest that the spectrum of primordial density perturbations is nearly flat on scales encompassing both the cluster and GAs,l=k ^–1(10,100)h ^–1 Mpc: _k ² k ³ P(k) k , =1 _–0.4 ^+0.6 , whereP(k) is the power spectrum of density perturbations.     

Latitudinal variation of solar wind velocity

Single station solar wind velocity measurements using the Ooty Radio Telescope (ORT) in India (operating at 327 MHz) are reported for the period August 1992 to August 1993. Interplanetary scintillation (IPS) observations on a large number of compact radio sources covering a latitudinal range of ±80° were used to derive solar wind velocities using the method of fitting a power law model to the observed IPS spectra. The data shows a velocity versus heliographic latitude pattern which is similar to that reported by Rickett and Coles (1991) for the 1981–1982 period. However, the average of the measured equatorial velocities are higher, being about 470 km s^–1 compared to their value of 400 km s^–1. The distribution of electron density variations (N _e ) between 50R and 90R was also determined and it was found that N _e was about 30% less at the poles as compared to the equator.     

Spectral observation of the composite supernova remnant G 29.7-0.3

The X-ray properties of the supernova remnant G 29.7-0.3 are discussed based on spectral data from the EXOSAT satellite. In the 2 to 10 keV range a featureless power-law spectrum is obtained, the best-fit parameters being: energy spectral index =-0.77, hydrogen column density on the line of sight N_H=2.3.10²² cm^–2. The incident X-ray flux from the source is (3.6±0.1) 10¹¹ erg cm^–2 s^–1 in the 2 to 10 keV range corresponding to an intrinsic luminosity of about 2. 10³⁶ erg s^–1 for a distance of 19 kpc. The source was not seen with the imaging instrument thus constraining the hydrogen column density to be N_H=(3.3 ±0.3) 10²² cm^–2 and the energy spectral index =1.0±0.15. This new observation is consistent with emission by a synchroton nebula presumably fed by an active pulsar. An upper limit of 1.5% for the pulsed fraction in the range of periods 32ms to 10⁴ s has been obtained.     

The synthesis of helium and CNO isotopes in massive stars

We present helium and CNO isotopic yields for massive mass-losing stars in the initial mass range 15M M _i 50M . We investigate their dependence on assumptions about mass loss rates, internal mixing processes, and metallicity, and specify the contributions from stellar winds and from supernova ejecta.     

Mass distribution in our Galaxy

This article summarizes recent work on the luminosity and mass distribution of the galactic bulge and disk, and on the mass of the Milky Way's dark halo. A new luminosity model consistent with the COBE NIR data and the apparent magnitude distributions of bulge clump giant stars has bulge/bar length of 3.5 kpc, axis ratios of 1:(0.3–0.4):0.3, and short disk scale-length (2.1 kpc). Gas-dynamical flows in the potential of this model with constant M/L fit the terminal velocities in 10° le|l|le50° very well. The luminous mass distribution with this M/L is consistent with the surface density of known matter near the Sun, but still underpredicts the microlensing optical depth towards the bulge. Together, these facts argue strongly for a massive, near-maximal disk in our L ^*, Sbc spiral galaxy. While the outer rotation curve and global mass distribution are not as readily measured as in similar spiral galaxies, the dark halo mass estimated from satellite velocities is consistent with a flat rotation curve continuing on from the luminous mass distribution.