ALBERT

Description: We compute optical and infrared light curves of the pulsating class of delayed detonation models for Type Ia supernovae (SN Ia's) using an elaborate treatment of the Local Thermodynamic Equilbrium (LTE) radiation transport, equation of state and ionization balance, expansion opacity including the cooling by CO, Co(+), and SiO, and a Monte Carlo gamma-ray deposition scheme. The models have an amount of Ni-56 in the range from approximately or equal to 0.1 solar mass up to 0.7 solar mass depending on the density at which the transition from a deflagration to a detonation occurs. Models with a large nickel production give light curves comparable to those of typical Type Ia supernovae. Subluminous supernovae can be explained by models with a low nickel production. Multiband light curves are presented in comparison with the normally bright event SN 1992bc and the subluminous events Sn 1991bg and SN 1992bo to establish the principle that the delayed detonation paradigm in Chandrasekhar mass models may give a common explosion mechanism accounting for both normal and subluminous SN Ia's. Secondary IR-maxima are formed in the models of normal SN Ia's as a photospheric effect if the photospheric radius continues to increase well after maximum light. Secondary maxima appear later and stronger in models with moderate expansion velocities and with radioactive material closer to the surface. Model light curves for subluminous SN Ia's tend to show only one 'late' IR-maximum. In some delayed detonation models shell-like envelopes form, which consist of unburned carbon and oxygen. The formation of molecules in these envelopes is addressed. If the model retains a C/O-envelope and is subluminous, strong vibration bands of CO may appear, typically several weeks past maximum light. CO should be very weak or absent in normal Sn Ia's.

Description: We present the results of a deep optical survey for distant solar system objects. An area of 1.2 sq deg of the ecliptic has been imaged to apparent red magnitude 25, resulting in the detection of seven trans-Neptunian objects. These are the first detected members of a trans-Neptunian disk that compries about 35 000 objects larger than 100 km in the 30-50 AU heliocentric distance range. We interpret the new measurements using a set of Monte Carlo models in which the effects of observational bias in the data are taken into account.

Description: We investigate observable effects of anisotropic turbulence on the velocity profiles and eclipse behavior of emission lines from accretion disks. Turbulence expands the local line broadening profile, enhancing the surface brightness of saturated emission lines. Anisotropic turbulence produces anisotropic emission in such lines. The effects become observable when the turbulence exceeds the thermal velocity. Each term in the velocity-velocity correlation matrix produces a distinctive azimuthal pattern of enhanced emission-line surface brightness on the face of the accretion disk. These patterns express themselves as changes in the observable shapes of the disk's emission lines. The best place to look for turbulence effects is in saturated emission lines of heavy elements such as Ca, Mg, and Fe, which have a smaller thermal velocity at a given sound speed and at moderate inclination (60-70 degrees), since the Keplerian shear broadening dominates at higher inclinations.

Description: We have calculated a grid of more than 700 model atmospheres valid for a wide range of parameters encompassing the coolest known M dwarfs, M subdwarfs, and brown dwarf candidates: 1500 less than or equal to T(sub eff) less than or equal to 4000 K, 3.5 less than or equal to log g less than or equal to 5.5, and -4.0 less than or equal to (M/H) less than or equal to +0.5. Our equation of state includes 105 molecules and up to 27 ionization stages of 39 elements. In the calculations of the base grid of model atmospheres presented here, we include over 300 molecular bands of four molecules (TiO, VO, CaH, FeH) in the JOLA approximation, the water opacity of Ludwig (1971), collision-induced opacities, b-f and f-f atomic processes, as well as about 2 million spectral lines selected from a list with more than 42 million atomic and 24 million molecular (H2, CH, NH, OH, MgH, SiH, C2, CN, CO, SiO) lines. High-resolution synthetic spectra are obtained using an opacity sampling method. The model atmospheres and spectra are calculated with the generalized stellar atmosphere code PHOENIX, assuming LTE, plane-parallel geometry, energy (radiative plus convective) conservation, and hydrostatic equilibrium. The model spectra give close agreement with observations of M dwarfs across a wide spectral range from the blue to the near-IR, with one notable exception: the fit to the water bands. We discuss several practical applications of our model grid, e.g., broadband colors derived from the synthetic spectra. In light of current efforts to identify genuine brown dwarfs, we also show how low-resolution spectra of cool dwarfs vary with surface gravity, and how the high-regulation line profile of the Li I resonance doublet depends on the Li abundance.

Description: The sensitivity of the solar g-mode oscillation spectrum to variability in the universal gravitational constant G is described. Solar models in varying G cosmologies were constructed by evolving a zero-age main-sequence stellar model to the Sun's current age, while allowing the value of G to change according to the power law G(t) proportional to t(exp -beta), where Beta approximately equals delta G/GH and H is the Hubble constant. All solar models were constrained to the observed luminosity and radius at the current age of the Sun by adjusting the helium abundance and the mixing-length parameter of the models in the usual way for standard stellar models. Low-l g-mode oscillation periods were calculated for each of the models and compared to the claimed observation of the solar g-mode oscillation spectrum by Hill & Gu (1990). If one accepts Hill & Gu's claims, then within the uncertainties of the physics of the solar model calculation, our models rule out all but (delta G/GH) less than approximately 0.05. In other words, we conclude that G could not have varied by more than 2% over the past 4.5 Gyr, the lifetime of the present-day Sun. This result lends independent support to the validity of the standard solar model.

Description: A total of 663 galaxies with known redshifts in a 12 deg x 12 deg field centered on A2634, including 211 new measurements, are used to study in detail the structure of the region. In it we find six main galaxy concentrations: the nearby clusters A2634 and A2666, two groups in the vicinity of A2634, and two distant clusters at approximately 18,000 (A2622) and approximately 37,000 km/s seen in projection near the core of A2634. For A2634, the most richly sampled of those concentrations, we are able to apply strict cluster membership criteria. Two samples - one containing 200 galaxies within 2 deg from the cluster center and a second, magnitude-limited, of 118 galaxies within the central half degree - are used to examine the structure, kinematics, dynamics, and morphological segregation of the cluster. We show that early type galaxies appear to be a relaxed system, while the spiral population eschews the center of the cluster and exhibits both a multimodal velocity distribution and a much larger velocity dispersion that the ellipticals. We propose that the spiral galaxies of A2634 represent a dynamically young cluster population. For the galaxy component of A2634, we find no evidence of significant substructure in the central regions. We also conclude that the adoption of lenient membership criteria that ignore the dynamical complexity of A2634 are unlikely to be responsible for the conflicting results reported on the motion of this cluster with respect ot the CMB. The kinematical and dynamical analysis is extended to A2634's close companion, A2666, and the two distant background clusters.

Description: The multifluid diffusive model of G. P. Zank et al. (1994), which describes the interaction of the solar wind with a cometary plasma in the outer coma, has been used to model the structure of the Halley bow shock. The theoretical results are compared to in situ observations made by Giotto. We compare the solar wind and cometary ion number densities and pressures upstream and through the quasi-perpendicular and quasi-parallel shocks (observed on the inbound and outbound legs of the encounter, respectively). In general, good agreement is found between theory and observations in terms of shock structure, strength, and location, especially for the quasi-parallel shock. The comparison between tha quasi-perpendicular shock observations and theory is complicated by the apparently nonstationary behavior of the shock, a feature which has been remarked upon by other investigators. The cometary bow shock appears to be an excellent example of an energetic-particle-mediated shock where the energetic particles comprise less than 10% of the total number density.

Description: The Saturn magnetosphere model of Richardson and Sittler (1990) is extended to include the outer magnetosphere. The inner magnetospheric portion of this model is updated based on a recent reanalysis of the plasma data near the Voyager 2 ring plane crossing. The result is an axially symmetric model of the plasma parameters which is designed to provide accurate input for models needing either in situ or line-of-sight data and to be a useful tool for Cassini planning.

Description: We present a detailed study of the highest-frequency component of smooth radio emission observed during the Voyager 2 encounter with Neptune in August 1989. This emission occurs during three distinct periods on August 24 and 25, 1989, in the frequency range of 550 to 900 kHz. By assuming straight-line propagation from sources of both fundamental and second harmonic gyroemission, we perform a detailed analysis of the observed polarization of the emission. The data are most consistent with an L-O mode source in the north magnetic polar region, around 50 deg W, 50 deg N. A second possible source is in the north magnetic polar region, around 270 deg W, 50 deg N. This source must emit in the R-X mode.

Description: We present a complete solution for a set of magnetohydrodynamic (MHD) Riemann problems in which the upstream and downstream states have the same total pressure, and in which the normal component of the magnetic field is very small. These solutions are pertinent to subfast flows in the earth's magnetic tail and near the magnetopause. In a coplanar situation a family of solutions exists that depend on two parameters as well as on dissapation mechanisms. In the parallel case the transverse magnetic field either does not change direction or changes the direction twice by involving two intermediate shocks. In the antiparallel case an intermediate shock is always required, except when the solution consists of two switch-off shocks. In a noncoplanar case the solution is not self-similar as a function of x/t, but continues to evolve. At early times the evolution is similar to the coplanar case. In general two time-dependent intermediate shocks are required to rotate the magnetic fields. The velocity shear has a strong effect on the Riemann solution. In some cases no Riemann solution can exist because of the cavitation caused by the slow refraction waves. The calculated magnetopause structure resembles the observed structure for northward interplanetary magnetic field (IMF). However, for southward IMF, the MHD result shows the existence of a depletion layer, which is not supported by observations. We also show that on the magnetosheath side, the Walen relation, which is exact for a rotational discontinuity, can also be well satisfied by a slow shock, an intermediate shock, or the head of a slow rarefaction wave.

Description: L1551NE is a very young protostar (class I or perhaps class 0), located very close to L1551NE IRS 5. It is the second brightest far-infrared source in the Taurus molecular cloud complex, but its proximity to the brightest source IRS 5 has prevented effective observations of any molecular outflow. We here present evidence that it does indeed process an outflow, that the optical/infrared reflection nebula is associated with the blueshifted outflow lobe, and that the L1551W outflow does not originate from L1551NE, as has been suggested.

Description: A small fraction of galaxies appear to reside in dense compact groups, whose inferred crossing times are much shorter than a Hubble time. These short crossing times have led to considerable disagreement among researchers attempting to deduce the dynamical state of these systems. In this paper, we suggest that many of the observed groups are not physically bound but are chance projections of galaxies well separated along the line of sight. Unlike earlier similar proposals, ours does not require that the galaxies in the compact group be members of a more diffuse, but physically bound entity. The probability of physically separated galaxies projecting into an apparent compact group is nonnegligible if most galaxies are distributed in thin filaments. We illustrate this general point with a specific example: a simulation of a cold dark matter universe, in which hydrodynamic effects are included to identify galaxies. The simulated galaxy distribution is filamentary and end-on views of these filaments produce apparent galaxy associations that have sizes and velocity dispersions similar to those of observed compact groups. The frequency of such projections is sufficient, in principle, to explain the observed space density of groups in the Hickson catalog. We discuss the implications of our proposal for the formation and evolution of groups and elliptical galaxies. The proposal can be tested by using redshift-independent distance estimators to measure the line-of-sight spatial extent of nearby compact groups.

Description: Laboratory studies of infrared emission from gas-phase naphthalene in the 3.3 micrometer region following ultraviolet laser excitation are used to interpret the unidentified infrared bands observed in many astronomical objects. A time-resolved Fourier transform infrared emission technique acquires the time and spectrally resolved data. Two excitation wavelengths are employed: 193 nm and 248 nm. The infrared emission features are strongly dependent on the initial excitation energy. Wavelength-resolved spectra recorded 6.8 microseconds after the laser pulse show a 45/cm redshift from the gas-phase absorption spectra for 193 nm excitation and 25/cm for 248 nm excitation. We hypothesize that a series of sequence bands originating from the highly vibrationally excited ensemble of molecules is responsible for the observed shift. As collisional and radiative deactivation removes energy from the highly vibrationally excited molecules, the maximum in the emission profile gradually approaches the customary absorption maximum. This indicates that the amount of redshift is strongly dependent on the amount of internal vibrational energy in the molecule at the time of the vibrational transition.

Description: We employ an approximate treatment of dissipative hydrodynamics in three dimensions to study the coalescence of binary neutron stars driven by the emission of gravitational waves. The stars are modeled as compressible ellipsoids obeying a polytropic equation of state; all internal fluid velocities are assumed to be linear functions of the coordinates. The hydrodynamics equations then reduce to a set of coupled ordinary differential equations for the evolution of the principal axes of the ellipsoids, the internal velocity parameters, and the binary orbital parameters. Gravitational radiation reaction and viscous dissipation are both incorporated. We set up exact initial binary equilibrium configurations and follow the transition from the quasi-static, secular decay of the orbit at large separation to the rapid dynamical evolution of the configurations just prior to contact. A hydrodynamical instability resulting from tidal interactions significantly accelerates the coalescence at small separation, leading to appreciable radial infall velocity and tidal lag angles near contact. This behavior is reflected in the gravitational waveforms and may be observable by gravitational wave detectors under construction. In cases where the neutron stars have spins which are not aligned with the orbital angular momentum, the spin-induced quadrupole moment can lead to precession of the orbital plane and therefore modulation of the gravitational wave amplitude even at large orbital radius. However, the amplitude of the modulation is small for typical neutron star binaries with spins much smaller than the orbital angular momentum.

Description: We develop numerical models of accretions disks in cataclysmic variables (CVs), including and emphasizing the boundary layer region where the accretion disk meets the accreting white dwarf. We confine ourselves to solutions where the boundary layer region is vertically optically thick, and find that these solutions share several common features. The angular and radial velocities of the accreting material drop rapidly in a dynamical boundary layer, which has a radial width approximately 1%-3% of the white dwarf radius. The energy dissipated in this region diffuses through the inner part of the disk and is radiated from the disk surface in a thermal boundary layer, which has a radial width comparable to the disk thickness, approximately 5%-15% of the white dwarf radius. We examine the dependence of the boundary layer structure on the mass accretion rate, the white dwarf mass and rotation rate, and the viscosity parameter alpha. We delineate the boundary between optically thick and optically thin boundary layer solutions as a function of these parameters and suggest that by means of a careful comparison with observations it may be possible to estimate alpha in CVs. We derive an expression for the total boundary layer luminosities as a function of the parameters and show that it agrees well with the luminosites of our numerical solutions. Finally, we calcuate simple blackbody continuum spectra of the boundary layer and disk emission for our solutions and compare these to soft X-ray, EUV, and He II emission-line observations of CVs. We show that, through such comparisons, it may be possible to determine the rotation rates of the accreting stars in CVs, and perhaps also the white dwarf masses and the accretion rates. The spectra are quite insensitive to alpha, so the uncertainty in this parameter does not affect such comparisons.

Description: We model the temperature distribution at the surface of a magnetized neutron star and study the effects on the observed X-ray spectra and light curves. Generalrelativistic effects, i.e., redshift and lensing, are fully taken into account. Atmospheric effects on the emitted spectral flux are not included: we consider only blackbody emission at the local effective temperature. In this first paper we restrict ourselves to dipole fields. General features are studied and compared with the ROSAT data from the pulsars 0833 - 45 (Vela), 0656 + 14, 0630 + 178 (Geminga), and 1055 - 52, the four cases for which there is strong evidence that thermal radiation from the stellar surface is detected. The composite spectra we obtain are not very different from a blackbody spectrum at the star's effective temperature. We conclude that, as far as blackbody spectra are considered, temperature estimates using single-temperature models give results practically identical to our composite models. The change of the (composite blackbody) spectrum with the star's rotational phase is also not very large and may be unobservable inmost cases. Gravitational lensing strongly suppresses the light curve pulsations. If a dipole field is assumed, pulsed fractions comparable to the observed ones can be obtained only with stellar radii larger than those which are predicted by current models of neutron star struture, or with low stellar masses. Moreover, the shapes of the theoretical light curves with dipole fields do not correspond to the observations. The use of magnetic spectra may raise the pulsed fraction sufficiently but will certainly make the discrepancy with the light curve shapes worse: dipole fields are not sufficient to interpret the data. Many neutron star models with a meson condensate or hypersons predict very small radii, and hence very strong lensing, which will require highly nondipolar fields to be able to reproduce the observed pulsed fractions, if possible at all: this may be a new tool to constrain the size of neutron stars. The pulsed fractions obtained in all our models increase with photon energy: the strong decrease observed in Geminga at energies 0.3-0.5 keV is definitely a genuine effect of the magnetic field on the spectrum in contradistinction to the magnetic effects on the surface temperature considered her. Thus, a detailed analysis of thermal emission from the four pulsars we consider will require both complex surface field configurations and the inclusion of magnetic effects in the atmosphere (i.e., on the emitted spectrum).

Description: We have systematically analyzed a sample of 13 new and archival ROSAT Position Sensitive Proportional Counter (PSPC) observations of compact groups of galaxies: 12 Hickson compact groups plus the NCG 2300 group. We find that approximately two-thirds of the groups have extended X-ray emission and, in four of these, the emission is resolved into diffuse emission from gas at a temperature of kT approximately 1 keV in the group potential. All but one of the groups with extended emission have a spiral fraction of less than 50%. The baryon fraction of groups with diffuse emission is 5%-19%, similar to the values in clusters of galaxies. However, with a single exception (HCG 62), the gas-to-stellar mass ratio in our groups has a median value near 5%, somewhat greater than the values for individual early-type galaxies and two orders of magnitude than in clusters of galaxies. The X-ray luminosities of individual group galaxies are comparable to those of similar field galaxies, although the L(sub X)-L(sub B) relation for early-type galaxies may be flatter in compact groups than in the field.

Description: Plasma waves observed in the VLF range upstream of planetary bow shocks not only modify the particle distributions, but also provide important information about the acceleration processes that occur at the bow shock. Electron plasma oscillations observed near the tangent field line in the electron foreshock are generated by electrons reflected at the bow shock through a process that has been referred to as Fast Fermi acceleration. Fast Fermi acceleration is the same as shock-drift acceleration, which is one of the mechanisms by which ions are energized at the shock. We have generated maps of the VLF emissions upstream of the Venus bow shock, using these maps to infer properties of the shock energization processes. We find that the plasma oscillations extend along the field line up to a distance that appears to be controlled by the shock scale size, implying that shock curvature restricsts the flux and energy of reflected electrons. We also find that the ion acoustic waves are observed in the ion foreshock, but at Venus these emissions are not detected near the ULF forshock boundary. Through analogy with terrestrial ion observations, this implies that the ion acoustic waves are not generated by ion beams, but are instead generated by diffuse ion distributions found deep within the ion foreshock. However, since the shock is much smaller at Venus, and there is no magnetosphere, we might expect ion distributions within the ion foreshock to be different than at the Earth. Mapping studies of the terrestrial foreshock similar to those carried out at Venus appear to be necessary to determine if the inferences drawn from Venus data are applicable to other foreshocks.

Description: Hydrodynamical calculations of the evolution of a collapsing, rotating axisymmetric 10 solar masses molecular clump, including the effects of radiative acceleration but without magnetic fields, are represented. The initial cloud is assumed to be uniformly rotating, centrally condensed sphere with rho is proportional to r(exp -2). Several cases are considered, in which both the overall clump size and the total amount of angular momentum are varied. The calculations show how a warm, quasi-hydrostatic disk surrounding a central unresolved core of only a few solar masses forms and grows in size and mass. The disk is encased in two distinct accretion shock fronts, both of which are several scale heights above the equatorial plane. At the end of the calculation of our standard case, the central unresolved region is found to have a mass of 2.7 solar masses and a ratio of rotational to gravitational energy of approximately 0.45, sufficiently large to be unstable to nonaxisymmetric perturbations. In addition, the inner portions of the disk containing most of the mass are unstable according to the local Toomre criterion, implying that also in this region nonaxisymmetric perturbations will lead to rapid evolution. Under the assumption that gravitational torques would transport angular momentum out of this region, a central core of less than or approximately 8 solar masses with a stable disk of greater than or approximately = 2 solar masses should result. Frequency-dependent radiative transfer calculations of the standard case at selected ages show how the continuum spectrum of the structure depends on the disk's orientation and age and how the observed isophotal contours vary with wavelength. Because of the strong dependence on viewing angle, continuum spectra alone should not be used to estimate the evolutionary stage of development of these objects. Comparable results were obtained for the other cases considered.

Description: We have obtained the first average 2-500 keV spectra of Seyfert galaxies, using the data from Ginga and Compton Gamma-Ray Observatory's (CGRO) Oriented Scintillation Spectrometer Experiment (OSSE). Our sample contains three classes of objects with markedly different spectra: radio-quiet Seyfert 1's and 2's, and radio-loud Seyfert 1's. The average radio-quiet Seyfert 1 spectrum is well-fitted by a power law continuum with the energy spectral index alpha approximately equals 0.9, a Compton reflection component corresponding to a approximately 2 pi covering solid angle, and ionized absorption. There is a high-energy cutoff in the incident power law continuum: the e-folding energy is E(sub c) approximately equals 0.6(sup +0.8 sub -0.3) MeV. The simplest model that describes this spectrum is Comptonization in a relativistic optically-thin thermal corona above the surface of an accretion disk. Radio-quiet Seyfert 2's show strong netural absorption, and there is an indication that their X-ray power laws are intrinsically harder. Finally, the radio-loud Seyfert spectrum has alpha approximately equals 0.7, moderate neutral absorption E(sub C) = 0.4(sup +0.7 sub -0.2) MeV, and no or little Compton reflection. This is incompatible with the radio-quiet Seyfert 1 spectrum, and probably indicating that the X-rays are beamed away from the accretion disk in these objects. The average spectra of Seyferts integrated over redshift with a power-law evolution can explain the hard X-ray spectrum of the cosmic background.

Description: This contribution reports the results of an infrared imaging survey aimed at characterizing the stellar populations associated with the three densest star-forming cores in the Ophiuchus molecular cloud complex. The survey has sufficient sensitivity at J, H, and K (at 5 sigma limits of 16.5, 15.4, and 14.2) to provide a complete census of embedded young stellar objects (YSOs) with masses greater than the hydrogen-burning limit, provided that their ages are less than 3 Myr and that they are obscured by no more than approximately 18 mag of visual extinction. Our data suggest (1) a large fraction (greater than 70%) of the sources located within the cores are still surrounded by circumstellar disks and/or envelopes; and (2) the shape of the initial mass function for masses, M less than 1 solar mass, appears to be consistent with that derived from the solar neighborhood. We also report the results of a deeper imaging survey of centimeter continuum sources (14 sources) in these star-forming cores and in the larger Ophiuchus complex (eight sources). A large fraction (11/14) of the radio sources associated with the cores appear to have infrared excesses diagnostic of circumstellar accretion disks and/or infalling circumstellar envelopes. In these cases, the centimeter continuum radiation most likely diagnoses the ionized component of energetic winds or jets which characterizes YSOs during the disk accretion phase. By contrast, of the eight radio sources located outside dense cores, only two show infrared excesses.

Description: For the bright z = 3.02 radio-quiet quasar, HS 1946+7658, we have obtained radio to X-ray data within the past year: 5 GHz and 1415 MHz data from the Very Large Array (VLA); IR photometry at J, H, K, L prime (3.4 micrometers and N; IR spectroscopy; UBVRI photometry; optical spectrophotometry and high-resolution spectra; and an X-ray spectrum from the ROSAT Position Sensitive Proportional Counter (PSPC). The spectral energy distribution (SED) constructed from these data is compared to the mean SED for a set of low-redshift quasars, and while they appear generally similar, there are several differences. In relation to the low-redshift mean, the SED of HS 1946+7658 shows (1) only an upper limit at 10 micrometers (a rest wavelength of 2.5 micrometers, indicating that HS 1946+7658 does not have a strong near-IR excess such as hot dust would produce; (2) relatively weak Fe II and Balmer continuum emission at approximately 3000 A; and (3) a steeper turndown shortward of Ly alpha, even after correction for the Ly alpha forest absorption lines.

Description: We present a calculation of the profiles of emission lines originating in a relativistic, eccentric disk, and show examples of the resulting model profiles. Our calculations are motivated by the fact that in about one-quarter of the double-peaked emission lines observed in radio-loud active galactic nuclei (and in the mildly active nucleus of NGC 1097), the red peak is stronger than the blue peak, which is contrary to the prediction of relativistic, circular disk models. Using the eccentric disk model we fit some of the observed profiles that cannot be fitted with a circular disk model. We propose two possible scenarios for the formation of an eccentric disk in an active galactic nucleus: (a) tidal perturbation of the disk around a supermassive black hole by a smaller binary companion, and (b) formation of an elliptical disk from the debris resulting from the tidal disruption of a star by the central black hole. In the former case we show that the eccentricity can be long-lived because of the presence of the binary companion. In the latter case, although the inner parts of the disk may circularize quickly, we estimate that the outer parts will maintain their eccentricity for times much longer than the local viscous time. We suggest that it may be possible to detect profile variability on much shorter timescales than those ranging from a decade to several centuries by comparing the evolution of the line profile with detailed model predictions. We argue that line-profile variability may also be the most promising discriminant among competing models for the origin of asymmetric, double-peaked emission lines.

Description: Voyager 2 crossed the inbound or upstream Neptunian bow shock at 1430 spacecraft event time on August 24 in 1989 (Belcher et al., 1989). The plasma and magnetic field measurements allow us to study the solar wind interaction with the outermost gas giant. To fully utilize all of the spacecraft observations, an improved nonlinear least squares, 'Rankine-Hugoniot' magnetohydrodynamic shock-fitting technique has been developed (Szabo, 1994). This technique is applied to the Neptunian data set. We find that the upstream bow shock normal points nearly exactly toward the Sun consistent with any reasonable large-scale model of the bow shock for a near subsolar crossing. The shock was moving outward with a speed of 14 +/- 12 km/s. The shock can be characterized as a low beta, high Mach number, strong quasi-perpendicular shock. Finally, the shock microstructure features are resolved and found to scale well with theoretical expectations.

Description: An analytical model is developed for the diffusive equilibrium plasma density distribution in the Io plasma torus. The model has been employed successfully to follow the ray path of plasma waves in the multi-ion Jovian magnetosphere; it would also be valuable for other studies of the Io torus that require a smooth and continuous description of the plasma density and its gradients. Validity of the analytical treatment requires that the temperature of thermal electrons be much lower than the ion temperature and that superthermal electrons be much less abundant than the thermal electrons; these two conditions are satisfied in the warm outer region of the Io torus from L = 6 to L = 10. The analytical solutions agree well with exact numerical calculations for the most dense portion of the Io torus within 30 deg of the equator.

Description: We consider the galactic evolutionary history of He-3 in models which deplete deuterium by as much as a factor of 2 to approximately 15 from its primordial value to its present-day observed value in the interstellar medium (ISM). We show that when He-3 production in low-mass stars (1-3 solar mass) is included over the history of the galaxy, He-3 is greatly overproduced and exceeds the inferred solar values and the abundances determined in galactic H II regions. Furthermore, the ISM abundances show a disturbing dispersion which is difficult to understand from the point of view of standard chemical evolution models. In principle, resolution of the problem may lie in either (1) the calculated He-3 production in low-mass stars; (2) the observations of the He-3 abundance; or (3) an observational bias twoard regions of depleted He-3. Since He-3 observations in planetary nebula support the calculated He-3 production in low-mass stars, option (1) is unlikely. We will argue for option (3) since the He-3 interstellar observations are indeed made in regions dominated by massive stars in which He-3 is destroyed. In conclusion, we note that the problem with He-3 seems to be galactic and not cosmological.

Description: We have computed 8, 10, and 12 Gyr isochrones and physically consistent models of zero-age red horizontal branch stars for stellar masses between 0.55 and 1.3 solar mass, all at (fe/H) = +0.15. Comparison to the NGC 6791 BVI photometry of Kaluzny & Udalski (1992) and Montgomery et al. (1994) yields an age of 10.0 +/- 0.5 Gyr at an apparent distance modulus 13.49 less than (m-M)(sub V) less than 13.70. The color offsets required to fit the isochrones, combined with the spectroscopic results of Friel & Janes (1993), imply that the foreground reddening to NGC 6791 lies in the range 0.24 greater than E(B -V) greater than 0.19 with +0.27 less than (Fe/H) less than + 0.44. These results are derived using a technique by which we predict color shifts and apply these to the isochrones to simulate progressively higher metallicities. The zero-age horizontal branch model suggest that the red horizontal branch stars of NGC 6791 have masses approximately less than 0.7 solar mass. The masses are similar to those found for M67 red horizontal branch stars by Tripicco et al. (1993) and for globular cluster red horizontal branch stars, even though the M67 progenitors are approximately 0.2 solar-mass more massive, while the progenitors of globular cluster horizontal branch stars are similarly less massive. This suggests the presence of a mechanism, not strongly dependent on metallicity, which reduces stellar envelopes on the zero-age horizontal branch to a given mass rather than by a given amount.

Description: We report on two pointed observations of UZ For carried out by the imaging photometers aboard the Extreme Ultraviolet Explorer (EUVE), one as part of the EUVE Right Angle Program and one as an off-axis source during a guest observation. Both observations lasted approximately 3 days and covered a total of 72 orbits of the UZ For binary providing multiple coverage of all the orbital phases of UZ For. The resulting high signal-to-noise phase-folded light curve strongly constrains the emission and absorption geometry of UZ For. We have detected a narrow absorption dip that we attribute to the accretion stream at the location of the stagnation region many white dwarf radii away from the accretion spot and have also detected a broad dip caused by absorption much closer to the white dwarf surface. Both absorption effects are variable in time and phase. Based on the timescales of M-star eclipse ingress and egress, the angular spot size is constrained to be less than 5 deg; thus the ratio of spot area to white dwarf surface area is less than or equal to 0.0005. To explain the light curve phase duration given this small angular spot size, the extreme-ultraviolet (EUV) accretion spot must be raised vertically by approximately 5% of the white dwarf radius.

Description: We extend the approximate radiative transfer analysis of Hershkowitz, Linder, & Wagoner (1986) to a more general class of supernova model atmospheres, using a simple fit to the effective continuum opacity produced by lines (Wagoner, Perez, & Vasu 1991). At the low densities considered, the populations of the excited states of hydrogen are governed mainly by photoionization and recombination, and scattering dominates absorptive opacity. We match the asymptotic expressions for the spectral energy density J(sub nu) at the photosphere, whose location at each frequency is determined by a first-order calculation of the deviation of J(sub nu) from the Planck function B(sub nu). The emergent spectral luminosity then assumes the form L(sub nu) = 4 pi(squared)r(squared)(sub *) zeta(squared)B(sub nu)(T(sub p)), where T(sub p)(nu) is the photospheric temperature zeta is the dilution factor, and r(sub *) is a fiducial radius (ultimately taken to be the photospheric radius r(sub p)(nu)). The atmosphere is characterized by an effective temperature T(sub e) (varies as L(sup 1/4)r(sup -1/2)(sub *)) and hydrogen density n(sub H) = dependence of zeta on frequency nu and the parameters T(sub p), r(sub p), and alpha. The resulting understanding of the dependence of the spectral luminosity on observable parameters which characterize the relevant physical conditions will be of particular use in assessing the reliability of the expanding photosphere method of distance determination. This is particularly important at cosmological distances, where no information about the progenitor star will be available. This technique can also be applied to other low-density photosphere.

Description: We calculated bolometric temperature (T(sub bol)) and luminosity (L(sub bol)) for 128 young stellar objects (YSOs) in Taurus, 74 in the Ophiuchus 'core', and 33 in the Ophiuchus 'off-core' region. We have constructed the bolometric luminosity-temperature (BLT) diagram, the log-log plot of L(sub bol) versus T(sub bol), for the three samples. T(sub bol) is defined as the temperature of a blackbody having the same frequency as the observed continuum spectrum. It measures the redness (or coldness) of an astronomical source. The BLT diagram is analogous to the H-R diagram and allows for a direct and quantitative comparison of YSOs at a wide variety of evolutionary states, ranging from the most deeply embedded stars to T Tauri stars nearly on the main sequence. We found (1) T(sub bol) increases monotonically from embedded sources (approximately 60-500 K) to classical T Tauri stars (approximately 1000-3000 K) to weak-line T Tauri stars (approximately 2000-5000 K); (2) T(sub bol) correlates reasonably well with the age inferred from the evolutionary models of pre-main-sequence stars and protostars for embedded 'protostars' and weak-line T Tauri stars. There is no significant correlation for the classical T Tauri stars. These results can be understood in terms of dissipation of circumstellar dust envelope and disk during the early stages of stellar evolution. Sources in the three regions have different distributions in the BLT diagram. The Ophiuchus core has the highest fraction of cold sources among the three regions. These cold sources are also more luminous than the YSOs in the other regions. The Ophiuchus off-core sample is dominated by the more evolved pre-main-sequence stars. The Taurus sources have distributions intermediate in L(sub bol), T(sub bol), and age between the Ophiuchus core and off-core distributions. These may suggest differences in the star formation history, and possibly in the stellar masses and mass accretion rates in these star-forming regions.

Description: We simulate the origin and evolution of the bipolar nebula surrounding Eta Car using numerical two-dimensional gasdynamic models. The generalized interacting stellar winds scenario, wherein a stellar wind interacts with an aspherical circumstellar environment, is adopted. The eruption wind of 1840-1860, which is taken to be spherically symmetric, interacts with a preeruption toroidal density environment. Using reasonable assumptions of initial conditions and eruption parameters based on archival data, we have performed over 30 simulations in an effort to bracket the initial parameters which produce models that best match observations. We find that models with high pole-to-equator density contrasts (greater than 100) and toroidal density configurations nicely account for the observed morphology and kinematics of the homunculus.

Description: The effects of the radioactive decay of Al-26 on the circumstellar envelopes of asymptotic giant branch stars are analyzed. The gamma-rays emitted by the product nucleus Mg-26 escape most envelopes, but the beta-decay positrons are stopped and can ionize and heat the gas. The ionization may produce observable effects in C-rich circumstellar envelopes, particularly if the photospheric Al-26 abundance is as large as inferred from measurements of live Al-26 in the primitive solar nebula or the observations of interstellar 1.8 MeV gamma-rays. For the nearby carbon star IRC +10216, the measured abundance of the molecular ion HCO(+) provides an upper limit of about 4 x 10(exp -3) for the photospheric Al-26/Al-27 ratio, consistent with presolar SiC grains with about the same C-12/C-13 ratio.

Description: We calculate the evolution of the magnetic fields in white dwarfs, taking into account the Hall effect. Because this effect depends nonlinearly upon the magnetic field strength B, the time dependences of the various multipole field components are coupled. The evolution of the field is thus significantly more complicated than has been indicated by previous investigations. Our calculations employ recent white dwarf evolutionary sequences computed for stars with masses 0.4, 0.6, 0.8, and 1.0 solar mass. We show that in the presence of a strong (up to approximately 10(exp 9) G) internal toroidal magnetic field; the evolution of even the lowest order poloidal modes can be substantially changed by the Hall effect. As an example, we compute the evolution of an initially weak quadrupole component, which we take arbitrarily to be approximately 0.1%-1% of the strength of a dominant dipole field. We find that coupling provided by the Hall effect can produce growth of the ratio of the quadrupole to the dipole component of the surface value of the magnetic field strength by more than a factor of 10 over the 10(exp 9) to 10(exp 10) year cooling lifetime of the white dwarf. Some consequences of these results for the process of magnetic-field evolution in white dwarfs are briefly discussed.

Description: Soft gamma-ray repeaters appear to be relatively young, strongly magnetized neutron stars embedded within plerions. The observationally inferred presence of teragauss fields and energetic outflows in these sources provides some support for a model in which the bursts results from crustal disturbances that load and energize the magnetosphere with the ensuing generation of sheared Alfven waves. We here determine accurately the structure of the particle efflux and calculate the spectrum of the resultant radiation, incorporating the key effects due to mildly relativistic magnetic bremsstrahlung, a frequency-dependent photospheric radius and angle-dependent boosting of the intensity for emitting elements at different inclinations. We find that the average photospheric radius during the transient event is approximately 4 stellar radii and that the particles attain a Lorentz factor approximately 13.2 by the time the gas becomes optically thin. The calculated spectrum is a weighted integral over the boosted intensity determined from the physical conditions (i.e., magnetic field, particle density, and temperature) at the photospheric radius corresponding to the relevant frequency. When applied to SGR 1806-20, this procedure yields a distance D approximately equals 10 kpc to this source, which appears to be consistent with the range (approximately 8-17 kpc) inferred from its apparent coincidence with the supernova remnant G10.0-0.3.

Description: We identify a new regime of time dependent helium burning for high accretion rate neutron stars and suggest that this burning is the origin of the low-level luminosity variations (on timescales of 10-10(exp 4) s, designated the 'very low-frequency noise'(VLFN) by van der Klis and collaborators) always detected in the brightest accreting X-ray sources. Only two nuclear burning regimes were previously recognized. At accretion rates in excess of the Eddington limit (dot-M approximately greater than (1-3) x 10(exp -8) solar mass/yr), the accreted matter fuses steadily. At very low dot-M, the star's entire surface is rapidly (approximately less than 10 s) burned by a fast propagating convective burning front at regular intervals, giving quasi-periodic Type I X-ray bursts. We show that for the observationally interesting range of 5 x 10(exp -10) solar mass/yr approximately less than dot-M approximately less than 10(exp -8) solar mass/yr, parts of the stellar surface burn slowly. At these accretion rates, a local thermonuclear instability starts a fire which propagates horizontally at v approximately 300 cm/s. The fire propagates around the flammable surface in roughly the same time it takes to accrete enough fuel for the next instability (approximately 10(exp 3)-10(exp 4), so that only a few fires are burning at once, giving rise to large luminosity flares. Nuclear burning is always time dependent for sub-Eddington local accretion rates: a local patch undergoes a recurrent cycle, accumulation fuel for hours until it becomes thermally unstable or is 'ignited' by a nearby burning region. The global pattern of burning and the resulting luminosity are thus very dependent on how fast nuclear fires spread around the star. The nuclear burning luminosity is not uniform over the stellar surface and so may provide a handle on measuring, or constraining, the spin periods of these neutron stars.

Description: Emission in the forbidden lines of ionized sulfur arises in both predominantly ionized and predominantly neutral regions. The separation of (S II) emission from supernova remnants into nearly orthogonal dependences on (O I) and a power of (N II)/H alpha is discussed. (S II) may be excited via photoionization into forbidden-line-emitting states in regions devoid of electronic excitation of (O I). Emission along 'background' lines of sight in the Milky Way and in extragalactic 'froth' is considered in terms of contributions of neutral gas.

Description: Hubble Space Telescope (HST) Faint Object Camera (FOC) f/48 images of M31, M32, and NGC 205 (field of view 23 sec x 23 sec with 0.45 sec pixel size) are analyzed as observed through the combined UV filters F150W and F130LP. The absolute calibration of the data and the internal disagreement between observed and expected count rates in the UV region lead us to suggest that the filter combination F150W+F130LP suffers from a 5 times degraded UV sensitivity. A corrected efficiency curve is construced using the UV/optical spectral energy distributions of these three galaxies, which is consistent with all of the data analyzed here. Eighty-one individual stars are detected in M31, 10 stars in M32, and 78 stars in NGC 205. Comparisons with other UV images and optical images indicates that these stars are hot, UV-bright stars, even though our corrected efficiency curve suggests that flux from 1200-2450 A contributes only 7% of the counts in M32, 19% in M31, and 60% in NGC 205. The complex nucleus of M31 as seen by Lauer et al. (1993) is confirmed; M32 has a generally smooth appearance and NGC 205 is dominated by a UV-bright, somewhat resolved nucleus. Analysis of these data is done through the new, extensive stellar isochrones of Bertelli et al. (1994) and the population synthesis models of Bressan, Chiosi, and Fagotto (1994). This analysis shows that high-metal stars (Z greater than 0.05) evolve into UV-bright stars (P-EAGB, H-HB, and AGB-manque stars) that are less luminous and cooler but are significantly longer lived than the P-AGB stars produced by stars with Z less than 0.05. Moreover, the proportion of P-EAGB, H-HB, and AGB-manque stars is also a function of age, with older stars of fixed mean metallicity having a higher proportion than younger stars. Hence, with either metallicity or age differences as an interpretation of the line-strength luminosity correlation for ellipticals, the high-metallicity 'tail' of the stellar content of a galaxy can produce far-UV flux in much greater proportion than its actual proportion of galaxy mass. The resulting model of the sources of far-UV flux is inherently composite, with the total UV flux from a stellar population both rapidly increasing and changing its mean spectrum with increasing mean metallicity (or mean age). This model is consistent with several pieces of observational evidence.

Description: The pure rotational spectrum of the magnesium monoacetylide radical, MgCCH, has been recorded in the laboratory using millimeter/sub-mm direct absorption spectroscopy. These measurements constitute the first time this molecule has been observed in the laboratory by any spectroscopic method. Seventeen rotational transitions were observed in the frequency range 210-370 GHz for MgCCH, which appears to be a linear molecule with a (sup 2)Sigma ground electronic state. Rotational and fine structure constants were determined for this radical from a nonlinear least squares fit to the data. The rotational rest frequencies measured here will enable astronomical observations to be carried out for MgCCH towards IRC +10216, where the magnesium compounds MgCN and MgNC, as well as many acetylide species, are present.

Description: This work addresses the problem of assigning confidence intervals to estimated photometry data obtained from astronomical observations. The proposed solution is to estimate the Cramer-Rao bound, which is an analytical expression that describes the minimum obtainable mean square error associated with a given estimate of a parameter. This Letter presents a compact and simple form for the bound associated with a linear estimator such as a Wiener filter estimator. A prescription for estimating the variance associated with each element in a restored object was developed using an analytical model for observed data corrupted by either Poisson or Gaussian noise. Both one- and two-dimensional examples are presented.

Description: Magnetohydrodynamic simulations have been made of the formation of outflows from a Keplerian disk threaded by a magnetic field. The disk is treated as a boundary condition, where matter is ejected with Keplerian azimuthal speed and poloidal speed less than the slow magnetosonic velocity, and where boundary conditions on the magnetic field correspond to a highly conducting disk. Initially, the space above the disk, the corona, is filled with high specific entropy plasma in thermal equilibrium in the gravitational potential of the central object. The initial magnetic field is poloidal and is represented by a superposition of monopoles located below the plane of the disk. The rotation of the disk twists the initial poloidal magnetic field, and this twist propagates into the corona pushing and collimating matter into jetlike outflow in a cylindrical region. Matter outflowing from the disk flows and accelerates in the z-direction owing to both the magnetic and pressure gradient forces. The flow accelerates through the slow magnetosonic and Alfven surfaces and at larger distances through the fast magnetosonic surface. The flow velocity of the jet is approximately parallel to the z-axis, and the collimation results from the pinching force of the toroidal magnetic field. For a nonrotating disk no collimation is observed.

Description: The Lyman continuum radiation from the brightest extreme ultraviolet (EUV) source, the B2 II star epsilon Canis Majoris (Adara), is so intense that it dominates the local stellar EUV radiation field at wavelengths longer than 450 A and therefore sets a lower limit to the ionization of hydrogen in the Local Cloud. Using the EUV (70-730 A) spectrum of epsilon CMa taken with the Extreme Ultraviolet Explorer Satellite (EUVE) and simple models that extrapolate this spectrum to the Lyman edge at 912 A, we have determined the local interstellar hydrogen photionizatin parameter Gamma solely from epsilon CMa to be 1.1 x 10(exp -15)/s. This fiugre is a factor of 7 greater than previous estimates of Gamma calculated for all nearby stars combined (Bruhweiler & Cheng 1988). Using measured values of the density and temperature of neutral interstellar hydrogen gas in the Local Cloud, we derive a particle density of ionized hydrogen n(H(+)) and electrons n(sub e) of 0.015-0.019/cu cm assuming ionization equilibrium and a helium ionization fraction of less than 20%. These values correspond to a hydrogen ionizatin fraction, chi(sub H) from 19% to 15%, respectively. The range of these derived quantities is due to the uncertainties in the local values of the neutral hydrogen and helium interstellar densities derived from both (1) solar backscatter measurements of Ly alpha lines of hydrogen and helium (1216 and 584 A), and (2) the average neutral densities along the line of sight to nearby stars. The local proton density produced by epsilon CMa is enough to allow the ionization mechanism of Ripken & Fahr (1983) to work at the heliopause and explain the discrepancy between the neutral hydrogen density derived from solar backscatter measurements and line-of-sight averages to nearby stars. A large value of electron density in the Local Cloud of n(sub e) is approximately 0.3-0.7/cu cm (T = 7000 K) has recently been reported by Lallement et al. (1994) using observations of Mg II and Mg I toward Sirius A. We show that if such a high value exists, it cannot result from the EUV stellar radiation field and, therefore, must be due to a strong diffuse source of EUV radiation.

Description: We determine the physical properties of the accretion region of the AM Her-type binary VV Puppis using extreme ultraviolet (EUV) medium-resolution spectroscopy and photometry obtained with the Extreme Ultraviolet Explorer (EUVE) observatory. The EUV continuum from VV Pup was detected in the wavelength range from 75 to 135 A and was simultaneously recorded with the Deep Survey/Spectrometer (DS/S) imaging telescope, allowing for the extraction of an accurate light curve. VV Pup appeared to have entered a high-accretion state just prior to the pointed EUVE observations. We use the EUV light curve to infer the diameter of the accretion region (d = 220 km) assuming a hemispherical geometry and a radius of 9000 km for the white dwarf. We perform a model atmosphere analysis and, based on the light curve properties and assuming a distance of 145 pc, we derive an effective temperature of the accretion region in the range 270,000 is less than T(sub eff) is less than 360,000 K and a neutral hydrogen column density in the local interstellar medium of n(sub H) = 1.9 - 3.7 x 10(exp 19)/sq cm. The total EUV/soft X-ray energy radiated by the accretion region is approximately 3.5 x 10(exp 32) ergs/s. Our results provide a first verification of past suggestions that deep heating of the white dwarf surface produces the soft X-ray flux from the polars. We present a possible detection of O VI absortion features, and we suggest that extensive EUVE observations targeting high-accretion events may result in oxygen and heavier element abundance determination in the accretion region.

Description: The observed magnetism of asteroids such as Gaspra and Ida (and other small bodies in the solar system including the Moon and meteorites) may have resulted from an impact-induced shock wave producing a thermodynamic state in which iron-nickel alloy, dispersed in a silicate matrix, is driven from the usual low-temperature, low-pressure, alpha, kaemacite, phase to the paramagnetic, epsilon (hcp), phase. The magnetization was acquired upon rarefaction and reentry into the ferromagnetic, alpha, structure. The degree of remagnetization depends on the strength of the ambient field, which may have been associated with a Solar-System-wide magnetic field. A transient field induced by the impact event itself may have resulted in a significant, or possibly, even a dominant contribution, as well. The scaling law of Housen et al. (Housen, K. R., R. M. Schmidt, and K. A. Holsapple 1991) for catastrophic asteroid impact disaggregation imposes a constraint on the degree to which small planetary bodies may be magnetized and yet survive fragmentation by the same event. Our modeling results show it is possible that Ida was magnetized when a large impact fractured a 125 +/- 22-km-radius protoasteroid to form the Koronis family. Similarly, we calculate that Gaspra could be a magnetized fragment of a 45 +/- 15 km-radius protoasteroid.

Description: The satellite observations at comet Halley have shown strong heating of solar wind alpha particles over an extended region dominated by high-intensity, low-frequency turbulence. These waves are excited by the water group pickup ions and can energize the solar wind plasma by different heating processes. The alpha particle heating by the Landau damping of kinetic Alfven waves and the transit time damping of low-frequency hydromagnetic waves in this region of high plasma beta are studied in this paper. The Alfven wave heating was shown to be the dominant mechanism for the observed proton heating, but it is found to be insufficient to account for the observed alpha particle heating. The transit time damping due to the interaction of the ions with the electric fields associated with the magnetic field compressions of magnetohydrodynamic waves is found to heat the alpha particles preferentially over the protons. Comparison of the calculated heating times for the transit time damping with the observations from comet Halley shows good agreement. These processes contribute to the thermalization of the solar wind by the conversion of its directed energy into the thermal energy in the transition region at comet-solar wind interaction.

Description: The complex molecules vinyl cyanide (CH2CHCN), methyl formate (HCOOCH3), and ethyl cyanide (CH3CH2CN) were observed in the Sgr B2 star-forming region with the BIMA millimeter wavelength array. A region with diameter less than 0.1 pc toward the Sgr B2(N) molecular core is found to be the major source of these molecules. Also, this source is coincident with continuum emission from dust and a center of H2O maser activity. Ultracompact (UC) H 11 regions are located within 0.1 pc. Strikingly, none of these molecules is detected toward Sgr B2(M), a core located 1 minute south of Sgr B2(N). The existence of complex molecules, a large mass of dust, high-velocity H2O masers, and UC H 11 regions strongly suggests that the Sgr B2(N) region has just begun to form stars, while the absence of strong dust emission and large molecules suggests Sgr B2(M) is more evolved. The detection of large molecules coincident with continuum emission from dust supports the idea found in current chemical models that grain chemistry is of crucial importance for the formation of these molecules.

Description: A new metal-containing molecule, MgCN, has been detected toward the late-type star IRC + 10216, using the NRAO 12 m and IRAM 30 m telescopes. The N = 11 approaches 10, 10 approaches 9, and 9 approaches 8 transtions of this species which has a (sup 2)Sigma(sup +) ground state, have been observed in the outer envelope of this object at 3 mm. For the N = 11 approaches 10 transitions, the two spin-rotation components are clearly resolved and conclusively identify this new radical. These measurements imply a column of density for MgCN of N(sub tot) approximately 10(exp 12)/sq cm in the outer shell, which corresponds to a fractional abundance of f approximately 7x10(exp -10). This molecule, the metastable isomer of MgNC, is the third metal-bearing species thus far identified in the outer shell of IRC + 10216, and its detection implies a ratio of MgNC/MgCN approximately 22/1. MgCN may be formed through a reaction scheme involving magnesium and HNC or CN, both prominent outer shell molecules, or through synthesis on grains.

Description: We derive the velocity and mass-loss rate of a steady state Wolf-Rayet (WR) wind, using a nonisotropic diffusion approximation applied to the transfer between strongly overlapping spectral lines. Following the approach of Friend & Castor (1983), the line list is assumed to approximate a statistically parameterized Poisson distribution in frequency, so that photon transport is controlled by an angle-dependent, effectively gray opacity. We show the nonisotropic diffusion approximation yields good agreement with more accurate numerical treatments of the radiative transfer, while providing analytic insight into wind driving by multiple scattering. We illustrate, in particular, that multiple radiative momentum deposition does not require that potons be repeatedly reflected across substantial distances within the spherical envelope, but indeed is greatest when photons undergo a nearly local diffusion, e.g., through scattering by many lines closely spaced in frequency. Our results reiterate the view that the so-called 'momentum problem' of Wolf-Rayet winds is better characterized as an 'opacity problem' of simply identfying enough lines. One way of increasing the number of thick lines in Wolf-Rayet winds is to transfer opacity from saturated to unsaturated lines, yielding a steeper opacity distribution than that found in OB winds. We discuss the implications of this perspective for extending our approach to W-R wind models that incorporate a more fundamental treatment of the ionization and excitation processes that determine the line opacity. In particular, we argue that developing statistical descriptions of the lines to allow an improved effective opacity for the line ensemble would offer several advantages for deriving such more fundamental W-R wind models.

Description: The magnetometer onboard the Giotto spacecraft observed a diamagnetic cavity surrounding the nucleus of comet Halley. A narrow transition layer with enhanced plasma density is formed at this shock-like boundary as a flux of cometary ions flowing into it from within the cavity is removed by electron-ion recombination. We examine the structure of this layer using both a two-dimensional magnetohydrodynamical model and a one- dimensional hybrid code.

Description: We study the secular evolution and gravitational wave signature of a newly formed, rapidly rotating neutron star. The neutron star may arise from core collapse in a massive star or from the accretion-induced collapse of a white dwarf. After a brief dynamical phase, the nascent neutron star settles into an axisymmetric, secularly unstable state. Gravitational radiation drives the star to a nonaxisymmetric, stationary equilibrium configuration via the bar-mode instability. The emitted quasi-periodic gravitational waves have a unique signature: the wave frequency sweeps downward from a few hundred Hertz to zero, while the wave amplitude increase from zero to a maximum and then decays back to zero. Such a wave signal could detected by broadband gravitational wave interferometers currently being constructed. We also characterize two other types of gravitational wave signals that could arise in principle from a rapidly rotating, secularly unstable neutron star: a high-frequency (f greater than or approximately = 1000 Hz) wave which increases the pattern-speed of the star, and a wave that actually increases the angular momentum of the star.

Description: Variations in the electron pressure can produce significant distortions of the cosmic microwave background radiation (CMBR) spectrum through the Sunyaev-Zel'dovich effect. In this paper, we compute the amplitude of these distortions in various cosmological scenarios using the results of hydrodynamical simulations combined with analytic theory. For reionized cold dark matter (CDM) with h = 0.5 probed with the OVRO experiment (Theta approximately = 0.1 deg), we find that the amplitude of these fluctuations is (Delta T/T)(sub rms) = 1 x 10(exp -5) (Omega(sub i)/0.06) if the experiment looks at a random part of the sky and 6 x 10(exp -6)(Omega(sub i)/(0.06) if bright areas are avoided. Here Omega is the density of ionized baryons in units of the critical density, rho(sub c) = 2.78 x 10(exp 11)h(exp 2)/cu Mpc. Primary fluctuations would generate (Delta T/T)(sub rms) approximately = 2 x 10(exp -6) for this model. These distortions could be a dominant source of CMBR anisotropics on small and intermediate angular scales. It is important to note that these signals are generated relatively recently, when nonlinear effects are improtant, and do not require early reionization. We also examine the contribution to anisotropies from the Doppler effect, in both the linear and nonlinear regimes. In this case, the bulk of the signal is generated before a redshift of 10, provided that the intercluster medium is ionized back to the surface of last scattering. On the smallest scales thermal effects are bigger, while, on larger scales the Doppler-induced fluctuations are greater with the exact boundary determined by the scenario. However, over a wide range of scales, the two signals are equal to within a factor of a few.

Description: A brief overview of the physics thought to determine the properties of the heliospheric termination shock is given and the central role of hot interstellar pick-up ions and cosmic rays is emphasized.

Description: A general theoretical overview for the sources, sinks, gas-surface interactions, and transport dynamics of sodium and potassium in the exospheric atmosphere of the Moon is given. These four factors, which control the spatial distribution of these two alkali-group gases about the Moon, are incorporated in numerical models. The spatial nature and relative importance of the initial source atoms atmosphere (which must be nonthermal to explain observational data) and the ambient (ballistic hopping) atom atmosphere are examined. The transport dynamics, atmospheric structure, and lunar escape of the nonthermal source atoms are time variable with season of the year and lunar phase because of their dependence on the radiation acceleration experienced by sodium and potassium atoms as they resonantly scatter solar photons. The dynamic transport time of fully thermally accomodated ambient atoms along the surface because of solar radiation acceleration (only several percent of surface gravity) is larger than the photoionization lifetimes and hence unimportant in determining the local density, although for potassium the situation is borderline. The sodium model was applied to analyze sodium observations of the sunward brightness profiles acquired near last quarter by Potter & Morgan (1988b) extending from the surface to an altitude of 1200 km, and near first quarter by Mendillo, Baumgardner, & Flynn (1991), extending in altitude from approximately 1430 to approximately 7000 km. The observations at larger altitudes could be fitted only for source atoms having a velocity distribution with a tail that is mildly nonthermal (like an approximately 1000 K Maxwell-Boltzmann distribution). Solar wind sputtering appears to a be a viable source atom mechanism for the sodium observations with photon-simulated desorption also possible but highly uncertain, although micrometeoroid impact vaporization appears to have a source that is too small and too hot, with likely an incorrect angular distribution about the Moon.

Description: Interstellar scintillation (ISS), fluctuations in the amplitude and phase of radio waves caused by scattering in the interstellar medium, is important as a diagnostic of interstellar plasma turbulence. ISS is also of interest because it is noise for other radio astronomical observations. The unifying concern is the power spectrum of the interstellar electron density. Here we use ISS observations through the nearby (less than or approximately =1 kpc) (ISM) to estimate the spectrum. From measurements of angular broadening of pulsars and extragalactic sources, decorrelation bandwidth of pulsars, refractive steering of features in pulsar dynamic spectra, dispersion measured fluctuations of pulsars, and refractive scintillation index measurements, we construct a composite structure function that is approximately power law over 2 x 10(exp 6) m less than scale less than 10(exp 13) m. The data are consistent with the structure function having a logarithmic slope versus baseline less than 2; thus there is a meaningful connection between scales in the radiowave fluctuation field and the scales in the electron density field causing the scattering. The data give an upper limit to the inner scale, l(sub o) less than or approximately 10(exp 8) m and are consistent with much smaller values. We construct a composite electron density spectrum that is approximately power law over at least the approximately = 5 decade wavenumber range 10(exp -13)/m less than wavenumber less than 10(exp -8)/m and that may extend to higher wavenumbers. The average spectral index of electron density over this wavenumber range is approximately = 3.7, very close to the value expected for a Kolmogorov process. The outer scale size, L(sub o), must be greater than or approximately = 10(exp 13) m (determined from dispersion measure fluctuations). When the ISS data are combined with measurements of differential Faraday rotation angle, and gradients in the average electron density, constraints can be put on the spectrum at much smaller wave numbers. The composite spectrum is consistent with a Kolmogorov-like power law over a huge range (10 or more decades) of spatial wavenumber with an infrared outer scale L(sub o) greater than or approximately 10(exp 18)m. This power-law subrange-expressed as ratio of outer to inner scales-is comparable to or larger than that of other naturally occurring turbulent fluids, such as the oceans or the solar wind. We outline some of the theories for generating and maintaining such a spectrum over this huge wavenumber range.

Description: We present non-local thermodynamic equilibrium (non-LTE) synthetic spectra for the Type Ia supernovae SN 1992A and SN 1981B, near maximum light. At this epoch both supernovae were observed from the UV through the optical. This wide spectral coverage is essential for determining the density structure of a SN Ia. Our fits are in good agreement with observation and provide some insight as to the differences between these supernovae. We also discuss the application of the expanding photosphere method to SNe Ia which gives a distance that is independent of those based on the decay of Ni-56 and Cepheid variable stars.

Description: The Cosmic Background Explorer Satellite (COBE) detection of microwave anisotropies provides the best way of fixing the amplitude of cosmological fluctuations on the largest scales. This normalization is usually given for an n = 1 spectrum, including only the anisotropy caused by the Sachs-Wolfe effect. This is certainly not a good approximation for a model containing any reasonable amount of baryonic matter. In fact, even tilted Sachs-Wolfe spectra are not a good fit to models like cold dark matter (CDM). Here, we normalize standard CDM (sCDM) to the two-year COBE data and quote the best amplitude in terms of the conventionally used measures of power. We also give normalizations for some specific variants of this standard model, and we indicate how the normalization depends on the assumed values on n, Omega(sub B) and H(sub 0). For sCDM we find the mean value of Q = 19.9 +/- 1.5 micro-K, corresponding to sigma(sub 8) = 1.34 +/- 0.10, with the normalization at large scales being B = (8.16 +/- 1.04) x 10(exp 5)(Mpc/h)(exp 4), and other numbers given in the table. The measured rms temperature fluctuation smoothed on 10 deg is a little low relative to this normalization. This is mainly due to the low quadrupole in the data: when the quadrupole is removed, the measured value of sigma(10 deg) is quite consistent with the best-fitting the mean value of Q. The use of the mean value of Q should be preferred over sigma(10 deg), when its value can be determined for a particular theory, since it makes full use of the data.

Description: We present the results of an extensive grid of evolutionary synthesis models for populations of massive stars. The parameter space has been chosen to correspond to conditions typically found in objects like giant H II regions, H II galaxies, blue compact dwarf galaxies, nuclear starbursts, and infrared luminous starburst galaxies. The models are based on the most up-to-date input physics for the theory of stellar atmospheres, stellar winds, and stellar evolution. A population of massive stars is not only important in terms of its output of radiation but also via its deposition of mechanical energy. The output of radiative and mechanical luminosity is compared at various starburst epochs. In a supernova dominated instantaneous starburst, the mechanical luminosity can be as large as almost 10% of the total radiative luminosity. This occurs when most massive O stars have disappeared, and the synthetic spectrum in the optical and near-ultraviolet is dominated by B and A stars. During this epoch, the output of ionizing radiation below 912 A becomes very small, as indicated by a very large Lyman discontinuity and a very small ratio of ionizing over mechanical luminosity. We discuss the relevance of these results for the interpretation of starburst galaxies, active galactic nuclei, and the energetics of the interstellar medium.

Description: Estimates indicate there may be greater than or approximately equal to 10(exp 3) low-luminosity X-ray pulsars (L less than or approximately equal to 10(exp 34) ergs/s) in the Galaxy undergoing 'low-state' wind accretion in Be/X-ray binary systems, and approximately 10(exp 8)-10(exp 9) isolated neutron stars which may be accreting directly from the interstellar medium. Despite their low effective temperatures (kT(sub e) less than or approximately equal to 300 eV), low-luminosity accreting neutron stars with magnetic fields B approximately (0.7-7) x 10(exp 12) G could emit a substantial fraction (0.5%-5%) of their total luminosity in a moderately broadened (Epsilon/delta Epsilon approximately 2-4) cyclotron emission line which peaks in the energy range approximately 5-20 keV. The bulk of the thermal emission from these stars will be in the extreme ultraviolet/soft X-ray regime. In sharp contrast, the nonthermal cyclotron component predicted here will not be strongly absorbed, and consequently it may be the only distinguishing signature for the bulk of these low-luminosity sources. We propose a search for this cyclotron emission feature in long pointed observations of the newly discovered candidate isolated neutron star MS 0317.7-6477, and the Be/X-ray transient pulsar 4U 0115+63 in its quiescent state. We note that an emission-like feature similar to the one we predict here has been reported in the energy spectrum of the unusual X-ray pulsar 1E 2259+586.

Description: We present a new numerical study of the equilibrium and stability properties of close binary systems. We use the smoothed-particle hydrodynamics (SPH) technique both to construct accurate equilibrium configurations in three dimensions and to follow their hydrodynamic evolution. We adopt a simple polytropic equation of state p = K(sub rho)(exp gamma) with gamma = 5/3 and K = constant within each star, applicable to low-mass degenerate dwarfs as well as low-mass main-sequence stars. For degenerate configurations, we set the two polytropic constants equal, K = K prime, independent of the mass ratio. For main-sequence stars, we adjust K and K prime so as to obtain a simple mass-radius relation of the form R/R prime = M/M prime, where R prime and M prime are the radius and mass of the secondary. Along a sequence of binary equilibrium configurations for two identical stars, we demonstrate the existence of both secular and dynamical instabilities, confirming directly the results of recent analytic work. We use the SPH method to calculate the nonlinear development of the dynamical instability and to determine the final fate of the system. We find that the two stars merge together into a single, rapidly rotating object in just a few orbital periods. Equilibrium sequences are also constructed for systems containing two nonidentical stars. These sequences terminate at a Roche limit, which we can determine very accurately using SPH. For two low-mass main-sequence stars with mass ratio q approximately less than 0.4 we find that the (synchronized) Roche limit configuration is secularly unstable. We discuss the implications of our results for the evolution of double white-dwarf systems and W Ursae Majoris binaries.

Description: We present some analyses of warping modes in inviscid near Keplerian disks taking the three-dimensional structure fully into account. The results of this investigation verify the validity of a vertical averaging approximation for thin disks when the radial wavelength is significantly longer than the disk thickness. They also indicate that long wavelength disturbances may persist for long times. Shorter wavelength disturbances in non-self gravitating disks are found to propagate with little dispersion at a speed related to the sound speed. When self-gravity becomes important, fast and slow waves are found which also propagate with little dispersion. When a small viscosity is included, the evolution of the disturbances becomes more diffusive in character.

Description: We continue to investigate the possibility that interstellar turbulence is caused by nonlinear interactions among shear Alfven waves. Here, we restrict attention to the symmetric case where the oppositely directed waves carry equal energy fluxes. This precludes application to the solar wind in which the outward flux significantly exceeds the ingoing one. All our detailed calculations are carried out for an incompressible magnetized fluid. In incompressible magnetohydrodynamics (MHD), nonlinear interactions only occur between oppositely direct waves. We present a theory for the strong turbulence of shear Alfven waves. It has the following main characteristics. (1) The inertial-stage energy spectrum exhibits a critical balance between linear wave periods and nonlinear turnover timescales. (2) The 'eddies' are elongated in the direction of the field on small spatial scales; the parallel and perpendicular components of the wave vector, k(sub z) and k(perpendicular) are related by k(sub z) approximately equals k(sub perpendicular to)(exp 2/3) L(exp -1/3), where L is the outer scale of the turbulence. (3) The 'one-dimensional' energy spectrum is proportional to k(sub perpendicular)(exp -5/3)-an anisotropic Kolmogorov energy spectrum. Shear Alfvenic turbulence mixes specific entropy as a passive contaminant. This gives rise to an electron density power spectrum whose form mimics the energy spectrum of the turbulence. Radio wave scattering by these electron density fluctuations produces anisotropic scatter-broadened images. Damping by ion-neutral collisions restricts Alfvenic turbulence to highly ionized regions of the interstellar medium.

Description: Two consecutive Hubble Space Telescope (HST) Faint Object Spectrograph (FOS) spectra of the exposed white dwarf in the ultrashort-period, high-amplitude, dwarf nova WZ Sge, reveal a rich absorption line spectrum of neutron carbon and ionized metals, the Stark-broadened Lyman-alpha absorption wing, the H2 quasi-molecular Lyman-alpha 'satellite' absorption line, and a double-peaked C IV emission line which is variable with orbital phase. A synthetic spectral analysis of the white dwarf yields T(sub eff) = 14,900 K +/- 250 K, log g = 8.0. In order to fit the strongest C I absorption lines and account for the weakness of the silicon absorption lines, the abundance of carbon in the photosphere must be approximately 0.5 solar, silicon abundance is 5 x 10(exp -3) solar, with all other metal species appearing to be 0.1-0.001 times solar. The H2 quasi-molecular absorption is fitted very successfully. The photospheric metals have diffusion timescales of fractions of a year, and thus they must have been accreted long after the 1978 December outburst. The source of the most abundance metal, carbon, is considered. If the time-averaged accretion rate during quiescence is low enough for diffusive equilibrium to prevail, then the equilibrium accretion rate pf neutron carbon is 7 x 10(exp -16) solar mass/yr. A convective dredge-up origin for the concentration of carbon is extremely unlikely, given that the white dwarf atmosphere is H-rich while in single degenerates showing carbon and hydrogen, the C and H are trace elements in a helium background. Additional implications are explored.

Description: A new munerical method for computing non-Local Thermodynamic Equilibrium (non-LTE) model stellar atmospheres is presented. The method, called the hybird complete linearization/accelerated lambda iretation (CL/ALI) method, combines advantages of both its constituents. Its rate of convergence is virtually as high as for the standard CL method, while the computer time per iteration is almost as low as for the standard ALI method. The method is formulated as the standard complete lineariation, the only difference being that the radiation intensity at selected frequency points is not explicity linearized; instead, it is treated by means of the ALI approach. The scheme offers a wide spectrum of options, ranging from the full CL to the full ALI method. We deonstrate that the method works optimally if the majority of frequency points are treated in the ALI mode, while the radiation intensity at a few (typically two to 30) frequency points is explicity linearized. We show how this method can be applied to calculate metal line-blanketed non-LTE model atmospheres, by using the idea of 'superlevels' and 'superlines' introduced originally by Anderson (1989). We calculate several illustrative models taking into accont several tens of thosands of lines of Fe III to Fe IV and show that the hybrid CL/ALI method provides a robust method for calculating non-LTE line-blanketed model atmospheres for a wide range of stellar parameters. The results for individual stellar types will be presented in subsequent papers in this series.

Description: The observed boron/beryllium ratio in extreme Population II stars has been interpreted as evidence of Be and B synthesis by early Galactic cosmic rays. However, a recent reanalysis of the boron abundance in the Population II halo star HD 140823 suggests that B/H may be larger than previously reported, by as much as a factor of 4. This would yield a B/Be ratio lying in the range 14 less than or approximately equal to B/Be less than or approximately equal to 50. The possibilty of a high Population II B/Be ratio stresses the importance of the upper limit to the B/Be ratio arising from cosmic-ray production. It is found that the limit to cosmic-ray-produced B/Be depends upon the assumed cosmic-ray spectrum. For any Population II cosmic-ray spectrum that is a single power law in either total energy per nucleon or in momentum, the B/Be ratio constrained to lie in the range 7.6 less than or approximately equal to B/Be less than or approximately equal to 14. Thus, if the new B/Be ratio is correct, it requires either a bimodal cosmic-ray flux with a large low-energy component, or, for another B source, possibly the proposed nu-process in supernovae, either of which may be helpful in explaining the observed B-11/B-10 ratio. Finally, it is noted that the boron reanalysis highlights the uncertainty in our knowledge of the B/Be ratio, and the need for additional data on Be and B abundances.

Description: We make a detailed study of the optical light curve of PSR 1957 + 20. We show that the deep, smooth, and symmetrical modulation can be successfully modeled by a higl irradiated secondary. Two types of models are consistent with the data: (! a secodary close to filling its Rochhe lobe with 10%-20% of the incident flux converted to optical emission, and (2) a secondary considerably underfilling its Roche lobe, with a high degree of beaming of the neutron star flux. owever, the second model canbe rejected on the basis of current estimates of the extinction toward 1957 + 20. The proximity of the photosphere of the secondary to its Roche lobe facilitates mass loss by irradiation. Our ignorance of the (nonirradiated) luminosity of the secondary allows us to place only weak constraints on the inclination -50-80 deg. There is no evidence in the light curve for optical emission from any source in PSR 1957 + 20 other than the secondary itself.

Description: Smoothed particle hydrodynamics (SPH) discretization techniques are generalized to develop a method, smoothed particle interpolation (SPI), for solving initial value problems of systems of a nonhydrodynamical nature. Under this approach, SPH is viewed as strickly as interpolation scheme and, as such, suitable for solving general hyperbolic and parabolic equations. The SPI method is tested on (1) the wave equation with inhomogenous sound speed and (2) Burger's equation. The efficiency of SPI is studied by comparing SPI solutions to those obtained with standard finite-difference methods. It is shown that the power of SPI arises when the smoothing particles are free to move.

Description: Molecular clouds are observed to be partially supported by turbulent pressure. The kinetic energy of the turbulence is directly measurable, but the potential energy, which consists of magnetic, thermal, and gravitational potential energy, is largly unseen. We have extended previous results on equipartition between kinetic and potential energy to show that it is likely to be a very good approximation in molecular clouds. We have used two separate approaches to demonstrate this result: For small-amplitude perturbations of a static equilibrium, we have used the energy principle analysis of Bernstein et al. (1958); this derivation applies to perturbations of arbitary wavelength. To treat perturbations of a nonstatic equilibrium, we have used the Lagrangian analysis of Dewar (1970); this analysis applies only to short-wavelength perturbations. Both analysis assume conservation of energy. Wave damping has only a small effect on equipartition if the wave frequency is small compared to the neutral-ion collision frequency; for the particular case we considered, radiative losses have no effect on equipartition. These results are then incorporated in a simple way into analyses of cloud equilibrium and global stability. We discuss the effect of Alfvenic turbulence on the Jeans mass and show that it has little effect on the magnetic critical mass.

Description: A photometric search for objects with H alpha emission in the globular cluster NGC 6397 has been carried out with the Hubble Space Telescope (HST) Planetary Camera. Images were obtained through the F675W (similar to Johnson R) and F656N (H alpha) filters, and photometry carried out on approximately 900 stars detected in the central approximately 38 x 35 arcseconds of the cluster. Limiting magnitudes of R approximately 21 and 19.5 were reached in the F675W and F656N images, respectively. Three H alpha-bright stars are found which have apparent magnitudes in the range R = 17.8-19.5. The corresponding absolute magnitudes (M(sub R) = 5.6-7.3) and inferred H alpha emission-line strengths (EW(H alpha) greater than or approximately equal to 15-26 A) make these objects good candidates to be cataclysmic variables (CVs) in this nearby, post-collapse cluster. All three H alpha-bright objects are found to be UV-bright in a photometric comparison of the F675W images with archival HST Faint Object Camera (FOC) images of the cluster center through the F220W and F346M filters. Such UV excesses further support the identification of these objects as candidates to be CVs in the cluster. The H alpha-bright stars are within the error circles of three X-ray sources previously found with the ROSAT High Resoltion Imager (HRI) and are likley to be the dominant source of the X-ray emission. The implied X-Ray to optical flux ratios (f(sub x)/f(sub v) approximately 0.8-4.2) are consistent with CVs known elsewhere in the Galaxy. Two additional stars with possible H alpha emission are also discussed, along the stars identified on the basis of UV emisssion. Preliminary completeness estimates suggest that these observations should be sensitive to approximately one-third to two-thirds of most varieties of CVs in the cluster, and that therefore greater than or approximately 5-10 CVs are likely to be present in the observed section of NGC 6397. Scaling these numbers to the cluster as a whole would increase these numbers by up to a factor of 2. A total CV population on the order of 5-20 is in reasonable with previous estimates made on the basis of X-ray observations, and with theoretical predictions scaled to this cluster.

Description: We present Goddard High Resolution Spectrograph intermediate-resolution measurements of the 1233-1256 A spectral region of HD 156396, a halo star at l = 328.7 deg, b = -14.5 deg in the inner Galaxy with a line-of sight distance of 11.1 kpc and a z-distance of -2.8 kpc. The data have a resolution of 18 km/s Full Width at Half Maximum (FWHM) and a signal-to-noise ratio of approximately 50:1. We detect interstellar lines of Mg II, S II, S II, Ge II, and N V and determine log N/(Mg II) = 15.78 +0.25, -0.27, log N(Si II) greater than 13.70, log N(S II) greater than 15.76, log N(Ge II) = 12.20 +0.09,-0.11, and log N(N v) = 14.06 +/- 0.02. Assuming solar reference abundances, the diffuse clouds containing Mg, S, and Ge along the sight line have average logarithmic depletions D(Mg) = -0.6 +/- 0.3 dex, D(S) greater than -0.2 dex, and D(Ge) = -0.2 +/- 0.2 dex. The Mg and Ge depletions are approximately 2 times smaller than is typical of diffuse clouds in the solar vicinity. Galactic rotational modeling of the N v profiles indicates that the highly ionized gas traced by this ion has a scale height of approximately 1 kpc if gas at large z-distances corotates with the underlying disk gas. Rotational modeling of the Si iv and C iv profiles measured by the IUE satellite yields similar scale height estimates. The scale height results contrast with previous studies of highly ionized gas in the outer Milky Way that reveal a more extended gas distribtion with h approximately equals 3-4 kpc. We detect a high-velocity feature in N v and Si II v(sub LSR) approximately equals + 125 km/s) that is probably created in an interface between warm and hot gas.

Description: The unusal supernova remnant CTB, 80 containing the 39.5 ms pulsar PSR 1951 + 32, has been observed with the Position Sensitive Proportional Counter (PSPC) and the High Resolution Imager (HRI) aboard ROSAT. The HRI image, centered on the pulsar, is composed of a bright compact core of approximately 1 arcminute radius containing the pulsar and a compact nebula, as well as a diffuse nebula extending approximately 5 arcminutes eastward of the pulsar. The PSPC allowed us to model the spectra of the point source, the compact nebula and the 5 arcminute diffuse nebula. For a power-law spectrum with photon index Gamma approximately 2 and an interstellar column density of N(sub H) approximately 3 x 10(exp 21)/cm(exp 21) the derived luminosities are approximately 2.3 x 10(exp 33) d(sub 2.5)(exp 2) ergs/s from the pointlike source, approximately 3.9 x 10(exp 33) d(sup 2.5)(exp 2) ergs/s from the compact nebula, and approximately 1.8 x 10(exp 33) d(sub 2.5)(exp 2) ergs/s from the 5 arcminutes diffuse nebula. In addition, the 2 deg diameter circular field of view of the PSPC reveals a hard emission feature southeast of the pulsar with a conical geometry extending out to the edge of the detector. The spectrum from this region is well described by a two-temperature Raymond-Smith thermal plasma with an average temperature of approximately 10(exp 7) K and a luminosity of approximately 10(exp 34)d(sub 2.5)(exp 2) ergs/s. Pulsations from the 39.5 ms puslar, PSR 1951 + 32, are detected at the 99% confidence level. The implied pulsed fration is approximately 35% with a complicated energy-dependent behavior. The compact core and the extended diffuse nebula can be explained as synchrotron radiation from the relativistic pulsar wind confined by the ram pressure of the surrounding inhomgeneous medium. The conelike feature detected southeast of PSR 1951 + 32 is consistent with emission from an optically thin SNR in the radiative cooling phase of its evolution.

Description: We present the results of an optical search for ionized gaseous nebulae surrounding luminous, 'supersoft' X-ray sources in the Large and Small Magellanic Clouds. This relatively new and mysterious X-ray class has characteristic luminosities approximately 10(exp 37)-10(exp 38) ergs/s with effective temperatures in the range of 2-6 x 10(exp 5) K. The presence of a large flux of UV and soft X-ray photons from these objects has led to predictions of bright optical emission lines from the local interstellar medium. One such object, CAL 83 in the LMC, was known to have an associated nebula, and we quantify here the asymmetry and luminosity of this remarkable nebula. Deep images were made using narowband filters to isolate the emission lines of H alpha and (O III) (lambda 5007). In these emission lines, the nebula is detected out to distances as far as 25 pc from the central object, and the integrated luminositu in each line is of order approximately 100 solar luminosity. Model calculations of such nebulae for chemical abundances characteristic of the LMC indicate that approximately 1% of the X-ray luminosity of the central source is reprocessed into the nebular H-alpha and (O III) lambda 5007 emission lines, from which we conclude that the time-averaged X-ray luminosity of the central source, CAL 83, is greater than 3 x 10(exp 37) ergs/s. The bright inner nebula contains approximately 150 solar mass within 7.5 pc of CAL 83, which clearly indicates that the nebular material has its origin in the interstelar medium. In sharp contrast, there were null detections for nebulae associated with nine other luminous, supersoft X-ray sources in the LMC and SMC, with upper limits for the (O III) luminosity that are a factor of approximately 10 below that for CAL 83. For eight of these latter sources, we conclude that either their time-averaged luminosity is substantially below that of CAL 83, or that the local interstellar medium is much less dense. The latter effect may be enhanced by expansion of the ionized nebula during the last several million years.

Description: We present broad-band J, H, and K images and K-band spectroscopy of the luminous starburst galaxy UGC 8387. The images show a disturbed morphology, tidal tails, and a single elognated nucleus. Near infrared color maps constructed from the images reveal that the nucleus region is highly reddened. Strong emission from the central 3 arcseconds in the 2.166 micrometer Brackett gamma, 2.122 micrometer H2 v = 1-0 S(1), and 2.058 micrometer He I lines is present in the K-band spectrum. From the Brackett gamma and published radio fluxes, we find an optical depth toward the nucleus of tau(sub V) approximately 24. The CO band heads produce strong absorption in the spectral region long-ward of 2.3 micrometers. We measure a 'raw' CO index of 0.17 +/- 0.02 mag, consistent with a population of K2 supergiants of K4 giants. The nuclear colors, however, are not consistent with an obscured population of evolved stars. Instead, the red colors are best explained by an obscured mixture of stellar and warm dust emission. The amount of dust emission predicted by the near-infrared colors exceeds that expected from comparisons to galactic H II regions. After correcting the spectrum of UGC 8387 for dust emission and extinction, we obtain a CO index of greater than or equal to 0.25 mag. This value suggests the stellar component of the 2.2 micrometer light is dominated by young supergiants. The infrared excess, L(sub IR)/L(sub Ly alpha) derived for UGC 8387 is lower than that observed in galactic H II regions and M82. This implies that either the lower or upper mass cutoff of the initial mass function must be higher than those of local star-forming regions and M82. The intense nuclear starburst in this galaxy is presumably the result of merger activity; and we estimate the starburst age to be at least a few times 10(exp 7) yr.

Description: We have measured the T(sub 90) and T(sub 50) durations of bright and dim gamma-ray bursts detected by the Compton Gamma-Ray Observatory's (CGRO) Burst and Transient Source Experiment (BASTE). The T(sub 90) T(sub 50) duration is defined as the interval over which 5% (25%) to 95% (75%) of the burst counts accumulate. Out of 775 bursts observed by BATSE 159 bursts were analyzed; bursts with durations shorter than 1.5 s were excluded. A Kolmogorov-Smirnov test yields a probability of 6 x 10(exp -5) that the T(sub 50) durations of the dim and bright samples are drawn from the same parent population. We find that the centroid and extent of the duration distribution for the dim sample are scaled by approximately a factor of 2 relative to those of the bright sample. The measured time-dilation factor is not sensitive to choice of energy band. These results are quantitatively consistent with previous tests for time dilation in a smaller sample of BATSE bursts. The sources of dimmer bursts, if cosmological, would lie at redshifts of order 2.

Description: The global nonlinear time-dependent evolution of the inertial-acoustic mode instability in accretion disks surrounding black holes has been investigated. The viscous stress is assumed to be proportional to the gas pressure only, i.e., tau = alphap(sub g). It is found that an oscillatory nonsteady behavior exists in the inner regions of disks (r is less than 10r(sub g) where r(sub g) is the Schwarzschild radius) for sufficiently large alpha(greater than or approximately equal to 0.2) and for mass accretion rates less than about 0.3 times the Eddington value. The variations of the integrated bolometric luminosity from the disk, Delta L/L, are less than 3%. A power spectrum analysis of these variations reveals a power spectrum which can be fitted to a power-law function of the frequency Pis proportional to f(exp -gamma), with index gamma = 1.4-2.3 and a low-frequency feature at about 4 Hz in one case. In addition, a narrow peak centered at a frequency corresponding to the maximum epicyclic frequency of the disk at approximately 100-130 Hz and its first harmonic is also seen. The low-frequency modulations are remarkably similar to those observed in black hole candidate systems. The possible existence of a scattering corona in the inner region of the disk and/or other processes contributing to the power at high frequencies in the inner region of the accretion disk may make the detection of the high-frequency component difficult.

Description: A model of the approximately 0.1-10 GeV gamma-ray jets observed by the EGRET instrument on the Compton Gamma Ray Observatory (CGRO) is developed. It is shown that the soft X-ray background in an active galactic nuclei (AGN) contributes an opacity to pair production and that a gamma-ray photosphere or 'gamma-sphere' can be defined whose radius increases with gamma-ray energy E(sub gamma). It is proposed that the observed gamma-ray emission is due to inverse Compton scattering of the ambient soft X-rays by relativistic pairs accelerated in situ by shock fronts in a relativistic jet. For a wide range of assumed physical conditions, the emission at a given E(sub gamma) originates from near the associated gamma-spheres; emission from below the gamma-sphere initiates a cascade down to the energy where the gamma-rays can escape freely. In this model, the slope of the emergent gamma-ray spectrum is determined by the scattered, soft X-ray spectrum and the variation of the particle acceleration rate with jet radius. In general it is expected that the variation in the gamma-ray flux will be either slower or later at higher energy. It is also shown that the efficiency of conversion of energy from injected high-energy pairs to 0.1-10 GeV gamma-rays is typically high so that the models are radiatively efficient. It is argued that the observed gamma-ray jets are likely to be particle-dominated, though magnetically confined. The gamma-ray spectrum should continue down to an energy approximately 5 MeV emitted from an annihilation radius within which the pair content of the jet is limited by annihilation. This is probably the site of the beamed hard X-ray emission. It is speculated that the relativistic jets associated with radio-loud AGNs are powered electromagnetically by a spinning black hole and that they are collimated by an encircling MHD wind leaving the accretion disk at a slower speed. Powerful FR2 radio sources are formed when the hole spins rapidly and the relativistic core accelerates the MHD sheath; low-power FR1 sources ensue when the opposite occurs. Finally, it is suggested that the key factor which determines whether or not a given active nucleus can form a jet and a radio to gamma-ray nonthermal continuum is the central density of mass-losing stars which, when large, precludes the formation of a super-Alfvenic, collimating wind.

Description: The origin of gamma-ray bursts continues to be a great mystery. Here we review some relevant observations and a number of recent models. While no clear solution exists at the present time, the parameter space for Galactic halo models is becoming very constrained. Cosmological models on the other hand require both enormous total energy and the concentration of that energy into a small mass. This implies compact objects, probably accreting black holes of stellar size. We review some of the physics of such accreting black holes and point out that the jet formed from accretion into a rapidly accreting black hole of stellar mass might precess. This precession, coupled to beaming, could impose additional time structure on the burst and its spectrum. In the event of a 'failed' supernova model, the wind of the Wolf-Rayet star prior to the event could provide the beam dump where the jet generates gamma-rays. Enduring emission that grows harder with time might be expected for several hours as the density in the vicinity of the black hole declines.

Description: We show that ejecta from the envelope of one asymptotic giant branch star of M is approximately 3 solar mass may account for many of the short-lived nuclei in the early solar system and also for the recent evidence of the presence of Ca-41 (bar-tau(sub 41) = 1.50 x 10(exp 5) yr) in early solar nebular condensates. This would require that the injection into the protosolar molecular cloud took place within a narrow time interval of (5-7) x 10(exp 5) yr before the formation of the solar system. If true, this places extremely tight constraints on the whole process of injection mixing and collapse. The timescales for both Ca-41 and Al-26 require that the placental medium be a dense molecular cloud (2 x 10(exp 3) -8 x 10(exp 3( H/cc). If the observed residual Ca-41 is instead produced by a proton bombardment mechanism within the early solar system, similar to what appears necessary to explain Mn-53, then the time interval is relaxed but would still be (1-2) x 10(exp 6) yr from consideration of Al-26.

Description: We present Hubble Space Telescope (HST) and ground-based data on the Z(sub abs) = 0.8596 metal-line absorption system along the line of sight to PKS 0454+0356. The system is a moderate-redshift damped Ly-alpha system, with N(H I) = (5.7 +/- 0.3) x 10(exp 20)/sq cm as measured from the Faint Object Spectrograph (FOS) spectrum. We also present ground-based images which we use to identify the galaxy which most probably gives rise to the damped system; the most likely candidate is relatively underluminous by QSO absorber standards M(sub B) approximately -19.0 for A(sub 0) = 0.5 and H(sub 0) = 50 km/s/Mpc) and lies approximately 8.5/h kpc in projection from the QSO sight line. Ground-based measurements of Zn II, Cr II, and Fe II absorption lines from this system allow us to infer abundances of (Zn/H) = -1.1, (Cr/H) = -1.2, and (Fe/H) = -1.2 indicating overall metallicity similar to damped systems at z is greater than 2, and that the depletion of Cr and Fe onto dust grains may be even less important than in many of the high-redshift systems of comparable metallicity. Limits previously placed on the 21 cm optical depth in the z = 0.8596 system, together with our new N(H I) measurement, suggest a very high spin temperature for the H I, T(sub s) is greater than 580 K.

Description: Enhanced chromospheric emission, which corresponds to an outwardly increasing semiempirical temperature structure, can be produced by wave motion without any increase in the mean gas temperature. Hence, the Sun may not have a classical chromosphere in magnetic field-free internetwork regions. Other significant differences between the properties of dynamic and static atmospheres should be considered when analyzing chromospheric observations.

Description: It is suggested that tidal interaction of an accreting planetary embryo with the gaseous preplanetary disk may provide a mechanism to breach the so-called runaway limit during the formation of the giant planet cores. The disk tidal torque converts a would-be shepherding object into a 'predator,' which can continue to cannibalize the planetesimal disk. This is more likely to occur in the giant planet region than in the terrestrial zone, providing a natural cause for Jupiter to predate the inner planets and form within the O(10(exp 7) yr) lifetime of the nebula.

Description: We compute the optical depth and duration distribution of microlensing events towrd Baade's window in a model composed of a Galactic disk and a bar. The bar model is a self-consistent dynamical model built out of individual orbits that has been populated to be consistent with the COBE maps of the Galaxy and kinematic observations of the Galactic bulge. We find that most of the lenses are in the bulge with a line-of-sight distance 6.25 kpc (adopting R(sub 0) = 8 kpc). The microlensing optical depth of a 2 x 10(exp 10) solar mass bar plus a truncated disk is (2.2 +/- 0.45) x 10(exp -6), consistent with the large optical depth (3.2 +/- 1.2) x 10(exp -6) found by Udalski et al. (1994). This model optical depth is enhanced over the predictions of axisymmetric models by Kiraga & Paczynski (1994) by slightly more than a factor of 2, since the bar is elongated along the line of sight. The large Einstein radius and small transverse velocity dispersion also predict a longer event duration in the self-consistent bar model than in the Kiraga-Paczynski model. The event rate and duration distribution also depend on the lower mass cutoff of the lens mass function. With a 0.1 solar mass cutoff, five to seven events (depending on the contribution of disk lenses) with a logarithmic mean duration of 20 days are expected for the Optical Gravitational Lensing Experiment (OGLE) according to our model, while Udalski et al. (1994) observed nine events with durations from 8 to 62 days. On the other hand, if most of the lenses are brown dwarfs, our model predicts too many short-duration events. A Kolmogorov-Smirnov test finds only 7% probability for the model with 0.01 solar mass cutoff to be consistent with current data.

Description: We develop an analytical theory of reconnection between colliding, twisted magnetic flux tubes. Our analysis is restricted to direct collisions between parallel tubes and is based on the collision dynamics worked out by Bogdan (1984). We show that there is a range of collision velocities for which neutral point reconnection of the Parker-Sweet type can occur, and a smaller range for which reconnection leads to coalescence. Mean velocities within the solar convection zone are probably significantly greater than the upper limit for coalescence. This suggests that the majority of flux tube collisions do not result in merging, unless the frictional coupling of the tubes to the background flow is extremely strong.

Description: Two theories of stellar envelope convection are considered here in the context of red giants and red supergiants of intermediate to high mass: Boehm-Vitense's standard mixing-length theory (MLT) and Canuto & Mazzitelli's new theory incorporating the full spectrum of turbulence (FST). Both theories assume incompressible convection. Two formulations of the convective mixing length are also evaluated: l proportional to the local pressure scale height (H(sub P)) and l proportional to the distance from the upper boundary of the convection zone (z). Applications to test both theories are made by calculating stellar evolutionary sequences into the red zone (z). Applications to test both theories are made by calculating stellar evolutionary sequences into the red phase of core helium burning. Since the theoretically predicted effective temperatures for cool stars are known to be sensitive to the assigned value of the mixing length, this quantity has been individually calibrated for each evolutionary sequence. The calibration is done in a composite Hertzsprung-Russell diagram for the red giant and red supergiant members of well-observed Galactic open clusters. The MLT model requires the constant of proportionality for the convective mixing length to vary by a small but statistically significant amount with stellar mass, whereas the FST model succeeds in all cases with the mixing lenghth simply set equal to z. The structure of the deep stellar interior, however, remains very nearly unaffected by the choices of convection theory and mixing lenghth. Inside the convective envelope itself, a density inversion always occurs, but is somewhat smaller for the convectively more efficient MLT model. On physical grounds the FST model is preferable, and seems to alleviate the problem of finding the proper mixing length.

Description: Modern stellar evolutionary tracks are used to calculate the evolution of a very large number of massive binary star systems (M(sub tot) greater than or = 15 solar mass) which cover a wide range of total masses, mass ratios, and starting separations. Each binary is evolved accounting for mass and angular momentum loss through the supernova of the primary to the X-ray binary phase. Using the observed rate of star formation in our Galaxy and the properties of massive binaries, we calculate the expected high-mass X-ray binary (HMXRB) population in the Galaxy. We test various massive binary evolutionary scenarios by comparing the resulting HMXRB predictions with the X-ray observations. A major goal of this study is the determination of the fraction of matter lost from the system during the Roche lobe overflow phase. Curiously, we find that the total numbers of observable HMXRBs are nearly independent of this assumed mass-loss fraction, with any of the values tested here giving acceptable agreement between predicted and observed numbers. However, comparison of the period distribution of our HMXRB models with the observed period distribution does reveal a distinction among the various models. As a result of this comparison, we conclude that approximately 70% of the overflow matter is lost from a massive binary system during mass transfer in the Roche lobe overflow phase. We compare models constructed assuming that all X-ray emission is due to accretion onto the compact object from the donor star's wind with models that incorporate a simplified disk accretion scheme. By comparing the results of these models with observations, we conclude that the formation of disks in HMXRBs must be relatively common. We also calculate the rate of formation of double degenerate binaries, high velocity detached compact objects, and Thorne-Zytkow objects.

Description: We calculate the radiative driving force for winds around rapidly rotating oblate B stars, and we estimate the impact these forces should have on the production of a wind compressed disk. The effects of limb darkening, gravity darkening, oblateness, and an arbitrary wind velocity field are included in the computation of vector 'oblate finite disk' (OFD) factors, which depend on both radius and colatitude in the wind. The impact of limb darkening alone, with or without rotation, can increase the mass loss by as much as 10% over values computed using the standard uniformly bright spherical finite disk factor. For rapidly rotating stars, limb darkening makes 'sub-stellar' gravity darkening the dominant effect in the radial and latitudinal OFD factors, and lessens the impact of gravity darkening at other visible latitudes (nearer to the oblate limb). Thus, the radial radiative driving is generally stronger over the poles and weaker over the equator, following the gravity darkening at these latitudes. The nonradial radiative driving is considerably smaller in magnitude than the radial component, but is directed both away from the equatorial plane and in a retrograde azimuthal direction, acting to decrease the effective stellar rotation velocity. These forces thus weaken the equatorward wind compression compared to wind models computed with nonrotating finite disk factors.

Description: We present a cosmological model with good physical properties which is invariant not only under changes of the space and time coordinates but also under changes of an extra (Kaluza-Klein) coordinate related to rest mass. In frames where the latter is chosen to be constant we recover standard cosmology. In frames where it is chosen to be variable we obtain new astrophysical effects and gain insight into the nature of the big bang.

Description: We present the analysis of 13 orbital phase-resolved IUE (large and small aperture) SWP spectra, sampling the quiescent intervals of the proptotype dwarf nova U Geminorum following two individual outbursts, the outburst of 1992, 29 August-14 September, and the outburst of 1993, 19 March-5 April. The quiescent interval following the outburst 29 August 1992 is the longest ever recorded. During quiescence, the photospheric radiation of the exposed white dwarf dominates the far ultrviolet. The variations in absorption line strengths and continum flux levels are analyzed as a function of both orbital phase and elapsed time, since the return to optical quiescence, by fitting the IUE spectra with a grid of high gravity, LTE, model atmospheres with solar composition constructed with TLUSTY and SYNSPEC (Hubeny, I. Tlusty and Synnspec: A Users Guide (1994)). We present evidence from both absorption line variations and continuum variations, as a function of time since the outbursts, that the white dwarf photosphere has cooled by several thousand degrees. Within the signal-to-noise limitations, we find no evidence of a difference in heating and cooling between the two outbursts.

Description: We explore the kinematics of 36 rich RGH89 groups identified from the first two complete slices of the CfA redshift survey. These groups have more than five members identified by a friends-of-friends algorithm at a number density contrast delta rho/rho greater than or equal to 80. To examine the stability of the determination of the velocity dispersion for these systems, we compare results for the original 232 members with results for a larger redshift sample, including 334 fainter members in the redshift neighborhoods. On average, we double the number of group members in each system. The observed distribution of velocity dispersions is stable. In fact, the velocity dispersion based on the original members identified in the CfA redshift survey is a reliable predictor of the value for the enlarged sample in an individual group. The velocity dispersion is thus a stable physical parameter for discrimination among systems galaxies. A larger sample of groups, particularly one selected from a distance limited catalog, should provide an interesting constraint on models for the formation of large-scale structure. We take H(sub 0) = km/s/Mpc.

Description: The unusual variable star AM CVn has puzzled astronomers for over 40 years. This object, both a photometric and spectroscopic variable, is believed to contain a pair of hydrogen-deficient white dwarfs of extreme mass ratio, transferring material via an accretion disk. We examine the photometric properties of AM CVn, analyzing 289 hours of high-speed photometric data spanning 1976 to 1992. The power spectrum displays significant peaks at 988.7, 1248.8, 1902.5, 2853.8, 3805.2, 4756.5, and 5707.8 microHz (1011.4, 800.8, 525.6, 350.4, 262.8, 210.2, and 175.2 s). We find no detectable power at 951.3 microHz (1051 s), the previously reported main frequency. The 1902.5, 2853.9, and 3805.2 microHz peaks are multiplets, with frequency splitting in each case of 20.77 +/- 0.05 microHz. The 1902.5 microHz seasonal pulse shapes are identical, within measurement noise, and maintain the same amplitude and phase as a function of color. We have determined the dominant frequency to be 1902.50902 +/- 0.00001 microHz with dot P = +1.71 (+/- 0.04) x 10(exp -11) s/s. We discuss the implications of these findings on a model for AM CVn.

Description: We present a three-dimensional hydrodynamic simulation of the disrupted stellar wind in the high-mass X-ray binary system SMC X-1. The three dominant processes that determine the geometry of the wind in high X-ray luminosity systems such as SMC X-1 are the X-ray suppression of the stellar wind from the X-ray irradiated face of the primary star, the focusing of the radiatively driven wind in the X-ray shadow by the effects of stellar rotation, and the rapid X-ray heating of gas in the vicinity of the X-ray source, including the X-ray illuminated surface of the primary star. The resulting distribution of circumstellar gas provides a successful explanation for the asymmetric, extended eclipse transitions and the intensity of the deep eclipse X-ray emission in SMC X-1, as well as a possible explanation for the X-ray dips seen near superior conjunction of the X-ray source in Cyg X-1.

Description: Most physical properties derived for quasars, as single entities or as a population, depend upon the cosmology assumed. In this paper, we calculate the quasar luminosity function and some related quantities for a flat universe dominated by a cosmological constant Lambda (Lambda = 0.9, Omega = 0.1) and compare them with those deduced for a flat universe with zero cosmological constant (Lambda = 0, Omega = 1). We use the ATT quasar survey data (Boyle et al. 1990) as input in both cases. The data are fitted well by a pure luminosity evolution model for both the cosmologies but with different evolutionary parameters. From the luminosity function, we predict (extrapolate) a greater number of quasars at faint apparent magnitudes (twice the number at B = 24, z is less than 2.2) for the Lambda-dominated universe. This population of faint quasars at high redshift would result in a higher incidence of gravitational lensing. The total luminosity of the quasar population and the total mass tied up in black hole remnants of quasars is not sensitive to the cosmology. However, for a Lambda cosmology, this mass is tied up in fewer but more massive black holes.

Description: This paper describes the formation of nonlinear structure in flat (zero curvature) Friedmann cosmological models. We consider models with two components: the usual nonrelativistic component that evolves under gravity and eventually forms the large-scale structure of the universe, and a uniform dark matter component that does not clump under gravity, and whose energy density varies with the scale factor a(t) like a(t)(sup -n), where n is a free parameter. Each model is characterized by two parameters: the exponent n and the present density parameter Omega(sub 0) of the nonrelativistic component. The linear perturbation equations are derived and solved for these models, for the three different cases n = 3, n is greater than 3, and n is less than 3. The case n = 3 is relevant to model with massive neutrinos. The presence of the uniform component strongly reduces the growth of the perturbation compared with the Einstein-de Sitter model. We show that the Meszaros effect (suppression of growth at high redshift) holds not only for n = 4, radiation-dominated models, but for all models with n is greater than 3. This essentially rules out any such model. For the case n is less than 3, we numerically integrate the perturbation equations from the big bang to the present, for 620 different models with various values of Omega(sub 0) and n. Using these solutions, we show that the function f(Omega(sub 0), n) = (a/delta(sub +))d(delta)(sub +)/da, which enters in the relationship between the present density contrast delta(sub 0) and peculiar velocity field u(sub 0) is essentially independent of n. We derive approximate solutions for the second-order perturbation equations. These second-order solutions are tested against the exact solutions and the Zel'dovich approximation for spherically symmetric perturbations in the marginally nonlinear regime (the absolute value of delta is less than or approximately 1). The second-order and Zel'dovich solutions have comparable accuracy, significantly higher than the accuracy of the linear solutions. We then investigate the dependence of the delta(sub 0) - u(sub 0) relationship upon the value of n in the nonlinear regime using the second-order solutions for marginally nonlinear, general perturbations, and the exact solutions for strongly nonlinear, spherically symmetric perturbations. In both cases, we find that the delta(sub 0) - u(sub 0) relationship remains independent of n. We speculate that this result extends to strongly nonlinear, general perturbations as well. This eliminates any hope to determine the presence of the uniform component or the value of n using dynamical methods. Finally, we compute the nonlinear evolution of the skewness of the distribution of values of delta, assuming Gaussian initial conditions. We find that the skewness is not only independent of n, but also of Omega(sub 0). Thus the skewness cannot be used to discriminate among various models with Gaussian initial conditions. However, it can be used for testing the Gaussianity of the initial conditions themselves. We conclude that the uniform component leaves no observable signature in the present large-scale structure of the universe. To determine its presence and nature, we must investigate the relationship between the past and present universe, using redshift-dependent tests.

Description: An upper bound to the energy density of infrared background radiation is derived from considering the effect of gamma-gamma interactions on the observed TeV gamma-ray spectrum of the active galaxy Markarian 421. This upper bound proves to be the most restrictive for the wavelength range of 10-12 micrometers. These constraints are presently limited by the uncertainty of extrapolating the source spectrum from 5 to 500 GeV. Observations in the regime less than 100 GeV would significantly improve these limits, allowing for a wide range of IR production models to be constrained.

Description: The intracluster medium (ICM) within merging clusters of galaxies is likely to be in a violent or turbulent dynamical state which may have a significant effect on the evolution of cluster radio sources. We present results from a recent gas + N-body simulation of a cluster merger, suggesting that mergers can result in long-lived, supersonic bulk flows, as well as shocks, within a few hundred kiloparsecs of the core of the dominant cluster. These results have motivated our new two-dimensional and three-dimensional simulations of jet propagation in such environments. The first set of simulations models the ISM/ICM transition as a contact discontinuity with a strong velocity shear. A supersonic (M(sub j) = 6) jet crossing this discontinuity into an ICM with a transverse, supersonic wind bends continuously, becomes 'naked' on the upwind side, and forms a distended cocoon on the downwind side. In the case of a mildly supersonic jet (M(sub j) = 3), however, a shock is driven into the ISM and ISM material is pulled along with the jet into the ICM. Instabilities excited at the ISM/ICM interface result in the jet repeatedly pinching off and reestablishing itself in a series of 'disconnection events.' The second set of simulations deals with a jet encountering a shock in the merging cluster environment. A series of relatively high-resolution two-dimensional calculations is used to confirm earlier analysis predicting that the jet will not disrupt when the jet Mach number is greater than the shock Mach number. A jet which survives the encounter with the shock will decrease in radius and disrupt shortly thereafter as a result of the growth of Kelvin-Helmholtz instabilities. We also find, in disagreement with predictions, that the jet flaring angle decreases with increasing jet density. Finally, a three-dimensional simulation of a jet crossing an oblique shock gives rise to a morphology which resembles a wide-angle tailed radio source with the jet flaring at the shock and disrupting to form a long, turbulent tail which is dragged downstream by the preshock wind.

Description: We reanalyzed the CN images of Comet P/Halley, in which jets have been discovered for the first time, in search of shell structures. Shells were actually detected at the outer edges of the frames on those dates for which shells with radii small enough to be covered by the limited field of view of the CCD were predicted. The C2 images of the same data set were subjected to an analogous investigation which led to the discovery of shell structures in C2 as well. The morphology of the CN and the C2 shells is essentially equal on the same observational date. They have the same radii and show almost identical asymmetries which suggests that CN and C2 in the shells originate from the same general source. The comparison of the jets in both species before and after a two-dimensional continuum subtraction supports this supposition. The similar morphology of the jets indicates that both species are produced from the same bulk of precursor material which has been ejected in the form of jets from the same active area. However, similarly located and oriented jets in CN and C2 do not show similar relative intensities in most cases. These differences in the intensity distribution imply that the production rates of CN and C2 follow different laws.

Description: Asymptotic giant branch (AGB) stars are known to be sites of dust formation and undergo significant mass loss. The outflow is believed to be driven by radiation pressure on grains and momentum coupling between the grains and gas. While the physics of shell dynamics and grain formation are closely coupled, most previous models of circumstellar shells have treated the problem separately. Studies of shell dynamics typically assume the existence of grains needed to drive the outflow, while most grain formation models assume a constant veolcity wind in which grains form. Furthermore, models of grain formation have relied primarily on classical nucleation theory instead of using a more realistic approach based on chemical kinetics. To model grain formation in carbon-rich AGB stars, we have coupled the kinetic equations governing small cluster growth to moment equations which determine the growth of large particles. Phenomenological models assuming stationary outflow are presented to demonstrate the differences between the classical nucleation approach and the kinetic equation method. It is found that classical nucleation theory predicts nucleation at a lower supersaturation ratio than is predicted by the kinetic equations, resulting in significant differences in grain properties. Coagulation of clusters larger than monomers is unimportant for grain formation in high mass-loss models but becomes more important to grain growth in low mass-loss situations. The properties of the dust grains are altered considerably if differential drift velocities are ignored in modeling grain formation. The effect of stellar temperature, stellar luminosity, and different outflow velocities are investigated. The models indicate that changing the stellar temperature while keeping the stellar luminosity constant has little effect on the physical parameters of the dust shell formed. Increasing the stellar luminosity while keeping the stellar temperature constant results in large differences in grain properties. For small outflow velocities, grains form at lower supersaturation ratios and close to the stellar photosphere, resulting in larger but fewer grains. The reverse is true when grains form under high outflow velocities, i.e., they form at higher supersaturation ratios, farther from the star, and are much smaller but at larger quantities.

Description: The limb to limb Lyman-alpha reflectivities observed with the Voyager ultraviolet spectrometer (UVS) instruments during the fly-by of Saturn are reanalyzed using a revised H Lyman-alpha sensitivity for the Voyager 1 instrument. The new sensitivity reconciles the measured intensities to those of Voyager 2 and gives a coherent set of data. To fit the UV airglow observations, four sources are considered: (1) H resonance and H2 Rayleigh scattering of solar Lyman-alpha radiation, (2) the interplanetary Lyman-alpha radiation, (3) a possible internal source of unknown origin, (4) the possibility of atmospheric turbulence recently proposed to explain the Lyman-alpha bulge of Jupiter. The analysis supports neither a dominant collisional excitation source for the UV emissions nor the presence of strong atmospheric turbulence. The best fit, in terms of brightness but also in terms of shape of the limb to limb profile (that is to say independent on the absolute calibrations), is obtained for pure resonance and Rayleigh scattering of solar and interstellar wind line in an atmosphere enriched in atomic hydrogen up to three times the standard model. Influx of water from the rings of Saturn may provide a means for producing such enhanced H densities in the upper atmosphere.

Description: The spacecraft and ground based observations of comet Halley inner coma showed a localized ion density depletion region whose origin is not well understood. Although it has been linked to a thermal instability associated with negative ions, the photodetachment lifetime of negative ions (approximately 1 sec) is too short compared to the electron attachment time scale (approximately 100 sec) for this process to have a significant effect. A mechanism for the ion density depletion based on the thermal instability of the cometary plasma due to the excitation of rotational and vibrational levels of water molecules is proposed. The electron energy losses due to these processes peak near 4000 K (0.36 eV) and at temperatures higher than this value a localized cooling leads to further cooling (thermal instability) due to the increased radiation loss. The resulting increase in recombination leads to an ion density depletion and the estimates for this depletion at comet Halley agree with the observations.

Description: Using the output from a grid of 60 Type II supernova models (Woosley & Weaver 1995) of varying mass (11 approx. less than (M/solar mass) approx. less than 40) and metallicity (0, 10(exp -4), 0.01, and 1 solar metallicity), the chemical evolution of 76 stable isotopes, from hydrogen to zinc, is calculated. The chemical evolution calculation employs a simple dynamical model for the Galaxy (infall with a 4 Gyr e-folding timescale onto a exponential dsk and 1/r(exp 2) bulge), and standard evolution parameters, such as a Salpeter initial mass function and a quadratic Schmidt star formation rate. The theoretical results are compared in detail with observed stellar abundances in stars with metallicities in the range -3.0 approx. less than (Fe/H) approx. less than 0.0 dex. While our discussion focuses on the solar neighborhood where there are the most observations, the supernova rates, an intrinsically Galactic quality, are also discussed.

Description: A recent ultraviolet snaphsot imaging survey of the nuclei of nearby galaxies detected a compact nuclear ultraviolet source in only five of the 26 LINERs (low-ionization nuclear emission-line regions) included in the observed sample. Motivated by this observational result, we examine the possibility that all LINERs are powered by photoionization from a nuclear source, which is, however, active only for 20% of the time. We show that decay times of low-ionization species can be of the order of one to a few centuries, and we demonstrate through time-dependent photoionization calculations that if the nuclear ionizing source is active for only a fraction of the time, this would not be readily noticeable in the emission-line spectrum. We suggest that the activity cycle is related to episodic accretion events which are associated with the tidal disruption of stars by a central black hole. The time interval between tidal disruptions is of the same order as the emission-line decay time, with the accretion episode following each disruption lasting a few decades. These estimates appear to support the duty cycle hypothesis. Some observational consequences of the proposed scenario are also discussed.