An expanding circumstellar cloud of Zeta Aurigae

Strong absorption satellite lines of CaI 6572 were found on spectrograms taken on three successive days just after the fourth contact of the 1971–72 eclipse of Zeta Aurigae. The radial velocities of the satellite lines are –88 km s^–1, –74 km s^–1, and –180 km^–1, respectively, relative to the K-type primary star (K4 Ib). These absorptions should be due to a circumstellar cloud in which the column density of neutral calcium atoms is 1×10¹⁷ cm^–2 and the turbulent velocities come to 20–50 km s^–1. It is suggested that the cloud may be formed by the rocket-effect of the Lyman quanta of the B-type component (B6 V). We estimate the density in the cloud to be 2×10¹¹ atoms cm^–3 fors=10R _K and 2×10¹⁰ atoms cm^–3 fors=10² R _K, wheres denotes the distance of the cloud from the K star andR _K the K star's radius. The mass loss rate of the K-type component is also estimated to be about 10^–7 M yr^–1, assuming that the expansion of the K star occurs isotropically.     

Emission line ratios for C III in the sun

Theoretical electron-density-sensitive C III emission line ratios are presented forR ₁ =I(2s2p ³ P – 2p ^{2 3} P)/I(2s2p ¹ P – 2p ^{2 1} S) =I(1176 Å)/I(1247 Å),R ₂ =I(2s2p ³ P – 2p ^{2 3} P)/I(2s ^{2 1} S – 2s2p ³ P ₁) =I(1176 Å)/I(1908 Å), andR ₃ =I(2s2p ¹ P – 2p ^{2 1} S)/I(2s ^{2 1} S – 2s2p ³ P ₁) =I(1247 Å)/I(1908 Å). These are significantly different from those deduced previously, principally due to the adoption of improved electron impact excitation rates in the present analysis. Electron densities deduced from the present theoretical line ratios, in conjunction with observed values ofR ₁,R ₂, andR ₃ measured from solar spectra obtained by the Naval Research Laboratory's S082B instrument on boardSkylab, are found to be generally compatible. In contrast, previous diagnostic calculations imply electron densities fromR ₁,R ₂, andR ₃ that differ by up to two orders of magnitude. These results provide observational support for the accuracy of the atomic physics adopted in the present calculations, and the methods employed in the derivation of the theoretical line ratios.     

The identification of Feix and Nixi in the solar corona

The levels of the configuration 3s ²3 _p ⁵3d of Fe ix have been experimentally determined from their combinations with 3s3 _p ⁶3d ³ D in the region 300–400 Å. Wavelengths can now be accurately predicted for all transitions within 3s ²3 _p 53^d, and eight of these can be identified with coronal lines from 2042 to 4585 Å. Also, identifications of solar lines from 171 to 245 Å with electric-dipole and magnetic-quadrupole transitions to the ground state, 3s ²3p ^{6 1} S, are confirmed and extended. Solar identifications with corresponding transitions in Ni xi, both within 3s ²3 _p ⁵3d and to the ground state, are proposed on the basis of a short extrapolation.     

Damping constant and turbulence in the solar atmosphere

In the region of the formation of weak and medium-strong lines, the microturbulence increases with height (V _ver=0.7–0.9 km s^-1, V _hor= 1.1–1.5 km s^-1), the macroturbulence decreases (V _ver=1.6–1.4 km s^-1, V _hor= 2.4–1.5 km s^-1), and the total velocity field (vertical component) is depth-independent (1.7 km s^-1). The empirical damping constants for Fe, Ti, Cr, Ni lines are equal 1.3₆, 1.7₆, 1.6₆, 1.6₆, respectively. The correlation length (the Kubo-Anderson process has been used) in the solar photosphere is 520–550 km.     

The photospheric velocity field of Procyon

BUSS observations of the profiles of two well observed spectral lines in the ultraviolet spectrum of CMi (Procyon; F5 IV–V) are analysed with a Fourier transform method in order to determine values of various parameters of the velocity field of the upper photosphere. We find a microturbulent line-of-sight velocity componentL = 0.9 ± 0.4 km s^–1, a macroturbulent velocity componentL _M = 5.3 ± 0.2 km s^–1, and a rotational velocity componentv _R sini=10.0±1.2 km s^–1. In these calculations a single-moded sinusoidal isotropic macroturbulent velocity function was assumed. The result appears to be sensitive to the assumed shape of the macroturbulence function: for an assumed Gaussian shape the observations can be described withv _R sini=4 km s^–1 andL _M = 11.6 ± 2.7 km s^–1. A comparison is made with other results and theoretical predictions.     

Fe XVIII EMISSION LINES IN SOLAR X-RAY SPECTRA

We have calculated intensity ratios for emission lines of Fexviii in the 13–94 Å wavelength range at electron temperatures characteristic of the solar corona, T _e = 2–10 x 10⁶ K. Our model ion includes data for transitions among the 2s ²2p ⁵ , 2s2p ⁶, 2s ²2p ⁴3l, and 2s2p ⁵3l (l = s, p, and d) states. Test calculations which omit the 2s2p ⁵3l levels show that cascades from these are important. We compare our results with observed ratios determined from four solar X-ray instruments, a rocket-borne spectrograph, and spectrometers on the P78–1, OV1–17 and Solar Maximum Mission (SMM) satellites. In addition, we have generated synthetic spectra which we compare directly with flare observations from SMM. Agreement between theory and observation is generally quite good, with differences that are mostly less than 30%, providing limited support for the accuracy of the atomic physics data used in our calculations. However, large discrepancies are found for ratios involving the 2s ²2p ^{5 2}P_3/2- 2s2p ^{6 2}S line at 93.84 Å, which currently remain unexplained. Our analysis indicates that the FeXVIII feature at 15.83 Å is the 2s ²2p ^{5 2}P_3/2 - 2s ²2p ⁴(³P)3s ⁴P_3/2 transition, rather than 2s ²2p ^{5 2}P_3/2 - 2s ²2p ⁴(³P)3s ²P_3/2, as suggested by some authors.     

Some requirements of a colliding comet source of gamma ray bursts

Colliding comets in the Solar System may be an important source of gamma ray bursts. The spherical gamma ray comet cloud required by the results of the Venera Satellites (Mazets and Golenetskii, 1987) and the BATSE detector on the Compton Satellite (Meeganet al., 1992a, b) is neither the Oort Cloud nor the Kuiper Belt. To satisfy observations ofN(>P _max) vsP _max for the maximum gamma ray fluxes,P _max > 10^–5 erg cm^–2 s^–1 (about 30 bursts yr^–1), the comet density,n, should increase asn a ¹ from about 40 to 100 AU wherea is the comet heliocentric distance. The turnover above 100 AU requiresn a ^–1/2 to 200 AU to fit the Venera results andn a ^1/4 to 400 AU to fit the BATSE data. Then the masses of comets in the 3 regions are from: 40–100 AU, about 9 earth masses,m _E; 100–200 AU about 25m _E; and 100–400 AU, about 900m _E. The flux of 10^–5 erg cm^–2 s^–1 corresponds to a luminosity at 100 AU of 3 × 10²⁶ erg s^–1. Two colliding spherical comets at a distance of 100 AU, each with nucleus of radiusR of 5 km, density of 0.5 g cm^–3 and Keplerian velocity 3 km s^–1 have a combined kinetic energy of 3 × 10²⁸ erg, a factor of about 100 greater than required by the burst maximum fluxes that last for one second. Betatron acceleration in the compressed magnetic fields between the colliding comets could accelerate electrons to energies sufficient to produce the observed high energy gamma rays. Many of the additional observed features of gamma ray bursts can be explained by the solar comet collision source.     

Theoretical emission line ratios for Si XIII compared to solar observations

The electron collision excitation rates recently calculated for transitions in Si xiii by Keenan et al. (1987) are used to derive the electron temperature sensitive ratio G(=(f + i)/r and the density sensitive ratio R(=f/i), where i, f, and r are the intercombination (1s ^{2 1} S – 1s2p ³ P _{1, 2}) forbidden (1s ^{2 1} S – 1s2s ³ S), and resonance (1s ^{2 1} S – 1s2p ¹ P), transitions respectively. Also estimated are the values of R in the low-density limit (R ₀) as a function of electron temperature. The theoretical G ratio at the temperature of maximum emissivity for Si xiii, G(T _m) = 0.70, is in much better agreement with the observed G for the 1985, May 5 flare determined by McKenzie et al. (G = 0.60 ± 0.07) than is the earlier calculation of Pradhan, who derived G(T _m) = 0.85. The error in the observed R ₀ ratio is so large that both our result and Pradhan's fall within the acceptable limits of uncertainty and hence one cannot estimate which of the two is the more accurate.     

Molecular cores of the high-latitude cloud MBM 40

Towards the high-latitude cloud MBM 40, we identify 3 dense molecular cores of M0.2–0.5 M, and sizes of 0.2 pc in diameter embedded in the H I cloud of 8 M which is observed to be extended along the northeast–southwest direction. The molecular cloud is located almost perpendicularly to the H I emission. We confirm the previous result of Magnani et al. that MBM 40 is not a site for new star formations. We found a very poor correlation between the H I and the IRAS 100 μm emissions, but the CO (1–0) and 100 μm emissions show a better correlation of W_CO/I₁₀₀=1±0.2 K km s⁻¹ (MJy sr⁻¹)⁻¹. This ratio is larger by a factor of ≥5 than in dense dark clouds, which may indicate that the CO is less depleted in MBM 40 than in dense dark clouds.     

Stellar winds and mass loss of a rotating star

The mass loss to be expected from the corona of a rotating F2-star is calculated. The rotation is supposed to be rigid up to a certain distances, as if it were maintained by a strong magnetic field. Dependent on the values of the rotational velocity the mass loss can increase to 26–40% forv _rot up to 200 km s^–1.     

The masers of E-type methanol in orion KL and SGR B2

Using a simplified model the statistical equilibrium and radiative transfer equations of E-type-CH₃OH are solved for Orion KL and SgrB2. According to our calculation results and the observation data taken by Matsakiset al. (1980) and Morimotoet al. (1985a, b), the physical conditions of both sources are estimated. In theJ ₂-J ₁ E methanol maser region of Orion KL, the density, kinetic temperature, dust temperature, and the fractional abundance are 0.8–2×10⁶ cm^–3, 150, 30–90 K, 0.8–8×10^–6. In the 4_–1-3₀ E and 5_–1-4₀ E methanol maser region of Sgr B2 the correspondance physical conditions above are 10⁴ cm³, 45, 23 K, and 7×10^–7, respectively.     

Extreme Ultraviolet Emission Lines of ca xv in Solar and Laboratory Spectra

New R-matrix calculations of electron impact excitation rates in Caxv are used to derive theoretical electron density diagnostic emission line intensity ratios involving 2s ²2p ²–2s2p ³ transitions, specifically R ₁=I(208.70 Å)/I(200.98 Å), R ₂=I(181.91 Å)/I(200.98 Å), and R ₃=I(215.38 Å)/I(200.98 Å), for a range of electron temperatures (T _e=10^6.4–10^6.8 K) and densities (N _e=10⁹–10¹³ cm^–3) appropriate to solar coronal plasmas. Electron densities deduced from the observed values of R ₁, R ₂, and R ₃ for several solar flares, measured from spectra obtained with the Naval Research Laboratory's S082A spectrograph on board Skylab, are found to be consistent. In addition, the derived electron densities are in excellent agreement with those determined from line ratios in Caxvi, which is formed at a similar electron temperature to Caxv. These results provide some experimental verification for the accuracy of the line ratio calculations, and hence the atomic data on which they are based. A set of eight theoretical Caxv line ratios involving 2s ²2p ²–2s2p ³ transitions in the wavelength range 140–216 Å are also found to be in good agreement with those measured from spectra of the TEXT tokamak plasma, for which the electron temperature and density have been independently determined. This provides additional support for the accuracy of the theoretical line ratios and atomic data.     

A deconvolved structure of the nebulosity associated with 3C 206

The published photographic profile of 3C 206 (reported in the low redshift sample of quasars by Wyckoffet al., 1981) has been deconvolved from the PSF by means of an effective restoration procedure. The deconvolved photometric structure of the quasar consists of a central point-like source, containing 68% of the integrated luminosity, an intermediate region of about 10 kpc radius (H ₀=60 km s^–1 Mpc^–1,q ₀=0) and an external region with nearly-linear slope and brightness level of the profile similar to those of the corresponding regions in giant elliptical and cD galaxies. The result confirms the previous findings in 3C 273, PKS 2135+147, and PKS 0812+020 obtained in the same way.     

V792 Her=HD 155638: A totally eclipsing RS CVn binary

Photometry of HD 155638=V792 Her has been analyzed to determine the elements of this totally eclipsing RS CVn binary. The light variation outside eclipse was found to have a period of 27^d.07±0^d.07, which is slightly different from the 27^d.5384±0^d.0045 orbital period. Analysis of the eclipses was achieved by a modification of the Russell-Merrill technique. With the aid of radial velocity measures, absolute elements were obtained for the hot and cool stars, respectively;R _h=2.58R ,R _c=12.28R ,M _h=1.40M ,M _c=1.46M ,i=80^o.61 and velocity semi-amplitudesK _c=48.36 km s^–1±0.79 km s^–1, andK _h=50.50 km s^–1±0.33 km s^–1. The apparent magnitudes areV _h=9 ^m .73 andV _c=8 ^m .48. The distance to HD 155638 was estimated to be 310 parsecs.     

The cosmic distance scale

Summary The status of the cosmic distance scale problem in early 1989 is reviewed. Internally consistent distances to Local Group galaxies are given in Tables 5 and 6. Within the Local Group the distance scale is found to be 11±5% smaller than that previously adopted by Sandage and Tammann. Distances to nearby galaxies are used as stepping stones to the Virgo cluster. The interpretation of the Tully-Fisher observations of Virgo spirals is found to be ambiguous because it is not yet clear which spirals are cluster members and which are background objects. Distance estimates of the Virgo cluster obtained by different techniques are listed in Table 11. The distance modulus of the Virgo cluster is found to be 31.5±0.2, corresponding to a distance of 20±2 Mpc. The elliptical galaxies in the core of the Virgo cluster haveV ₀=1200±46 kms^–1, which corresponds toV _LG=1082±48 km s^–1. With an infall velocity of 250±50 km s^–1 this yields a cosmological redshiftV=1332±69 km s^–1, from which a Hubble parameter H₀=67±8 km s^–1 Mpc^–1 is obtained. Space Telescope observations of distant Cepheids, Tully-Fisher observations of spirals in the Hercules eluster, and interference filter observations of Virgo planetary nebulae in the light of [O_III], should soon result in a major improvement in the accuracy with which H₀ is known.     

Theoretical intensities of Fexiv in the solar EUV spectrum

Equilibrium population of Fexiv levels in coronal conditions was calculated including configurations 3s ²3p, 3s3p ², 3s ²3d, 3p ³, 3s3p3d, 3s ²4s, 3s ²4p, 3s ²4d, 3s ²4f. Relative populations of selected levels are given in Table VII. Figure 1 shows the dependence of relative intensities of the strongest lines on electron density. Certain line ratios can be used for the determination of N _e .E.g., at T=2 × 10⁶ K and with a dilution factor 0.4, the intensity ratio of λ211.3 and λ219.0 changes by a factor of 65 if N _e increases from 10⁷ to 10¹¹ (Table VIII). Cascades from the 3s3p3d and 3p ³ configurations are important in the population of some levels of 3s3p ² (Table VI). A possibility of identification of additional lines in the solar spectrum is indicated. NAS-NRC Resident Research Associate.     

Equilibrium Structures of Differentially Rotating Primary Components of Binary Stars

In this paper a method is proposed for computing the equilibrium structures and various other observable physical parameters of the primary components of stars in binary systems assuming that the primary is more massive than the secondary and is rotating differentially about its axis. Kippenhahn and Thomas averaging approach (1970) is used in a manner earlier used by Mohan, Saxena and Agarwal (1990) to incorporate the rotational and tidal effects in the equations of stellar structure. Explicit expressions for the distortional terms appearing in the stellar structure equations have been obtained by assuming a general law of differential rotation of the typeω² = b ₀+b ₁ s ²+b ₂ s ⁴, where ω is the angular velocity of rotation of a fluid element in the star at a distance s from the axis of rotation, and b ₀, b ₁, b ₂ are suitably chosen numerical constants. The expressions incorporate the effects of differential rotation and tidal distortions up to second order terms. The use of the proposed method has been illustrated by applying it to obtain the structures and observable parameters of certain differentially rotating primary components of the binary stars assuming the primary components to have polytropic structures. This revised version was published online in July 2006 with corrections to the Cover Date.     

Extreme ultraviolet spectra of solar active regions and their analysis

Extreme ultraviolet spectra of several active regions are presented and analyzed. Spectral intensities of 3 active regions observed with the NRL Skylab XUV spectroheliograph (170–630 Å) are derived. From this data density sensitive line ratios of Mg viii, Si x, S xii, Fe ix, Fe x, Fe xi, Fe xii, Fe xiii, Fe xiv, and Fe xv are examined and typically yield, to within a factor of 2, electron pressures of 1 dyne cm^–2 (n _e T = 6 × 10¹⁵ cm^–3 K). The differential emission measure of the brightest 35 × 35 portion of an active region is obtained between 1.4 × 10⁴ K and 5 × 10⁶ K from HCO OSO-VI XUV (280–1370 Å) spectra published by Dupree et al. (1973). Stigmatic EUV spectra (1170–1710 Å) obtained by the NRL High Resolution Telescope and Spectrograph (HRTS) are also presented. Doppler velocities as a function of position along the slit are derived in an active region plage and sunspot. The velocities are based on an absolute wavelength scale derived from neutral chromospheric lines and are accurate to ±2 km s^–1. Downflows at 10⁵ K are found throughout the plage with typical velocities of 10 km s^–1. In the sunspot, downflows are typically 5 to 20 km s^–1 over the umbra and zero over the penumbra. In addition localized 90 and 150 km s^–1 downflows are found in the umbra in the same 1 × 1 resolution elements which contain the lower velocity downflows. Spectral intensities and velocities in a typical plage 1 resolution element are derived. The velocities are greatest ( 10 km s^–1) at 10⁵ K with lower velocities at higher and lower temperatures. The differential emission measure between 1.3 × 10⁴ K and 2 × 10⁶ K is derived and is found to be comparable to that derived from the OSO-VI data. An electron pressure of 1.4 dynes cm^–2 (n _e T = 1.0 × 10¹⁶ cm^–3 K) is determined from pressure sensitive line ratios of Si iii, O iv, and N iv. From the data presented it is shown that convection plays a major role in determining the structure and dynamics of the active region transition zone and corona.     

Fourier analysis of the light curves of eclipsing variables,VIII

The main aim of this paper will be to develop explicit form of the moments of the light curvesA _2m(r ₁,r ₂,i) required for the solution for the geometrical elementsr _1,2 andi of eclipsing systems exhibiting annular eclipses (Sections 2 and 3), as well as partial eclipses (Section 4).In the concluding Section 5 we shall demonstrate that — regardless of the type of eclipse and distribution of brightness on the apparent disc of the eclipsed star, or indeed of the shape of the eclipsing as well as eclipsed components — the momentsA _2m satisfy certain simple functional equations — a fact which relates them to other classes of functions previously studied in applied mathematics.     

High-resolution analysis of solar photospheric oscillations

The coherent 5-min photospheric pressure oscillations with spherical harmonic degrees in the range 100 <l< 1000 were directly imaged over the photosphere with the monochromatic solar telescope FPSS at Meudon Observatory. Movie films were obtained with images spatially filtered to select sizes of increasing wave numbers (or l). Areas with ephemeral concentrations of coherent waves evolve in shape and may move horizontally with velocities of several tenths of km s^–1. When a large number of waves are interacting, the maximum vertical velocity V _max of the pulsation reaches around 1000 m s^–1, irrespective of the size. Extrapolation to the ideal case of a single isolated wave gives V _max proportional to size. For the areas of the smallest scale measured (l = 1000), when about 100 waves are interacting, V _max is found to be 260 + 25 m s^–1 at an altitude of 210 km above the reference level ₅₀₀₀ = 1 and increases vertically with a scale height of 750 ± 400 km.