title: A Search for X-Ray Bright Distant Clusters of Galaxies

abstract: We present the results of a search for X--ray luminous distant clusters of galaxies. We found extended X--ray emission characteristic of a cluster towards two of our candidate clusters of galaxies. They both have a luminosity in the ROSAT bandpass of $\simeq10^{44}{\rm \,erg\,s^{-1}}$ and a redshift of $>0.5$; thus making them two of the most distant X--ray clusters ever observed. Furthermore, we show that both clusters are optically rich and have a known radio source associated with them. We compare our result with other recent searches for distant X--ray luminous clusters and present a lower limit of $1.2\times10^{-7}\,{\rm Mpc^{-3}}$ for the number density of such high redshift clusters. This limit is consistent with the expected abundance of such clusters in a standard (b=2) Cold Dark Matter Universe. Finally, our clusters provide important high redshift targets for further study into the origin and evolution of massive clusters of galaxies. Accepted for publication in the 10th September 1994 issue of ApJ.


title: "The Contribution of Active Galactic Nuclei to the Diffuse Gamma Ray Background"

abstract: We consider the emission of high energy to very high energy $\gamma$-rays in radio-quiet active galactic nuclei (AGN) or the central regions of radio-loud AGN. We use our results to estimate the $\gamma$-ray flux from the central regions of nearby AGN, and then to calculate the contribution to the diffuse $\gamma$-ray flux from unresolved AGN.


title: Kurtosis as a Non--Gaussian Signature of the Large--Scale Velocity Field

abstract: We discuss the non--linear growth of the kurtosis of the smoothed peculiar velocity field (along an arbitrary direction), in an Einstein--de Sitter universe, induced by Gaussian primordial density fluctuations. Applying the perturbative theory, we show that, for different cosmological models, a departure from the original Gaussian distribution is gravitationally induced only on small scales ($\simlt 20-30$ Mpc). Models with scale--free power spectrum $P(k)\propto k^{\,n}$, with $-1< n \leq 1$, are also considered. When the fluid particles move according to the Zel'dovich approximation the probability distribution of the peculiar velocity field remains unaltered during the evolution.


title: Radio and Infrared Emission from a [WC]-Type Planetary Nebula in the LMC

abstract: Radio continuum emission has been detected from a planetary nebula in the Large Magellanic Cloud: this is the first radio continuum detection for any planetary nebula outside our galaxy. The radio flux density is a factor of two lower than predicted from the \hbeta\ flux. This could be due either to a two-component morphology or to a stellar contribution to the \hbeta\ emission. We have modeled the optical and infrared spectrum using the photo-ionisation code Cloudy: a very good fit is obtained if we assume the nebula is carbon rich. The derived diameter implies an evolutionary age of the nebula of $<1000\yr$, similar to what is derived from the IRAS colours alone. The central star, which is of [WC] type, has a much higher stellar temperature than Galactic [WC] stars showing similar circumstellar IRAS colours. An explanation could be that the expansion velocity of the nebula is lower than those of its galactic counterparts. This radio detection indicates that accurate nebular luminosities could be determined from their radio emission for many LMC planetary nebulae.


title: Confrontation of Double-Inflationary Models with Observations

abstract: We consider double-inflationary models with two noninteracting scalar fields, a light scalar field $\phi_l$ with potential $\frac{1}{2}m_l^2\phi_l^2$ and a heavy scalar field $\phi_h$ with potential $\frac{\lambda}{n}\phi_h^n$ with $n=2,4$. CDM with the initial spectrum of adiabatic perturbations produced in these models is compared with observations. These models contain two more free parameters than the standard CDM model with an initial scale-invariant spectrum. We normalize our spectra to COBE DMR and compare the predictions with observations on the biasing factor, large-scale peculiar velocities, quasar and galaxy formation and the Stromlo-APM counts-in-cells analysis. The model with $n=4$ is excluded by the data while for the $n=2$ model, taking cosmic variance into account, a small window of parameters compatible with observations is found.


title: An Axisymmetric Distribution Function for the Galactic Bulge

abstract: We describe a method for parameterizing two-integral distribution functions, based on triangular tesselations of the integral plane. We apply the method to the axisymmetric isotropic rotator model for the Galactic bulge of Kent~(1992), and compare the results with observations of proper motions in Baade's Window, and with radial velocity surveys. In spite of mounting evidence from surface photometry and from study of the gas kinematics that the Galactic bulge is not axisymmetric, the stellar kinematics in Baade's Window are very similar to those of an isotropic oblate rotator. Another field at large radius does not fit this model, though. In any case, the edge-on kinematics of a hot stellar population are a poor handle on the existence or otherwise of a bar.


title: On the Origin of Highest Energy Cosmic Rays

abstract: In this paper we show that the conventional diffusive shock acceleration mechanism for cosmic rays associated with relativistic astrophysical shocks in active galactic nuclei (AGNs) has severe difficulties to explain the highest energy cosmic ray events. We show that protons above around $2\times10^{20}\eV$ could have marginally been produced by this mechanism in an AGN or a rich galaxy cluster not further away than around $100\Mpc$. However, for the highest energy Fly's Eye and Yakutsk events this is inconsistent with the observed arrival directions. Galactic and intergalactic magnetic fields appear unable to alter the direction of such energetic particles by more than a few degrees. We also discuss some other options for these events associated with relativistic particles including pulsar acceleration of high $Z$ nuclei. At the present stage of knowledge the concept of topological defects left over from the early universe as the source for such events appears to be a promising option. Such sources are discussed and possible tests of this hypothesis are proposed.


title: Anisotropies of the Cosmic Background Radiation in a Reionized Universe

abstract: We trace the evolution of cosmic microwave background photons propagating through a reionized model universe. The reionization of the intergalactic medium is achieved by UV photons emitted from the decaying `hot' dark matter neutrinos. The model universe is represented by a comoving cube through which we propagate photons and subject them to the Sachs-Wolfe effect and Thomson scattering off the ionized fraction of the IGM, beginning at z=1900 until the present z=0. The simulation follows the evolving matter inhomogeneities and microwave background photons into the present nonlinear regime and yields temperature maps of the microwave sky as a comoving observer would detect. Our temperature maps display anisotropies on scales ranging from 2 to 8 degrees with rms amplitudes ranging from DT/T_{rms} =2.8 to 3.4 times 10^{-5} without reionization. If we allow the neutrinos to decay and subsequently reionize the IGM, the level of the temperature anisotropies is reduced to DT/T=0.8 to 1.1 times 10^{-5}. This includes the Doppler shifts suffered by the CMB photons as they scatter off the moving ionized fraction of the IGM. If, however, reionization is caused by a late and sudden energy input into the IGM, we find that some reduction in $\Delta T/T$ occurs, but not enough to make this scenario consistent with observations. We find that the decaying neutrino hypothesis is more efficient at reducing primordial anisotropies and generates smaller Doppler anisotropies than a model reionized via a late and sudden energy input into the IGM.


title: The Fate of the Magellanic Stream

abstract: We show that HI clouds in the Magellanic Clouds are stripped by the ram-pressure due to the halo and disk gases of the Galaxy. Molecular clouds are swept to the edge of the LMC, showing an eccentric distribution. The stripped HI clouds form a narrow band on the sky, and mimics the Magellanic Stream, when the LMC takes a polar orbit. We point out that the Magellanic Stream will fall into the Galaxy, and will be finally accreted by the Galactic disk. The accretion may cause warping of the inner gas disk of the Galaxy. Some clouds are accreted by the nuclear disk, which could explain the peculiar distribution and kinematics in the central region. The stripping of the interstellar gas from the Magellanic Clouds will rapidly change their morphological type into dwarf ellipticals.


title: The Clustering Properties of IRAS Galaxies

abstract: We investigate the clustering of galaxies in the QDOT redshift survey of IRAS galaxies. We find that the redshift space two-point correlation function is well approximated by a power-law of slope -1.11+\-0.09, with clustering length 3.87+\-0.32 Mpc/h out to pair separations of 25 Mpc\h. On scales larger than 40 Mpc/h, the correlation function is consistent with zero. The r.m.s. fluctuation in the count of QDOT galaxies above the Poisson level in spheres of radius 8 Mpc/h is \sigma_8^{IRAS}=0.58+\-0.14, showing that fluctuations in the distribution of IRAS galaxies on these scales are smaller than those of optical galaxies by a factor of about 0.65, with an uncertainty of 25%. We find no detectable difference between the correlation functions measured in redshift space and in real space, leading to a $2\sigma$ limit of b_{IRAS}/\Omega^{0.6} >1.05, where b_{IRAS} is the bias factor for IRAS galaxies and \Omega is the cosmological density parameter. The QDOT autocorrelation function calculated in concentric shells increases significantly with shell radius. This difference is more likely due to sampling fluctuations than to an increase of the clustering strength with galaxy luminosity, but the two effects are difficult to disentangle; our data allow at most an increase of \sim 20% in clustering strength for each decade in luminosity. For separations greater than \sim 3 Mpc/h, the cross-correlation function of Abell clusters (with Richness R\ge 1) and QDOT galaxies is well approximated by a power-law of slope -1.81+\-0.10, with clustering length 10.10+\-0.45 Mpc/h, and no significant signal beyond \gsim 50 Mpc\h. This cross-correlation depends only weakly on cluster richness.


title: $^7$Li Abundances in Halo Stars: Testing Stellar Evolution Models and the Primordial $^7$Li Abundance

abstract: A large number of stellar evolution models with [Fe/H] = -2.3 and -3.3 have been calculated in order to determine the primordial $^7$Li abundance and to test current stellar evolution models by a comparison to the extensive database of Li abundances in extremely metal poor halo stars observed by Thorburn (1994). Standard models do a good job of fitting the observed Li abundances in stars hotter than 5600 K. They predict a primordial $^7$Li abundance of Log N(Li) = 2.24\pm 0.03$. Models which include microscopic diffusion predict a downward curvature in the $^7$Li destruction isochrones at hot temperatures which is not present in the observations. Thus, the observations clearly rule out models which include uninhibited microscopic diffusion of $^7$Li from the surface of the star. The [Fe/H] = -2.3 stellar models which include both diffusion and rotational mixing provide an excellent match to the observations. Both the plateau stars and the heavily depleted cool stars are well fit by these models. The rotational mixing leads to considerable $^7$Li depletion in these models and the primordial $^7$Li abundance inferred from these models is Log N(Li) = $3.08\pm 0.1$.


title: The Evolution of Voids in the Adhesion Approximation

abstract: We apply the adhesion approximation to study the formation and evolution of voids in the Universe. Our simulations -- carried out using 128$^3$ particles in a cubical box with side 128 Mpc -- indicate that the void spectrum evolves with time and that the mean void size in the standard COBE-normalised Cold Dark Matter (hereafter CDM) model with $h_{50} = 1,$ scales approximately as $\bar D(z) = {\bar D_0\over \sqrt {1+z}},$ where $\bar D_0 \simeq 10.5$ Mpc. Interestingly, we find a strong correlation between the sizes of voids and the value of the primordial gravitational potential at void centers. This observation could in principle, pave the way towards reconstructing the form of the primordial potential from a knowledge of the observed void spectrum. Studying the void spectrum at different cosmological epochs, for spectra with a built in $k$-space cutoff we find that, the number of voids in a representative volume evolves with time. The mean number of voids first increases until a maximum value is reached (indicating that the formation of cellular structure is complete), and then begins to decrease as clumps and filaments merge leading to hierarchical clustering and the subsequent elimination of small voids. The cosmological epoch characterizing the completion of cellular structure occurs when the length scale going nonlinear approaches the mean distance between peaks of the gravitational potential. A central result of this paper is that


title: Null Result in Gamma-Ray Burst Lensed Echo Search

abstract: We have searched for gravitational-lens induced echoes between gamma-ray bursts (GRBs) in NASA's orbiting {\it Compton} Gamma Ray Observatory's Burst and Transient Source Experiment (BATSE) data. The search was conducted in two phases. In the first phase we compared all GRBs in a brightness complete sample of the first 260 GRBs with recorded angular positions having at least 5 \% chance of being coincident from their combined positional error. In the second phase, we compared all GRB light curves of the first 611 GRBs with recorded angular positions having at least 55 \% chance of being coincident from their combined positional error. No unambiguous gravitational lens candidate pairs were found in either phase, although a ``library of close calls" was accumulated for future reference. This result neither excludes nor significantly constrains a cosmological origin for GRBs.


title: Modelling of Shallow and Inefficient Convection in the Outer Layers of the Sun Using Realistic Physics

abstract: In an attempt to understand the properties of convective energy transport in the solar convection zone, a numerical model has been constructed for turbulent flows in a compressible, radiation-coupled, non-magnetic, gravitationally stratified medium using a realistic equation of state and realistic opacities. The time-dependent, three-dimensional hydrodynamic equations are solved with minimal simplifications. The statistical information obtained from the present simulation provides an improved understanding of solar photospheric convection. The characteristics of solar convection in shallow regions is parameterized and compared with the results of Chan and Sofia's simulations of deep and efficient convection (Chan and Sofia 1989). We assess the importance of the zones of partial ionization in the simulation, and confirm that the radiative energy transfer is negligible throughout the region except in the uppermost scale heights of the convection zone, a region of very high super-adiabaticity. When the effects of partial ionization are included, the dynamics of flows are altered significantly. However, we confirm the Chan and Sofia result that kinetic energy flux is non-negligible and can have a negative value in the convection zone.


title: Chaos and the Dynamical Evolution of Tidal Capture Binaries in Globular Clusters

abstract: Globular clusters cores harbour many low-mass X-ray binaries, cataclysmic variables, and millisecond pulsar binaries which are likely to have been formed via the process of tidal capture. Tidal capture binaries were originally thought to be responsible for halting core collapse and for subsequent re-expansion via the process of ``binary heating''. The standard model, now suggests that these binaries are no longer viable as a direct energy source. We present the results of a study which suggests that tidal capture binaries may indeed be a significant direct source of energy for the cores of globular clusters. We show that following capture, these binaries go through a short, violent chaotic phase, with the orbital eccentricity (and hence the tidal energy) suffering large changes on a very short timescale. This is followed by a long quiescent phase in which the tidal energy remains small while the orbit circularizes only via normal dissipative effects. This quiescent phase allows captured main sequence stars to evolve in the absence of large tides, a fact important for the production of low-mass X-ray binaries and cataclysmic variables.


title: The Hamburg Quasar Monitoring Program (HQM) at Calar Alto. I. Low amplitude variability in quasars

abstract: HQM is an optical broad-band photometric monitoring program carried out since Sept.~1988. Our main intention is to search for indications of microlensing in a sample of $\sim$\,100 selected quasars; however, we also want to study the intrinsic variability. We use a CCD camera equipped to the MPIA 1.2$\,$m telescope. Fully automatic photometric reduction relative to stars in the frame is done within a few minutes after each exposure, thus interesting brightness changes can be followed in detail. The typical photometric error is 1--2\,\% for a 17.5\,mag quasar, making HQM the most accurate long-term quasar monitoring program yet carried out. The main results of HQM which we discuss here are: (1) Concerning variability, quasars form two, clearly distinct classes, optically violent variables (OVVs) and non-OVVs. (2) All OVVs are radio loud and probably belong to the blazar class. (3) Non-OVVs have lightcurve gradients of at most several 0.1\,mag\,yr$^{-1}$ in the quasar restframe and can be well fitted by polynomials of low order. (4) Although our data cover only a relatively short timespan, we conclude that there is a large fraction of quasars which would be undetectable in photographic surveys using optical variability alone as the selection criterium. (5) A broad class of flat-spectrum radio quasars are no blazars, they are even less variable in the optical than radio quiet objects. (6) There is some statistical evidence for microlensing in our sample; if present it does, however, lead only to low lightcurve gradients.


title: The Hamburg Quasar Monitoring Program (HQM) at Calar Alto. II. Lightcurves of weakly variable objects

abstract: HQM is an optical broad-band photometric monitoring program carried out since Sept.~1988. We use a CCD camera equipped to the MPIA 1.2$\,$m telescope. Fully automatic photometric reduction relative to stars in the frames is done within a few minutes after each exposure, thus interesting brightness changes can be followed in detail. The typical photometric error is 1--2\,\% for a 17.5\,mag quasar. We here present lightcurves already evaluated but not shown in Paper I. We also discuss existing literature data.


title: The Hamburg Quasar Monitoring Program (HQM) at Calar Alto III. Lightcurves of optically violent variable sources

abstract: HQM is an optical broad-band photometric monitoring program carried out since September 1988. We use a CCD camera at the MPIA 1.2$\,$m telescope. Fully automatic photometric reduction relative to stars in the frames is done within a few minutes after each exposure, thus interesting brightness changes can be followed in detail. The typical photometric error is 1--2\,\% for a 17.5\,mag quasar. We here present lightcurves of 14 known violently variable sources and compare them with literature data. For two BL\,Lac objects, 1E\,1229+645 and 4C\,56.27 (1823+568), this paper is the first variability study. We have also carried out POSS photometry to obtain indications for variability on a longer timescale.


title: The Orion OB1 Association I. Stellar content

abstract: Walraven photometry of established and probable members of the Orion OB1 association is presented. Effective temperature, surface gravity, luminosity and mass are derived for all stars, using atmosphere models by Kurucz (1979). Absolute magnitudes are calculated using the Straizys and Kuriliene (1981) tables. Distance moduli and visual extinctions are determined. A comparison of the visual extinctions to IRAS $100\um$ data shows that the near edge of the Orion A and B clouds lies at a distance of $\sim 320\pc$, while the far edge is at $\sim 500\pc$. A method for deriving the ages of the subgroups by comparing theoretical isochrones to the observations in the log g, log T plane is presented. The derived ages suggest, contrary to earlier studies, that subgroup 1b is younger than 1c, which can possibly be explained by past geometries of the system of stars and gas. The initial mass function for Orion OB1 is derived with the aid of the Kolmogorov-Smirnoff test. Through extensive simulations, we show that it is very difficult to derive accurately the IMF from the available data. To within somewhat weak limits the IMF is found to be of the form $\xi (\log M)=AM^{-1.7\pm 0.2}$ for all subgroups. The energy output of the subgroups in the form of stellar winds and supernovae is calculated and compared to the observed size and expansion velocity of the Orion-Eridanus bubble. It is shown that the energy output of the association can account for the morphology and kinematics of the ISM.


title: Advection-Dominated Accretion: A Self-Similar Solution

abstract: We consider viscous rotating accretion flows in which most of the viscously dissipated energy is stored as entropy rather than being radiated. Such advection-dominated flows may occur when the optical depth is either very small or very large. We obtain a family of self-similar solutions where the temperature of the accreting gas is nearly virial and the flow is quasi-spherical. The gas rotates at much less than the Keplerian angular velocity; therefore, the central stars in such flows will cease to spin up long before they reach the break-up limit. Further, the Bernoulli parameter is positive, implying that advection-dominated flows are susceptible to producing outflows. Convection is likely in many of these flows and, if present, will tend to enhance the above effects. We suggest that advection-dominated accretion may provide an explanation for the slow spin rates of accreting stars and the widespread occurrence of outflows and jets in accreting systems.


title: The Zel'dovich Approximation and the Relativistic Hamilton-Jacobi Equation

abstract: Beginning with a relativistic action principle for the irrotational flow of collisionless matter, we compute higher order corrections to the Zel'dovich approximation by deriving a nonlinear Hamilton-Jacobi equation for the velocity potential. It is shown that the velocity of the field may always be derived from a potential which however may be a multi-valued function of the space-time coordinates. In the Newtonian limit, the results are nonlocal because one must solve the Newton-Poisson equation. By considering the Hamilton-Jacobi equation for general relativity, we set up gauge-invariant equations which respect causality. A spatial gradient expansion leads to simple and useful results which are local --- they require only derivatives of the initial gravitational potential.


title: Relaxing the Big-bang Bound to the Baryon Density

abstract: In the standard picture of big-bang nucleosynthesis the yields of D, $^3$He, $^4$He, and $^7$Li only agree with their inferred primordial abundances if the fraction of critical density contributed by baryons is between $0.01h^{-2}$ and $0.02h^{-2}$ ($h$ is the present value of the Hubble constant in units of $100\kms\Mpc^{-1}$). This is the basis of the very convincing and important argument that baryons can contribute at most 10\% of critical density and thus cannot close the Universe. Nonstandard scenarios involving decaying particles,$^1$ inhomogeneities in the baryon density,$^2$ and even more exotic ideas$^3$ put forth to evade this bound have been largely unsuccessful.$^4$ We suggest a new way of relaxing the bound: If the tau neutrino has a mass of $20\MeV-30\MeV$ and lifetime of $200\sec -1000\sec$, and its decay products include electron neutrinos, the bound to the baryon mass density can be loosened by a about factor of $10$. The key is the decay-generated electron antineutrinos: around the time of nucleosynthesis they are captured by protons to produce neutrons, thereby changing the outcome of nucleosynthesis. Experiments at $e^\pm$ colliders should soon be sensitive to a tau-neutrino mass in the required range.


title: Diffuse X-Ray Emission Around High-Redshift, Radio-Loud Qsos

abstract: We announce the detection of correlations on angular scales of $\ga10'$ between optically bright, high-redshift, radio-loud QSOs with diffuse X-ray emission seen by ROSAT in the {\it All-Sky Survey}. These correlations reach significance levels of up to $99.8\%$. A comparison of the results with a sample of control fields, bootstrapping analyses, and Kolmogorov-Smirnov tests provide unambiguous support for the statistical significance of the correlations found. We argue that the detected enhanced diffuse X-ray emission is in the foreground of the QSOs, and that it is probably due to galaxy clusters which magnify the QSOs by their gravitational lensing effect, thereby giving rise to a magnification bias in the background source sample. A comparison of the results presented below with correlations previously found between the same QSO sample and either Lick or IRAS galaxies provides further evidence for this interpretation, and identifies positions in the sky where weak gravitational lensing may be detected by searching for coherent distortions of background galaxy images.


title: Light Propagation in Arbitrary Spacetimes and the Gravitational Lens Approximation

abstract: We reformulate the transport equation which determines the size, shape and orientation of infinitesimal light beams in arbitrary spacetimes. The behaviour of such light beams near vertices and conjugate points is investigated, with special attention to the singular behaviour of the optical scalars. We then specialize the general transport equation to the case of an approximate metric of an inhomogeneous universe, which is a Friedmann metric `on average' with superposed isolated weak matter inhomogeneities. In a series of well-defined approximations, the equations of gravitational lens theory are derived. Finally, we derive a relative optical focusing equation which describes the focusing of light beams relative to the case that the beam is unaffected by matter inhomogeneities in the universe, from which it follows immediately that no beam can be focused less than one which is unaffected by matter clumps, before it propagates through its first conjugate point.


title: Long-term planetary integration with individual time steps

abstract: We describe an algorithm for long-term planetary orbit integrations, including the dominant post-Newtonian effects, that employs individual timesteps for each planet. The algorithm is symplectic and exhibits short-term errors that are $O(\epsilon\Omega^2\tau^2)$ where $\tau$ is the timestep, $\Omega$ is a typical orbital frequency, and $\epsilon\ll1$ is a typical planetary mass in solar units. By a special starting procedure long-term errors over an integration interval $T$ can be reduced to $O(\epsilon^2\Omega^3\tau^2T)$. A sample 0.8 Myr integration of the nine planets illustrates that Pluto can have a timestep more than 100 times Mercury's, without dominating the positional error. Our algorithm is applicable to other $N$-body systems.


title: The Structure and Evolution of Lyman-alpha Forest Clouds in the Minihalo Model

abstract: Results are presented for one-dimensional numerical hydrodynamics computations of the structure and evolution of Lya forest clouds gravitationally confined by dark matter minihalos. The clouds are developed from linear perturbations at high redshift and exposed to either a QSO- or galaxy-dominated metagalactic radiation field at moderate redshifts. While the emphasis is on spherical systems, slab symmetry is also considered. Three zones may be identified in a collapsed cloud: (1) a quasi-hydrostatic core in thermal equilibrium, (2) a nonhydrostatic intermediate zone out of thermal equilibrium, and (3) a cosmological accretion layer joining onto the Hubble expansion. Most of the measured Lya forest column densities arise in the intermediate zone. The development of the core would result in a flattening in the column density distribution near an HI column density of 10**15-10**16 cm**-2. The cloud diameters corresponding to an HI column density of 10**14 cm**-2 lie in the range 10-60 kpc, while systems with column densities exceeding 10**15 cm**-2 have diameters smaller than 10 kpc. Systems with circular velocities exceeding 50 km/s result in clouds which contract until they become Jeans unstable and collapse. The critical column density for collapse is 10**17-10**18 cm**-2. A mild correlation of Doppler parameter with neutral hydrogen column density is found in several models for systems with column densities less than 10**13.5 cm**-2, with Doppler parameters occurring as low as b > 20 km/s for 10**13 cm**-2 column density systems. No lines with b < 15 km/s are found.


title: A Search for Extragalactic Methanol Masers

abstract: A sensitive search for 6.7--GHz methanol maser emission has been made towards 10 galaxies that have yielded detectable microwave molecular--line transitions. These include several which show OH megamaser or superluminous \water\/ maser emission. Within the Galaxy, \methanol\/ and OH masers often occur in the same star formation regions and, in most cases, the \methanol\/ masers have a greater peak flux density than their OH counterparts. Thus we might expect \methanol\/ masers to be associated with extragalactic OH maser sources. We failed to detect any emission or absorption above our 60--mJy detection limit. We conclude that if the physical conditions exist to produce \methanol\/ megamaser emission, they are incompatible with the conditions which produce OH megamaser emission.


title: On the Mass of the Dark Compact Objects in the Galactic Disk

abstract: Recently the Polish-American collaboration OGLE has reported the observation of four possible microlensing events by monitoring, over several months, the brightness of millions of stars in the region of the galactic bulge. If these events are due to microlensing, the most accurate way to get information on the mass of the dark compact objects, that acted as gravitational lenses, is to use the method of the mass moments. Here I apply this method to the analysis of the events detected by OGLE. The average mass turns out to be 0.28$M_{\odot}$, suggesting that the lens objects are faint disk stars. The same method applied to the five microlensing events detected so far by the EROS and MACHO collaborations, which monitor stars in the Large Magellanic Cloud, leads to an average value of 0.08$M_{\odot}$ for the dark compact halo objects.


title: Gravitino Production in the Inflationary Universe and the Effects on Big-Bang Nucleosynthesis

abstract: Gravitino production and decay in the inflationary universe are reexanimed. Assuming that gravitino mainly decays into photon and photino, we have calculated the upperbound of the reheating temperature. Compared to previous works, we have essentially improved the following two points; (i) the helicity $\pm\frac{3}{2}$ gravitino production cross sections are calculated by using the full relevant terms in the supergravity lagrangian, and (ii) high energy photon spectrum is obtained by solving the Boltzmann equations numerically. Photo-dissociation of the light elements (D, T, $^3$He, $^4$He) leads to the most stringent upperbound of the reheating temperature, which is given by ($10^{6}$--$10^{9}$)GeV for the gravitino mass 100GeV--1TeV. On the other hand, requiring that the present mass density of photino should be smaller than the critical density, we find that the reheating temperature have to be smaller than ($10^{11}$--$10^{12}$)GeV for the photino mass (10--100)GeV, irrespectively of the gravitino mass. The effect of other decay channel is also considered.


title: VLA Measurements of a sample of Planetary Nebulae

abstract: We report on new radio measurements of Galactic planetary nebulae, aimed at resolving the controversies on the reliability of older VLA flux densities and the suggested deviations from the standard Galactic extinction law found for planetary nebulae. We show that for faint ($<$10 mJy) objects observed at high angular resolution, previous determinations are indeed too low. For the bright objects we find no significant differences. The new values are the most accurate flux determinations yet for planetary nebulae, reaching 1\% for the brightest objects in the sample. Based on the new data, we confirm that there is a systematic difference between the extinction derived from the radio/H$\beta$ flux ratio and derived from the Balmer decrement, which led to the suggestion of deviations from the standard extinction law. However, final confirmation of this has to await the availibility of more, accurate measurements of the Balmer (and/or Paschen) lines.


title: On the Luminosity Function of Early--Type Galaxies

abstract: In a recent paper Loveday et al. (1992) have presented new results on the luminosity function for a sample of galaxies with $b_J \le 17.15$. After having morphologically classified each galaxy (early--type, late--type, merged or uncertain), they have estimated the parameters of a Schechter luminosity function for early-- and late--type galaxies. However, in their sample there is a bias against identifying early--type galaxies at large distances and/or faint magnitudes: in fact, many of the early--type galaxies at faint magnitudes have probably been classified as ``uncertain". As discussed in Loveday et al., the existence of such a bias is indicated by the fact that for these galaxies $<V/V_{max}>=0.32$. In this paper we show, both theoretically and through the use of simulated samples, that this incompleteness strongly biases the derived parameters of the luminosity function for early--type galaxies. If no correction for such incompleteness is applied to the data (as done by Loveday et al.), one obtains a flatter slope $\alpha$ and a brighter $M^*$ with respect to the real parameters.


title: Possible Systematic Decreases in the Age of Globular Clusters

abstract: The ages of globular clusters inferred from observations depends sensitively on assumptions like the initial helium abundance and the mass loss rate. A high helium abundance (e.g., $Y\approx$0.28) or a mass loss rate of $\sim$10$^{-11}M_\odot$ yr$^{-1}$ near the main sequence turn-off region lowers the current age estimate from 14 Gyr to about 10--12 Gyr, significantly relaxing the constraints on the Hubble constant, allowing values as high as 60km/sec/Mpc for a universe with the critical density and 90km/sec/Mpc for a baryon-only universe. Possible mechanisms for the helium enhancement in globular clusters are discussed, as are arguments for an instability strip induced mass loss near the turn-off. Ages lower than 10 Gyr are not possible even with the operation of both of these mechanisms unless the initial helium abundance in globular clusters is $>0.30$, which would conflict with indirect measurements of helium abundances in globular clusters.


title: Testing the Velocity Field in Non--Scale Invariant Cold Dark Matter Models

abstract: We analyze the cosmic peculiar velocity field as traced by a sample of 1184 spiral, elliptical and S0 galaxies, grouped in 704 objects. We perform a statistical analysis, by calculating the bulk flow, Cosmic Mach Number and velocity correlation function for this sample and for mock catalogs extracted from a set of N--body simulations. We run four cold dark matter (CDM) simulations: two tilted models (with spectral index $n=0.6$ and $n=0.8$), the standard model ($n=1$) and a ``blue" one ($n=1.2$), with different values of the linear bias parameter $b$. By means of a Maximum Likelihood analysis we estimate the ability of our models to fit the observations, as measured by the above statistics, and to reproduce the Local Group properties. On the basis of this analysis we conclude that the best model is the unbiased standard model $(n,~b)=(1,~1)$, even though the overall flatness of the joint likelihood function implies that one cannot strongly discriminate models in the range $0.8 \le n \le 1$, and $1 \le b \le 1.5$. Models with $b \geq 2.5$ are rejected at the $95\%$ confidence level. For $n=0.8$ the values of $b$ preferred by the present analysis, together with the {\em COBE} data, require a negligible contribution to $\Delta T/T$ by gravitational waves. Finally, the blue model, normalized to {\em COBE}, does not provide a good fit to the velocity data.


title: On Determining the Spectrum of Primordial Inhomogeneity from the COBE DMR Sky Maps: I. Method

abstract: The natural approach to a spectral analysis of data distributed on the sky employs spherical harmonic decomposition. A common problem encountered in practical astronomy is the lack of full sky coverage in the available data. For example, the removal of the galactic plane data from the COBE DMR sky maps compromises Fourier analysis of the cosmic microwave background (CMB) temperature distribution due to the loss of orthogonality of the spherical harmonics. An explicit method for constructing orthonormal functions on an incomplete (e.g. Galaxy-cut) sphere is presented. These functions should be used in the proper Fourier analysis of the COBE DMR sky maps to provide the correct input for the determination of the spectrum of primordial inhomogeneity. The results of such an analysis are presented in an accompanying Letter. A similar algebraic construction of appropriate functions can be devised for other astronomical applications.


title: On Determining the Spectrum of Primordial Inhomogeneity from the COBE DMR Sky Maps: II. Results of Two Year Data Analysis

abstract: A new technique of Fourier analysis on a cut sky (Gorski, 1994) has been applied to the two year COBE DMR sky maps. The Bayesian power spectrum estimation results are consistent with the Harrison-Zel'dovich n=1 model. The maximum likelihood estimates of the parameters of the power spectrum of primordial perturbations are n=1.22 (1.02) and Q_{rms-PS}=17 (20) uK including (excluding) the quadrupole anisotropy. The marginal likelihood function on n renders n=1.10 \pm 0.32 (0.87 \pm 0.36).


title: An Optical Interferometer with Wavelength Dispersion

abstract: MAPPIT is an optical interferometer installed at the coude focus of the 3.9-m Anglo-Australian Telescope. The instrument combines non-redundant masking with wavelength dispersion and is able to record fringes simultaneously over a wide bandwidth. For typical observations centred near 600 nm, the bandwidth is 55 nm and the spectral resolution is 3 nm. This paper describes the instrument and the data processing methods and presents some results. We find the star sigma Sgr to be a close binary; the system is only partially resolved, with a separation of (11.5 +/- 2) milliarcsec (assuming the components to have equal magnitudes). We also give angular diameter measurements of two red giant stars, alpha Sco and beta Gru. The observations of beta Gru (spectral type M5 III) resolve the star for the first time and give an equivalent uniform-disk diameter of (27 +/- 3) milliarcsec.


title: Has A Standard Model Solution to the Solar Neutrino Problem Been Found?

abstract: The claim by Dar and Shaviv that they have found a standard model solution to the solar neutrino problem is based upon an incorrect assumption made in extrapolating nuclear cross sections and the selective use of a small fraction of the nuclear physics and of the neutrino data. In addition, five different solar model codes show that the rate obtained for the chlorine experiment using the Dar-Shaviv stated parameters differs by at least $14 \sigma$ from the observed rate.


title: The Settling of Warped Disks in Oblate Dark Halos

abstract: When a galaxy forms, the disk may initially be tilted with respect to a flattened dark halo. The misalignment between the disk and the halo is a common explanation for galactic disk warps, since in this state, disks have precessing bending modes which resemble real warps. A dubious assumption in this theory is that the gravitational response of the halo is negligible. We therefore calculate the response of an oblate halo to a precessing inclined exponential disk using a variety of techniques. We construct models with a rigid exponential disk precessing in a particle halo, a particle disk precessing inside a static bulge/halo potential, and a self-consistent model with a particle disk, bulge and halo. When the disk:halo mass ratio is small ($\sim$10\%) within 5 exponential scale radii, the disk settles to the equatorial plane of the halo within 5 orbital times. When the disk and halo mass are comparable, the halo rapidly aligns with the disk within a few orbital times, while the disk inclination drops. The rapid response of the halo to a inclined disk suggests that the warps seen in galactic disks are not likely the bending modes of a precessing disk inside of triaxial halos. If a galaxy forms inclined to the principal plane of a dark halo, either the disk will settle to a principal plane or the inner halo will twist to align with the disk. The outer halo will remain misaligned for a much longer time and therefore may still exert a torque. Warped bending modes may still exist if the misalignment of the outer halo persists for a Hubble time.


title: Contribution to the local cosmic-ray flux from the Geminga supernova

abstract: The contribution to the local cosmic-ray flux from the Geminga supernova is calculated assuming shock acceleration to 10^14 eV in a remnant which was formed several 10^5 years ago along with the Geminga pulsar. The particles are propagated to Earth using a simple diffusion model. In the region below the knee in the spectrum, it is found the supernova may contribute 10% of the local cosmic-ray flux, assuming plausible explosion parameters. The contribution to the amplitude of the anisotropy is not in conflict with the data in this energy region.


title: Velocity Bias in Clusters

abstract: Hierarchical gravitational clustering creates galaxies that usually do not fully share the dynamical history of an average particle in the mass field. In particular, galaxy tracers identified in numerical simulations can have individual velocity dispersions in virialized regions a factor \bvs\ lower than the dark matter. The field average of the pairwise velocity dispersion depends on the statistical weighting of collapsed regions, so that the tracer pairwise dispersion is a different factor, \bvp, times the density field value. A model of a cool equilibrium tracer population demonstrates that mass to light segregation is very sensitive to single particle velocity bias. For the $\phi=-GM/(r+a)$ potential a $\bvs=0.9$ tracer indicates a virial mass about a factor of 5 low. The likely value of \bvs\ is estimated and the simple equilibrium model is tested for applicability using a $10^6$ particle simulation of the formation of a single cluster from cosmological initial conditions. From this simulation, the value of \bvs\ is estimated as $0.8\pm0.1$. The pairwise velocity dispersion bias, \bvp, which is equal to \bvs\ augmented with any anti-bias of galaxies against high velocity dispersion clusters. The value of \bvp\ is estimated to be 15\% less than \bvs, consequently $\bvp=0.6^{+0.2}_{-0.1}$. Comparison of cluster mass and luminosity profiles at large radii is a test for the existence of single particle velocity bias. [compressed postscript file cluster.ps.Z via anonymous ftp from orca.astro.washington.edu in ~ftp/pub/hpcc, or carlberg@astro.washington.edu]


title: On Solar Model Solutions to the Solar Neutrino Problem

abstract: Without assuming any solar models and neutrino flavor conversions, $^8$B neutrinos seen by the Kamiokande experiment should contribute 2.6$\pm 0.45$ SNU to the chlorine experiment. When this rate is compared with the total event rate of 2.3$\pm 0.2$ SNU observed by the Homestake experiment which should inlude a 0.2 SNU contribution from uncertainty-free pep neutrinos, there may still be a possible evidence that $^7$Be neutrinos are more severely suppressed than the $^8$B neutrinos with respect to the predictions of standard solar models, which cannot be explained by any known astrophysics solution. Given a Kamiokande event rate of larger than 36$\%$ (2$\sigma$) of the prediciton of Bahcall and Pinsonneault's standard solar model, variations of standard solar models yield minimal rates of 3.6 SNU for the Homestake experiment and 114 SNU for GALLEX and SAGE. Therefore, variations of standard solar models as solutions to the solar neutrino problem are so far inconsistent with the Homestake experiment and only marginally allowed by the gallium experiments. If the gallium experiments confirm a flux significantly below 114 SNU, it would imply unconventional physics.


title: Wiener Filtering of the COBE DMR Data

abstract: We derive an optimal linear filter to suppress the noise from the COBE DMR sky maps for a given power spectrum. We then apply the filter to the first-year DMR data, after removing pixels within $20^\circ$ of the Galactic plane from the data. The filtered data have uncertainties 12 times smaller than the noise level of the raw data. We use the formalism of constrained realizations of Gaussian random fields to assess the uncertainty in the filtered sky maps. In addition to improving the signal-to-noise ratio of the map as a whole, these techniques allow us to recover some information about the CMB anisotropy in the missing Galactic plane region. From these maps we are able to determine which hot and cold spots in the data are statistically significant, and which may have been produced by noise. In addition, the filtered maps can be used for comparison with other experiments on similar angular scales.


title: Geminga: A Cooling Superfluid Neutron Star

abstract: We show that Geminga's surface temperature can be understood by both types of neutrino cooling scenarios, i.e, slow neutrino cooling by the modified Urca process or fast neutrino cooling by the direct Urca process or by some exotic matter, and thus does not allow us to discriminate between these two competing schemes. However, for both types of scenarios, agreement with the observed temperature can only be obtained if BARYON PAIRING IS PRESENT IN MOST, IF NOT ALL, OF THE CORE OF THE STAR. We also comment on the recent temperature estimates for PSR 0656+14 and PSR 1055-52, which pertain to the same photon cooling era. We argue though that observational evidence for the slow neutrino cooling model (the ``standard'' model) is in fact very dim and that the interpretation of the surface temperature of all neutron stars could be done with a reasonable theoretical a priori within the fast neutrino cooling scenarios only. In this case, Geminga, PSR 0656+14, and PSR 1055-52 all show evidence of baryon pairing down to their very centers. xxx----To be published in Ap.J. (June 10 isue).---xxx


title: The Amplitude and Spectral Index of the Large Angular Scale Anisotropy in the Cosmic Microwave Background Radiation

abstract: In many cosmological models, the large angular scale anisotropy in the cosmic microwave background is parameterized by a spectral index, $n$, and a quadrupolar amplitude, $Q$. For a Peebles-Harrison-Zel'dovich spectrum, $n=1$. Using data from the Far Infra-Red Survey (FIRS) and a new statistical measure, a contour plot of the likelihood for cosmological models for which $-1 < n < 3$ and $0 \le Q \le 50~\mu{\rm K}$ is obtained. We find that the likelihood is maximum at $(n, Q) = (1.0, 19 uK)$. For constant $n$ the likelihood falls to half its maximum at $Q \approx 14 uK$ and $25 uK$ and for constant $Q$ the likelihood falls to half its maximum at $n \approx 0.5$ and $1.4$. Regardless of $Q$, the likelihood is always less than half its maximum for $n < -0.4$ and for $n > 2.2$, as it is for $Q < 8 uK$ and $Q > 44$.


title: The Correlation Function of Flux-Limited X-Ray Clusters

abstract: We show that the spatial correlation function of a flux-limited sample of X-ray selected clusters of galaxies will exhibit a correlation scale that is smaller than the correlation scale of a volume-limited, richness-limited sample of comparable apparent spatial density. The flux-limited sample contains clusters of different richnesses at different distances: poor groups are found nearby and rich clusters at greater distances. Since the cluster correlation strength is known to increase with richness, the flux-limited sample averages over the correlations of poor and rich clusters. On the other hand, a volume-limited, richness-limited sample has a minimum richness threshold, and a constant mixture of richnesses with redshift. Using the observed correlation scale for rich ($R\ge 1$) clusters, $r_o (R\ge 1) = 21\pm 2 h^{-1}$Mpc that was determined from previous volume-limited studies, we derive for the ROSAT flux-limited X-ray cluster sample $r_o({\rm flux}{\rm -}{\rm limited}) \approx 14h^{-1}$Mpc, in agreement with the recently observed value of $13.7\pm 2.3 h^{-1}$Mpc.


title: A Hydrodynamic Approach to Cosmology: The Mixed Dark Matter Cosmological Scenario

abstract: We compute the evolution of spatially flat, mixed cold and hot dark matter (``MDM") models containing both baryonic matter and two kinds of dark matter. The mean final Zeldovich-Sunyaev $y$ parameter is estimated to be $\bar y=(5.4\pm 2.7)\times 10^{-7}$, with a rms fluctuation of approximately $\bar{\delta y}=(6.0\pm 3.0)\times 10^{-7}$ on arc minute scales. The rate of galaxy formation peaks at an even later epoch ($z\sim 0.3$) than in the standard ($\Omega=1$, $\sigma_8=0.67$) CDM model ($z\sim 0.5$) and, at a redshift of z=4 is nearly a factor of a hundred lower than for the CDM model with the same value of $\sigma_8$. The model yields a cluster mass function which is about a factor of 4 higher than observed but a cluster-cluster correlation length lower by a factor of 2 than what is observed but both are closer to observations than in the COBE normalized CDM model. The one dimensional pairwise velocity dispersion is $605\pm 8$km/s at $1h^{-1}$ separation, lower than that of the CDM model normalized to COBE, but still significantly higher than observations (Davis \& Peebles 1983). A plausible velocity bias $b_v=0.8\pm 0.1$ on this scale will reduce but not remove the discrepancy. Overall, the model is similar both on large and intermediate scales to the standard CDM model normalized to the same value of $\sigma_8$, but the problem with regard to late formation of galaxies is more severe in this model than in that CDM model. Adding hot dark matter significantly improves the ability of COBE normalized CDM scenario to fit existing observations, but the model is in fact not as good as the CDM model with the same $\sigma_8$ and is still probably unsatisfactory with regard to several critical tests.


title: X-Ray Clusters in a CDM$+\Lambda$ Universe: A Direct, Large-Scale, High Resolution, Hydrodynamic Simulation

abstract: A new, three-dimensional, shock capturing, TVD hydrodynamic code is utilized to determine the distribution of hot gas in a CDM$+\Lambda$ model universe. This open model succeeds in matching local observations of clusters in contrast to the standard $\Omega=1$, CDM model, which fails. It predicts an order of magnitude decline in the number density of bright ($h\nu = 2-10$keV) clusters from $z=0$ to $z=2$ in contrast to a slight increase in the number density for standard $\Omega=1$, CDM model. This COBE-normalized CDM$+\Lambda$ model produces approximately the same number of X-ray clusters having $L_x>10^{43}$erg/s as observed. The number density of bright clusters increases to $z\sim 0.2-0.5$ and declines, but the luminosity per typical cluster decreases monotonically with redshift, with the result that the number density of bright clusters shows a broad peak near $z=0.5$, and then a rapid decline as $z\rightarrow 3$. The most interesting point which we find is that the temperatures of clusters in this model freeze out at later times ($z\le 0.3$), while previously we found in the CDM model that there was a steep increase during the same interval of redshift. The trend should be detectable even with a relatively ``soft" X-ray instrument such as ROSAT, providing a powerful discriminant between $\Omega=1$ and $\Omega<1$ models. Examining the ratio of gas-to-total mass in the clusters, we find a slight antibias ($b=0.9$ or $({\Omega_{gas}\over\Omega_{tot}})_{cl} =0.083\pm 0.007$), which is consistent with observations [$({\Omega_{gas}\over\Omega_{tot}})_{obs}=0.097\pm 0.019$ for the Coma cluster for the given value of $h$, \cf, White 1991].


title: Hot Gas in the CDM Scenario: X-Ray Clusters from a High Resolution Numerical Simulation

abstract: We exmaine the distribution of hot gas in a standard CDM model of the universe using high resolution hydrodynamic simulations. Adopting standard parameters determined from COBE and light element nucleosynthesis, $\sigma_8=1.05$, $\Omega_b=0.06$ and assuming $h=0.5$, we find the X-ray emitting clusters, compute the luminosity function at several wavelengths, the temperature distribution and estimated sizes as well as the evolution of these quantities with redshift. This standard CDM model, normalized to COBE, produces approximately 5 times too much emission from clusters having $L_x>10^{43}$erg/s, a not unexpected result. If all other parameters were unchanged, we would expect adequate agreement for $\sigma_8=0.6$. This provides a new and independent argument for lower small scale power than standard CDM at the $8h^{-1}$Mpc scale. The background radiation field at 1keV due to clusters in this model is approximately $1/3$ of the observed background which, after correction for numerical effects, again indicates approximately 5 times too much emission and the appropriateness of $\sigma_8=0.6$. If we had used the observed ratio of gas to total mass in clusters, rather than basing the mean density on light element nucleosynthesis, then the computed luminosity of each cluster would have increased still further, by a factor of approximately ten. Examining the ratio of gas to total mass in the clusters (which we find to be anti-biased by a factor of approximately 0.6), normalized to $\Omega_b h^2=0.015$, and comparing to observations, we conclude, in agreement with S. White, that the cluster observations argue for an open universe.


title: A Comparison of Cosmological Hydrodynamic Codes

abstract: We present a detailed comparison of the simulation results of various cosmological hydrodynamic codes. Starting with identical initial conditions based on the Cold Dark Matter scenario for the growth of structure, we integrate from redshift $z=20$ to $z=0$ to determine the physical state within a representative volume of size $L^3$ where $L=64 h^{-1} {\rm Mpc}$. Five independent codes are compared: three of them Eulerian mesh based and two variants of the Smooth Particle Hydrodynamics "SPH" Lagrangian approach. The Eulerian codes were run at $N^3=(32^3,~64^3,~128^3,~{\rm and},~256^3)$ cells, the SPH codes at $N^3= 32^3$ and $64^3$ particles. Results were then rebinned to a $16^3$ grid with the expectation that the rebinned data should converge, by all techniques, to a common and correct result as $N \rightarrow \infty$. We find that global averages of various physical quantities do, as expected, tend to converge in the rebinned model, but that uncertainties in even primitive quantities such as $\langle T \rangle$, $\langle \rho^2\rangle^{1/2}$ persists at the 3\%-17\% level after completion of very large simulations. The two SPH codes and the two shock capturing Eulerian codes achieve comparable and satisfactory accuracy for comparable computer time in their treatment of the high density, high temperature regions as measured in the rebinned data; the variance among the five codes (at highest resolution) for the mean temperature (as weighted by $\rho^2$) is only 4.5\%. Overall the comparison allows us to better estimate errors, it points to ways of improving this current generation of hydrodynamic codes and of suiting their use to problems which exploit their individually best features.


title: Adiabatic Invariants in Stellar Dynamics: I. Basic concepts

abstract: The adiabatic criterion, widely used in astronomical dynamics, is based on the harmonic oscillator. It asserts that the change in action under a slowly varying perturbation is exponentially small. Recent mathematical results precisely define the conditions for invariance show that this model does not apply in general. In particular, a slowly varying perturbation may cause significant evolution stellar dynamical systems even if its time scale is longer than any internal orbital time scale. This additional `heating' may have serious implications for the evolution of star clusters and dwarf galaxies which are subject to long-term environmental forces. The mathematical developments leading to these results are reviewed, and the conditions for applicability to and further implications for stellar systems are discussed. Companion papers present a computational method for a general time-dependent disturbance and detailed example.


title: Adiabatic Invariants in Stellar Dynamics: II. Gravitational shocking

abstract: A new theory of gravitational shocking based on time-dependent perturbation theory shows that the changes in energy and angular momentum due to a slowly varying disturbance are not exponentially small for stellar dynamical systems in general. It predicts significant shock heating by slowly varying perturbations previously thought to be negligible according to the adiabatic criterion. The theory extends the scenarios traditionally computed only with the impulse approximation and is applicable to a wide class of disturbances. The approach is applied specifically to the problem of disk shocking of star clusters.


title: Adiabatic Invariants in Stellar Dynamics: III. Application to Globular Cluster Evolution

abstract: The previous two companion papers demonstrate that slowly varying perturbations do not result in adiabatic cutoffs and provide a formalism for computing the long-term effects of time-dependent perturbations on stellar systems. Here, the theory is implemented in a Fokker-Planck code and a suite of runs illustrating the effects of shock heating on globular cluster evolution are described. Shock heating alone results in considerable mass loss for clusters with $R_g\lta8\kpc$: a concentration $c=1.5$ cluster with $R_g=8\kpc$ loses up to $95\%$ of its initial mass in $15\Gyr$. Only those with concentration $c\lta1.3$ survive disk shocks inside of this radius. Other effects, such as mass loss by stellar evolution, will increase this survival bound. Loss of the initial halo together with mass segregation leads to mass spectral indices, $x$, which may be considerably larger than their initial values.


title: Testing Higher--Order Lagrangian Perturbation Theory Against Numerical Simulations -- 2. Hierarchical Models

abstract: We present results showing an improvement of the accuracy of perturbation theory as applied to cosmological structure formation for a useful range of scales. The Lagrangian theory of gravitational instability of Friedmann--Lema\^\i tre cosmogonies investigated and solved up to the third order in the series of papers by Buchert (1989, 1992, 1993), Buchert \& Ehlers (1993), Buchert (1994), Ehlers \& Buchert (1994), is compared with numerical simulations. In this paper we study the dynamics of hierarchical models as a second step. In the first step (Buchert, Melott and Wei{\ss} 1994) we analyzed the performance of the Lagrangian schemes for pancake models, i.e., models which initially have a truncated power spectrum. We here explore whether the results found for pancake models carry over to hierarchical models which are evolved deeply into the non--linear regime............We find that for spectra with negative power--index the second--order scheme performs considerably better than TZA in terms of statistics which probe the dynamics, and slightly better in terms of low--order statistics like the power--spectrum. In cases with much small--scale power the gain from the higher--order schemes is small, but still measurable. However, in contrast to the results found for pancake models, where the higher--order schemes get worse than TZA at late non--linear stages and on small scales, we here find that the second--order model is as robust as TZA, retaining the improvement at later stages and on smaller scales. In view of these


title: An Accelerated Lambda Iteration Method for Multilevel Radiative Transfer III. Noncoherent Electron Scattering

abstract: Since the mass of the electron is very small relative to atomic masses, Thomson scattering of low-energy photons ($h\nu \ll m_ec^2$) produces thermal Doppler frequency shifts that are much larger than atomic Doppler widths. A method is developed here to evaluate the electron scattering emissivity from a given radiation field which is considerably faster than previous methods based on straightforward evaluation of the scattering integral. This procedure is implemented in our multilevel radiative code (MALI), which now takes full account of the effects of noncoherent electron scattering on level populations, as well as on the emergent spectrum. Calculations using model atmospheres of hot, low-gravity stars display not only the expected broad wings of strong emission lines but also effects arising from the scattering of photons across continuum edges. In extreme cases this leads to significant shifts of the ionization equilibrium of helium.


title: The Motions of Clusters and Groups of Galaxies

abstract: The distributions of peculiar velocities of rich clusters and of groups of galaxies are investigated for different cosmological models and are compared with observations. Four cosmological models are studied: standard ($\Omega=1$) CDM, low-density CDM, HDM ($\Omega=1$), and PBI. We find that rich clusters of galaxies exhibit a Maxwellian distribution of peculiar velocities in all models, as expected from a Gaussian initial density fluctuation field. The cluster 3-D velocity distribution is generally similar in the models: it peaks at $v \sim 500$ km s$^{-1}$, and extends to high cluster velocities of $v \sim 1500$ km s$^{-1}$. Approximately 10\% of all model rich clusters move with high peculiar velocities of $v \ge 10^3$ km s$^{-1}$. The highest velocity clusters frequently originate in dense superclusters. The group velocity distribution is, in general, similar to the velocity distribution of the rich clusters. In all but the low-density CDM model, the mass exhibits a longer tail of high velocities than do the clusters. This high-velocity tail originates mostly from the high velocities that exist within rich clusters. The model velocity distributions of groups and clusters of galaxies are compared with observations. The data are generally consistent with the models, but exhibit a somewhat larger high-velocity tail, to $v_r \sim 3000$ km s$^{-1}$. While this high-velocity tail is similar to the HDM model predictions, the data are consistent with the other models studied, including the low-density CDM model, which best fits most other large-scale structure observations. The observed velocity


title: Microlensing and Halo Mass in Form of Stars

abstract: Paczy\'nski (1986) suggested that ``dark" objects in the halo of our Galaxy could enhance the luminosity of foreground stars, acting as gravitational microlenses. Such events has been recently reported by different collaborations. We assume that these microlensing events are produced by baryonic objects in the halo of our Galaxy. Rather than adopting a mean mass for the lensing objects we consider them distributed along a mass function similar to the disk one. Even if the number of microlensing events so far detected is not large, we show that it is possible to constrain the fraction of mass density in the halo of our Galaxy in form of stars. We estimate that this fraction varies between 0.1 and 0.8, if halo stars trace the total halo mass density; in the opposite case this range becomes narrower. If lensing objects have not been formed apart from the other halo stars, but in a more general context of star formation, neutron stars must also be present. We discuss these results in relation to galactic halo models for gamma-ray bursts.


title: Active Optics and Coronography with the Hubble Space Telescope

abstract: In the field of planet and proto-planetary disk detection, achieving high angular resolution and high dynamic range is a necessity. Coronography coupled with adaptive optics on Hubble Space Telescope is a way to get both good spatial resolution and high dynamic range. However, because of the residual figure errors on the primary and on the secondary, HST has a scattered light level that prevents it from detecting extra-solar planets. Our simulations show that by using an active mirror (400-1000 actuators) in the optical path of HST with adaptive optics, we can correct the mirror errors and decrease the scattering level by a factor $10^{2}$ (from $10^{-4}$ to $10^{-6}$ fainter than the star). Furthermore, by controlling the spatial frequencies of the active mirror with a {\em dark hole\/} algorithm we can decrease the scattering level in image zones where planet detection is likely. Using this technique, we have succeeded in decreasing the scattering level to $3\times 10^{-8}$ of the star intensity within 1 arcsec from the central star. This will allow the detection of a Jupiter-like planet $10^{-9}$ times dimmer than the central star located 10 pc away in 1 hour of integration time with a signal-to-noise ratio of 5. This paper describes the method used to determine the actuator strokes applied to a deformable mirror to achieve planet detection and the design of a coronograph which implements this novel technique.


title: Gauge-Invariant Cosmological Perturbations in Generalized Einstein Theories

abstract: Using the covariant approach and conformal transformations, we present a gauge-invariant formalism for cosmological perturbations in generalized Einstein theories (GETs), including the Brans-Dicke theory, theories with a non-minimally coupled scalar field and certain curvature-squared theories. We find an enhancement in the growth rate of density perturbations in the Brans-Dicke theory, and discuss attractive features of GETs in the structure formation process.


title: Inside the Supernova: A Powerful Convective Engine

abstract: Condensed Abstract: We present an extensive study of the inception of supernova explosions by following the evolution of the cores of two massive stars (15 Msun and 25 Msun) in two dimensions. Our calculations begin at the onset of core collapse and stop several 100 ms after the bounce, at which time successful explosions of the appropriate magnitude have been obtained. (...) Guided by our numerical results, we have developed a paradigm for the supernova explosion mechanism. We view a supernova as an open cycle thermodynamic engine in which a reservoir of low-entropy matter (the envelope) is thermally coupled and physically connected to a hot bath (the protoneutron star) by a neutrino flux, and by hydrodynamic instabilities. (...) In essence, a Carnot cycle is established in which convection allows out-of-equilibrium heat transfer mediated by neutrinos to drive low entropy matter to higher entropy and therefore extracts mechanical energy from the heat generated by gravitational collapse. We argue that supernova explosions are nearly guaranteed and self-regulated by the high efficiency of the thermodynamic engine. (...) Convection continues to accumulate energy exterior to the neutron star until a successful explosion has occurred. At this time, the envelope is expelled and therefore uncoupled from the heat source (the neutron star) and the energy input ceases. This paradigm does not invoke new or modified physics over previous treatments, but relies on compellingly straightforward thermodynamic arguments. It provides a robust and self-regulated explosion mechanism to power supernovae which is effective under a wide range of physical parameters.


title: Gamma-Rays from Grazing Incidence Cosmic Rays in the Earth's Atmosphere

abstract: Interactions of grazing incidence, ultra high energy cosmic rays with the earth's atmosphere may provide a new method of studying energetic cosmic rays with gamma-ray satellites. It is found that these cosmic ray interactions may produce gamma-rays on millisecond time scales which may be detectable by satellites. An extremely low gamma-ray background for transient gamma-ray events and a large area of interaction, the earth's surface, make the scheme plausible. The effective cross section of detection of interactions for cosmic rays above 10^20 eV is found to be more than two orders of magnitude higher than earth based detection techniques. This method may eventually offer an efficient way of probing this region of the cosmic ray energy spectrum where events are scarce. In this paper, a conceptual model is presented for the production of short bursts of gamma-rays based on these grazing incidence encounters with the earth's atmosphere.


title: Color-Magnitude Diagram Distribution of the Bulge Red Clump Stars - Evidence for the Galactic Bar

abstract: The color-magnitude diagrams of $\sim 5\times 10^5$ stars obtained for 13 fields towards the Galactic bulge with the OGLE project reveal a well-defined population of bulge red clump stars. We find that the distributions of the extinction-adjusted apparent magnitudes of the red clump stars in fields lying at $l=\pm5\deg$ in galactic longitude differ by $0.37\pm0.025\;mag$. Assuming that the intrinsic luminosity distribution of the red clump stars is the same on both sides of the Galactic center, this implies that the distances to the red clump stars in the two fields differ by a factor of $1.185\pm0.015$. A plausible explanation of the observed difference in the luminosity distribution is that the Galactic bulge is a triaxial structure, or bar, which is inclined to the line of sight by no more than $45\deg$, with the part of the bar at the positive galactic longitude being closer to us. This agrees rather well with other studies indicating the presence of the bar in the center of the Galaxy. Color-magnitude diagram data are accessible over the computer network with anonymous {\tt ftp}.


title: Observational Prospects for Extra-Galactic Microlensing Events

abstract: We consider the feasibility of directly observing gravitational microlensing in extra-galactic sources, whose stars are not generally resolved. This precludes use of the simple optical depth to microlensing formulation, which is applicable only to resolved stars. We, instead, extend this method to consider observational constraints, such as seeing, sky background and minimum detectable change in magnitude. Our analysis provides quantitative relations between these constraints and the expected observational results, event duration and number of detections per year. We find that an ambitious ground-based observer should detect several events per year in M31. We also consider detection of microlensing in visual binary galaxies. We find that, although these objects may present hundreds of events per year, extremely short durations would yield poor prospects for observation.


title: The Vertical Equilibrium of Molecular Gas in the Galactic Disk

abstract: We examine the vertical structure and equilibrium of the molecular gas layer in the galactic disk, measuring its scale height and velocity dispersion as a function of Galactic radius by modeling the CO emission at the tangent points. The model takes into account emission from a large path length along the line of sight, corresponding to an interval \Delta R ~ 200-400 pc; and is parametrized by the scale height of the gas, the centroid in z, the rotation velocity and the velocity dispersion. This model is then fit to the $^{12}CO$ survey of Knapp et al.(1985) to determine the best fit parameters. The reduced \chi^2 from fitting the models (considering only the photon statistics) range from 2 to 19. The main source of error is the `shot noise' due to the small number of clouds. Simulations using discrete molecular clouds are carried out to estimate the errors. The scale height of the molecular gas increases with radius. The velocity dispersion varies from ~ 2 km/s to 11 km/s with a typical uncertainty of ~3 km/s and shows a monotonic increase with galactic radius. The midplane mass density of the disk calculated from the scale-height and velocity dispersion is consistent with the local value \rho_0(R_0) of 0.2 \msun pc^{-3}.


title: High Altitude Molecular Clouds

abstract: A population of molecular clouds with a significantly greater scale height than that of Giant Molecular Clouds has been identified by examining maps of the latitude distribution of the $^{12}CO(1-0)$ emission in the first quadrant of the Galaxy. These clouds are found by identifying emission more than 2.6 times the scale-height away from the galactic midplane (centroid of CO emission) at the tangent points. Since the distance to the tangent points is known, we know the height and the sizes of these clouds. They are smaller and fainter than the GMCs and do not seem to be gravitationally bound. These clouds have properties similar to the high latitude clouds in the solar neighborhood. Although they lie outside the molecular cloud layer, the high altitude clouds are well within the HI layer in the Galaxy and coincide with distinct peaks in the HI distribution. These clouds represent a galaxy wide population of small molecular clouds having a larger scale height. They may be clouds in transition between molecular and atomic phases.


title: The Cosmological Dependence of Cluster Density Profiles

abstract: We use N-body simulations to study the shape of mean cluster density and velocity profiles formed via gravitational instability. The dependence of the final structure on both cosmology and initial density field is examined, using a grid of cosmologies and scale-free initial power spectra P\propto k^n. For each model, we stack clusters to define an average density profile in the non-linear regime. The profiles are well fit by a power law over 99% of the cluster volume, with a clear trend toward steeper slopes with both increasing n and decreasing Omega_o. For models with a Omega_o = 0.2, the profile slopes are consistently higher than those for Omega-1.0. Cluster density profiles are thus potentially useful cosmological diagnostics. We find no evidence for a constant density core in any of the models, although the density profiles do tend to flatten at small radii. Much of the flattening is due to the force softening required by the simulations, and an attempt is made to recover the unsoftened profiles assuming angular momentum invariance. The recovered profiles in the Omega=1 cosmologies are consistent with a pure power law up to the highest density contrasts (10^6) accessible with our resolution. The low density models show significant deviations from a power law above density contrasts \sim 10^5. We interpret this curvature as reflecting the non scale-invariant nature of the background cosmology in these models.


title: Hydrodynamics of Rotating Stars and Close Binary Interactions: Compressible Ellipsoid Models

abstract: We develop a new formalism to study the dynamics of fluid polytropes in three dimensions. The stars are modeled as compressible ellipsoids and the hydrodynamic equations are reduced to a set of ordinary differential equations for the evolution of the principal axes and other global quantities. Both viscous dissipation and the gravitational radiation reaction are incorporated. We establish the validity of our approximations and demonstrate the simplicity and power of the method by rederiving a number of known results concerning the stability and dynamical oscillations of rapidly rotating polytropes. In particular, we present a generalization to compressible fluids of Chandrasekhar's classical results for the secular and dynamical instabilities of incompressible Maclaurin spheroids. We also present several applications of our method to astrophysical problems of great current interest such as the tidal disruption of a star by a massive black hole, the coalescence of compact binaries driven by the emission of gravitational waves, and the development of instabilities in close binary systems.


title: Momentum Transfer by Astrophysical Jets

abstract: We have used 3-D smoothed particle hydrodynamical simulations to study the basic properties of the outflow that is created by a protostellar jet in a dense molecular cloud. The dynamics of the jet/cloud interaction is strongly affected by the cooling in the shocked gas behind the bow shock at the head of the jet. We show that this cooling is very rapid, with the cooling distance of the gas much less than the jet radius. Thus, although ambient gas is initially driven away from the jet axis by the high thermal pressure odf the post-shock gas, rapid cooling reduces the pressure and the outflow subsequently evolves in a momentum-conserving snowplow fashion. The velocity of the ambient gas is high in the vicinity of the jet head, but decreases rapidly as more material is swept up. Thus, this type of outflow produces extremely high velocity clumps of post shock gas which resemble the features seen in outflows. We have investigated the transfer of momentum from the jet to the ambient medium as a function of the jet parameters. We show that a low Mach number (<6) jet slows down rapidly because it entrains ambient material along its sides. On the other hand, the beam of a high Mach number jet is separated from the ambient gas by a low density cocoon of post-shock gas, and this jet transfers momentum to the ambient medium principally at the bow-shock. In high Mach number jets, as those from young stellar objects, the dominant interaction is therefore at the bow shock at the head of the jet.


title: Internal Velocity and Mass Distributions in Clusters of Galaxies for a Variety of Cosmogonic Models

abstract: The mass and velocity distributions in the outskirts ($0.5-3.0h^{-1}$Mpc) of clusters of galaxies are examined for a suite of cosmogonic models utilizing large-scale Particle-Mesh (PM) simulations ($500^3$ cells, $250^3$ particles and box size of $100h^{-1}$Mpc, giving a nominal resolution of $0.2h^{-1}$Mpc with the true resolution about $0.5h^{-1}$Mpc). Through a series of model computations, designed to isolate the different effects, we find that both $\Omega_0$ and $P_k$ ($\lambda\le 16h^{-1}$Mpc) are important to the mass and velocity distributions in clusters of galaxies. More details in the full abstract.


title: A Connection between Submillimeter Continuum Flux and Separation in Young Binaries

abstract: We have made sensitive 800-micron continuum observations of low-mass, pre-main sequence (PMS) binary stars with projected separations less than 25 AU in Taurus-Auriga to study disks in the young binary environment. We did not detect any of the observed binaries, with typical 3-sigma upper limits of about 30 mJy. Combining our observations with previous 1300-micron observations of PMS Taurus binaries by Beckwith et al. (1990) and others, we find that the submillimeter fluxes from binaries with projected separations between 1 AU and 50 AU are significantly lower than fluxes from binaries with projected separations > 50 AU. The submillimeter fluxes from the wider binaries are consistent with those of PMS single stars. This may indicate lower disk surface densities and masses in the close binaries. Alternatively, dynamical clearing of gaps by close binaries is marginally sufficient to lower their submillimeter fluxes to the observed levels, even without reduction of surface densities elsewhere in the disks.


title: Has a Standard Physics Solution to the Solar Neutrino Problem Been Found? - a Response

abstract: In a recent paper Dar and Shaviv (1994) presented results of calculations with an improved standard solar model (SSM) that suggest a standard physics solution to the solar neutrino problem. In a subsequent publication with a similar title Bahcall et al. (1994) have claimed to refute the results obtained by Dar and Shaviv. Here we show that the criticism and conclusions of Bahcall et al. are unjustified. In particular, their attempts to reproduce the results of Dar and Shaviv failed because they did not use the same nuclear cross sections and did not include in their calculations major physics improvements which were included in the improved SSM code used by Dar and Shaviv.


title: Measurements of Anisotropy in the Cosmic Microwave Background Radiation at 0.5 Degree Angular Scales Near the Star Gamma Ursae Minoris

abstract: We present results from a four frequency observation of a 6 x 0.6 degree strip of the sky centered near the star Gamma Ursae Minoris during the fourth flight of the Millimeter-wave Anisotropy eXperiment (MAX). The observation was made with a 1.4 degree peak-to-peak sinusoidal chop in all bands. The FWHM beam sizes were 0.55 +/- 0.05 degrees at 3.5 cm-1 and 0.75 +/-0.05 degrees at 6, 9, and 14 cm-1. During this observation significant correlated structure was observed at 3.5, 6 and 9 cm-1 with amplitudes similar to those observed in the GUM region during the second and third flights of MAX. The frequency spectrum is consistent with CMB and inconsistent with thermal emission from interstellar dust. The extrapolated amplitudes of synchrotron and free-free emission are too small to account for the amplitude of the observed structure. If all of the structure is attributed to CMB anisotropy with a Gaussian autocorrelation function and a coherence angle of 25', then the most probable values of DeltaT/TCMB in the 3.5, 6, and 9 cm-1 bands are 4.3 (+2.7, -1.6) x 10-5, 2.8 (+4.3, -1.1) x 10-5, and 3.5 (+3.0, -1.6) x 10-5 (95% confidence upper and lower limits), respectively.


title: Biased Power Spectrum and Bispectrum for an Ensemble of Three-Dimensional Scale Free Numerical Simulations

abstract: We examine the effect of a threshold bias on the power spectrum and the bispectrum in an ensemble of numerical simulations (Gaussian initial perturbations with power law spectra P(k) \sim k^n, n=+1, 0, -1, -2) and compare our results with theoretical predictions. Our simulations are evolved sufficiently that on the scale where we apply the threshold the rms fluctuation has developed significantly into the nonlinear regime. Thus, predictions based on perturbation theory do not necessarily apply. Nevertheless, we find our results for the power spectrum, biased power simply amplified by a numerical factor, follow predicted trends, far beyond the regime where perturbation theory is expected to be valid. We find that the biased bispectrum continues to follow the so-called hierarchical form, with reduced three-point amplitude Q \approx 1 in the strongly nonlinear regime, independent of initial spectrum. In the quasi-linear perturbative regime the three-point amplitude depends on configuration shape, a behavior that is found to give useful information about the amount of bias without information about the unbiased matter distribution.


title: Unsteady Outflow Models for Cosmological Gamma-Ray Bursts

abstract: The 'event' that triggers a gamma ray burst cannot last for more than a few seconds. This is, however, long compared with the dynamical timescale of a compact stellar-mass object ($\sim 10^{-3}$ seconds). Energy is assumed to be released as an outflow with high mean lorentz factor $\Gamma$. But a compact stellar-mass collapse or merger is, realistically, likely to generate a mass (or energy) flux that is unsteady on some timescales in the range $10^{-3}$ - 10 seconds. If $\Gamma$ fluctuates by a factor of $\sim 2$ around its mean value, relative motions within the outflowing material will themselves (in the comoving frame) be relativistic, and can give rise to internal shocks. For $\Gamma \sim 10^2$, the resultant dissipation occurs outside the 'photosphere' and can convert a substantial fraction of the overall outflow energy into non-thermal radiation. This suggests a mechanism for cosmological bursts that demands less extreme assumptions (in respect of $\Gamma$-values, freedom from baryonic contamination, etc) than earlier proposals.


title: Yet Another Model of Soft Gamma Repeaters

abstract: We develop a model of SGR in which a supernova leaves planets orbiting a neutron star in intersecting eccentric orbits. These planets will collide in $\sim 10^4$ years if their orbits are coplanar. Some fragments of debris lose their angular momentum in the collision and fall onto the neutron star, producing a SGR. The initial accretion of matter left by the collision with essentially no angular momentum may produce a superburst like that of March 5, 1979, while debris fragments which later lose their angular momentum produce an irregular pattern of smaller bursts.


title: The Standard Cosmology

abstract: These lectures provide an introductory review of big bang cosmology. I discuss the expanding Friedmann-Robertson-Walker universe, summarizing the observational evidence which has led to its adoption as the `standard' cosmological model and reviewing its basic properties. Subsequent lectures provide an overview of the early universe. The final lectures give an introduction to the inflationary universe, beginning with the motivating puzzles of the standard cosmology (the horizon and flatness problems) and ending with the inflationary production of quantum field fluctuations and their possible role in seeding the large-scale structure of the Universe.


title: Numerical analysis of the dynamics of a cosmic string loop as a vortex

abstract: Time evolution of a circular cosmic string loop is investigated by numerically solving the field equations for the scalar and the gauge fields consisting of the vortex. It is shown that the result agrees with an analytic estimate based on the Nambu-Goto action, which supports its validity in analyzing nonstraight and rapidly moving strings.


title: Synthetic and Observed Photometric Indices for Globular Clusters in the Galaxy and M31

abstract: Buzzoni's (1989) grid of synthetic spectral energy distributions, representative of old stellar populations, was used to derive colours in different photometric systems, and to compare the theoretical predictions with the observational data referring to about $120$ globular clusters in the Galaxy and to $159$ objects of the globular cluster system of M31. Synthetic and observed indices display an overall agreement in the composite planes of two-colour diagrams, thus in agreement with the standard evolutionary scenario leading, for globular clusters, to old stellar populations consistent with an age of $15$ Gyr and a Salpeter initial mass function (IMF). The two main parameters modulating the cluster colour distributions are, as known, metallicity and horizontal branch morphology, while IMF slope and mass loss rate from stars in the red-giant branch and asymptotic-giant branch evolutionary stages produce only minor, although not negligible, effects on the integrated colours. The M31 and Galactic cluster populations are found to be substantially similar, at least as far as the spectral energy distribution characteristics in the $U \to V$ spectral range are concerned.


title: A new Contribution to Cosmological Perturbations of some Inflationary Models

abstract: We show that there are inflationary models for which perturbations in the energy momentum tensor, which are of second order in the scalar field, cannot be neglected. We first specify the conditions under which the usual first order perturbations are absent. We then analyze classically, the growth and decay of our new type of perturbations for one mode of fluctuations $\de\phi_k$ in the scalar field. We generalize this analysis, considering the contribution from the whole spectrum of $\de\phi$ to a given wavelength of geometrical perturbations. Finally, we discuss the evolution of the perturbations during the subsequent radiation dominated era and discuss the resulting spectrum of density fluctuations. In the case of a massless scalar field we find a spectral index $n=4$. For massive scalar fields we obtain $n=0$ but the resulting amplitude of fluctuations for inflation around a GUT scale are by far too high. Hence, 'conventional' inflationary models must not be influenced by this new type of perturbations, in order to lead to acceptable matter and CMB fluctuations.


title: Optimized Lagrangian Approximations for Modelling Large--Scale Structure at Non--Linear Stages

abstract: We report on a series of tests of Newtonian Lagrangian perturbation schemes using N--body simulations for various power--spectra with scale--independent indices in the range $-3$ to $+1$. The models have been evolved deeply into the non--linear regime of structure formation in order to probe the dynamical and statistical performance of the Lagrangian perturbation schemes (whose first--order solution contains as a subset the celebrated ``Zel'dovich--approximation'' (hereafter ZA). These tests reveal properties of the approximations at stages beyond the obvious validity of perturbation theory. Recently, another series of tests of different analytical and semi--numerical approximations for large--scale structure was conducted with the result that ZA displays the best dynamical performance in comparison with the N--body simulations, if the initial data were smoothed before evolving the model, i.e., a truncated form of ZA (TZA). We show in this Letter that the excellent performance of TZA can be further improved by going to second order in the Lagrangian perturbation approach. The truncated second--order Lagrangian scheme provides a useful improvement over TZA especially for negative power indices, which suggests it will be very useful for modelling standard scenarios such as ``Cold--'', ``Hot--'' and ``Mixed--Dark--Matter''.


title: Formal Results Regarding Metric-Space Techniques for the Study of Astrophysical Maps

abstract: We extend a newly developed formal system for the description of astrophysical maps. In this formalism, we consider the difference between maps to be the distance between elements of a pseudometric space (the space of all such maps). This ansatz allows us to measure quantitatively the difference between any two maps and to order the space of all maps. For each physical characteristic of interest, this technique assigns an ``output'' function to each map; the difference between the maps is then determined from the difference between their corresponding output functions. In this present study, we show that the results of this procedure are invariant under a class of transformations of the maps and the domains of the maps. In addition, we study the propagation of errors (observational uncertainties) through this formalism. We show that the uncertainties in the output functions can be controlled provided that the signal to noise ratios in the original astrophysical maps are sufficiently high. The results of this paper thus increase the effectiveness of this formal system for the description, classification, and analysis of astrophysical maps.


title: A Quantitative Analysis of IRAS Maps of Molecular Clouds

abstract: We present an analysis of IRAS maps of five molecular clouds: Orion, Ophiuchus, Perseus, Taurus, and Lupus. For the classification and description of these astrophysical maps, we use a newly developed technique which considers all maps of a given type to be elements of a pseudometric space. For each physical characteristic of interest, this formal system assigns a distance function (a pseudometric) to the space of all maps; this procedure allows us to measure quantitatively the difference between any two maps and to order the space of all maps. We thus obtain a quantitative classification scheme for molecular clouds. In this present study we use the IRAS continuum maps at 100$\mu$m and 60$\mu$m to produce column density (or optical depth) maps for the five molecular cloud regions given above. For this sample of clouds, we compute the ``output'' functions which measure the distribution of density, the distribution of topological components, the self-gravity, and the filamentary nature of the clouds. The results of this work provide a quantitative description of the structure in these molecular cloud regions. We then order the clouds according to the overall environmental ``complexity'' of these star forming regions. Finally, we compare our results with the observed populations of young stellar objects in these clouds and discuss the possible environmental effects on the star formation process. Our results are consistent with the recently stated conjecture that more massive stars tend to form in more ``complex'' environments.


title: Gamma Ray Bursts from Interaction of Relativistic Flows with Radiation Fields

abstract: Relativistic flows resulting from sudden explosive events upscatter ambient interstellar photons of local radiation fields. For Lorentz factor $ > 100$ and dense optical - UV radiation fields the emergent signal is a typical gamma ray burst. Presumably the explosions occur in dense globular clusters or in galactic nuclei, at cosmological distances.


title: The Age of the Universe, the Hubble Constant and QSOs in a Locally Inhomogeneous Universe

abstract: A local void in the globally Friedmann-Robertson-Walker cosmological model with the critical density ($\Omega_{0}=1$) is studied. The inhomogeneity is described using a Lema\^{\i}tre-Tolman-Bondi solution for a spherically symmetric distribution of matter. The scale of the central underdense region is $\sim 150$ Mpc. We investigate the effects this has on the cosmological time scale, the measurement of the Hubble constant and the redshift--luminosity distance for moderately and very distant objects ($z \sim 0.1$ and more). The results indicate that if we happened to live in such a void, but insisted on interpreting cosmological observations through the FRW model, we could go wrong in a few instances. For example, the Hubble constant measurement could give results depending on the separation of the source and the observer, the quasars could be younger than we think and also less distant (less energetic).


title: A New Multiplet Table for Fe

abstract: We have recorded spectra of iron-neon and iron-argon hollow cathode lamps in the region 1700A -- 5um (59000 -- 2000cm-1), with Fourier transform (FT) spectrometers at the National Solar Observatory, Tucson, Arizona, U.S.A. and Imperial College, London, U.K., and with a high resolution grating spectrograph at the National Institute of Standards and Technology, Gaithersburg, U.S.A. The uncertainty of the strongest lines in the FT spectra is <0.002cm-1 (0.2mA at 3000A; 8mA at 2um). Pressure and current-dependent shifts are <0.001cm-1 for transitions between low lying levels, increasing to 0.006 cm-1 for transitions between the most highly excited levels. We report 28 new energy levels of Fe I and revised values of another 818 levels. We have identified 9501 lines as due to 9759 transitions in Fe I, and these are presented in the form of a new multiplet table and finding list. This compares with the 5500 lines due to 467 energy levels in the multiplet tables of Moore (1950, 1959). The biggest increase is in the near ultraviolet and near infra-red, and many of the new lines are present in the solar spectrum. Experimental log(gf) values are included where available. A further 125 unidentified lines due to Fe I are given.


title: Analysis of the 3d6 4s(6D)4f-5g Supermultiplet of Fe I in Laboratory and Solar Infrared Spectra

abstract: The combined laboratory and solar analysis of the highly-excited subconfigurations 4f and 5g of Fe I has allowed us to classify 87 lines of the 4f-5g supermultiplet in the spectral region 2545-2585 cm-1. The level structure of these JK-coupled configurations is predicted by semiempirical calculations and the quadrupolic approximation. Semiempirical gf-values have been calculated and are compared to gf values derived from the solar spectrum. The solar analysis has shown that these lines, which should be much less sensitive than lower excitation lines to departures from LTE and to temperature uncertainties, lead to a solar abundance of iron which is consistent with the meteoritic value (A_Fe = 7.51).


title: Hierarchical Clustering and the Butcher-Oemler Effect

abstract: We show that the rapid evolution in the fraction of blue, star-forming galaxies seen in clusters as a function of redshift (the Butcher-Oemler effect) can be explained very simply if structure formation in the universe proceeeds hierarchically. We show that a rich cluster observed at high redshift has had a significantly different evolutionary history to a cluster of the same richness observed today. High redshift clusters take longer to assemble and thus under- go more merging at small lookback times. We have investigated two models of star formation in cluster galaxies: 1) a model in which star formation is induced by galaxy-galaxy mergers and interactions and 2) a model in which star formation is regulated by the infall of galaxies onto larger systems such as groups and clusters. Both models produce trends consistent with the Butcher- Oemler effect. Our models of cluster formation and evolution allow us to make predictions about trends in the observed properties of clusters with redshift. We find that there should be a correlation between the mass of the cluster or group and the strength of the observed Butcher-Oemler effect, with more massive systems exhibiting more evolution than less massive systems. We also predict that both the blue galaxy fraction and the incidence of interacting galaxies in rich clusters should rise continuously with redshift. Finally, we have explored the influence of cosmological parameters on our predictions for cluster evolution. We find that models in which structure formation occurs at very early epochs, such as low $\Omega$ models, predict rather little recent star formation and merging activity in clusters at redshifts of around 0.4.