Characteristics and interpretation of the photometric variability of eta carinae and its nebula

A short review is given on the history of the peculiar variable object Car and on a number of relevant references describing and discussing its physical characteristics and behaviour, based on different types of observational techniques. The star is known to be variable since the 17th century. The excessive mass loss to which it was subject during the 19th century is now visible as an ellipsoidal reflection nebula of 15 diameter: the so-called homunculus. The remainder of the paper is spent on different kinds of problems partly based on the results of a decade of photometric monitoring in the VBLUW photometric system of Walraven. Foreground reddening and reddening by dust in the homunculus are determined and amount to E(B - V) _J = 0 50 and < 6, respectively. Scanning of the homunculus revealed an estimate for the photometric characteristics of the central object, which presumably consists of a massive hot star surrounded by a cooler gas envelope. The total luminosity is derived using fluxes of various sources in the wavelength region 0.15 < < 175 n resulting in M _bol = - 12 3 ± 0 2. The total observed flux corrected for foreground extinction corresponds to a star with R 96 R if T _eff 30 000 K. The mass may be near 150 M . The excess luminosity in 1843, when the star was presumably bolometrically at least 2 5 brighter than at present, may have been caused by envelope-energized pulsations when the star's luminosity was close to its Eddington limit. The temperature should then have been 50 000 K. The mass loss rate, during the excess luminosity phase lasting 30 yr, is estimated to amount to M 4 × 10^-3 M yr^-1. At present the mass loss may be M 10^-4M yr^-1. Since the homunculus is mainly built up from material expelled in the 30 yr interval (from 1830 to 1860), its total mass amounts to M _hom 0.15 M . The historical observations of the colours of Car and a comparison with the characteristics of S Dor type stars, suggest that Car was subject to a number of S Dor type phases similar to those of P Cyg (in the 17th century), S Dor and others. A satisfactory explanation is found for the complete historical light curve. The decrease in light after the 1843 maximum by 9 ^m , presumably consists of a fading of the luminosity excess and the S Dor effect by 2 5 and 3^m, respectively, and a 3 5 extinction by circumstellar dust. The small amplitude light variations which Car showed during the last decade, were studied with the aid of the variations of the Balmer jump. They are presumably caused by temperature variations of the central star.Percy and Welch (1983) (Publ. Astron. Soc. Pacific 95, 491) have observed P Cyg on a number of nights in 1982 and found for the photometric variations a time scale of 30 to 50 days and an amplitude of 0 _. ^m 15.Based partly on observations collected at the ESO, La Silla, Chile.     

Nonlinear theory of the two-point correlation function

We consider the influence of the nonlinear stage of gravitational instability on the two-point correlation functions of gravitationally bound objects. Based on the theory of nonlinear gravitational contraction of a single density peak of dissipationless matter (Gurevich and Zybin, 1988a,b; 1990) we develop a method for calculating the two-point correlation functions of different objects of any mass. The method works good in the region of strong correlations and can be easily extended to calculate higher correlation functions. We show that the main contribution to the correlation function _i in the region of strong correlations _i 1 is made by pair systems located outside large clusters of objects. In this region the shape of _i is determined only by the nonlinear dynamics of gravitational contraction of dissipationless matter and has the form _i C ^–, where 1.8 is a universal parameter.     

A dissection of 30 Doradus and a discussion of other Magellanic Cloud OB associations

I Present the results of ground-based and Hubble Space Telescope photometry and spectroscopy of the stars in the central region (roughly 7×7 arcmin) of 30 Doradus in the Large Magellanic Cloud (LMC). Using photometric data for over 2400 stars (complete toV18 mag), and spectroscopic observations of over 150 stars in the region, the best estimate of the initial mass function (IMF) yields a slope of =–1.5±0.2 for masses > 12M, where the Salpeter slope is =–1.35. I compare these results to other measurements of the IMF for OB associations in the Magellanic Clouds.     

Some results of ionospheric investigations by means of coherent radiowaves emitted by satellites

In this paper we discuss theoretical expressions, determining the difference of Doppler shifts of various coherent radiowave frequencies emitted by a radiator moving in the ionosphere or interplanetary medium. The rotating Doppler effect (Faraday effect) caused by the Doppler shifts ±H of the ordinary and extraordinary waves is also considered. In a three-dimensional inhomogeneous ionosphere, stationary in time (N/t = 0), is determined in the general case, by an equation with three variables. The equation for proper depends only on the local value of the electron concentration N _c around the radiator and on integral values, determining, by means of additional calculations, the angle of refraction or its components, the horizontal gradients of electron concentration N/x and N/y, and in some cases, the integral electron concentration ₀ ^zcN dz. We describe the analysis of the measurements, made with the satellites Cosmos I, II and partially XI, assuming that N/t = N/y = 0, with a two variables equation. The expected errors are considered. The results coincide well for different points (Moscow, The Crimea, Sverdlovsk) and thus agree with the measurements of _H and with height-frequency ionospheric characteristics. The curve giving electron concentration versus height N (z) in the outer ionosphere (above the maximum of F2), shows a new maximum higher than the main maximum of the ionosphere N _MF2 at 120–140 km. At this maximum the value of N (z) is (0.9–0.95) N _MF2. The new data on the large-scale horizontal inhomogeneities of the ionosphere, exceed the previous ones by about a factor 10. By means of the irregular variations of the spectrum W() of the inhomogenous formation is determined. Three unknown constant maxima with values 16 to 18 km, 28 to 32 km and 100 to 120 km are found. The spectrum W () mainly characterizes the local properties of the ionosphere along the orbit of the satellite.     

Fourier parameters of heliospheric current sheet and their significance

If the path of the neutral line on the coronal source surface is expressible as a singlevalued function (colatitude vs longitude ), then Fourier analysis of ctn with respect to leads to a simple algorithm for realistically mapping the neutral line outward to model the heliospheric current sheet (HCS) at distancesr1 AU. To be compatible with MHD, the source surface used for this mapping should be prolate (aligned with dipole axis) rather than spherical. Orientation of the Sun's magnetic-dipole moment is indicated by them=1 Fourier amplitude (a ₁ sin +b ₁ cos ) of ctn on the source surface. Physical features (including the neutral line) on a prolate source surface intrinsically map to lower dipole latitudes atr1 AU in the heliosphere, and Ulysses observations of a unipolar field at latitudes beyond 30°S (when the neutral line on the source surface still reached 39°S) confirm the expected geometry.     

Pulsars as gamma ray sources: Nebular shocks and magnetospheric gaps

I summarize the results of recent research on the structure and particle acceleration properties of relativistic shock waves in which the magnetic field is transverse to the flow direction in the upstream medium, and whose composition is primarily electrons and positrons with an admixture of heavy ions. Shocks which contain heavy ions that are a minority constituent by number but which carry most of the energy density in the upstream medium put 20% of the flow energy into a nonthermal population of pairs downstream, whose distribution in energy space is N(E) E ^-2, where N(E)dE is the number of particles with energy between E and E+dE. Synchrotron maser activity in the shock front, stimulated by the quasi-coherent gyration of the whole particle population as the plasma flowing into the shock reflects from the magnetic field in the shock front, provides the mechanism of thermalization and non-thermal particle acceleration. The maximum energy achievable by the pairs is _± m _± c ² = m _i c ² ₁/Z _i, where ₁ is the Lorentz factor of the upstream flow and Z _i is the atomic number of the ions. The shock's spatial structure contains a series of overshoots in the magnetic field, regions where the gyrating heavy ions compress the magnetic field to levels in excess of the eventual downstream value. These overshoots provide a new interpretation of the structure of the inner regions of the Crab Nebula, in particular of the wisps, surface brightness enhancements near the pulsar. The wisps appear brighter because the small Larmor radius pairs are compressed and radiate more efficiently in the regions of more intense magnetic field. This interpretation suggests that the structure of the shock terminating the pulsar's wind in the Crab Nebula is spatially resolved, and allows one to measure ₁ 4 × 10⁶, the upstream magnetic field B ₁ to be 3 × 10^-5 Gauss, as well as to show that the total ion flow is 3 × 10³⁴ elementary charges/sec, in good agreement with the total current flow predicted by the early Goldreich and Julian (1969) model. The total pair outflow is shown to be about 5 × 10³⁷ pairs per second, in good agreement with the particle flux required to explain the nebular X—ray source.The energetics of particle acceleration within the magnetospheres of rotation powered pulsars and the consequences for pulsed gamma ray emission are also briefly discussed. The gamma ray luminosity above 100 MeV is shown to scale in proportion to _R ^1/2 , as is in accord with some of the simplest ideas about polar cap models. Models based on acceleration in the outer magnetosphere are also briefly discussed.     

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.     

NGC 4151 and MCG 8-11-11: Two X-ray seyfert galaxies with strong soft γ-ray emission

During a balloon flight of the MISO telescope on the 30th September 1979, the Seyfert galaxies NGC 4151 and MGC 8-11-11 were studied in the hard X-ray range (E_X > 20 keV) and low-energy -ray range up to 19 MeV. An emission at the 4.5 level above 20 keV (4 above 260 keV) was detected in the direction of NGC 4151. -ray emission at the 3.9 level above 90 keV was also observed from the direction of MCG 8-11-11. The emission photon spectrum shows a high-energy cutoff at about 3 MeV. A large amount of the observed low-energy -ray diffuse background could be produced by a few percent of the X-ray emitting Seyfert galaxies having a -ray luminosity comparable to that observed from the regions of NGC 4151 or MCG 8-11-11.     

Determination of the characteristics of the gamma-ray telescope Gamma-1

The Gamma-1 telescope has been developed through a collaboration of scientists in the USSR and France in order to conduct -ray astronomical observations within the energy range from 50 to 5000 MeV. The major characteristics of the telescope were established by Monte-Carlo simulations and calibrations made with the aid of electron and tagged -ray beams produced by an accelerator, and these have been found to be as follows: the effective area for photons coming along the instrument's axis varies from about 50 cm² at E = 50 MeV to approximately 230 cm² at E 300 MeV; the angular resolution (half opening of the cone embracing 68% events) is equal to 2.7° at E = 100 MeV, and 1.8° at E = 300 MeV; the energy resolution (FWHM) varies from 70% to 35% as the energy of the detected photons increases from 100 to 550 MeV; the telescope's field-of-view at the half-sensitivity level is 300–450 square degrees depending upon the spectrum of the detected radiation, and the event selection logic. Proceeding from the thus obtained characteristics it is demonstrated that a point source producing a photon flux J (E 100 MeV) = 3 × 10^-7 cm^-2 s^-1, can be detected with a 5 significance by observing it during 10⁶ s at the level of the Cygnus background, and a source having intensity J (E 100 MeV) = 10^-6 cm^-2 s^-1 can be detected to within a mean square positional accuracy of about 15.     

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.     

Electron coronal density irregularity\overline {n^2 } /(\bar n)^2 from measurements of K-coronal and lyman lines brightnesses

We propose a technique to derive the coronal density irregularity factor , wheren is the electron density. The absolute photometric comparison between the intensity of UV lines and the white-light K-coronal polarized brightness (pB) provides an unique constraint on the inhomogeneity of the corona. The ratio of the measured H I Lyman (Ly-) line intensity to the resonant-scattering dominated H I Lyman (Ly-) intensity can be used to extract the collisonal component of the Ly-. This component yields an estimate of . The quantity is then obtained from white-light K-coronal measurements. The use of lines of the same atomic species minimizes the effects due to outflow velocities (i.e., Doppler dimming), and reduces the errors introduced by the uncertainties in the ionization balance, the atomic parameters, and the solar abundances. The UVCS/SOHO unique capability of performing cotemporal and cospatial measurements of the Ly- and Ly- lines, and ofpB makes this instrument ideal for implementing this technique.     

The importance of anisotropic interstellar turbulence and molecular-cloud magnetic mirrors for galactic cosmic-ray propagation

Recent studies suggest that when magnetohydrodynamic (MHD) turbulence is excited by stirring a plasma at large scales, the cascade of energy from large to small scales is anisotropic, in the sense that small-scale fluctuations satisfy the inequality k k , where k and k are, respectively, the components of a fluctuations wave vector and to the background magnetic field. Such anisotropic fluctuations are very inefficient at scattering cosmic rays. Results based on the quasilinear approximation for scattering of cosmic rays by anisotropic MHD turbulence are presented and explained. The important role played by molecular-cloud magnetic mirrors in confining and isotropizing cosmic rays when scattering is weak is also discussed.     

Some problems of gamma-astronomy

The modern state of gamma-ray astronomy is reviewed, the paper being mainly devoted to the theoretical models that describe generation of gamma-ray emission under astrophysical conditions. Basic information on the processes of generation and absorption of gamma-rays, as well as the results of observations for various gamma-ray photon energies are reported.In the region of soft gamma-ray emission (i.e., for energies less than tens of MeV), where emission in gamma-ray lines dominates, we also discuss the nature of gamma-ray bursts, the origin of gamma-ray emission from the galactic centre, etc.Discrete sources and, in particular, the mysterious source Cyg X-3 are discussed in the region of very high (E > 10¹² eV) and ultra-high (E > 10¹⁵ eV) energy gamma-ray emission.A larger portion of the review is devoted to the analysis of cosmic-ray origin on the basis of the available gamma-ray data in the region from several tens of MeV to several GeV. The peculiarity of this energy range is, in particular, in the fact that the diffuse galactic emission was observed mainly there. We also discuss the problem of determination of the cosmic-ray density gradient from the gamma-ray data.The origin of high-latitude gamma-ray emission, the problem of galactic gamma-ray halo, etc., are discussed.The theoretical models explaining the nature of unidentified gamma-ray sources, as well as the results of measurements and theoretical estimations of a gamma-ray flux from SN1987A are analysed.List of Notations m electron mass, m = 9.108 × 10^–28 g, - M proton mass, M = 1.672 × 10^–24 g, - e electron charge, e = 4.803 × 10^–10 CGS - c velocity of light, c = 2.9979 × 10¹⁰ cm s^–1, - k Boltzmann constant, k = 1.380 × 10^–16 erg grad^–1, - e electron - p proton - gamma-ray photon - p antiproton - ^{0 0}-meson - -lepton - e ⁺ positron - r, , x radio-frequency, gamma-ray, and X-ray emission bands - E total energy of a particle - E _k kinetic energy - p particle momentum - spectral index for particles - spectral index for emission - n particle density (concentration) - H magnetic field strength - T temperature - _ph energy of low-energy photons - emission frequency - r _H Larmor radius of relativistic particles - k wave number - , z cylindric coordinates, in this case the coordinate (radius) along the galactic disk, z perpendicular to the galactic disk - M solar mass, M = 1.99 × 10³³ g.     

Tests for coalescing binary neutron stars as cosmological origin of gamma-ray bursts

Theoretical logN-logS distributions and (V/V _max) tests of gamma-ray bursts in the model of coalescence of neutron star (NS+NS) and/or NS+black hole (NS+BH) binaries are calculated for a flat Universe (=1) with different values of the cosmological constant and under various assumptions about the star formation history. The observed logN-logS distribution and value of (V/V _max)=0.33 for 411 bursts with knownC _max/C _lim from the 2d BATSE catalogue are best fitted with a model for which = 0.2 and primary star formation occurs at redshiftsz5–6.     

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.     

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     

Voyager energetic particle observations at interplanetary shocks and upstream of planetary bow shocks: 1977–1990

The Voyager 1 and 2 spacecraft include instrumentation that makes comprehensive ion (E 28 keV) and electron (E 22 keV) measurements in several energy channels with good temporal, energy, and compositional resolution. Data collected over the past decade (1977–1988), including observations upstream and downstream of four planetary bow shocks (Earth, Jupiter, Saturn, Uranus) and numerous interplanetary shocks to 30 AU, are reviewed and analyzed in the context of the Fermi and shock drift acceleration (SDA) models. Principal findings upstream of planetary bow shocks include the simultaneous presence of ions and electrons, detection of tracer ions characteristic of the parent magnetosphere (O, S, O⁺), power-law energy spectra extending to 5 MeV, and large (up to 100:1) anisotropies. Results from interplanetary shocks include observation of acceleration to the highest energies ever seen in a shock ( 22 MeV for protons, 220 MeV for oxygen), the saturation in energy gain to 300 keV at quasi-parallel shocks, the observation of shock-accelerated relativistic electrons, and separation of high-energy (upstream) from low-energy (downstream) populations to within 1 particle gyroradius in a near-perpendicular shock. The overall results suggest that ions and electrons observed upstream of planetary bow shocks have their source inside the parent magnetosphere, with first order Fermi acceleration playing a secondary role at best. Further, that quasi-perpendicular interplanetary shocks accelerate ions and electrons most efficiently to high energies through the shock-drift process. These findings suggest that great care must be exercised in the application of concepts developed for heliosphere shocks to cosmic ray acceleration through shocks at supernova remnants.     

Near Earth Current Meander (Necm) Model of Substorms

We propose that the appropriate instability to trigger a substorm is a tailward meander (in the equatorial plane) of the strong current filament that develops during the growth phase. From this single assumption follows the entire sequence of events for a substorm. The main particle acceleration mechanism in the plasma sheet is curvature drift with a dawn-dusk electric field, leading to the production of auroral arcs. Eventually the curvature becomes so high that the ions cannot negotiate the sharp turn at the field-reversal region, locally, at a certain time. The particle motion becomes chaotic, causing a local outward meander of the cross-tail current. An induction electric field is produced by Lenz's law, E ^ind=–A/t. An outward meander with B _z>0 will cause E×B flow everywhere out from the disturbance; this reaction is a macroscopic instability which we designate the electromotive instability. The response of the plasma is through charge separation and a scalar potential, E ^es=–. Both types of electric fields have components parallel to B in a realistic magnetic field. For MHD theory to hold the net E must be small; this usually seems to happen (because MHD often does hold), but not always. Part of the response is the formation of field-aligned currents producing the well-known substorm current diversion. This is a direct result of a strong E ^ind (the cause) needed to overcome the mirror force of the current carriers; this enables charge separation to produce an opposing electrostatic field E ^es (the effect). Satellite data confirm the reality of a strong E in the plasma sheet by counter-streaming of electrons and ions, and by the inverse ion time dispersion, up to several 100 keV. The electron precipitation is associated with the westward traveling surge (WTS) and the ion with omega () bands, respectively. However, with zero curl, E ^es cannot modify the emf =Edl=–d^M/dt of the inductive electric field E ^ind (a property of vector fields); the charge separation that produces a reduction of E must enhance the transverse component E . The new plasma flow becomes a switch for access to the free energy of the stressed magnetotail. On the tailward side the dusk-dawn electric field with EJ<0 will cause tailward motion of the plasma and a plasmoid may be created; it will move in the direction of least magnetic pressure, tailward. On the earthward side the enhanced dawn-dusk induction electric field with EJ>0 will cause injection into the inner plasma sheet, repeatedly observed at moderate energies of 1–50 keV. This same electric field near the emerging X-line will accelerate particles non-adiabatically to moderate energies. With high magnetic moments in a weak magnetic field, electrons (ions) can benefit from gradient and curvature drift to attain high energies (by the ratio of the magnetic field magnitude) in seconds (minutes).     

Non-BBN Constraints on the Key Cosmological Parameters

Since the baryon-to-photon ratio 10 is in some doubt at present, we ignore the constraints on 10 from big bang nucleosynthesis (BBN) and fit the three key cosmological parameters (h, M, 10) to four other observational constraints: Hubble parameter (ho), age of the universe (to), cluster gas (baryon) fraction (fo fGh3/2), and effective shape parameter (o). We consider open and flat CDM models and flat CDM models, testing goodness of fit and drawing confidence regions by the 2 method. CDM models with M = 1 (SCDM models) are accepted only because we allow a large error on ho, permitting h < 0.5. Open CDM models are accepted only for _M 0.4. CDM models give similar results. In all of these models, large 10 ( 6) is favored strongly over small 10 ( 2), supporting reports of low deuterium abundances on some QSO lines of sight, and suggesting that observational determinations of primordial 4He may be contaminated by systematic errors. Only if we drop the crucial o constraint are much lower values of M and 10 permitted.     

β 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).