title: Damped Lyman Alpha Systems and Galaxy Formation

abstract: We examine the constraints on theories of galaxy formation that are obtained from observations of damped $\lya$ (DL) systems, assuming they are gaseous protodisks in dark matter halos. Using the Press-Schechter formalism, we find that the mixed dark matter model, with $\ohdm = 0.3$, $\ocdm = 0.65$, $\obaryon = 0.05$, and $h=0.5$, is ruled out because the number of galactic halos at $z\simeq 3$ is too small to account for the total gaseous mass in DL systems, even under the assumption that all the gas in collapsed halos has settled into disks of neutral gas. The standard CDM model can account for the gas in DL systems if the bias is $b\lsim 2$; the same is true for the CDM model with a cosmological constant, if $b\lsim 1.5$ for $\Lambda = 0.8$. However, one still needs to assume that a fraction $\gsim 0.4$ of the baryons in collapsed halos at $z\simeq 3$ is in the form of neutral gas in disks. We also calculate the column density distribution $f(\nhi)$ of the DL systems, in terms of the surface density profiles of disks and the distribution of their central column densities. It is shown that the form of $f(\nhi)$ at the high end of column density is a diagnostic for the nature of DL systems.


title: Constraints on Models of Galaxy Formation from the Evolution of Damped Ly$\alpha$ Absorption Systems

abstract: There is accumulating observational evidence suggesting that damped Ly$\alpha$ absorption systems systems are the progenitors of present-day spiral galaxies. We use the observed properties of these systems to place constraints on the history of star formation in galactic disks, and on cosmological theories of structure formation in the universe. We show that the observed increase in $\Omega_{HI}$ contributed by damped Ly$\alpha$ systems at high redshift implies that star formation must have been considerably less efficient in the past. We also show that the data can constrain cosmological models in which structure forms at late epochs. A mixed dark matter (MDM) model with $\Omega_{\nu}=0.3$ is unable to reproduce the mass densities of cold gas seen at high redshift, even in the absence of any star formation. We show that at redshifts greater than 3, this model predicts that the total baryonic mass contained in dark matter halos with circular velocities $V_c > 35$ km s$^{-1}$ is less than the observed mass of HI in damped systems. At these redshifts, the photo-ionizing background would prevent gas from dissipating and collapsing to form high column density systems in halos smaller than 35 km s$^{-1}$. MDM models are thus ruled out by the observations.


title: Photometric and H$\alpha$ observations of LSI+61303: detection of a $\sim$26 day V and JHK band modulation

abstract: We present new optical and infrared photometric observations and high resolution H$\alpha$ spectra of the periodic radio star \lsi. The optical photometric data set covers the time interval 1985-1993 and amounts to about a hundred nights. A period of $\sim$26 days is found in the V band. The infrared data also present evidence for a similar periodicity, but with higher amplitude of variation (0\rmag 2). The spectroscopic observations include 16 intermediate and high dispersion spectra of \lsi\ collected between January 1989 and February 1993. The H$\alpha$ emission line profile and its variations are analyzed. Several emission line parameters -- among them the H$\alpha$ EW and the width of the H$\alpha$ red hump -- change strongly at or close to radio maximum, and may exhibit periodic variability. We also observe a significant change in the peak separation. The H$\alpha$ profile of \lsi\ does not seem peculiar for a Be star. However, several of the observed variations of the H$\alpha$ profile can probably be associated with the presence of the compact, secondary star.


title: Dark Matter in ms1224 from Distortion of Background Galaxies

abstract: We explore the dark matter distribution in ms1224.7+2007 using the gravitational distortion of the images of faint background galaxies. Projected mass image reconstruction reveals a highly significant concentration coincident with the X-ray and optical location. The concentration is seen repeatably in reconstructions from independent subsamples and the azimuthally averaged tangential shear pattern is also clearly seen in the data. The projected mass within a $2.76'$ radius aperture is $\simeq 3.5 \times 10^{14} h^{-1} M_{\sun}$. This is $\simeq 3$ times larger than that predicted if mass traces light with $M/L = 275 h$ as derived from virial analysis. It is very hard to attribute the discrepancy to a statistical fluctuation, and a further indication of a significant difference between the mass and the light comes from a second mass concentration which is again seen in independent subsamples but which is not seen at all in the cluster light. We find a mass per galaxy visible to $I = 22$ of $\simeq 8 \times 10^{12}h^{-1}M_{\sun}$ which, if representative of the universe, implies a density parameter $\Omega \sim 2$. We find a null detection of any net shear from large-scale structure with a precision of 0.9\% per component. This is much smaller than the possible detection in a recent comparable study, and the precision here is comparable to to the minimum level of rms shear fluctuations implied by observed large-scale structure.


title: Anisotropic Scattering of OH/IR Stars Towards the Galactic Center

abstract: We have made angular broadening measurements in the satellite line of the OH molecule at 1612 MHz of 8 OH/IR stars located near the Galactic Center with the Very Large Array. Early observations with the Very Long Baseline Array had revealed a region of pronounced scattering approximately centered on SgrA*. These new observations show that the compact line emission from these stars is scattered by an amount that is two to three times larger than had been previously seen. Additionally, we measure non-circular scattering disks for these objects, evidence that the electron density variations are distributed anisotropically. The present data suggest that the enhanced turbulence is located near the Galactic Center. However, we cannot completely rule out that the scattering is caused by an unrelated region along the line of sight.


title: Computational aspects of the gravitational instability problem for a multicomponent cosmological medium

abstract: The paper presents results for deriving closed-form analytic solutions of the non-relativistic linear perturbation equations, which govern the evolution of inhomogeneities in a homogeneous spatially flat multicomponent cosmological model. Mathematical methods to derive computable forms of the perturbations are outlined.


title: Astrophysical thermonuclear functions

abstract: As theoretical knowledge and experimental verification of nuclear cross sections increases it becomes possible to refine analytic representations for nuclear reaction rates. In this paper mathematical/statistical techniques for deriving closed-form representations of thermonuclear functions are summarized and numerical results for them are given.The purpose of the paper is also to compare numerical results for approximate and closed-form representations of thermonuclear functions.


title: Second-Order Reconstruction of the Inflationary Potential

abstract: To first order in the deviation from scale invariance the inflationary potential and its first two derivatives can be expressed in terms of the spectral indices of the scalar and tensor perturbations, $n$ and $n_T$, and their contributions to the variance of the quadrupole CBR temperature anisotropy, $S$ and $T$. In addition, there is a ``consistency relation'' between these quantities: $n_T= -{1\over 7}{T\over S}$. We discuss the overall strategy of perturbative reconstruction and derive the second-order expressions for the inflationary potential and its first two derivatives and the first-order expression for its third derivative, all in terms of $n$, $n_T$, $S$, $T$, and $dn/d\ln k$. We also obtain the second-order consistency relation, $n_T =-{1\over 7}{T\over S}[1+0.11{T\over S} + 0.15 (n-1)]$. As an example we consider the exponential potential, the only known case where exact analytic solutions for the perturbation spectra exist. We reconstruct the potential via Taylor expansion (with coefficients calculated at both first and second order), and introduce the Pad\'{e} approximant as a greatly improved alternative.


title: High Resolution Observations of Cepheus a

abstract: New high resolution 6 cm observations have been made on Cepheus A using MERLIN, and combined with new VLA observations at 3.7 cm. Angular resolution with the latter was 200 mas, and with MERLIN was 60 mas, except for isolated unresolved sources where 33 mas was achieved. Unresolved objects at 60 mas were observed in Sources 2, 3, and in particular 9 which also was not resolved at 33 mas. There is no evidence for any other object as small as this with any significant flux density, although Source 8 was quiescent at the time. The upper limit to the size of Source 9 sets a minimum brightness temperature of 4.3 $\times$ 10$^5$K, and adds credence to a previous suggestion that it is a gyrosynchrotron source. The compact objects of Sources 2 and 3 are thought to be produced by mass outflow from stars, which could be of spectral type B0 - B1, but this is uncertain. There is a discussion regarding the powerhouse for the molecular outflow. Sources 8 and 9, which are the highly time dependent objects, appear at the centre of the disruption of the high density gas, and their estimated high temperature of 10$^7$ - 10$^8$K indicates that they could produce high velocity winds. On the other hand, the OH masers surrounding Source 2(ii) show an outward velocity of about 10 km s${-1}$, which is small, but higher velocity winds could tunnel through adjacent spaces, and even be responsible for the 300 km s${-1}$ bullet of Source 7. An extrapolation of the orthogonals to the IR polarization vectors are not accurate enough to pinpoint the source of the IR radiation, but it is estimated that Sources 2(ii) and 3(d)(ii) could


title: Signatures of Cosmic Strings on the Microwave Background

abstract: ABSRACT: We review recent progress on testing the hypothesis of the existence of cosmic string perturbations in microwave background maps. Using an analytical model for the string network we show that the predicted amplitude and spectrum of MBR fluctuations are consistent with the COBE data for a reasonable value of the single free parameter of the string model (the mass per unit length of the string $\mu$). To {\it distinguish} the predictions of cosmic strings from those of Gaussian models it is necessary to apply specific {\it statistical tests} which are sensitive to non-random phases of the MBR temperature maps. We discuss two such statistical tests: First, the probability distribution and the moments of fluctuations ${{\delta T}\over T}$ and second the corresponding statistical quantities for the {\it gradient} of fluctuations along a fixed direction. We show that the second statistic can detect cosmic string specific signatures on MBR temperature maps with resolution of a few arcminutes or smaller.


title: Spectral Analysis of Microwave Background Perturbations Induced by Cosmic Strings

abstract: Using a simple analytic model based on seed superposition we obtain the spectrum of microwave background perturbations induced by cosmic strings on all angular scales larger than about 2 armin. We assume standard recombination in an Einstein de Sitter universe with $h=1/2$ and study the fluctuation spectrum along a great circle in the sky. Doppler and potential perturbations on the last scattering surface (Sachs-Wolfe effect) are shown to dominate over post-recombination perturbations on scales smaller than about 2 degrees. Using a filter function corresponding to the COBE experiment we obtain an effective power spectrum index $n_{eff}\simeq 1.35$ in good agreement with the recently announced second year COBE data showing $n_{eff}\simeq 1.5$. The only free parameter ($G\mu$) of the string model is fixed by normalizing on the COBE detection leading to $G\mu=1.6$. Other parameters (e.g. the {\it rms} string velocity) are fixed by comparing with string simulations. Using these parameter values we compare the {\it rms} temperature fluctuations $({{\Delta T}\over T})_{\it rms}$ predicted for ongoing experiments (Tenerife, SP91, MAX, OVRO etc) with detections and with the corresponding predictions of inflationary models.


title: The Distribution of Galactic Disk Stars in Baade's Window

abstract: The color-magnitude diagram of $ \sim 3 \times 10^5 $ stars obtained for Baade's Window towards the Galactic bulge with the OGLE project reveals a surprisingly narrow main sequence due to galactic disk stars at a distance of $d \sim 2 $ kpc, i.e. at the location of the Sagittarius spiral arm. A more careful analysis indicates there is an excess in the number of disk stars by a factor $\sim 2$ between us and $d \sim 2.5$ kpc and a rapid drop by a factor $\sim 10$ beyond that distance. It is unlikely that the observed structure is an artifact of the interstellar extinction, but careful determination of the extinction is needed before the structure is firmly established. The narrow main sequence extends down to stars as faint as $M_{_V} \approx 7$, i.e. it is composed of old stars. This is not expected in a conventional disk model, or in a conventional model of a spiral structure. However, a strong concentration of old stars towards spiral arms has been noticed in some other galaxies, like M51 (Rix \& Rieke 1993), with the near infrared surface photometry. We have also found that the relative distribution of red clump and red giant stars in the Galactic bulge implies that there is a relatively young stellar population present there.


title: An Alternative View of Flat Rotation Curves. II. The Observations

abstract: The rotation curves of 20 spiral galaxies are examined in the light of a toy model (Soares 1992) which has as the main feature the assignment of a high M/L ratio (=30; Ho=50 km/s/Mpc) to the visible matter. The observed rotation of all galaxies can be accommodated without the assumption of a dark halo. Moreover, the suggestion is made that the fact that almost all available rotational velocity measurements are derived from emission lines emitted by galaxian gas (either neutral or ionized) makes them inappropriate as tracers of the galaxy gravitational potential. To account for that, the model introduces an effective potential meant to describe the hydrodynamics inside a gaseous disk. The general morphology of the curves (i.e., the presence of a plateau in V(r) X r) is interpreted in this framework as a consequence of the hydrodynamical characteristics of galaxian disks. The Tully-Fisher relation is expressed in terms of model parameters and used as an additional constraint in the process of fitting the model to the observed rotation of the galaxies.


title: Gravitational Microlensing by Dark Clusters in the Galactic Halo

abstract: The dark matter in Galactic halos, or some fraction of it, may be in the form of dark clusters which consist of small mass objects. Carr & Lacey (1987) have derived the permissible properties of such systems, and proposed the existence of dark clusters with mass of order $10^6\solarmass$ to explain some of the observed dynamical properties of the stellar disk of the Galaxy. A population of bound systems with mass of $\sim 10^5-10^6\solarmass$ is also an attractive possibility since it is close to the baryon Jeans mass at recombination, which may be the preferred mass scale for the first bound objects to form in the universe. At the present, the existence of dark clusters which consist of brown dwarfs, Jupiters, or black hole remnants of an early generation of stars, is not indicated, nor can be excluded on observational grounds. We describe how dark clusters can be discovered in a sample of gravitational microlensing events in LMC stars. Alternatively, it could provide strict bounds on the fraction of halo mass which resides in such systems. If MACHOs are clustered, the implied degeneracy in their spatial and velocity distributions would result in a strong autocorrelation in the sky position of microlensing events on an angular scale $\lesssim 20$ arcsec, along with a correlation in the event duration. We argue that a small number of events could be enough to indicate the existence of clusters, and demonstrate that a sample of $\simeq 10$ events would be sufficient to reject the proposal of Carr & Lacey (1987) at the $95\%$ confidence level.


title: Is a Massive Tau Neutrino Just What Cold Dark Matter Needs?

abstract: The cold dark matter (CDM) scenario for structure formation in the Universe is very attractive and has many successes; however, when its spectrum of density perturbations is normalized to the COBE anisotropy measurement the level of inhomogeneity predicted on small scales is too large. This can be remedied by a tau neutrino of mass $1\MeV - 10\MeV$ and lifetime $0.1\sec - 100\sec$ whose decay products include electron neutrinos because it allows the total energy density in relativistic particles to be doubled without interfering with nucleosynthesis. The anisotropies predicted on the degree scale for ``$\tau$CDM'' are larger than standard CDM. Experiments at $e^\pm$ colliders may be able to probe such a mass range.


title: Microlensing by brown dwarfs: the case of dark clusters

abstract: Experiments are looking for diffuse brown dwarfs in the dark galactic halo through the gravitational lens effect. If brown dwarfs are clumped in dark clusters, the event rate is not changed, but events are spatially clustered, and stars nearby a micro-lensed one are likely to be micro-lensed in a near future. Therefore, an intensive survey of the region where a micro-lensing occurred should reveal many other events.


title: Multiple High-Velocity Emission-Line Systems in the E+S Pair CPG 29

abstract: A detailed study of the mixed-morphology galaxy pair CPG 29 (Arp 119, VV 347) shows spectacular spectroscopic peculiarities in the southern (spiral) component (Mkn 984) including a spatially resolved region, roughly aligned along the minor axis of the galaxy, with multiple emission-line components redshifted by as much as 1300 km s$^{-1}$ with respect to the LINER nucleus. Strong [OI]\l6300 and [SII]\l\l6716,6731 emission suggest shock ionization. The rest of the galaxy's disk is spectroscopically undistinguished, with emission lines characteristic of disk H II regions and, in large part, ordered rotation. We discuss the following four possible explanations for the morphological and spectroscopic peculiarities in the spiral: (1) a super--wind scenario, supported by the high FIR luminosity of the spiral, and by emission line activity similar to that observed in FIR strong galaxies; (2) motion of the spiral through a hot, relatively dense intergalactic medium, possibly associated with the elliptical, that could affect the star formation rate and create a brightened disk rim; (3) a faint companion in direct collision with the disk, at $\Delta v_r \sim 1000$ \kms, that would straightforwardly explain the existence of the different redshift systems; (4) a pole on collision by the elliptical, that could produce the disturbed morphology and other properties of the spiral. The elliptical/spiral pole--on collision explains the high velocity line emitting regions as tidally stripped gas, and accounts for the presence of gas ionized by moterate velocity shocks. The morphology of Mrk 984 is consistent


title: Monopoles as Big as a Universe

abstract: We show that, contrary to the standard belief, primordial monopoles expand exponentially during inflation in the new inflationary universe scenario. Moreover, inflation of monopoles continues without an end even when inflation ends in the surrounding space. Therefore primordial monopoles (as well as other topological defects produced during inflation) can serve as seeds for the process of eternal self-reproduction of inflationary universe.


title: How Cosmic Background Correlations at Large Angles Relate to Mass Autocorrelations in Space

abstract: The Sachs-Wolfe effect is known to produce large angular scale fluctuations in the Cosmic Microwave Background Radiation (CMBR) due to gravitational potential fluctuations. We show how the angular correlation function of the CMBR can be expressed explicitly in terms of the mass autocorrelation function $\xi (r)$ in the Universe. We derive analytic expressions for the angular correlation function and its multipole moments in terms of integrals over $\xi (r)$ or its second moment, $J_3 (r)$, which does not need to satisfy the sort of integral constraint that $\xi (r)$ must. We derive similar expressions for bulk flow velocity in terms of $\xi$ and $J_3$. One interesting result that emerges directly from this analysis is that, for angles $\theta$, there is a substantial contribution to the correlation function from a wide range of distance $r$ and that the radial shape of this contribution does not vary greatly with angle.


title: A Singularity-Free Cosmological Model with a Conformally Coupled Scalar Field

abstract: We explore the possibility of describing our universe with a singularity--free, closed, spatially homogeneous and isotropic cosmological model, using only general relativity and a suitable equation of state which produces an inflationary era. A phase transition to a radiation--dominated era occurs as a consequence of boundary conditions expressing the assumption that the temperature cannot exceed the Planck value. We find that over a broad range of initial conditions, the predicted value of the Hubble parameter is approximately $47$ km$\cdot$~s$^{-1}\cdot$~Mpc$^{-1}$. Inflation is driven by a scalar field, which must be conformally coupled to the curvature if the Einstein equivalence principle has to be satisfied. The form of the scalar field potential is derived, instead of being assumed a priori.


title: Giant Explosion at the Galactic Center and Huge Shocked Shells in the Halo

abstract: We simulate the propagation of a shock front through the galactic halo induced by an explosion and/or a starburst at the galactic center. A huge dumbbell ($\Omega$)-shaped shock front produced by an explosion of energy $3\times 10^{56}$ ergs about fifteen million years ago can mimic the radio and X-ray North Polar Spur as well the southern large X-ray spur. The post-shock high-temperature gas in the corona will explain the observed X-ray bulge around the galactic center. Key workds Galaxy: center -- Galaxy: halo -- Radio continuum: general -- Shock wave -- X-rays: general


title: The Gunn-Peterson Effect in the Spectrum of the Z=4.7 QSO 1202-0725: The Intergalactic Medium at Very High Redshift

abstract: A measure of the average depression between Lyman absorption lines in the spectrum of the faint quasar BR1202-0725 ($z_{em}=4.695$) is presented. The relatively high resolution of the spectrum ($\sim 40$ km s$^{-1}$) allows the selection of regions free of strong absorption lines in the Lyman alpha forest. A reliable evaluation of the continuum shape is based on the careful flux calibration and on the large wavelength interval covered (5000--9300 ~\AA). A best estimate of $\tau _{_{GP}}\leq 0.02\pm 0.03$ has been found for the Gunn-Peterson optical depth at the highest absorption redshift observed at this resolution, $z \simeq 4.3$. The derived baryon density of the intergalactic medium is $\Omega _{IGM}\mincir 0.01$ if the observed quasars are the major contributor to the ionizing UV background flux. This limit, when compared with the total baryon density deduced from the nucleosynthesis, could imply that most of the baryons are already in bound systems at $z\sim 5$.


title: Lensing from the Light-Traces-Mass Map of MS1224+20

abstract: The MS1224+20 cluster of galaxies is a high luminosity X-ray source at z=0.325. To compare the lensing mass to the virial mass we have completed a uniform, high precision redshift survey over a field of $7\arcm\times9\arcm$, obtaining 75 redshifts. The velocity dispersion of 30 cluster galaxies is $775$ \kms\ and the projected harmonic radius is 0.32 \hmpc. The virial mass is $2.1\times 10^{14} h^{-1} \msun$. Correcting for faint cluster galaxies without redshifts and allowing for a modest evolution in the galaxy luminosity function gives a current epoch mass-to-light ratio $M/L_V(0)=255h \msun/\lsun$. The same field contains a $z=0.225$ cluster with $\sigma_v\simeq 500$ \kms, and a group at $z=0.412$ with $\sigma_v\simeq 400$ \kms. The clusters' gravitational field induces image ellipticities that are calculated from the light-traces-mass density distribution and compared to the observed average tangential distortions from Fahlman \et (1994). Between 1.5 and 3 virial radii, the observed lensing distortions are $4\pm2$ times stronger than the light-traces-mass model and the virial M/L predict.


title: Window Functions for CMB Experiments

abstract: We discuss the applicability and derivation of window functions for cosmic microwave background experiments on large and intermediate angular scales. These window functions describe the response of the experiment to power in a particular mode of the fluctuation spectrum. We give general formulae, illustrated with specific examples, for the most common observing strategies.


title: Reconstructing the CMB power spectrum

abstract: We discuss a method for a model independent reconstruction of the CMB temperature fluctuation power spectrum on small and intermediate scales that is geared to individual experiments. The importance of off-diagonal correlations for determining the shape of the power spectrum is emphasized and some examples of a reconstruction method are given. By using this method to ``map'' the power spectrum on the scales to which they are sensitive several experiments could be combined to map out the full power spectrum, and test consistency of observed features. For example we find that the GUM scan of the 3rd flight of the MAX experiment prefers a positive slope to the power spectrum near $\ell\simeq160$, which provides weak evidence for the presence of a Doppler peak.


title: A General Class of Self-similar Self-gravitating Fluids

abstract: (accepted for publication in the Ap.J.) I present a general classification of self-similar solutions to the equations of gravitational hydrodynamics that contain many previous results as special cases. For cold flows with spherical symmetry, the solution space can be classified into several regions of behavior similar to the Bondi solutions for steady flow. A full description of these solution is possible, which serves as the asymptotic limit for the general problem. By applying a shock jump condition, exact general solutions can be constructed. The isothermal case allows an extra exact integral, and can be asymptotically analyzed in the presence of finite pressure. These solutions serve as analytic models for problems such as spherical accretion for star formation, infall or outflow of gas into galaxies, Lyman alpha cloud dynamics, etc. Most previous self-similar results are obtained as special cases. The critical values for a cosmological flow with $\Omega=1$ and $\gamma=4/3$ turn out to play a special role.


title: The contribution of galaxies to OMEGA_0

abstract: We estimate the contribution of galaxies to the cosmological density Omega_0 including also their extended dark halos. We find Omega_gal< 0.1, so implying that the matter in galaxies cannot by itself close the Universe but marginally could be baryonic. (latex file)two figures available upon request).


title: How Well Can Cosmological Parameters Be Estimated from CMB Observations?

abstract: The CMB anisotropy depends sensitively upon the slope and amplitude of primordial density and gravitational wave fluctuations, the baryon density, the Hubble constant, the cosmological constant, the ionization history, {\it etc.} We report on recent work showing how well multi-scale measurements of the anisotropy power spectrum can resolve these factors. We identify a hypersurface in cosmic parameter space that can be accurately localized by observations, but along which the likelihood will hardly vary. Other cosmic observations will be needed to break this degeneracy.


title: Interactions of Cosmic Ray Nuclei

abstract: We present convenient formulae for the energy losses of energetic atomic nuclei over the entire energy range relevant to the physics of cosmic rays. Results are applied to a leaky-box equation with a complete loss term. Thereby we derive the equilibrium spectrum of cosmic rays in various types of galaxies. We emphasize a spectral break energy at 450 MeV independent of the matter density, resulting from the transition from Coulomb and ionization losses to pion production losses as the relevant cooling process for the cosmic ray nuclei . We comment on the possible cosmic ray origin of the cosmic gamma ray background.


title: CMB Broad-Band Power Spectrum Estimation

abstract: The natural outcome of theoretical calculations of microwave background anisotropy is the angular power spectrum ${\cal C}_\ell$ as a function of multipole number $\ell$. Experimental ${\cal C}_\ell$'s are needed for direct comparison. Estimation procedures using statistics linear in the pixel amplitudes as well as the conventional but less useful quadratic combinations are described. For most current experiments, a single broad-band power amplitude is all that one can get with accuracy. Results are given for the Capri-meeting detections. Mapping experiments, sensitive to many base-lines, can also give spectral ``colour'' information, either with a series of contiguous narrow-band powers or as parameterized by a local ``colour'' index $n_{\Delta T}$ (scale invariant is -2, white noise is 0). Bayesian analyses of the full first year DMR and FIRS maps give very similar band-powers (\eg $Q_{rms,PS}=17.9 \pm 2.9 \mu K$ {\it c.f.} $18.6 \pm 4.7 \mu K$ for $n_{\Delta T}=-2$) and colour indices (with 1 and 2 sigma error bars) $n_{\Delta T}+3=2.0^{+0.4;+0.7}_{-0.4;-1.0}$ and $ 1.8^{+0.6;+0.9}_{-0.8;-1.3}$ ({\it c.f.} the value 1.15 for a ``standard'' scale invariant CDM model). The 53 and 90 GHz DMR maps, as well as the FIRS map, have residual short-distance noise which steepens $n_{\Delta T}$. Allowing the pixel error bars to increase absorbs much of the residual, but further exploration is needed to see if a second residual evident in the data which is responsible for the high $n_{\Delta T}$ is from systematic errors or is physical.


title: Topological Defects in an Open Universe

abstract: (To appear in Nuclear Physics B Supplements Proceedings section) This talk will explore the evolution of topological defects in an open universe. The rapid expansion of the universe in an open model slows defects and suppresses the generation of CBR fluctuations at large angular scale as does the altered relationship between angle and length in an open universe. Defect models, when normalized to COBE in an open universe, predict a galaxy power spectrum consistent with the galaxy power spectrum inferred from the galaxy surveys and do not require an extreme bias. Neither defect models in a flat universe nor standard inflationary models can fit either the multipole spectrum or the power spectrum inferred from galaxy surveys.


title: Power Spectrum Constraints from Spectral Distortions in the Cosmic Microwave Background

abstract: %The content of this replacement paper is identical to the original. %We have attempted to fix the postscript so that it will print out on %a larger number of printers. Using recent experimental limits on $\mu$ distortions from COBE FIRAS, and the large lever-arm spanning the damping of sub-Jeans scale fluctuations to the scale of the COBE DMR fluctuations, we set a constraint on the slope of the primordial power spectrum $n$. It is possible to analytically calculate the contribution over the full range of scales and redshifts, correctly taking into account fluctuation growth and damping as well as thermalization processes. We find that the 95\% upper limit is weakly dependent on cosmological parameters, e.g. $n<1.54 (h=0.5)$ and $n<1.56 (h=1.0)$ for $\Omega_0=1$ with marginally weaker constraints for $\Omega_0<1$ in a flat $\Omega_0 +\Omega_\Lambda=1$ universe.


title: Are Analogues of Hot Subdwarf Stars Responsible for the UVX Phenomenon in Galaxy Nuclei?

abstract: We present the case that populations of sdB/sdOB/sdO-type stars may be a common constituent of galactic stellar populations, responsible for the UV upturn (`UVX') observed in the spectra of spiral bulges and normal galaxy nuclei. Extreme Horizontal Branch stars with $\log g > 5$ and $\log \teff > 20,000K$ have emerged in the last few years as the most likely candidate for the origin of the UVX. The magnitude of this far-UV flux in some systems (e.g. NGC~1399, NGC~4649) indicates that galactic nuclear regions must contain larger numbers of these subdwarfs than does the solar neighbourhood. This paper summarizes the results of a quantitative study of the UV radiation from evolved stellar populations. We have computed a large grid of stellar models in advanced stages of evolution, as well as a set of isochrones for ages 2-20 Gyr, for a wide range in composition. We use these calculations to derive synthetic UV colour indices for stellar populations with hot components.


title: Is the High-Energy Emission from Centaurus A Compton-Scattered Jet Radiation?

abstract: We consider whether the hard X-ray and soft gamma-ray emission from Centaurus A is beamed radiation from the active nucleus which is Compton-scattered into our line-of-sight. We derive the spectrum and degree of polarization of scattered radiation when incident beamed radiation is scattered from a cold ($kT<<m_ec^2$) electron cloud moving with bulk relativistic motion along the jet axis, and calculate results for an unpolarized, highly-beamed incident power-law photon source. We fit the OSSE data from Centaurus A with this model and find that if the scatterers are not moving relativistically, then the angle the jet makes with respect to our line-of-sight is $61^\circ\pm 5^\circ$. We predict a high degree of polarization of the scattered radiation below $\sim300$ keV. Future measurements with X-ray and gamma-ray polarimeters could be used to constrain or rule out such a scenario.


title: Gravitational Lensing of Distant Field Galaxies by Rich Clusters: I. -- Faint Galaxy Redshift Distributions

abstract: {}From deep optical images of three clusters selected by virtue of their X-ray luminosity and/or optical richness (1455+22; $z=0.26$, 0016+16; $z=0.55$ and 1603+43; $z=0.89$), we construct statistically-complete samples of faint field galaxies ($I \leq 25$) suitable for probing the effects of gravitational lensing. By selecting clusters across a wide redshift range we separate the effects of the mean redshift distribution of the faint field population well beyond spectroscopic limits and the distribution of dark matter in the lensing clusters. A significant lensing signature is seen in the two lower redshift clusters whose X-ray properties are well-constrained. Based on these and dynamical data, it is straightforward to rule out field redshift distributions for $I \leq 25$ which have a significant low redshift excess compared to the no evolution prediction, such as would be expected if the number counts at faint limits were dominated by low-$z$ dwarf systems. The degree to which we can constrain any high redshift tail to the no evolution redshift distribution depends on the distribution of dark matter in the most distant lensing cluster. In the second paper in this series, we use the lensing signal to reconstruct the full two-dimensional mass distribution in the clusters and, together with high resolution X-ray images, demonstrate that their structural properties are well-understood. The principal result is therefore the absence of a dominant low-$z$ dwarf population to $I \leq25$.


title: Gravitational Lensing of Distant Field Galaxies by Rich Clusters: II. -- Cluster Mass Distributions

abstract: Using a non-parametric procedure developed by Kaiser \& Squires (1993), we analyse the statistical image distortions of faint field galaxies to $I\ls25.5$ in two distant X-ray luminous clusters 1455+22 ($z_{cl}=0.26$) and 0016+16 ($z_{cl}=0.55$) to derive two dimensional projected mass distributions for the clusters. The mass maps of 1455+22 and 0016+16 are presented at effective resolutions of 135 kpc and 200 kpc respectively (for $H_o$=50 kms sec$^{-1}$ Mpc$^{-1}$, $q_o=0.5$) with a mean signal to noise per resolution element of 17 and 14. We compare our 2-D mass distributions on scales up to $\sim$1 Mpc with those defined by the spatial distribution of colour-selected cluster members and from deep high resolution X-ray images of the hot intracluster gas. Despite the different cluster morphologies, one being cD-dominated and the other not, in both cases the form of the mass distribution derived from the lensing signal is strikingly similar to that traced by both the cluster galaxies and the hot X-ray gas. We find some evidence for a greater central concentration of dark matter with respect to the galaxies. The overall similarity between the distribution of total mass and that defined by the baryonic components presents a significant new observational constraint on the nature of dark matter and the evolutionary history of rich clusters.


title: The Formation of Galaxies: A Challenge for Supercomputers -- a Simple Task for GRAPE ?

abstract: We present numerical simulations of galaxy formation, one of the most challenging problems in computational astrophysics. The key point in such simulations is the efficient solution of the N--body problem. If the gas of a galaxy is treated by means of smoothed particle hydrodynamics (SPH), the hydrodynamic equations can be reduced to a form similar to that of the N--body problem. A straightforward implementation requires a computational effort $\propto N^2$, making it prohibitively expensive to simulate systems larger than $10^5$ particles even on the largest available supercomputers. After a description of the physical and numerical problems, we shortly review the standard numerical methods to tackle these problems and discuss their advantages and drawbacks. We also present a completely different approach to perform such simulations using a workstation in combination with the special purpose hardware {\sc Grape}. After a discussion of the main features of {\sc Grape}, we present a simple implementation of a {\sc SPH}--N--body code on such a configuration. Comparing results and performance of these two approaches, we show, that with an investment of US \$ 50000, the problem can be solved up to 5 times faster than on a {\sc Cray} YMP.


title: Spherical Accretion onto Neutron Stars Revisited: Are Hot Solutions Possible ?

abstract: Stationary, spherical accretion onto an unmagnetized neutron star is here reconsidered on the wake of the seminal paper by Zel'dovich \& Shakura (1969). It is found that new ``hot'' solutions may exist for a wide range of luminosities. These solutions are characterized by a high temperature, $10^{9}\div 10^{11}$ K, and arise from a stationary equilibrium model where the dominant radiative mechanisms are multiple Compton scattering and bremsstrahlung emission. For low luminosities, $\mincir 10^{-2} \ L_{E}$, only the ``cold'' (\`a la Zel'dovich and Shakura) solution is present.


title: Typical Scales in the Spatial Distribution of QSOs

abstract: We present results of searching for the possible typical scales in the spatial distribution of QSOs. Our method is based on the second derivative of the two-point correlation function. This statistic is sensitive to the scale of the maximum in the spectrum $P(k)$ of the density perturbation in the universe. This maximum or bend scale can be detected as the wavelengths of the periodic component in the second derivative of the integral correlation function. For various QSO samples compiled from surveys of pencil-beam and bright QSOs, a typical scale of about 93 $\pm$ 10 h$^{-1}$Mpc for $q_0=0.5$ has been detected. This typical scale is in good agreement with that found in the spatial distributions of galaxies, clusters of galaxies, and CIV absorption systems of QSOs if $q_0$ is taken to be $\sim 0.2$. Therefore, it is likely a common or universal scale in the large scale structure traced by these objects. This result is consistent with the assumption that the typical scale comes from a characteristic scale in the spectrum of the density perturbation in the universe.


title: Noise Correlation in Cosmic Microwave Background Experiments

abstract: Many analyses of microwave background experiments neglect the correlation of noise in different frequency or polarization channels. We show that these correlations, should they be present, can lead to severe misinterpretation of an experiment. In particular, correlated noise arising from either electronics or atmosphere may mimic a cosmic signal. We quantify how the likelihood function for a given experiment varies with noise correlation, using both simple analytic models and actual data. For a typical microwave background anisotropy experiment, noise correlations at the level of 1\% of the overall noise can seriously {\it reduce} the significance of a given detection.


title: Evidence for Hco+ Infall toward T Tauri?

abstract: High spectral and spatial resolution observations of the HCO+ molecule toward the T Tau binary system obtained with the Owens Valley Millimeter Array reveal a broad emission line profile (V_LSR = 7.2 km/s, width = 2.0 km/s) upon which a narrow, redshifted absorption feature (V_LSR = 8.5 km/s, width +/- 0.5 km/s) is superposed. One possible interpretation of the absorption is that it arises from infall of molecular cloud envelope material onto the immediate circumbinary region, which would imply that large scale dynamic accretion processes are still active well into the optically visible stages of star formation. Provided the mass is deposited at radii sufficiently close to the young stars, estimates of the accretion rates indicate that infall may add enough material to the system to account for the FU~Orionis type flares recently observed in T~Tau in =< 100 yr.


title: Discovery of hydroxyl and water masers in R Aquarii and H1-36 Arae

abstract: We present the first results from an all-sky maser-line survey of symbiotic Miras. Interferometric spectral-line observations of R Aqr and H1-36 Arae have revealed a 22-GHz water maser in the former and 1612-MHz hydroxyl and weak 22-GHz water maser emission from the latter. H1-36 has thus become the first known symbiotic OH/IR star. We have also detected weak OH line emission from the vicinity of R Aqr, but we note that there are small discrepencies between the OH- and H2O-line velocities and positions. These detections demonstrate unequivocally that dust can shield some circumstellar hydroxyl and water molecules from dissociation, even in systems which possess intense local sources of UV. Finally, we discuss some of the implications of these observations. The narrow profile of the water maser in R Aqr means that there may finally be an opportunity to determine the system's orbital parameters. We also point out that high resolution synthesis observations may trace the distribution of dust in H1-36 and R Aqr, possibly throwing light on the mass-loss process in symbiotic Miras and placing constraints on the amount of collimation experienced by UV radiation from their hot, compact companions.


title: A search for symbiotic behaviour amongst OH/IR colour mimics

abstract: Recent maser surveys have shown that many potential OH/IR stars have no OH masers in their circumstellar envelopes, despite the modest requirements which should be implicitly met by IRAS colour-selected candidates. It has been suggested that these OH/IR colour mimics must have a degenerate companion which dissociates OH molecules and disrupts the masing action, ie. that they are related to symbiotic Miras. Coincidentally, there is a paucity of long-period symbiotic Miras and symbiotic OH/IR stars. Phenomonologically, those that are known seem to cluster in the zone where field Miras transform into OH/IR stars. If it could be proven that OH/IR colour mimics contain a degenerate star, that observable evidence of this star is hidden from view by CS dust whilst it slowly accretes from the wind of its Mira companion, then we have an excellent explanation for not only the existence of OH/IR colour mimics, but also for the low observed frequency of symbiotic OH/IR stars and the common occurrence of very slow novae in long-period symbiotic Miras. Here, we employ radio continuum radiation (which should escape unhindered from within the dust shells) as a simple probe of the postulated hot degenerate companions which would inevitably ionize a region of their surrounding gas. We compare the radio and infrared properties of the colour mimics with those of normal symbiotic Miras, using the strong correlation between radio and mid-IR emission in symbiotic stars. We show that if a hot companion exists then, unlike their symbiotic counterparts, they must produce radiation-bounded nebulae. Our observations provide no support for the above scenario for the lack of observed masers, but neither do they permit a rejection of this scenario.


title: Comparison of Dynamical Approximation Schemes for Non-Linear Gravitational Clustering

abstract: I report on controlled comparison of gravitational approximation schemes linear/lognormal/adhesion/frozen-flow/Zel'dovich(ZA) and ZA's second--order generalization. In the last two cases we also created new versions of the approximation by truncation, i.e., by finding an optimum smoothing window (see text) for the initial conditions. The Zel'dovich approximation, with optimized initial smoothing, worked extremely well. Its second-order generalization was slightly better. The success of our best-choice was a result of the treatment of the phases of nonlinear Fourier components. The adhesion approximation produced the most accurate nonlinear power spectrum and density distribution, but its phase errors suggest mass condensations were moved somewhat incorrectly. Due to its better reproduction of the mass density distribution function and power spectrum, adhesion might be preferred for some uses. We recommend either n-body simulations or our modified versions of ZA, depending on the purpose. Modified ZA can rapidly generate large numbers of realizations of model universes with good accuracy down to galaxy group (or smaller) mass scales.


title: Substructure in Clusters and Central Galaxy Peculiar Velocities

abstract: Formation theories for central dominant galaxies in clusters require them to be located at the minimum of the cluster gravitational potential. However, 32\% (8 out of 25) of the clusters with more than 50 measured redshifts have central galaxies with significant velocity offsets (with respect to other cluster members). By studying their velocity distributions and correlations between velocity and position, I show that the presence of a large peculiar velocity is strongly correlated with the presence of substructure in these massive systems. About 85\% (21 of 25) of all well-studied clusters show some evidence for substructure, in contrast to the 30-40\% found when using only galaxy or gas distributions. The correlation between substructure and central galaxy location verifies the hypothesis of Merritt (1985) and Tremaine (1990) that high peculiar velocities are indicative of recent merger events between less-massive systems of galaxies. Dynamical friction should act quickly to pull the central galaxy, the most massive discrete object in a cluster, to the minimum of the potential. The less-massive galaxies retain information about their primordial subclusters for a longer period of time. I use an objective partitioning algorithm to assign cluster galaxies to their host subclumps. When galaxies are allocated in this fashion to their subclusters, 75\% of the significant velocity offsets are eliminated. Only 2 out of the 25 clusters have central galaxies which are not centrally-located when substructure is considered in the analysis.


title: Non-Gaussian statistics of pencil beam surveys

abstract: We study the effect of the non-Gaussian clustering of galaxies on the statistics of pencil beam surveys. We find that the higher order moments of the galaxy distribution play an important role in the probability distribution for the power spectrum peaks. Taking into account the observed values for the kurtosis of galaxy distribution we derive the general probability distribution for the power spectrum modes in non-Gaussian models and show that the probability to obtain the 128$\hm$ periodicity found in pencil beam surveys is raised by roughly one order of magnitude. The non-Gaussianity of the galaxy distribution is however still insufficient to explain the reported peak-to-noise ratio of the periodicity, so that extra power on large scales seems required.


title: Microwave Background Fluctuations in an Open Universe

abstract: ( to appear in: Proceedings of the Nishonomiya Yukawa Memorial Symposium Edited by M. Sasaki) There are several models for generating fluctuations in an open universe that are compatible with the microwave background fluctuations detected by COBE {\it and} observations of large scale structure. Topological defects, such as strings and textures, appear to be more successful in a low-$\Omega$ model than in an $\Omega = 1$ model. These models predict that the cosmic microwave background (CMB) fluctuations are non-Gaussian with distinctive signatures that may lead to their confirmation. If there exists a mechanism for generating scale invariant mass fluctuations, $(\delta M/M)^2 \propto M^{-4/3}$, then these models are also compatible with COBE. Both models predict mass fluctuations on the 8/h Mpc scale, $\sigma_8$, $\sim 0.5 - 1$ for $\Omega h \sim 0.2 - 0.3$. A variety of observations of large-scale structure, clusters and galaxies, as well as the age problem, suggest that the low-$\Omega$ models are attractive.


title: Short time scale monitoring of SiO sources

abstract: We present the results of a short time scale monitoring of SiO maser emission (v=1 J=1-0 transition) in four known strong sources. These sources were monitored nightly for a period of about a month. The aim of these observations is to investigate the possible presence of variations in the maser lines on time scales of a few days to weeks, due to sound waves propagating out from the central star. If sound waves are responsible for the mass loss of certain cool giants, as suggested by Pijpers and Hearn (1989) and Pijpers and Habing (1989), local variations in density and relative velocity are expected just above the stellar photosphere. These could give rise to variations in any narrow spectral line formed in this region, and therefore in particular in the SiO maser lines. Our observations indicate that variations in the line shape (leading to relative changes in the intensity of about 20\%) occur in the SiO emission of Mira type stars, within short time scales of 10-20 days. The main component of the profile variability is consistent with a displacement of the velocity centroid of the dominant maser peaks, by about 1 km s$^{-1}$ in the average. Apparent variations in the total line flux were also found, but could be partially due to calibration uncertainties.


title: Smoothed Particle Interpolation

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


title: Evolution of Galaxies in Clusters

abstract: We use new models of stellar population synthesis to estimate the fraction of stars formed during the last major bursts of star formation in E/S0 galaxies in low-redshift clusters ($z\simlt0.4$) from the spectral signatures of intermediate-age stars. We find that the mass fraction of stars formed in late bursts in early-type galaxies in clusters must have decreased smoothly with redshift, from about 8\% at $z\approx1$ to less than 1\% at $z\approx0$. This result, which we interpret as a constraint on stellar mass added in mergers, is nearly independent of the assumed ages and morphological types of the progenitor galaxies prior to the last major bursts of star formation. We then compute the implied color and 4000 {\AA} break evolution for progenitors of E/S0 galaxies in clusters at redshifts $0\simlt z\simlt1$. We investigate a conservative model in which all present-day E/S0 galaxies are assumed to initially be elliptical galaxies and to then undergo bursts of star formation at the rate estimated from the signatures of intermediate-age stars at low redshifts. This model reproduces well the observed spread of colors and 4000 {\AA} breaks of galaxies in high-redshift clusters, but underestimates the fraction of galaxies in the blue tail of the distribution. Such a discrepancy may be interpreted as an increasing fraction of spiral galaxies in clusters at high redshift, as suggested independently by recent {\it HST} observations of intermediate-redshift clusters. The current data do not seem to require morphological evolution of cluster galaxies out to $z\sim0.4$, but suggest that either morphological or strong luminosity evolution might have played a major role at $z\simgt0.7$.


title: Fabry-Perot Measurements of the Dynamics of Globular Cluster Cores: M15 (NGC~7078)

abstract: We report the first use of the Rutgers Imaging Fabry-Perot Spectrophotometer to study the dynamics of the cores of globular clusters. We have obtained velocities for cluster stars by tuning the Fabry-Perot to take a series of narrow-band images at different wavelengths across one of the Na D (5890 AA) absorption lines. Measuring the flux in every frame yields a short portion of the spectrum for each star simultaneously. This proves to be a very efficient method for obtaining accurate stellar velocities; in crowded regions we are able to measure hundreds of velocities in 3-4 hours of observing time. We have measured velocities with uncertainties of less than 5 km/s for 216 stars within 1.5' of the center of the globular cluster M15 (NGC 7078). The paper is a uuencoded compressed postscript file.


title: Non-Linear Approximations to Gravitational Instability: A Comparison in the Quasi-Linear Regime

abstract: We compare different nonlinear approximations to gravitational clustering in the weakly nonlinear regime, using as a comparative statistic the evolution of non-Gaussianity which can be characterised by a set of numbers $S_p$ describing connected moments of the density field at the lowest order in $<\delta^2>$: $<\delta^n>_c \simeq S_n<\delta^2>^{n-1}$. Generalizing earlier work by Bernardeau (1992) we develop an ansatz to evaluate all $S_p$ in a given approximation by means of a generating function which can be shown to satisfy the equations of motion of a homogeneous spherical density enhancement in that approximation. On the basis of the values of we show that approximations formulated in Lagrangian space (such as the Zeldovich approximation and its extensions) are considerably more accurate than those formulated in Eulerian space such as the Frozen Flow and Linear Potential approximations. In particular we find that the $n$th order Lagrangian perturbation approximation correctly reproduces the first $n+1$ parameters $S_n$. We also evaluate the density probability distribution function for the different approximations in the quasi-linear regime and compare our results with an exact analytic treatment in the case of the Zeldovich approximation.


title: Collapse of Axionic Domain Wall and Axion Emission

abstract: We examine the collapse of an axion domain wall bounded by an axionic string. It is found that the collapse proceeds quickly and axion domain walls disappear. However axions are emitted in the collapse and its energy density increases during radiation dominated era and contributes significantly to the present mass density of the universe. In particular the axion emitted from the wall can account for the dark matter in the universe for $F_a\gsim 10^{10}\GeV$.


title: Lick Optical Spectra of Quasar HS 1946+7658 at 10 km/sec Resolution: Lyman-Alpha Forest and Metal Absorption Systems

abstract: We present optical spectra of QSO HS~1946+7658 with either high resolution (FWHM=10 km/s) or high signal to noise ratio (SNR=40-100). We find 113 Lyman alpha and six metal line systems. The metal systems at Zabs=2.844 and 3.050 have complex velocity structures. We find that the system at 2.844 is a damped Ly-a absorption system, with neutral hydrogen column density of logN(HI)=20.2+/-0.4, and it is the cause of the Lyman limit break at 3520\AA. The metal abundance in the gas phase of the system is [M/H]=-2.6+/-0.3, with a best estimate of [M/H]= -2.8, with ionization parameter Gamma=-2.75, from a photoionization model. We show that the abundance ratios rule out appreciable dust in this DLA system. We do not see the enhancement of O over C reported in Nature by Reimers last years. We see CII*(1335) offset by 15 km/s with respect to CII(1334), presumably because the gas density varies from 2 - 8 cm(-3) with changing velocity in the DLA system. These densities imply that the damped component is 6 - 25 pc thick, which is reasonable for a single cloud in a cold spiral disk. We see a 2.6sigma lack of Ly-a forest lines well away from the QSO redshift, which may be a chance fluctuation. We also see a correlation between column density N(HI) and Doppler parameter b for 96 unsaturated Ly-a forest absorption lines, and although this correlation persists in the 36 Ly-a lines which lie in regions where SNR=8-16, we agree with Rauch etal (1993) that it is probably a bogus effect of low SNR. We show that there are few low-b Ly-a forest lines in moderate SNR data.


title: Ohm's Law for a Relativistic Pair Plasma

abstract: We derive the fully relativistic Ohm's law for an electron-positron plasma. The absence of non-resistive terms in Ohm's law and the natural substitution of the 4-velocity for the velocity flux in the relativistic bulk plasma equations do not require the field gradient length scale to be much larger than the lepton inertial lengths, or the existence of a frame in which the distribution functions are isotropic.


title: Merger rates in hierarchical models of galaxy formation. II: Comparison with N-body simulations

abstract: We have made a detailed comparison of the results of large N-body simulations with the analytical description of the merging histories of dark matter halos presented in Lacey & Cole 1993, which is based on an extension of the Press- Schechter method (Bond etal 1991,Bower 1991). We find the analytical predictions for the halo mass function, merger rates and formation times to be remarkably accurate. The analytical model is however expected to apply for arbitrary Omega and general power spectra. Dark matter halos were identified in the simulations using two different methods. For halos selected at mean overdensities 100-200, the analytical mass function was found to provide a good fit to the simulations with a collapse threshold close to that predicted by the spherical collapse model, with a typical error of <30% over a range of 1000 in mass, which is the full dynamic range of our N-body simulations. This was insensitive to the type of filtering used. Over a range of 100 - 1000 in mass, there was also good agreement with the analytical predictions for merger rates, including their dependence on the masses of the two halos involved and the time interval being considered, and for formation times, including the dependence on halo mass and formation epoch. The analytical Press-Schechter mass function and its extension to halo lifetimes and merger rates thus provide a very useful description of the growth of dark matter halos through hierarchical clustering, and should providea valuable tool in studies of the formation and evolution of galaxies and galaxy clusters.


title: Simulating Self-gravitating Hydrodynamic Flows

abstract: An efficient algorithm for solving Poisson's equation in two and three spatial dimensions is discussed. The algorithm, which is described in detail, is based on the integral form of Poisson's equation and utilizes spherical coordinates and an expansion into spherical harmonics. The solver can be applied to and works well for all problems for which the use of spherical coordinates is appropriate. We also briefly discuss the implementation of the algorithm into hydrodynamic codes which are based on a conservative finite--difference scheme.


title: A ``Skewed'' lognormal approximation to the probability distribution function of the large-scale density field

abstract: I propose a method to fit the probability distribution function (hereafter PDF) of the large scale density field rho, motivated by a Lagrangian version of the continuity equation. It consists in applying the Edgeworth expansion to the quantity Phi=log rho - < log rho >. The method is tested on the matter particle distribution in two cold dark matter N-body simulations of different physical sizes to cover a large dynamic range. It is seen to be very efficient, even in the non-linear regime, and may thus be used as an analytical tool to study the effect on the PDF of the transition between the weakly non-linear regime and the highly non-linear regime.


title: Relating Spectral Indices to Tensor and Scalar Amplitudes in Inflation

abstract: Within an expansion in slow-roll inflation parameters, we derive the complete second-order expressions relating the ratio of tensor to scalar density perturbations and the spectral index of the scalar spectrum. We find that ``corrections'' to previously derived formulae can dominate if the tensor to scalar ratio is small. For instance, if $V V''/(V')^2\neq 1$ or if $m_{Pl}^2/(4\pi) ~|V'''/V'|\ga 1$, where $V(\phi)$ is the inflaton potential and $m_{Pl}$ is the Planck mass, then the previously used simple relations between the indices and the tensor to scalar ratio fails. This failure occurs in particular for natural inflation, Coleman--Weinberg inflation, and ``chaotic'' inflation.


title: Boron in the Extreme Pop II Star HD 140283 and the Production of Light Elements in the Early Galaxy

abstract: Using observations of the 2496.7 A B I line with the HST GHRS at a nominal resolution of 90,000, we have found the abundance of boron of HD 140283 to be log epsilon_B (= 12 + log (N_B/N_H)) = 0.34 \pm 0.20. This result is found when a significant non-LTE effect in the formation of the B I line is taken into account. The resulting N_B/N_Be ratio is about 17 (in the range 9 -- 34), which is in very good agreement with what is expected from spallation by cosmic rays. We conclude that this origin of Be and B in the Early Galaxy is the most probable of recently suggested formation mechanisms. (Submitted to A&A. Postscript figures available via annymous ftp to "chiron.astro.uu.se" in directory pub/Preprints/P87)


title: X-Ray and Optical Properties of Groups of Galaxies

abstract: We examine the relationship between the group x-ray luminosity in the 0.3-3.5 keV band and the measured velocity dispersion, and galaxy surface number density. We find definite correlations. Richer groups follow the same relation as rich clusters (cf. Quintana & Melnick 1982) with $L_x \propto \sigma^{4.0\pm 0.6}$, but the relation flattens for lower luminosity systems which have velocity dispersions below 300 km/s). We suggest that the $L_x-\sigma$ relation arises from a combination of extended cluster emission and emission associated with individual galaxies. The x-ray emission for the richer groups is dominated by emission from the intra-group medium, as for the richer clusters; emission from the poorer clusters is dominated by less extended emission associated with the individual group galaxies. This is a corrected version of the paper which appeared in the Astronomical Journal (February 1994). We include a new version of the affected figure.


title: Testing Inflation with the Cosmic Microwave Background

abstract: Cosmic microwave background (CMB) anisotropy may result from both scalar and tensor perturbations. For a sufficiently narrow range of angular scales, CMB perturbations can be characterized by four parameters. Results from the Cosmic Background Explorer fix one combination of the parameters, reducing the parameters to three. If CMB perturbations are from inflation, there is an additional relation, reducing the parameters to two. An appropriate combination of a medium-angle and a small-angle CMB observation can test the inflation hypothesis because inflation cannot explain a high signal in one experiment and a low signal in the other.


title: Texture-Induced Microwave Background Anisotropies

abstract: We use numerical simulations to calculate the cosmic microwave background anisotropy induced by the evolution of a global texture field, with special emphasis on individual textures. Both spherically symmetric and general configurations are analysed, and in the latter case we consider field configurations which exhibit unwinding events and also ones which do not. We compare the results given by evolving the field numerically under both the expanded core (XCORE) and non-linear sigma model (NLSM) approximations with the analytic predictions of the NLSM exact solution for a spherically symmetric self-similar (SSSS) unwinding. We find that the random unwinding configuration spots' typical peak height is 60--75\% and angular size typically only 10\% of those of the SSSS unwinding, and that random configurations without an unwinding event nonetheless may generate indistinguishable hot and cold spots. The influence of these results on analytic estimates of texture induced microwave anisotropies is examined, and comparison made with other numerical work.


title: Numerical Methods in Cosmological Global Texture Simulations

abstract: Numerical simulations of the evolution of a global topological defect field have two characteristic length scales --- one macrophysical, of order the field correlation length, and the other microphysical, of order the field width. The situation currently of most interest to particle cosmologists involves the behaviour of a GUT-scale defect field at the epoch of decoupling, where the ratio of these scales is typically of order $10^{50}$. Such a ratio is unrealisable in numerical work, and we consider the approximations which may be employed to deal with this. Focusing on the case of global texture we outline the implementation of the associated algorithms, and in particular note the subtleties involved in handling texture unwinding events. Comparing the results in each approach then establishes that, subject to certain constraints on the minimum grid resolution, the methods described are both robust and consistent with one another.


title: Particle Physics and Cosmology

abstract: Lectures presented at the 42nd Scottish Universities Summer School in Physics, St. Andrews, Scotland, August 1993.


title: Diffuse Molecular Clouds and the Molecular Interstellar Medium from 13CO Observations of M33

abstract: We have obtained 12CO and 13CO J=1-0 observations of the nearby spiral galaxy M33 to try to resolve the long-standing discrepancy between 12CO/13CO line ratios measured in Galactic giant molecular clouds and external galaxies. Interferometer maps of the molecular cloud MC20 give a 12CO/13CO line ratio of 7.5+/-2.1, which agrees reasonably well with the line ratio measured in Galactic giant molecular clouds. In contrast, the 12CO/13CO line ratio obtained from single dish data is 10.0+/-0.9, significantly higher than Galactic values but in good agreement with line ratios measured in other galaxies. The interferometer map of MC20 reveals that the cloud has similar spatial and velocity extents in the two lines, and thus the high single dish line ratio cannot be due to different filling factors in the two lines. In addition, the single dish data show no evidence for significant variations in the line ratio with metallicity, which eliminates abundance changes as the explanation for the high single dish line ratio. We conclude that the high 12CO/13CO line ratio s observed in M33, and in the disks of spiral galaxies in general, are due to the presence of a population of diffuse molecular clouds. The lower limit to the fraction of the total 12CO emission from M33 that originates in the diffuse clouds is 30+/-30%, while the upper limit is ~60%.


title: On the Formation of Elliptical Galaxies

abstract: It is shown that the violent relaxation of dissipationless stellar systems leads to universal de Vaucouleurs profiles only outside 1.5 effective radii $R_e$. Inside $1.5 R_e$ the surface density profiles depend strongly on the initial conditions and are in general not in agreement with the de Vaucouleurs law. This result is in contradiction to the observations which show that all elliptical galaxies have de Vaucouleurs profiles inside $1.5 R_e$ and show strong deviations outside $1.5 R_e$. The $r^{1/4}$-profiles of E-galaxies are therefore not a direct result of collisionless violent relaxation and angular momentum redistribution. It is rather a universal initial density profile which is required in order to explain the universal structure of E-galaxies. A new model for the formation of E-galaxies is presented which assumes that ellipticals and bulges formed from isothermal, gaseous spheres which accumulated in the core regions of dark matter halos and condensed into stars through a cooling instability. This scenario predicts that early type galaxies trace directly the core radii and core masses of their dark matter halos, providing important information on the nature of the dark matter component in galaxies.


title: Is there any Scaling in the Cluster Distribution ?

abstract: We apply fractal analysis methods to investigate the scaling properties in the Abell and ACO catalogs of rich galaxy clusters. We also discuss different technical aspects of the method when applied to data sets with small number of points as the cluster catalogs. Results are compared with simulations based on the Zel'dovich approximation. We limit our analysis to scales less than 100 $\hm$. The cluster distribution show a scale invariant multifractal behavior in a limited scale range. For the Abell catalog this range is 15--60$\hm$, while for the ACO sample it extends to smaller scales. Despite this difference in the extension of the scale--range where scale--invariant clustering takes place, both samples are characterized by remarkably similar multifractal spectra in the corresponding scaling regime. In particular, the correlation dimension turns out to be $D_2\simeq 2.2$ for both Abell and ACO clusters. Although it is difficult to point out the scale at which homogeneity is reached with the present size of these redshift surveys, our results indicate that the cluster distribution shows a tendency to homogeneity at large scales, disproving the picture of a pure scale invariant fractal structure extending to arbitrarily large distances.


title: Large-amplitude isothermal fluctuations and high-density dark-matter clumps

abstract: Large-amplitude isothermal fluctuations in the dark matter energy density, parameterized by $\Phi\equiv\delta\rhodm/\rhodm$, are studied within the framework of a spherical collapse model. For $\Phi \ga 1$, a fluctuation collapses in the radiation-dominated epoch and produces a dense dark-matter object. The final density of the virialized object is found to be $\rho_F \approx 140\, \Phi^3 (\Phi+1) \rhoeq$, where $\rhoeq$ is the matter density at equal matter and radiation energy density. This expression is valid for the entire range of possible values of $\Phi$, both for $\Phi \gg 1$ and $\Phi \ll 1$. Some astrophysical consequences of high-density dark-matter clumps are discussed.


title: Abundance of moderate-redshift clusters in the Cold + Hot dark matter model

abstract: Using a set of $\pppm$ simulation which accurately treats the density evolution of two components of dark matter, we study the evolution of clusters in the Cold + Hot dark matter (CHDM) model. The mass function, the velocity dispersion function and the temperature function of clusters are calculated for four different epochs of $z\le 0.5$. We also use the simulation data to test the Press-Schechter expression of the halo abundance as a function of the velocity dispersion $\sigma_v$. The model predictions are in good agreement with the observational data of local cluster abundances ($z=0$). We also tentatively compare the model with the Gunn and his collaborators' observation of rich clusters at $z\approx 0.8$ and with the x-ray luminous clusters at $z\approx 0.5$ of the {\it Einstein} Extended Medium Sensitivity Survey. The important feature of the model is the rapid formation of clusters in the near past: the abundances of clusters of $\sigma_v\ge 700\kms$ and of $\sigma_v\ge 1200 \kms$ at $z=0.5$ are only 1/4 and 1/10 respectively of the present values ($z=0$). Ongoing ROSAT and AXAF surveys of distant clusters will provide sensitive tests to the model. The abundance of clusters at $z\approx 0.5$ would also be a good discriminator between the CHDM model and a low-density flat CDM model both of which show very similar clustering properties at $z=0$.


title: Magnification Bias in Galactic Microlensing Searches

abstract: It is shown that a significant amount of detectable gravitational microlensing events that could potentially be found by MAssively Parallel Photometry (MAPP) projects (such as the MACHO, EROS and OGLE collaborations) will occur for stars too dim to be easily noticed individually by these projects. This is the result of a large magnification bias effect, a bias of including high magnification events in any flux-limited sample. The probability of detecting these events may be as high as 2.3 times the lensing probability of stars currently being monitored by MAPP collaborations.


title: The Near-IR-Optical-UV Emission of BL Lacertae Objects

abstract: Near--infrared, optical and ultraviolet quasi--simultaneous observations of 11 BL Lacertae objects are reported. For all but one source the dereddened spectral flux distribution in the $8\cdot10^{13}-2\cdot10^{15}$ Hz frequency range can be described by a single power law f$_\nu \propto \nu^{-\alpha}$ with average spectral index $<\alpha>$ = 0.88 $\pm$ 0.42 (standard deviation) plus, where relevant, the contribution of the host galaxy. In most cases the non simultaneous soft X--ray fluxes obtained by the {\it Einstein Observatory} lie on or below the extrapolation of the power law. The results are compared with the average spectral properties of other samples of BL Lacs studied separately in the IR--optical and in the UV bands. The implications for existing models of the objects are shortly discussed.


title: Amplitudes of Stellar Oscillations: The Implications for Asteroseismology

abstract: There are no good predictions for the amplitudes expected from solar-like oscillations in other stars. In the absence of a definitive model for convection, which is thought to be the mechanism that excites these oscillations, the amplitudes for both velocity and luminosity measurements must be estimated by scaling from the Sun. In the case of luminosity measurements, even this is difficult because of disagreement over the solar amplitude. This last point has lead us to investigate whether the luminosity amplitude of oscillations (dL/L) can be derived from the velocity amplitude v_osc. Using linear theory and observational data, we show that p-mode oscillations in a large sample of pulsating stars satisfy (dL/L)_bol proportional to v_osc/T_eff. Using this relationship, together with the best estimate of v_osc(Sun) = (23.4 +/- 1.4) cm/s, we estimate the luminosity amplitude of solar oscillations at 550 nm to be dL/L = (4.7 +/- 0.3) ppm. Next we discuss how to scale the amplitude of solar-like (i.e., convectively-powered) oscillations from the Sun to other stars. The only predictions come from model calculations by (Christensen-Dalsgaard & Frandsen, Sol. Phys. 82, 469). However, their grid of stellar models is not dense enough to allow amplitude predictions for an arbitrary star. Nevertheless, although convective theory is complicated, we might expect that the general properties of convection -- including oscillation amplitudes -- should change smoothly through the colour-magnitude diagram. Indeed, we find that the velocity amplitudes predicted by the model


title: Lagrangian Evolution of the Weyl Tensor

abstract: We derive the evolution equations for the electric and magnetic parts of the Weyl tensor for cold dust from both general relativity and Newtonian gravity. In a locally inertial frame at rest in the fluid frame, the Newtonian equations agree with those of general relativity. We give explicit expressions for the electric and magnetic parts of the Weyl tensor in the Newtonian limit. In general, the magnetic part does not vanish, implying that the Lagrangian evolution of the fluid is not purely local.


title: On the spin parameter of dark-matter haloes

abstract: The study by White (1984) on the growth of angular momentum in dark haloes is extended towards a more detailed investigation of the spin parameter $\lambda\equiv L\sqrt{E}/{G M^{2.5}}$. Starting from the Zel'dovich approximation to structure formation, a dark halo is approximated by a homogeneous ellipsoid with the inertial tensor of the (highly irregular) Lagrangian region $\Upsilon$ from which the dark halo forms. Within this approximation, an expression for the spin parameter can be derived, which depends on the geometry of $\Upsilon$, the cosmological density parameter $\Omega_0$, the overdensity of the dark halo, and the tidal torque exerted on it. For Gaussian random fields, this expression can be evaluated statistically. As a result, we derive a probability distribution of the spin parameter which gives $\lambda\simeq0.07^{+0.04}_{-0.05}$, consistent with numerical investigations. This probability distribution steeply rises with increasing spin parameter, reaching its maximum at $\lambda\simeq0.025$. The 10 (50,90) percentile values are $\lambda=0.02$ (0.05,0.11, respectively). There is a weak anticorrelation of the spin parameter with the peak height $\nu$ of the density fluctuation field $\lambda\propto \nu^{-0.29}$. The dependence on $\Omega_0$ and the variance $\sigma$ of the density-contrast field is very weak; there is only a marginal tendency for the spin parameter to be slightly larger for late-forming objects in an open universe. Due to the weak dependence on $\sigma$, our results should be quite generally applicable and independent on


title: Axion String Constraints

abstract: We study the production of cosmological axions in the standard scenario in which a global string network forms at the Peccei-Quinn phase transition. We make detailed calculations of the axions produced by string loops (a previously overlooked source), comparing these with estimates of other contributions from long strings and domain walls. We dilineate key uncertainities in the string network evolution, the chief of which is the magnitude of the radiative backreaction. We conclude that axions produces by the these topological contributions provide the significant cosmological constraint on the symmetry breaking scale $\fa$ and the mass $\ma$


title: The Parker instability in disks with differential rotation

abstract: We present a detailed study of the growth of the Parker instability in a differentially rotating disk embedded in an azimuthal equilibrium magnetic field, such as the interstellar gas or an accretion disk. Basic properties of the instability without shear are first recalled. Differential rotation is modeled in the shearing sheet approximation, classical in the theory of spiral density waves. The action of differential rotation is reduced to two different effects, (i) a linear time-dependence of the radial wavenumber, and (ii) a radial differential force. We present both exact numerical solutions, and approximate analytical ones based on the WKB approximation in the limit of weak differential rotation. Most important are (i) a transient natural stabilization of the Parker mode due to the radial differential force (ii) the generation of magnetosonic and Alfvenic waves, and (iii) in a certain parameter range a possible ``turn-over'' of the perturbation whereby, quite surprisingly, matter which had started being elevated by the instability may end up dropping towards the disk midplane. A simplified model shows the possible observable effects of this turn-over.


title: The Effects of Smoothing on the Statistical Properties of the Large-Scale Cosmic Fileds

abstract: It has been shown that the large--scale correlation functions of the density field (and velocity divergence field) follow a specific hierarchy in the quasilinear regime and for Gaussian initial conditions (Bernardeau 1992). The exact relationships between the cumulants of the probability distribution functions (the so-called $S_p$ parameters) are however sensitive to the smoothing window function applied to the fields. In this paper, I present a method to derive the whole series of the $S_p$ parameters when the density field is smoothed with a top--hat window function. The results are valid for any power spectrum and any cosmological parameters. Similar calculations are presented for the velocity divergence field. The resulting shapes of the one--point probability distribution functions of the cosmic density and the velocity divergence fields are given as a function of the power spectrum and $\Omega$. Simple analytical fits are proposed when the index of the power spectrum is -1. Comparisons with numerical simulations prove these analytical results to be extremely accurate.


title: Correlated Errors in the COBE DMR Sky Maps

abstract: The {\it COBE} DMR sky maps contain low-level correlated noise. We obtain estimates of the amplitude and pattern of the correlated noise from three techniques: angular averages of the covariance matrix, Monte Carlo simulations of two-point correlation functions, and direct analysis of the DMR maps. The results from the three methods are mutually consistent. The noise covariance matrix of a DMR sky map is diagonal to an accuracy of better than 1\%. For a given sky pixel, the dominant noise covariance occurs with the ring of pixels at an angular separation of $60 \deg$ due to the $60 \deg$ separation of the DMR horns. The mean covariance at $60 \deg$ is $0.45\% ^{+0.18}_{-0.14}$ of the mean variance. Additionally, the variance in a given pixel is $0.7\%$ greater than would be expected from a single beam experiment with the same noise properties. Auto-correlation functions suffer from a $\sim 1.5\; \sigma$ positive bias at $60 \deg$ while cross-correlations have no bias. Published {\it COBE} DMR results are not significantly affected by correlated noise. COBE pre-print 94-


title: Microlensing by Stars in the Disk of M31

abstract: The optical depth to microlensing toward M31 due to known stars in the disk of M31 itself is $\tau\sim 2\times 10^{-7}e^{-r/d}$ where $d$ is the disk scale length and $r$ is the distance along the major axis. Thus, there can be significant lensing toward the M31 disk even if M31 contains no dark compact objects. The optical depth has a strong dependence on azimuthal angle: at fixed radius $\tau\propto[1+ (h/d)\tan i\cos\phi]^{-2}$ where $h$ is the scale height of the disk, $i=75^\circ$ is the inclination of M31, and $\phi$ is the azimuthal angle relative to the near minor axis. By measuring the optical depth as a function of radial and azimuthal position, it is possible to estimate $h$ and $d$ for the mass of the M31 disk, and so determine whether the disk light traces disk mass. Ground-based observations in $0.\hskip-2pt''5$ seeing of 0.8 $\rm deg^2$ once per week could yield $\sim 3$ events per year. With an ambitious space-based project, it would be possible to to observe $\sim 80$ events per year. If lensing events were dominated by the M31 spheroid rather than the disk, the event rate would be higher and the spheroid's parameters could be measured. The pattern of optical depths differs substantially for a disk and a spheroid.


title: Higher--Order Lagrangian Perturbation Theory

abstract: Fundamental assumptions which form the basis of models for large-scale structure in the Universe are sketched in light of a Lagrangian description of inhomogeneities. This description is introduced for Newtonian self-gravitating flows. On its basis a Lagrangian perturbation approach is discussed and compared with the standard Eulerian theory of gravitational instability. The performance of Lagrangian perturbation solutions up to the third order is demonstrated in comparison with numerical N-body simulations. First results of this comparison are presented for large scales (PM-code) and for small scales (tree-code).


title: Statistics of N-Body Simulations. II. Equal Masses after Core Collapse

abstract: This paper presents and analyses statistical results from a large number of $N$-body simulations of isolated systems with equal masses, in which $250\le N\le 2000$. It concentrates on the phase starting around the end of core collapse, when binaries start to play a crucial role. Interactions of hard binaries are consistent with the Spitzer (1987) cross section. The evolution of the half-mass radius after core collapse nearly follows classical theory, The evolution of the inner parts of the system around the time of core bounce is consistent with these simplified models provided that the continuous production of energy, as usually assumed, is replaced by a model of {\sl stochastic} energy production. Similarly, post-collapse evolution of the core requires a modest recalibration of the coefficient of energy generation, especially for small $N$. The distributions of velocity dispersion and anisotropy become remarkably homologous soon after core bounce. The bound mass of the systems very nearly follows a power law with time. A small number of escapers, presumed to be those associated with binary activity, dominate the energy which is carried off. The ``internal" energy of escaping binaries is consistent with theoretical expectations.


title: Polarization Profiles of Scattered Emission Lines. I. General Formalism for Optically Thin Rayleigh Scattering

abstract: A general theoretical framework is developed for interpreting spectropolarimetric observations of optically thin emission line scattering from small dust particles. Spatially integrated and spatially resolved line profiles of both scattered intensity and polarization are calculated analytically from a variety of simple kinematic models. These calculations will provide a foundation for further studies of emission line scattering from dust and electrons in such diverse astrophysical environments as Herbig-Haro objects, symbiotic stars, starburst galaxies and active galactic nuclei.


title: Polarization Profiles of Scattered Emission Lines. II. Upstream Dust Scattering in the HH 1 Jet

abstract: Detailed comparisons are made between observations of scattered light upstream of the head of the HH~1 jet and predictions of simple scattering models. It is shown that, in order to unambiguously determine the velocity of the head of the jet (bow shock) with respect to the upstream dust, existing spectroscopic observations are insufficient and that spectropolarimetric observations of the scattered light are necessary. Such an independent measure of the bow shock velocity is important in order to test ``multiple outflow'' theories of Herbig-Haro jets. It is also shown 2that the scattering dust must have a very forward-throwing scattering phase function ($\langle\cos\theta\rangle\msim 0.7$) and slight evidence is found for a dust-gas ratio that is higher than average.


title: Doppler Ghosts in a Hall of Mirrors: Polarisation Profiles of Scattered Emission Lines

abstract: The Doppler shifts of optical emission lines which have been scattered by surrounding dust and electrons can provide useful information about the kinematics, geometry and physical conditions of astrophysical flows. In principle, the scatterers can provide views of the line-emitting gas from different directions, allowing the 3-d velocity field of the emitting gas to be determined and revealing sources which are hidden from direct view. Unfortunately, as in a Hall of Mirrors, the resultant multiple images can be confusing and hard to interpret. In general, the scattering geometry will be unknown, which makes it difficult to disentangle the effect of the motion of the scatterers from that of the motion of the emitting source. In this situation, spectropolarimetric observations can be a great help, since the scattered light will be partially polarised to a degree dependent on the angle of scattering and on the details of the scattering process. Analysis of the polarisation profiles of the scattered emission lines can then determine the relative orientation and velocity of the source and the scatterers. Such techniques are applied to the upstream scattered emission lines in the HH 1 jet, where it is shown that spectropolarimetric measurements would allow an independent determination of the bowshock velocity with respect to the upstream matter, which would test ``multiple outburst'' models of the jet.


title: Properties of the Cosmological Density Distribution Function

abstract: The properties of the probability distribution function of the cosmological continuous density field are studied. We present further developments and compare dynamically motivated methods to derive the PDF. One of them is based on the Zel'dovich approximation (ZA). We extend this method for arbitrary initial conditions, regardless of whether they are Gaussian or not. The other approach is based on perturbation theory with Gaussian initial fluctuations. We include the smoothing effects in the PDFs. We examine the relationships between the shapes of the PDFs and the moments. It is found that formally there are no moments in the ZA, but a way to resolve this issue is proposed, based on the regularization of integrals. A closed form for the generating function of the moments in the ZA is also presented, including the smoothing effects. We suggest the methods to build PDFs out of the whole series of the moments, or out of a limited number of moments -- the Edgeworth expansion. The last approach gives us an alternative method to evaluate the skewness and kurtosis by measuring the PDF around its peak. We note a general connection between the generating function of moments for small r.m.s $\sigma$ and the non-linear evolution of the overdense spherical fluctuation in the dynamical models. All these approaches have been applied in 1D case where the ZA is exact, and simple analytical results are obtained. The 3D case is analyzed in the same manner and we found a mutual agreement in the PDFs derived by different methods in the the quasi-linear regime. Numerical CDM simulation was used to validate the accuracy of considered approximations. We explain the successful log-normal fit of the PDF from that simulation at moderate $\sigma$ as mere fortune, but not as a universal form of density PDF in general.


title: Equations of Gravitational Instability are Non-Local

abstract: Few recent generations of cosmologists have solved non-local newtonian equations of the gravitational instability in an expanding universe. In this approach pancaking is the predominant form of first collapsing objects. Relativistic counterparts of these equations contain the electric and magnetic parts of the Weyl tensor. In the linear theory the magnetic part is associated with gravitational waves. If the magnetic part is ignored, then the newtonian limit of the relativistic equations is reduced to the closed set of the local Lagrangian equations. Recently this fact drew much attention since the gravitational instability in that form would greatly simplify the study of cosmic structure formation. In particular, the filamentary structure of collapsing is predicted. In this paper we resolve the contradiction between the newtonian theory and relativistic version adopted in some recent papers. We show that dropping the magnetic part from the basic relativistic equations is {\it incorrect}. The correct newtonian limit is derived by the $1/c$-expansion of the GR equations and the Bianchi identities for the Weyl tensor. The last ones begin with $\sim 1/c^3$ order, therefore one {\it must} take into account the magnetic part in the post newtonian order $\sim 1/c^3$, which contains non-local terms, related to the non-local gravitational interaction. For the first time we rigorously show that the GR equations with the magnetic part are reduced precisely to the canonic newtonian non-local equations. Thus, the correct treatment of the relativistic version of the gravitational instability resurrects the canonic picture of the structure formation.


title: High-Resolution Simulation of Deep Pencil Beam Surveys

abstract: We report on a recent paper (Wei\ss\ \& Buchert 1993), where we carry out pencil beam constructions in a high-resolution simulation of the large-scale structure of galaxies. As an example we present the results for the case of ``Hot-Dark-Matter" (HDM) initial conditions (with scale-free $n = 1$ power index on large scales and $\Omega = 1$) as a representative of models with sufficient large-scale power. We use an analytic approximation for particle trajectories of a self-gravitating dust continuum and apply a local dynamical biasing of volume elements to identify luminous matter in the model. Using this method, we are able to resolve formally a simulation box of 1200h^-1 Mpc (e.g. for HDM initial conditions) down to the scale of galactic halos using $2160^3$ particles. Pencil beam probes are taken for a given epoch using the parameters of observed beams. In particular, our analysis concentrates on the detection of a quasi-periodicity in the beam probes. The simulation is designed for application to parameter studies which prepare future observational projects. We find that a large percentage of the beams shows quasi-periodicities with periods which cluster at a certain length scale. The periods found range between one and eight times the cutoff length in the initial fluctuation spectrum.


title: Thermal History Constraints on the Isocurvature Baryon Model

abstract: We present a thorough and detailed investigation of baryon isocurvature models including realistic thermal histories, \ie, late or partial reionization of the universe and compact object formation after standard recombination. Constraints on these models are imposed from spectral distortion in the cosmic microwave background (CMB), number fluctuations in galaxy counts $\sigma_8$, and recent CMB anisotropy observations. Since the status of degree scale anisotropies is as yet unclear, the lack of spectral distortions is the most serious constraint on the spectral index of initial fluctuations: $-1.2 \simlt n \simlt -0.5$ in partially ionized models and $-1.2 \simlt n \simlt -0.9$ in compact object dominated universes in order to avoid significant bias in galaxy formation. Full ionization of the universe before $z = 300$ is also forbidden by this constraint. Intermediate scale CMB fluctuations which are significantly larger than standard CDM is a {\it robust} prediction of these models. Many models are consequently in conflict with the recent degree scale detections by Python, ARGO, and MSAM, \eg, compact object dominated models with $\Omega_0 \simlt 0.2$. Moreover, {\it all} models will be ruled out if the low CMB fluctuations on degree scales as detected by the SP91 13 point scan are confirmed. On the other hand, most models fare well compared with the high Tenerife and MAX detections.


title: A Relativistic Approch to Gravitational Instability in the Expanding Universe: Second Order Lagrangian Solutions

abstract: A Lagrangian relativistic approach to the non--linear dynamics of cosmological perturbations of an irrotational collisionless fluid is considered. Solutions are given at second order in perturbation theory for the relevant fluid and geometric quantities and compared with the corresponding ones in the Newtonian approximation. Specifically, we compute the density, the volume expansion scalar, the shear, the ``electric" part, or tide, and the ``magnetic" part of the Weyl tensor. The evolution of the shear and the tide beyond the linear regime strongly depends on the ratio of the characteristic size of the perturbation to the cosmological horizon distance. For perturbations on sub--horizon scales the usual Newtonian approximation applies, at least at the considered perturbative order; on super--horizon scales, instead, a new picture emerges, which we call ``silent universe", as each fluid element evolves independently of the environment, being unable to exchange signals with the surrounding matter through either sound waves or gravitational radiation. For perturbations inside the Hubble radius particular attention is paid in singling out non--local effects during the non--linear evolution of fluid elements. These non--local effects are shown to be carried by a traceless and divergenceless tensor, contained in the magnetic part of the Weyl tensor, which is dynamically generated as soon as the system evolves away from the linear regime.