title: Detection of Lens Candidates for the Double QSO Q2345+007

abstract: Luminous objects associated with the 7.06\arcsec\ separation double quasar Q2345+007 have, until now, escaped detection. In this letter we present the results of the deepest known imaging of the region surrounding the quasar. The total exposure times were 47600 seconds in B$_j$ (101 frames) and 32160 seconds in R (82 frames). The frames came from CFHT, CTIO, and KPNO. We detect a B$_j$ = 25.0 mag galaxy (B$_j$ -- R = 0.5 mag) in close proximity to the fainter QSO image, possibly at $z = 1.49$ given by several absorption features in the QSO spectra. Furthermore, there is a 2.9$\sigma$ enhancement in the number density of faint galaxies ($24 \le$ B$_j$ $\le 28$, B$_j-$ R $ \le 1.5$) near the quasar and another 3.1$\sigma$ enhancement further away. These discoveries support the hypothesis that Q2345+007 is being lensed by one or more distant mass concentrations and may imply that compact ``clusters'' already exist at early epochs. We discuss several lens models for the system and the cosmological implications.


title: Gravitational Lensing in Clusters of Galaxies: New Clues Regarding the Dynamics of Intracluster Gas

abstract: Long arcs in clusters of galaxies, produced by gravitational lensing, can be used to estimate the mass interior to the arcs and therefore, constrain the cluster mass distribution. The radial density distribution of the intracluster gas (ICM) can be extracted from the X-ray surface brightness observations. If the gas temperature is also known, it is then possible to probe the dynamical state of the gas and test whether the ICM is in hydro- static equilibrium within the gravitational potential of the cluster as a result of thermal pressure support. We analyze three clusters that exhibit large arcs, whose X-ray surface brightness profiles have been observed, and whose gas temperatures have been determined. In two of the clusters, A2218 and A1689, the central mass implied by lensing is a factor of $2$--$2.5$ too large for the gas at the observed temperature to be in hydrostatic equilibrium solely due to thermal pressure support. In other words, if we accept the mass estimate derived from the lensing analysis and demand that the X-ray surface brightness profile be consistent with the observations, the required gas temperature is a factor of $2$--$2.5$ higher than observed. The results for the third cluster, A2163 (the most luminous and the hottest cluster known), are more ambiguous. The discrepancy between the X-ray and the lensing mass estimates arise because the presence of arcs imply a highly concentrated cluster mass distribution whereas the observed X-ray profiles imply a more extended mass distribution. The large X-ray core radii are not the result of the limited resolution of the X-ray detectors. We consider various possibilities that could account for the discrepancy.


title: Neutrinos from the Sun: experimental results confronted with solar models

abstract: For standard neutrinos, recent solar neutrino results together with the assumption of a nuclearly powered Sun imply severe constraints on the individual components of the total neutrino flux: \Phi_{Be}<0.7*10^9cm^-2 s^-1, \Phi_{CNO}< 0.6*10^9 cm^-2 s^-1, and $64*10^9 cm^-2 s^-1< \Phi_{pp+pep} < 65*10^9 cm^-2 s^-1 (at 1 \sigma level). The bound on \Phi_{Be} is in strong disagreement with the standard solar model prediction \Phi_{Be}^{SSM}\approx 5*10^9 cm^-2 s^-1. We study a large variety of non-standard solar models with low inner temperature, finding that the temperature profiles T(m) follow the homology relationship: T(m)=kT^{SSM}(m), so that they are specified just by the central temperature T_c. There is no value of T_c which can account for all the available experimental results. Even if we only consider the Gallium and Kamiokande results, they remain incompatible. Lowering the cross section p+Be7 \to \gamma+B8 is not a remedy. The shift of the nuclear fusion chain towards the pp-I termination could be induced by a hypothetical low energy resonance in the He3+He3 reaction. This mechanism gives a somehow better, but still bad fit to the combined experimental data. We also discuss what can be learnt from new generation experiments about the properties of neutrinos and of the Sun.


title: Winds from Neutron Stars and Strong Type I X--Ray Bursts

abstract: A model for stationary, radiatively driven winds from X--ray bursting neutron stars is presented. General relativistic hydrodynamical and radiative transfer equations are integrated from the neutron star surface outwards, taking into account for helium nuclear burning in the inner, dense, nearly hydrostatic shells. Radiative processes include both bremsstrahlung emission--absorption and Compton scattering; only the frequency--integrated transport is considered here. It is shown that each solution is characterized by just one parameter: the mass loss rate $\Mdot$, or, equivalently, the envelope mass $\Menv$. We found that, owing to the effects of Comptonization, steady, supersonic winds can exist only for $\Mdot$ larger than a limiting value $\Mdot_{min} \approx\Mdot_{E}$. Several models, covering about two decades in mass loss rate, have been computed for given neutron star parameters. We discuss how the sequence of our solutions with decreasing $\Menv$ can be used to follow the time evolution of a strong X--ray burst during the expansion/contraction phase near to the luminosity maximum. The comparison between our numerical results and the observational data of Haberl {\it et al.\/} (1987) for the bursts in 4U/MXB 1820-30 gives an estimate for both the spectral hardening factor and the accretion rate in this source.


title: Accretion Rates in X--Ray Bursting Sources

abstract: We present estimates for the accretion rates in 13 X--ray bursting sources which exhibit photospheric expansion, basing on theoretical models of stationary, radiatively driven winds from neutron stars. The relatively high values obtained, $\Mdot_{acc}\magcir 10^{-9} \MS$, are in accordance with theoretical limits for unstable helium burning, and, at the same time, almost never exceed the ``dynamical'' limit for stationary accretion, $\sim 10 \Mdot_{Edd}$. The only exceptions are 1820-30, already known to be a very peculiar object, and 1608-522; there are indications, however, that in both sources, accretion could be non--stationary.


title: Hydrodynamics and High-Energy Physics of WR Colliding Winds

abstract: One of the main properties of Wolf-Rayet (WR) stars is a very intense outflow of gas. No less than 40\% \ of WR stars belong to binary systems. Young massive O and B stars are the secondary components of such systems. OB stars also have an intense stellar wind. If the intensities of the stellar winds of WR and OB stars are more or less comparable or if the distance between the components of the binary is large enough, the winds flowing out of WR and OB stars can collide and the shock waves are formed. In the shock the gas is heated to temperature $\sim 10^7$ K and generates X-ray emission. Stellar wind collision may be responsible not only for the X-ray emission of WR + OB binaries and for their radio, IR and $\gamma$-ray emision as well. Stellar wind collision, gas heating, particle acceleration, and generation of X-ray, $\gamma$-ray, radio and IR emission in WR + OB binaries are discussed.


title: Abundance Analysis and Origin of the Zeta Sculptoris Open Cluster

abstract: We have determined chemical abundances and radial velocities for stars in the field of the zeta Sculptoris cluster. We find that the cluster metal deficiency previously found from UBV photometry is too high; the cluster overall metallicity, [Fe/H]=+0.24, is about 70% higher than the solar value. The chemical abundance pattern is unusual: the Ni/Fe ratio is significantly lower than typical for field stars with similar Fe and Ca abundances, and the cluster may also be deficient in the alpha element Si. For its age, the cluster is unusually far from the galactic disk, approximately 240 pc. The adopted heliocentric cluster radial velocity, +3.9 km/s, shows that it is close to its maximum distance from the galactic plane. We suggest that the zeta Sculptoris cluster was formed in the galactic disk 45 Myr ago by the interaction of a high velocity cloud with the interstellar medium, and that its formation may be connected with that of Gould's Belt. (Table3, postscript figures and text available by anonymous ftp from "chiron.astro.uu.se" in directory pub/articles/atmos/P90)


title: The EMSS catalog of X-ray-selected clusters of galaxies. I. An atlas of CCD images of 41 distant clusters

abstract: An atlas of deep, wide-field R-band CCD images of a complete sample of distant, X-ray-selected clusters of galaxies is presented. These clusters are the 41 most distant (z>= 0.15) and most X-ray luminous (Lx >= 2x10**44 erg\s-1) clusters in the Einstein Observatory Extended Medium Sensitivity Survey (EMSS) catalog that are observable from Mauna Kea (delta > -40 deg). The sample spans a redshift range of 0.15 <= z < 0.81 and includes at least 2 and possibly as many as 6 rich clusters with z>0.5. For the most part, the data are of superior quality, with a median seeing of 0''.8 FWHM and coverage of at least 1 Mpc X 1 Mpc in the cluster frame (Ho=50; q0=1/2). In addition, we update the available optical, X-ray and radio data on the entire EMSS sample of 104 clusters. We outline the cluster selection criteria in detail, and emphasize that X-ray-selected cluster samples may prove to be more useful for cosmological studies than optically selected samples. The EMSS cluster sample in particular can be exploited for diverse cosmological investigations, as demonstrated by the detection of evolution in the X-ray luminosity function previously reported, and more recently by the discovery of a large number of gravitationally lensed images in these clusters.


title: Observational Support for the Gurzadyan-Kocharyan Relation in Clusters of Galaxies

abstract: We show that observational data for four Abell clusters of galaxies support the Gur\-za\-dyan-Kocharyan relation between the Hausdorff dimension and the dynamical properties of a galaxy system. The Hausdorff dimension is calculated using the two-point correlation function, while the dynamical parameters are estimated using available data and reasonable assumptions on the mass function of galaxies. This result can have essential consequences in the understanding of the dynamical mechanisms that determine the fractal distribution of galaxies.


title: Ejection of Fragments in Supernova Explosions

abstract: Recent observations by the ROSAT X-ray satellite of the Vela supernova remnant have revealed, in addition to the previously identified compact nebula, a nearly circular emitting region with a radius of about 4 degrees. The Vela pulsar is slightly off the center of this circular region, consistent with its measured proper motion of about 100km/s and an age of about 10^4 yr. The emitting region is bounded by the main supernova shock. Just outside the shock, the X-ray image reveals several well-defined V-shaped features extending radially outwards. These features are most likely wakes produced by objects moving supersonically through the outside medium. The shapes and orientations of the wakes suggest that these objects have been ejected from the center of the supernova explosion. Their present positions indicate that they have been moving with a mean velocity of a few thousand km/s. We show that pre-existing objects, such as planets in orbit around the progenitor star, could not have been accelerated to sufficiently high velocities or would have been destroyed. Instead, we propose that the observed objects are fragments ejected during the formation of the neutron star. Fragmentation during gravitational collapse is a natural consequence of both convective and rotational instabilities.


title: Reionization and the Cosmic Microwave Background

abstract: If the COBE detection of CMB fluctuations is used to normalize the power spectrum of primordial density fluctuations in a cold dark matter--dominated universe, early reionization is likely to result in a substantial diminution of primordial temperature fluctuations on degree scales. I argue that the reionization may be non--Gaussian because of feedback effects. Secondary fluctuations on arc--minute scales provide an important probe of the efficiency of rescattering.


title: The Primordial Abundance of $^6$Li and $^9$be

abstract: Light element ($^6$Li, $^7$Li and $^9$Be) depletion isochrones for halo stars have been calculated with standard stellar evolution models. These models include the latest available opacities and are computed through the sub-giant branch. If $^6$Li is not produced in appreciable amounts by stellar flares, then the detection of $^6$Li in HD 84937 by Smith, Lambert \& Nissen (1993) is compatible with standard stellar evolution and standard big bang nucleosynthesis only if HD 84937 is a sub-giant. The present parallax is inconsistent with HD 84937 being a sub-giant star at the $2.5\, \sigma$ level. The most metal poor star with a measured $^9$Be abundance is HD 140283, which is a relatively cool sub-giant. Standard stellar evolution predict that $^9$Be will have been depleted in this star by $\sim 0.3$ dex (for ${\rm T_{eff}} = 5640$ K). Revising the abundance upward changes the oxygen-beryllium relation, suggesting incompatible with standard comic ray production models, and hence, standard big bang nucleosynthesis. However, an increase in the derived temperature of HD 140283 to 5740 K would result in no depletion of $^9$Be and agreement with standard big bang nucleosynthesis.


title: Evolution in the Clustering of Galaxies to r = 26 (TAPIR GRP-390)

abstract: We present results for the two-point angular correlation function of galaxies to a limiting magnitude of r=26. The final sample is 97% complete to r=26.0, yielding 5730 galaxies over a 90.1 sq. arcmin field. The correlation function for our faint galaxies can be parameterised by a power law, $A \theta^{-0.8}$, in agreement with the clustering statistics of shallower catalogues. The derived amplitude, $A$, is small, but non-zero. We combine this measurement with the latest statistical constraints on faint galaxy redshifts from gravitational lensing studies, which imply that the bulk of the r<26 field galaxies should be at redshifts of order 1. Our derived $A$ is significantly lower than that predicted from the local bright galaxy correlation function using the lensing-determined galaxy redshift distribution and modest growth of clustering. This simplistic model does not include the variation in observed morphological mix as a function of redshift and apparent magnitude in our sample. At our faintest limits we reach sufficiently high redshifts that differential $K$ corrections will result in the observed galaxy mix being dominated by star bursting dwarf and low surface brightness irregulars, rather than the early-type systems used to define the local bright galaxy correlation function. Adopting the correlation function measured locally for these low surface brightness galaxies and assuming modest clustering evolution, we obtain reasonable agreement between our model and observations. This model supports the scenario in which the high number density of faint galaxies is produced by normally clustered star forming dwarf galaxies at modest redshifts.


title: Probing the Large-Scale Velocity Field with Clusters of Galaxies

abstract: What is the role of clusters of galaxies in probing the large-scale velocity field of the universe? We investigate the distribution of peculiar velocities of clusters of galaxies in the popular low-density ($\Omega=0.3$) flat Cold-Dark-Matter (CDM) cosmological model, which best fits many large-scale structure observations. An $\Omega=1$ CDM model is also studied for comparison. We find that clusters of galaxies are efficient tracers of the large-scale velocity field. The clusters exhibit a Maxwellian distribution of peculiar velocities, as expected from Gaussian initial density fluctuations. The cluster 3-D velocity distribution for the $\Omega=0.3$ model peaks at $v \sim 400$ km s$^{-1}$, and extends to high velocities of $v \sim 1200$ km s$^{-1}$. The rms peculiar velocity of the clusters is $440$ km s$^{-1}$. Approximately 10\% of all model clusters move with high peculiar velocities of $v \ge 700$ km s$^{-1}$. The observed velocity distribution of clusters of galaxies is compared with the predictions from cosmological models. The observed data exhibit a larger velocity tail than seen in the model simulations; however, due to the large observational uncertainties, the data are consistent at a $\sim 3\sigma$ level with the model predictions, and with a Gaussian initial density field. The large peculiar velocities reported for some clusters of galaxies ($v \geq 3000$ km s$^{-1}$) are likely to be overestimated, if the current model is viable.


title: Imprints of Galaxy Clsutering Evolution on Delta T / T

abstract: The nonlinear evolution of matter clustering induces temperature anisotropies in the cosmic microwave background. The fluctuations in temperature are produced at a relatively low redshift, $z\simlt 10$, and consequently are not affected by the uncertain reionization proccesses that occurred at earlier epochs. The amplitude of the effect depends on the evolution model assumed for galaxy clustering. For plausible parameter choices, the temperature fluctuations are in the range $10^{-6}\simlt \Delta T/T \simlt 10^{-5}.$ Observed limits on $\Delta T/T $ would be violated if nonlinear clustering on present cluster scales commenced prior to $z\sim 10.$


title: Compton Scattering in Jets: A Mechanism for $\sim$ 0.4 and $\lsim$ 0.2 Mev Line Production

abstract: We show that gamma ray line emission at $\sim$0.4 MeV and $\lsim 0.2$ MeV can be produced by Compton scattering of beamed radiation in the jets of Galactic black hole candidates. This mechanism has the novel feature of not invoking the presence of $e^+$--$e^-$ pairs. To produce the two lines we employ a symmetric double sided jet with bulk flow velocity of about 0.5 c and incident beam radiation with a hard energy spectrum. We show that the two lines can be seen at viewing angle cosines relative to the jet ranging from 0.2 to 0.6. This comprises 40\% of the total solid angle. In addition, the line radiation is approximately 10\% polarized. Depending on the bulk flow and viewing angle the model can produce lines at other energies as well. In particular a broad feature near 1 MeV can be seen by viewing the jet close to its axis. Our model can also accommodate single line spectra if the beamed gamma ray emission or the jets themselves are asymmetric.


title: The Companion of M4a: A Planet or a Star?

abstract: The $\ddot P$ (Backer 1992, Backer et al. 1993, Thorsett et al. 1993) observed for PSR~B1620-26 in the globular cluster M4 is most likely due to external gravitational jerk from a bound companion. Solving for the $\ddot P$ observed, constrained by the observed $\dot P$, two natural solutions are a jovian mass in a $\gtorder 10 \ AU, e=0.3-0.5$ orbit, or a solar mass companion in a $\sim 50 \ AU$ highly eccentric orbit (Michel 1994, Rasio 1994). If the companion is solar mass, the system must be young, and its location outside the cluster core is somewhat of a mystery. If the companion is planetary, it was either exchanged with the secondary, or angular momentum transport in excretion disks and planet formation is extraordinarily efficient. The implications for planet formation around stellar systems are potentially profound.


title: Modeling the Radial Distribution of Blue Stragglers in M3

abstract: Recent surveys of the blue straggler (BS) population in the galactic globular cluster M3 (NGC 5272) give the first complete characterization of the number density of BSs as a function of radius over an entire globular cluster. The BSs in M3 are overabundant at large radii, and significantly underabundant at intermediate radii. Here we present the result of a simulation of the dynamical evolution of a population of BSs in a multi--mass model of M3. Assuming the BSs were formed in the core through binary interactions (Hut \&\ Verbunt 1983, Leonard 1989, Sigurdsson \&\ Phinney 1993, Hut \etal 1992, Davies, Benz and Hills 1994), and given some very general assumptions about the recoil that occurs during stellar mergers in interacting binaries; we find an excellent fit to the observed radial distribution of BSs, suggesting strongly that most of the BSs in M3 were formed through binary collisions in the core.


title: Effects of Photoionization on Galaxy Formation

abstract: We investigate the effects of the observed UV background radiation on galaxy formation. Photoionization by UV radiation decreases the cooling rate of the gas in halos, so that objects with only large density contrasts can self-shield against the background radiation, thereby allowing the shielded neutral cores to cool and form stars. In the context of the CDM model, we use the criterion that self-shielding is essential for star formation to calculate the mass function of galaxies based on both Press-Schechter and the peaks formalisms. The ionizing UV radiation results in the inhibition of galaxy formation --- we show the decrease in the number density of the galaxies quantitatively. We also find that the merging in general is made inefficient by the UV photons through photoionization of the gas in the bigger system into which small objects are incorporated. The latter means that, in a merging dominated region, where the number density at the low-mass end ($M_b \la 10^{10}M_{\odot}$) is usually expected to decrease with time, the trend is {\it reversed} (the number of low-mass galaxies is more at lower redshifts $z$) due to the decreasing UV flux with time after $z \sim 2$. We further discuss the implication of our results for the number counts of galaxies and possible evolution of the luminosity function of the galaxies.


title: On the Nature of Photospheric Oscillations in Strong X--Ray Bursts

abstract: A possible sound origin for the photospheric oscillations in the X--ray bursting sources 1608-522 and 2127+119 is suggested. It is shown that standing sound waves in an expanding spherical envelope can have periods very close to the observed ones. The quite large ratio, $\sim$ 10, of the periods in the two sources is explained in terms of different wave regimes. The relevance of sound oscillations to the observed QPO in type II bursts of the Rapid Burster is also discussed.


title: On the nature of gamma-ray burst time dilations

abstract: The recent discovery that faint gamma-ray bursts are stretched in time relative to bright ones has been interpreted as support for cosmological distances: faint bursts have their durations redshifted relative to bright ones. It was pointed out, however, that the relative time stretching can also be produced by an intrinsic correlation between duration and luminosity of gamma-ray bursts in a nearby, bounded distribution. While both models can explain the average amount of time stretching, we find a generic difference between them in the way the duration distribution of faint bursts deviates from that of bright ones. This allows us to distinguish between these two broad classes of model on the basis of the duration distributions of gamma-ray bursts, leading perhaps to an unambiguous determination of the distance scale of gamma-ray bursts. We apply our proposed test to the second BATSE catalog and conclude, with some caution, that the data favor a cosmological interpretation of the time dilation.


title: Cosmological Simulations of Galaxy Formation Including Hydrodynamics (hyper-abridged)

abstract: The formation of galaxies in hierarchical cosmogonies is studied using high resolution N-body plus SPH hydrodynamics simulations. The collapse of structure is followed self-consistently from Mpc scale filamentary structures to kpc scale galactic objects. The characteristics and formation processes of the galaxy like objects are studied in detail, along with the aggregation into a poor cluster. Related studies consider the effects of modelling star formation, the reliability of tracing galaxies in simulations, and tests of SPH methods. This submission serves first to notify that the full text and figures of my thesis are available in compressed PostScript form via anonymous ftp from astro.princeton.edu in the directory /summers/thesis (122 files, 19 MB compressed, 65 MB uncompressed). See the README file first. Second, this submission contains the title page, abstract, table of contents, introductory chapter, summary chapter, and references for my thesis. Those who are curious about the work may scan these pages to identify which chapters may be interesting to get via ftp.


title: Wiener Reconstruction of Density, Velocity, and Potential Fields from All-Sky Galaxy Redshift Surveys

abstract: We present a new method for recovering the cosmological density, velocity, and potential fields from all-sky redshift catalogues. The method is based on an expansion of the fields in orthogonal radial (Bessel) and angular (spherical harmonic) functions. In this coordinate system, peculiar velocities introduce a coupling of the radial harmonics describing the density field in redshift space but leave the angular modes unaffected. In the harmonic transform space, this radial coupling is described by a distortion matrix which can be computed analytically within the context of linear theory; the redshift space harmonics can then be converted to their real space values by inversion of this matrix. Statistical noise is mitigated by regularizing the matrix inversion with a Wiener filter. The method yields a minimum variance estimate of the density field in real space. In this coordinate system, the minimum variance harmonics of the peculiar velocity and potential fields are related to those of the density field by simple linear transformations. Tests of the method with simulations of a CDM universe and comparison with previously proposed methods demonstrate it to be a very promising new reconstruction method for the local density and velocity field. A first application to the 1.2 Jy IRAS redshift survey is presented.


title: Introduction: Why and How to Probe the Zone of Avoidance ?

abstract: The motivation and major ways for probing the Zone of Avoidance (ZOA) are reviewed. Galaxies hidden behind the ZOA may have important implications for the internal dynamics of the Local Group, for the origin of its motion relative to the Microwave Background, and for the connectivity of the large scale structure. Current direct (`observational') methods for exploring the ZOA include eye-balling of plates, source identification in the IRAS data base and pointed and blind-search observations in 21 cm. Interesting regions identified so far include the two crossing points of the Supergalactic Plane by the Galactic Plane (at Galactic longitude $l \sim 135^o$, near Perseus-Pisces, and $l \sim 315^o$, near the Great Attractor), the Puppis cluster (at $l \sim 240^o, cz \sim 1500$ km/sec) and the Ophiuchus cluster (at $l \sim 0^o, cz \sim 8400$ km/sec). New promising wavelengths are the $2 \mu$ and the X-ray band. Indirect (`theoretical') approaches include `Wiener reconstruction' from incomplete and noisy data, and using the peculiar velocity field as a probe of the mass distribution hidden behind the ZOA. The problem of source confusion at low Galactic latitude can be addressed by novel statistical methods, e.g. Artificial Neural Networks.


title: Wiener Reconstruction of All-Sky Spherical Harmonic Maps of the Large Scale Structure

abstract: A statistical method for reconstructing large scale structure behind the Zone of Avoidance is presented. It also corrects for shot-noise and for redshift distortion in galaxy surveys. The galaxy distribution is expanded in an orthogonal set of spherical harmonics. We show that in the framework of Bayesian statistics and Gaussian random fields the $4 \pi$ harmonics can be recovered and the shot-noise can be suppressed, giving the optimal picture of the underlying density field. The correction factor from observed to reconstructed harmonics turns out to be the well-known Wiener filter (the ratio of signal to signal+noise), which is also derived by requiring minimum variance. We apply the method to the 1.2 Jy IRAS survey. A reconstruction of the projected galaxy distribution confirms the connectivity of the Supergalactic Plane across the Galactic Plane (at Galactic longitude $l \sim 135^o$ and $l \sim 315^o$) and the Puppis cluster behind the Galactic Plane ($ l \sim 240^o$). The method is extended to 3-D, and is used to recover from the 1.2 Jy redshift survey the density, velocity and potential fields in the local universe.


title: New Evidence for the Cosmological Origin of Gama-Ray Bursts

abstract: We compare the burst distribution of the new (2B) BATSE catalogue to a cosmological distribution. We find that the distribution is insensitive to cosmological parameters such as Omega and Lambda and to the width of the bursts luminosity function. The maximal red shift of the long bursts is ~2.1 (assuming no evolution) while Zm(long) of the short bursts is significantly lower Zm(short) =~ 0.5 In agreement with this relatively nearby origin of the short burst we find an indication that these bursts are correlated ( >=2 sigma level at 10 degrees) with Abell clusters. This is the first known correlation of the bursts with any other astrophysical population and if confirmed by further observations it will provides additional evidence for the cosmological origin of those bursts.


title: Perturbative Lagrangian Approach to Gravitational Instability

abstract: This paper deals with the time evolution in the matter era of perturbations in Friedman-Lemaitre models with arbitrary density parameter $\Omega$, with either a zero cosmological constant, $\Lambda = 0$, or with a non-zero cosmological constant in a spatially flat Universe. Unlike the classical Eulerian approach where the density contrast is expanded in a perturbative series, this analysis relies instead on a perturbative expansion of particles trajectories in Lagrangian coordinates. This brings a number of advantages over the classical analysis. In particular, it enables the description of stronger density contrasts. Indeed the linear term is the famous Zel'dovich approximate solution (1970). We present here a systematic and detailed account of this approach. We give analytical results (or fits to numerical results) up to the third order. We then proceed to explore the link between the lagrangian description and statistical measures. We show in particular that Lagrangian perturbation theory provides a natural framework to compute the effect of redshift distortions, using the skewness of the density distribution function as an example. Finally, we show how well the second order theory does as compared to other approximat- ions in the case of spherically symmetric perturbations. We also compare this second order approximation and Zel'dovich solution to N-body simulations in the description of large-scale structure formation starting from a power law (n=-2) power spectrum of Gaussian perturbation. We find that second order theory is both simple and powerful.


title: Non--Gaussian Effects in the Cosmic Microwave Background from Inflation

abstract: The presence of non--Gaussian features in the Cosmic Microwave Background (CMB) radiation maps represents one of the most long--awaited clues in the search for the actual structure of the primordial radiation, still needing confirmation. These features could shed some light on the non trivial task of distinguishing the real source of the primeval perturbations leading to large scale structure. One of the simplest non--Gaussian signals to search is the (dimensionless) skewness ${\cal S}$. Explicit computations for ${\cal S}$ are presented in the frame of physically motivated inflationary models (natural, intermediate and polynomial potential inflation) in the hope of finding values in agreement with estimated quantities from large angle scale (e.g., {\em COBE} DMR) maps. In all the cases considered the non--Gaussian effects turn out to lie below the level of theoretical uncertainty (cosmic variance). The possibility of unveiling the signal for ${\cal S}$ with multiple--field models is also discussed.


title: The Flattened Dark Halo of Polar Ring Galaxy NGC 4650A: A Conspiracy of Shapes?

abstract: Kinematics and photometry of the polar ring galaxy NGC 4650A, including new observations of the rotation and velocity dispersion of its central stellar disk, are used to infer the presence of a dark matter halo and to measure its shape. Fits to the observed disk and polar ring rotation curves from detailed mass and photometric modeling rule out a spherical dark halo. The best fit models have halos with isodensity surfaces that are flattened to a shape between E6 and E7 (axis ratios between 0.4 and 0.3); the asymptotic equatorial speeds of these models are in excellent agreement with the I-band Tully-Fisher relation. This degree of dark halo flattening is larger than that expected from N-body collapse simulations of dissipationless dark matter. Since the kinematics and surface brightness profile of the central luminous body indicate that its light has an intrinsic axis ratio $c/a \lsim 0.4$, in NGC 4650A the radial ``conspiracy" between the dark and luminous components that leads to flat rotation curves may extend to the {\it shape} of the mass distribution as well.


title: Lagrangian Dynamics in Non-Flat Universes and Non-Linear Gravitational Evolution

abstract: We present a new formalism which allows to derive the general Lagrangian dynamical equations for the motion of gravitating particles in a non--flat Friedmann universe with arbitrary density parameter $\Omega$ and no cosmological constant. We treat the set of particles as a Newtonian collisionless fluid. The non--linear dynamical evolution of the fluid element trajectories is then described up to the third--order expansion in Lagrangian coordinates. This work generalizes recent investigations carried out by Bouchet \etal (1992) and Buchert (1994).


title: Constraints to the Masses of Brown Dwarf Candidates from the Lithium Test

abstract: We present intermediate dispersion (0.7-2.2 \AA ~pix$^{-1}$) optical spectroscopic observations aimed at applying the ``Lithium Test'' to a sample of ten brown dwarf candidates located in the general field, two in young open clusters, and two in close binaries. We find evidence for strong Li depletion in all of them, and thus infer lower mass limits of 0.065~M$_\odot$, depending only slightly ($\pm$0.005~M$_\odot$) on the interior models. None of the field brown dwarf candidates in our sample appears to be a very young (age $<$~10$^8$~yr) substellar object. For one of the faintest proper motion Pleiades members known (V=20.7) the Li test implies a mass greater than $\sim$0.08~M$_\odot$, and therefore it is not a brown dwarf. From our spectra we estimate spectral types for some objects and present measurements of Halpha emission strengths and radial velocities. Finally, we compare the positions in the H-R diagram of our sample of brown dwarf candidates with the theoretical region where Li is expected to be preserved (Substellar Lithium Region). We find that certain combinations of temperature calibrations and evolutionary tracks are consistent with the constraints imposed by the observed Li depletion in brown dwarf candidates, while others are not.


title: Abundances of four very metal-poor stars of the BPS survey

abstract: High resolution CASPEC and EMMI spectra (R$\approx$10$^4$) of the faint metal- deficient candidates CS 22891--209, CS 22897--8, CS 22948--66 and CS 22968--14 taken from the survey of Beers, Preston and Shectman (1985, BPS) are analyzed for abundances. The fine analysis with LTE model atmospheres gives an iron metallicity for the four stars between $-$3.15 and $-$3.45 in good agreement with the abundances inferred from CaII K index by Beers et al. (1992). These four objects considerably augment the number of spectroscopically studied stars at the low metallicity extreme. We found the $\alpha$-elements Mg, Ca, Si, and Ti enhanced, and the iron-peak elements Sc, Mn and Cr tracking closely iron with the typical pattern of the more metallic halo stars. Aluminum results deficient with respect both to iron and to magnesium showing a strong {\it odd-even\/} effect, but for [Fe/H] $\leq$ $-$3.0 [Al/Mg] levels out into a plateau around $-$1.0, showing a primary behaviour. [Na/Fe] is $\approx$ $-$0.2 and only compared to magnesium it shows a mild odd-even effect. The heavy elements Sr and Ba show a large scatter which is not observed for other elements and cannot be due to observational errors only. In the case of CS 22891--209 they are particularly overabundant with [Sr,Ba/Fe] $\approx$ 0.9. The spread might reflect the presence of strong inhomogeneities in the interstellar gas at the earliest stages of Galaxy formation. This might be an important signature of the still elusive processes of the early formation of these elements observed in old stars.


title: M Dwarfs, Microlensing, and the Mass Budget of the Galaxy

abstract: We show that faint red stars do not contribute significantly to the mass budget of the Galaxy or to microlensing statistics. Our results are obtained by analyzing two long exposures of a high-latitude field taken with the Wide Field Camera on the newly repaired {\it Hubble Space Telescope}. Stars are easily distinguished from galaxies essentially to the limiting magnitudes of the images. We find five stars with $2.0<V-I<3.0$ and $I<25.3$ and no stars with $V-I>3.0$. Therefore, main-sequence stars with $M_I>10$ that are above the hydrogen-burning limit in the dark halo or the spheroid contribute $<6%$ of the unseen matter. Faint red disk stars, M-dwarfs, contribute at most $ 15%$ to the mass of the disk. We parameterize the faint end of the cumulative distribution of stars, $\Phi$, as a function of luminosity $L_V$, $d\Phi/d \ln L_V \propto L_V^{-\gamma}$. For spheroid stars, $\gamma<0.32$ over the range $6<M_V<17$, with 98% confidence. The disk luminosity function falls, $\gamma<0$, for $15\lsim M_V\lsim 19$. Faint red stars in the disk or thick disk, and stars with $M_V<16$ in the spheroid contribute $\tau< 10^{-8}$ to the optical depth to microlensing toward the Large Magellanic Cloud.


title: Power Spectrum of Cosmic String Perturbations on the Microwave Background

abstract: I review recent progress towards a detailed understanding of the power spectrum of cosmic string perturbations induced on the Microwave Background (CMB). The use of a simple analytic model allows to include important effects that have not been included in previous studies. These effects are potential fluctuations on the last scattering surface (the Sachs-Wolfe effect) and Doppler CMB perturbations. These previously neglected fluctuations are shown to dominate over post-recombination fluctuations (Gott-Kaiser-Stebbins effect) on angular scales less than about $2^\circ$ assuming that no reionization occurs. The effective power spectrum index $n_{eff}$ is calculated on COBE scales and is shown to be somewhat larger than 1 ($1.35 \lsim n_{eff} \lsim 1.5$). After normalizing the model on COBE data and string simulations, I derive the predicted values of $({{\delta T}\over T})_{rms}$ for several ongoing CMB experiments on medium and small angular scales. A comparison is also made with current observations and predictions of models based on inflation.


title: $\gamma$-Ray Pulsars: Emission Zones and Viewing Geometries, A Computer Animation

abstract: The computer animation illustrates the geometries described in a paper by the same authors. The preprint is available as number 9401045. The opening scene shows dipole field lines emanating from the polar caps of a rotating neutron star. The dipole axis is inclined along the green rods. The field lines shown are defined from the condition that they be tangent to the light cylinder (the cylindrical radius at which the tangential velocity of rotation reaches the speed of light). The static dipole field lines are smoothly morphed into the correct retarted-potential vacuum solutions. A red surface spanning these field lines is painted. In the next scene the blue surfaces represent the outer gaps above the surface of last closed field lines. High energy emission (blue) is produced in these outer gaps, and is beamed tangentially along the field lines. The radio emission (green) originates close to the surface of the star and is beamed along the dipole axes. The inclination angle $\alpha$ of the dipole and the viewing angle $\zeta$ are chosen to match the Crab parameters; $\alpha$ = 70, $\zeta$ = 65. The corresponding light curve is computed and shown for these angles, and the red dot traces rotation phase. The next scene shows the situation for angles appropriate to PSR1706-44; $\alpha$ = 45, $\zeta$ = 65. The final scene is a possibility for Geminga; $\alpha$ = 20, $\zeta$ = 75. These angles are poorly constrained as there is no radio emission.


title: Secular Evolution of the Eccentricity in the PSR 1620-26 Triple System

abstract: A simple analytic calculation is presented for the secular evolution of the eccentricities in a hierarchical triple system such as the one containing the millisecond pulsar PSR 1620-26 in M4. If the second companion of PSR 1620-26 is of stellar mass ($m_2\ga0.1\,M_\odot$), an eccentricity as large as that observed today for the inner binary ($e_1\simeq0.03$) could very well have been induced by secular perturbations. In contrast, such a large eccentricity cannot be induced by a second companion of planetary mass.


title: Environmental Effects on Local Active Galactic Nuclei

abstract: Using an extensive sample of nearby galaxies (the Nearby Galaxies Catalog, by Tully), we investigate the environment of the galaxies hosting low-luminosity AGNs (Seyferts and LINERs). We define the local galaxy density, adopting a new correction for the incompleteness of the galaxy sample at large distances. We consider both a complete sample of bright and nearby AGNs, identified from the nuclear spectra obtained in available wide optical spectroscopic surveys, and a complete sample of nearby Seyferts. Basically, we compare the local galaxy density distributions of the AGNs with those of non-AGN samples, chosen in order to match the magnitude and morphological type distributions of the AGN samples. We find, only for the early-type spirals more luminous than $\sim M^*$, that both LINERs and Seyferts tend to reside in denser environments on all the scales tested, from tenths of Mpc to a few Mpc; moreover Seyferts show an enhanced small-scale density segregation with respect to LINERs. This gives support to the idea that AGNs can be stimulated by interactions. On larger scales, tens of Mpc, we find that the AGNs hosted in luminous early-type spirals show a tendency to stay near the center of the Local Supercluster. Finally we discuss the interpretations of our findings and their consequences for some possible scenarios of AGN formation and evolution and for the problem of how AGNs trace the large-scale structures.


title: Stability of Fireballs and $\gamma$-Ray Bursts

abstract: Fireballs are an essential part of any cosmological $\gamma$-ray burst. We drive here a stability criterion for fireballs and show that fireballs are Rayleigh-Taylor unstable in any region in which the entropy decreases outward. The instability begins to operate when the fireball becomes matter dominated. Among the possible implication of the instability are: (i) Conversion of a fraction of the radiation energy to a convective energy expressed in the motion of bubbles relative to each other. (ii) Penetration of fast bubbles through slower ones and creation of high $\gamma$ regimes which are essential for efficient conversion of the energy to $\gamma$-rays. (iii) Formation of rapid variation (of the scale of the bubbles) in the observed $\gamma$-rays.


title: Integrated photometric properties of open clusters

abstract: Galactic open clusters provide an abundant sample of stellar aggregates of various sizes, ages and metal abundances, apt to constitute a template for comparison with star systems in other galaxies. In this paper we present and discuss a standard methodology to synthesize {\sl U,B,V} fluxes and colours, and apply it to a set of 138 open clusters. Results are compared with previous available ones in the literature. We were able to calibrate a mass--luminosity relation by which we evaluated the mass of $\simeq 400$ open clusters, leading to a well defined present--day mass function. The number--complete sample of galactic open clusters presented in Battinelli \& Capuzzo--Dolcetta (1991) is enlarged of a $15\%$.}


title: The Imprint of $\Omega$ on the Cosmic Microwave Background

abstract: We investigate the dependence of large angular scale cosmic microwave background anisotropies on various initial conditions, including both adiabatic and isocurvature perturbations and the initial power law index $n$, in a variety of low-$\Omega$ cosmological models. Cosmological constant- dominated flat models and inflationary open adiabatic models, with $n_{\rm eff}<1$, and open isocurvature models, with $n_{\rm eff} \simeq 2$, are significantly constrained by the present observations.


title: New Automatic Identification Technique for OB Associations in Unresolved Galaxies

abstract: We present a new automatic technique based on Principal Component Analysis and Cluster Analysis, with the aim of its application to the identification of those clumps in unresolved galaxies which likely represent regions of star formation. We test the method by applying it to the galaxy M~31, for which there are already several sets of identifications of OB associations based on multi--colour images of resolved stars. We use small--scale digital images of M~31 and compare the associations that we detect from these unresolved data with previously--published large--scale data, finding a rather good agreement. We obtain a strict agreement of our identification with the most compact associations of the original van den Bergh (1964) identification. We then apply the technique to CCD images of the more distant spiral galaxy NGC~2903 and identify 68 OB association candidates.


title: Black Holes from Blue Spectra

abstract: Blue primordial power spectra with a spectral index $n>1$ can lead to a significant production of primordial black holes in the very early Universe. The evaporation of these objects leads to a number of observational consequences and a model independent upper limit of $n \approx 1.4$. In some cases this limit is strengthened to $n=1.3$. Such limits may be employed to define the boundary to the region of parameter space consistent with generalized inflationary predictions. [To appear in Proceedings of the CASE WESTERN CMB WORKSHOP, April 22-24 1994. Figures available on request from J.H.Gilbert@qmw.ac.uk]


title: The Mass Function of Cosmic Structures with Non-Spherical Collapse

abstract: Non-spherical dynamical approximations and models for the gravitational collapse are used to extend the well-known Press \& Schechter (PS) approach, in order to determine analytical expressions for the mass function of cosmic structures. The problem is rigorously set up by considering the intrinsic Lagrangian nature of the mass function. The Lagrangian equations of motion of a cold and irrotational fluid in single-stream regime show that the shear, which is non-locally determined by all the matter field, is the quantity which characterizes non-spherical perturbations. The Zel'dovich approximation, being a self-consistent first-order Lagrangian and local one, is used as a suitable guide to develop realistic estimates of the collapse time of a mass clump, starting from the local initial values of density and shear. Both Zel'dovich-based \an\ and models and the homogeneous ellipsoidal model predict that more large-mass objects are expected to form than the usual PS relation. In particular, the homogeneous ellipsoid model is consistent at large masses with a Press \& Schechter mass function with a lower value of the \dc\ parameter, in the range 1.4$\div$1.6. This gives a dynamical explanation of why lower \dc\ values have been found to fit the results of several N-body simulations. When more small-scale structure is present, highly non-linear dynamical effects can effectively slow down the collapse rate of a perturbation, increasing the effective value of \dc. This may have interesting consequences on the abundance of large-mass high-redshift objects.


title: Evidence from Gravitational Lensing for a Non-Thermal Pressure Support in the Cluster of Galaxies A2218

abstract: The central mass distribution of clusters of galaxies can be inferred from gravitationally lensed arcs with known redshifts. For the Abell cluster 2218, this method yields a core mass which is larger by a factor of $2.5 \pm 0.5$ than the value deduced from X-ray observations, under the assumptions that the gas is supported by thermal pressure and that the cluster is spherical. We show that a non-thermal pressure support is the most plausible explanation for this discrepancy. Such a pressure can be naturally provided by strong turbulence and equipartition magnetic fields ($\sim 50\mu$G) that are tangled on small spatial scales (<10 kpc). The turbulent and magnetic pressures do not affect the measured Sunyaev-Zel'dovich effect for this cluster. Intracluster magnetic fields with a comparable magnitude ($\sim 10^{1-2} \mu$G) have already been detected by Faraday rotation in other clusters. If generic, a small-scale equipartition magnetic field should affect the structure of cooling flows and must be included in X-ray determinations of cluster masses.


title: Implications of Neutrino Balls as the Source of Gamma-Ray Bursts

abstract: (To appear in the Astrophysical Journal) Holdom and Malaney (1994) have suggested a mechanism for gamma-ray bursts which requires that stars be captured by a neutrino ball. Neutrino balls would be, for the most part, denser than main sequence stars, but their density would decrease as their mass increased. We show that small neutrino balls would subject stars to tidal forces sufficient to disrupt them. We thus argue that if neutrino balls existed at the centres of galaxies, only the largest would be able to act as a source of gamma-ray bursts. Such neutrino balls would have a mass of order $10^7\Msun$. Tidal capture of stars by a neutrino ball would not be important, but dynamical friction against the neutrinos or star-disc interactions could both be important capture mechanisms. We find that a gamma-ray burst would occur in a galaxy containing such a neutrino ball roughly every $10^2\y$, and the fraction of all galaxies contributing to the gamma-ray burst flux would be $\sim 10^{-4}$, assuming that this was the mechanism of all gamma-ray bursts. These numbers have implications for neutrino ball models of active galaxies, assuming that all gamma-ray bursts and all AGN come from neutrino balls. Either a small fraction $\sim 10^{-2}$ of the lifetime of such an object could be spent as an AGN, or that the probability of a neutrino ball becoming an AGN would be $10^{-2}$. It is not possible to rule out the possibility that neutrino balls might exist at the centres of galaxies through direct ground-based observation of stellar kinematics.


title: Hydrodynamics of Binary Coalescence. II. Polytropes with Gamma=5/3

abstract: We present a new numerical study of the equilibrium and stability properties of close binary systems using the smoothed-particle hydrodynamics (SPH) technique. We adopt a simple polytropic equation of state $p=K\rho^\gam$ with $\gam=5/3$ and $K=\,$constant within each star, applicable to low-mass degenerate dwarfs as well as low-mass main-sequence stars. Along a sequence of binary equilibrium configurations for two identical stars, we demonstrate the existence of both secular and dynamical instabilities, confirming directly the results of recent analytic work. We use the SPH method to calculate the nonlinear development of the dynamical instability and to determine the final fate of the system. We find that the two stars merge together into a single, rapidly rotating object in just a few orbital periods. Equilibrium sequences are also constructed for systems containing two nonidentical stars. These sequences terminate at a Roche limit, which we can determine very accurately using SPH. For two low-mass main-sequence stars with mass ratio $q\lo0.4$ we find that the (synchronized) Roche limit configuration is secularly unstable. Degenerate binary configurations remain hydrodynamically stable all the way to the Roche limit for all mass ratios $q\ne1$. Dynamically unstable mass transfer can also lead to the rapid coalescence of a binary system, but the details of the hydrodynamic evolution are quite different. We discuss the implications of our results for double white-dwarf and W Ursae Majoris systems.


title: Dynamical Model of Pulsar-Planet Systems

abstract: About two years ago, Wolszczan and Frail announced the detection of a possible planetary system consisting of two Earth-mass planets around a millisecond pulsar. It was pointed out shortly thereafter that the mutual gravitational interaction of the planets leads to predictable evolution of their orbits. Because millisecond pulsars are very stable ``clocks'' and allow pulse timing measurements with microsecond precision, the orbit evolution would lead to a measurable modulation of the arrival times of the pulses. Detection of this effect would serve to corroborate the planetary interpretation of the data and provide additional constraints on the system. The theory of planetary perturbations in the context of pulsar timing observations is summarized here.


title: Simulations of Clusters of Galaxies

abstract: The degree of complexity and, to a somewhat lesser degree, realism in simulations has advanced rapidly in the past few years. The simplest approach - modeling a cluster as collisionless dark matter and collisonal, non--radiative gas is now fairly well established. One of the most fruitful results of this approach is the {\sl morphology--cosmology connection} for X-ray clusters. Simulations have provided the means to make concrete predictions for the X-ray morphologies of clusters in cosmologies with different $\Omega_o$, with the result that low $\Omega_o$ cosmologies fair rather poorly when compared to observations. Another result concerns the accuracy of \xray binding mass estimates. The standard, hydrostatic, isothermal model estimator is found to be accurate to typically better than $50\%$ at radii where the density contrast is between $10^2$ and $10^3$. More complicated approaches, which attempt to explicitly follow galaxy formation within the proto--cluster environment are slowly being realized. The key issue of {\sl dynamical biasing} of the galaxy population within a cluster is being probed, but conclusive answers are lacking. The dynamics of multi--phase gas, including conversion of cold, dense gas into stars and the feedback therefrom, is the largest obstacle hindering progress. An example demonstrating the state--of--the--art in this area is presented.


title: Viscosity Prescriptions in Accretion Disks with Shock Waves

abstract: We study the evolution of viscous, isothermal, rotating, thin, axisymmetric accretion disks around a compact object using Smoothed Particle Hydrodynamics. We emphasize the effects of different choices of viscosity prescriptions on the evolution of angular momentum as well as other physical quantities of the disk. We show that a flow with the Shakura-Sunyaev viscosity prescription may produce only {\it shear} shock, where the angular momentum changes significantly across the shock. We present and study the effects of other prescriptions of viscosity which render the angular momentum continuous across the shock waves. In general, it is observed that for flows with a small viscosity, the shocks are weaker, form farther away and are wider as compared to the shocks in inviscid flows. If the viscosity is high, shocks do not form at all. The flow remains subsonic and Keplerian throughout the disk and becomes supersonic only very close to the horizon.


title: Determining Gravitational Masses of Galaxy Clusters With (1)Virial Equilibrium And (2)Arc-like Images

abstract: The mass derived from gravitational lensing reflects the total mass contained in the lensing system, independent of the specific matter contents and states. A comparison of the dynamical masses from hydrostatic equilibrium with the gravitational masses from arc-like images of background galaxies is made for four clusters of galaxies at intermediate redshits. It is found that virial analysis has underestimated the total cluster masses (from lensing) by a factor of $3\sim6$ within a radius of $\sim0.3$ Mpc $h_{50}^{-1}$ around the cluster centers, indicating that clusters of galaxies might not be regarded as the well relaxed virialized systems. The increase of the total cluster masses obtained from lensing leads to the decrease of the baryon fractions of clusters of galaxies, which provides a crue for solving the ``$\Omega_0$ disprepancy puzzle" in cosmology.


title: VLBI Observations of a Complete Sample of Radio Galaxies. IV-The Radio Galaxies NGC2484, 3C109 and 3C382

abstract: We present here new VLBI observations of one FR-I radio galaxy (NGC2484) and two Broad Line FR-II radio galaxies (3C109 and 3C382). For 3C109 new VLA maps are also shown. These sources belong to a complete sample of radio galaxies under study for a better knowledge of their structures at parsec resolution. The parsec structure of these 3 objects is very similar: asymmetric emission which we interpret as the core plus a one-sided jet. The parsec scale jet is always on the same side of the main kpc-scale jet. The limit on the jet to counterjet brightness ratio, the ratio of the core radio power to the total radio power and the Synchrotron-self Compton model allow us to derive some constraints on the jet velocity and orientation with respect to the line of sight. From these data and from those published on 2 other sources of our sample, we suggest that parsec scale jets are relativistic in both FR-I and FR-II radio galaxies and that pc scale properties in FR-I and FR-II radio galaxies are very similar despite of the large difference between these two classes of radio galaxies on the kpc scale.


title: Can Standard Cosmological Models Explain the Observed Abell Cluster Bulk Flow?

abstract: Lauer \& Postman (LP) observe that all Abell clusters with redshifts less than 15,000\kms\ appear to be participating in a bulk flow of 689 km s$^{-1}$ with respect to the Cosmic Microwave Background. We find this result difficult to reconcile with all popular models for large-scale structure formation that assume Gaussian initial conditions. This conclusion is based on Monte-Carlo realizations of the LP data, drawn from large Particle-Mesh $N$-body simulations. We have taken special care to treat properly the longest-wavelength components of the power spectra. Bulk flows with amplitude as large as that reported by LP are not uncommon in the Monte-Carlo datasets. However, the $\chi^2$ of the observed bulk flow, taking into account the anisotropy of the error ellipsoid, is much more difficult to match in the simulations. The models examined are ruled out at confidence levels between 94\% and 98\%. Any model that has {\it intrinsic} flows of less than 480\kms\ on the scales probed by LP scales can be ruled out at a similar level.


title: Solar Neutrinos: Sensitivity to Pre-Main Sequence Evolution and to the Depth of the Convective Zone

abstract: Predicted rates for solar neutrino experiments that are obtained with a modified stellar evolution code originally developed to study the advanced stages of stellar evolution are shown to be in agreement with other recently-calculated precise solar models to about 2% (0.2 SNU for the chlorine experiment). Different scenarios for pre-main sequence evolution are explored and are found to change the predicted rates for solar neutrino experiments by less than or of order 1%. The influence of the depth of the solar convective zone on the predicted solar neutrino fluxes is established by direct calculation. It is shown that a change in the calculated depth of the convective zone that is five times larger than the quoted helioseismological measurement uncertainty determined by Christensen-Dalsgaard, Gough, & Thompson causes a change in the predicted $^8$B neutrino flux of less than 7% and a change in the $^7$Be neutrino flux of less than 4%. Additionally, it is shown that the radiative opacities near the depth of the convective zone cannot differ from the standard OPAL values by more than about 12% without causing the calculated and measured depths of the convective zone to differ by more than five times the helioseismological measurement uncertainty.


title: Interpreting the candidate Galactic microlensing events

abstract: Four microlensing collaborations are presently searching for compact matter in the Galaxy and all have detected possible candidates. Using the detection efficiencies recently published by the MACHO and OGLE collaborations, we present Monte-Carlo calculations of the expected optical depth, rates and timescales, along the LMC and Galactic bulge lines of sight, for dark matter in a four-component `standard Galaxy' model with a spherically-symmetric halo and spheroid. Using the typically observed event durations we show that, whilst the halo fraction comprised of compact matter is likely to be $f_{h} < 0.4$, a `no halo compact matter' hypothesis is ruled out at greater than the 80\% confidence level, unless the LMC itself has a substantial halo of such objects. On the basis of the timescales observed by OGLE towards the bulge we find the rate predicted by the model to be in good agreement with the number of OGLE detections. We compute lens mass probability distributions for the various components and compare these estimates with current observational and theoretical constraints on the mass scale of baryonic dark matter. We assess the uniformity of the amplification distributions for the published EROS, MACHO and OGLE events and find that they are quite consistent with the microlensing hypothesis, although the OGLE candidate selection criteria mean that its data are particularly sensitive to photometric selection effects. The EROS team has recently placed strong limits on the density contribution of very low mass halo objects from their short timescale CCD search. On the basis of a recent study of the flux amplification of finite-size sources by Simmons, Newsam \& Willis \shortcite{simm95} we suggest that EROS may have detected up to 5 low-amplification events due to halo lenses with mass $m\sim 10^{-7}~\sm$.


title: Simulations of Relativistic Extragalactic Jets

abstract: We describe a method for the numerical solution of the relativistic Euler equations which we have found to be both robust and efficient, and which has enabled us to simulate relativistic jets. The technique employs a solver of the Godunov-type, with approximate solution of the local Riemann problems, applied to laboratory frame variables. Lorentz transformations provide the rest frame quantities needed for the estimation of wave speeds, etc. This is applied within the framework of an adaptive mesh refinement algorithm, allowing us to perform high-resolution, 2-D simulations with modest computing resources. We present the results of nonrelativistic, and relativistic ($\gamma=5$ and $10$) runs, for adiabatic indices of $5/3$ and $4/3$. We find the same gross morphology in all cases, but the relativistic runs exhibit little instability and less well-defined structure internal to the jet: this might explain the difference between (relatively slow) BL~Lacs and (faster) QSOs. We find that the choice of adiabatic index makes a small but discernible difference to the structure of the shocked jet and ambient media.


title: Femtolensing: Beyond the Semi-Classical Approximation

abstract: Femtolensing is a gravitational lensing effect in which the magnification is a function not only of the positions and sizes of the source and lens, but also of the wavelength of light. Femtolensing is the only known effect of $(10^{-13}-10^{-16} M_{\sun})$ dark-matter objects and may possibly be detectable in cosmological gamma-ray burst spectra. We present a new and efficient algorithm for femtolensing calculations in general potentials. The physical-optics results presented here differ at low frequencies from the semi-classical approximation, in which the flux is attributed to a finite number of mutually coherent images. At higher frequencies, our results agree well with the semi-classical predictions. Applying our method to a point-mass lens with external shear, we find complex events that have structure at both large and small spectral resolution. In this way, we show that femtolensing may be observable for lenses up to $10^{-11}$ solar masses, much larger than previously believed. Additionally, we discuss the possibility of a search for femtolensing of white dwarfs in the LMC at optical wavelengths.


title: The Galaxy Group/Cosmology Connections

abstract: \def\spose#1{\hbox to 0pt{#1\hss}} \def\lta{\mathrel{\spose{\lower 3pt\hbox{$\mathchar"218$}} \raise 2.0pt\hbox{$\mathchar"13C$}}} Groups of galaxies are highly linked to cosmology: 1) groups are tidally destroyed by the tidal field of the cluster they fall into; 2) spherical infall leads to the young cosmo-dynamical state of loose groups, a fundamental surface for groups, $\Omega_0 \simeq 0.3$, and the mixed nature of compact groups (virialized groups, groups at full collapse, and chance alignments within collapsing loose groups, for decreasing compact group velocity dispersion) ; 3) X-ray analyses lead to $\Omega_0 \lta 0.5$.


title: Star Forming Galaxy Models: Blending Star Formation into TREESPH

abstract: We have incorporated star formation algorithms into a hybrid N-body/smoothed particle hydrodynamics code (TREESPH) in order to describe the star forming properties of disk galaxies over timescales of a few billion years. The models employ a Schmidt law of index $n\sim 1.5$ to calculate star formation rates, and explicitly include the energy and metallicity feedback into the ISM. Modeling the newly formed stellar population is achieved through the use of hybrid SPH/young star particles which gradually convert from gaseous to collisionless particles, avoiding the computational difficulties involved in creating new particles. The models are shown to reproduce well the star forming properties of disk galaxies, such as the morphology, rate of star formation, and evolution of the global star formation rate and disk gas content. As an example of the technique, we model an encounter between a disk galaxy and a small companion which gives rise to a ring galaxy reminiscent of the Cartwheel (AM 0035-35). The primary galaxy in this encounter experiences two phases of star forming activity: an initial period during the expansion of the ring, and a delayed phase as shocked material in the ring falls back into the central regions.


title: Stochastic Inflation: New Results

abstract: This is an old article which has never been posted for public use and which only appeard in a relatively hard-to-get Proceedings of the Sakharov Conference in Moscow (May, 1991). The subject of this article has received a lot of attention in last two years. Although we acknowledge that there has been a very big progress in the field of relating the spectrum of fluctuations to scalar field potentials of inflation models, which have outdated some of our results (particularly because we never took into account the tensor fluctuations which have been shown to be crucial for reconstructing the inflaton potential) we think that some of the other ideas and approaches which we presented in this article maybe still of interest. We did not rewrite the article to include any of the newer approaches and we did not add any of the newer references. The reader is advised to refer to recent papers (for example, see the paper astro-ph/9403001 by Edward Kolb and Sharon Vadas) for state of the art review and up-to-date references. ==The original abstract of the article==: We prove that, in stochastic approach, there exists an equivalence relation between different inflationary models under some redefinition of field and time variables. The postinflationary physics is insensitive to it and one can say that related theories are in fact indistinguishable from the viewpoint


title: Polarization and Anisotropy of the Microwave Sky

abstract: We study the polarization-polarization and polarization-temperature correlations in standard adiabatic scenarios for structure formation. Temperature anisotropies due to gravitational potential wells and oscillations in the photon-baryon-electron fluid on the surface of last scattering are each associated with a correlated polarization pattern. While the `correlated part' of the polarization has an r.m.s. of only a third of the total signal, it may still be measurable by mapping a large area on the sky. We calculate the expected signal to noise ratio for various measures of the polarization in a hypothetical mapping experiment such as those now being planned.


title: Statistical Constraints on the Inflation Effective Potential from the COBE DMR Results

abstract: We explore constraints on various forms for the effective potential during inflation based upon a statistical comparison between inflation-generated fluctuations in the cosmic microwave background temperature and the COBE DMR results. Fits to the first year 53 A + B x 90 A + B cross correlation function using an effective potential of the form V(phi) = lambda*phi^x/x!, yield upper limits of x le 107 and x le 397 at the 1 sigma and 2 sigma confidence levels, while the same analysis produces a power law of index of n = 1.10 ^{+0.47}_{-0.64} which would imply an underestimate of x le 63 from the simple analytic relation between x and n . We also quantify new limits on the parameters for polynomial effective potentials. This work highlights the importance of a careful statistical treatment when seeking constraints on the inflation-generating effective potential.


title: Lyman alpha absorption at low redshifts and hot gas in galactic haloes

abstract: Motivated by recent observation of Lanzetta et al. that most luminous galaxies at low redshifts produce $\lya$ absorption at impact parameter $l\lsim 160 \kpch$, we propose that these absorbers are clouds confined by the pressure of ambient hot gas in galactic haloes. We determine the properties of this hot gas and of the absorption systems on the basis of observational and theoretical constraints. The absorbing clouds need to be replenished on about one orbital time ($\sim 10^9$ yrs) in the galactic halo. The pressure and temperature of the gas at radius $r\sim 100\kpc$ are $P=(10-100){\rm cm^{-3} K}$, $T=10^{(5.5-6.5)}{\rm K}$. The model requires that about 10 per cent of the gas in low-redshift galactic haloes is in the hot phase. Such gas in galactic haloes emits x-ray with bolometric luminosity of the order $10^{37-40}{\rm erg\,s^{-1}}$. The plausibility for such gas to exist in current models of galaxy formation is discussed.


title: The Optical Redshift Survey: Sample Selection and the Galaxy Distribution

abstract: This is the first in a series of papers describing the {\it Optical Redshift Survey} (ORS), a redshift survey of optically selected galaxies covering 98\% of the sky above $|b| = 20^\circ$ (8.09 ster). The survey is drawn from the UGC, ESO, and ESGC galaxy catalogues, and contains two sub-samples, one complete to a $B$ magnitude of 14.5, the other complete to a $B$ major axis diameter of $1.9^\prime$. The entire sample consists of 8457 objects, of which redshifts are now available for 8286; 171 objects remain without measured redshifts. The ORS provides the most detailed and homogeneous sampling of the large-scale galaxy distribution to date in these areas. The density field of bright optical galaxies is well-defined to $8000 \kms$, and is dominated by the Virgo, Telescopium-Pavo-Indus, Hydra-Centaurus, Pisces-Perseus, and Coma-A1367 Superclusters. The dense sampling provided by ORS allows a detailed analysis of the galaxy density field, and will be used to test its dependence on morphology and other galaxy parameters.


title: A Two-Fluid Approximation for Calculating the Cosmic Microwave Background Anisotropies

abstract: We present a simple, yet accurate approximation for calculating the cosmic microwave background anisotropy power spectrum in adiabatic models. It consists of solving for the evolution of a two-fluid model until the epoch of recombination and then integrating over the sources to obtain the CMB anisotropy power spectrum. The approximation is useful both for a physical understanding of CMB anisotropies, as well as for a quantitative analysis of cosmological models. Comparison with exact calculations shows that the accuracy is typically better than 20 percent over a large range of angles and cosmological models, including those with curvature and cosmological constant. Using this approximation we investigate the dependence of the CMB anisotropies on the cosmological parameters. We identify six dimensionless parameters that uniquely determine the anisotropy power spectrum within our approximation. CMB experiments on different angular scales could in principle provide information on all these parameters. In particular, mapping of the Doppler peaks would allow an independent determination of baryon mass density, matter mass density and Hubble constant.


title: The gravitational deflection of light in MOND

abstract: The deflection angle $\Delta\phi$ of light rays by the gravitational field of a spherical system $M(r)$ is calculated using the MOdified Newtonian Dynamics (MOND). It is shown that $\Delta\phi$ with an impact parameter $r_0$ can be expressed by the measured rotation velocity $v(r)$ as $\Delta\phi(r_0)=2\int_{r_0}^{\infty} \frac{v^2(r)}{c^2}\frac{r_0dr}{r\sqrt{r^2-r_0^2}}$, where $v(r)=\left\{\left(Ga_0M(r)\right)^{1/4}, & r_0>r_c;\left(\frac{GM(r)}{r}\right)^{1/2},& r_0\leq r_c,\right$, and $r_c$ is the critical radius that is determined by the critical acceleration $a_0$. In the Newtonian limit of the gravitational acceleration $a\gg a_0$, $\Delta\phi$ approaches $\Delta\phi=2Gm(r_0)/c^2r_0$ with the projected surface mass $m(r_0)$. Whilst the asymptotic value of $\Delta\phi$ reaches a constant $\pi(v_{\infty}/c)^2$ in the low-acceleration limit of $a\ll a_0$. Taking the empirical correction of a factor of 2 from the theory of general relativity into account and utilizing the relation between rotation velocity $v$ and velocity dispersion $\sigma$, MOND results naturally in a constant deflection angle, $4\pi(\sigma/c)^2$, which has been widely used in the present-day study of gravitational lensing by galaxies and clusters of galaxies, implying that without introducing the massive halos acting as $r^{-2}$ for dark matter MOND has no difficulty in reproducing the known cases of gravitational lensing associated with galaxies and clusters of galaxies.


title: Mapping Cluster Mass Distributions via Gravitational Lensing of Background Galaxies

abstract: We present a new method for measuring the projected mass distributions of galaxy clusters. The gravitational amplification is measured by comparing the joint distribution in redshift and magnitude of galaxies behind the cluster with that of field galaxies. We show that the total amplification is directly related to the surface mass density in the weak field limit, and so it is possible to map the mass distribution of the cluster. The method is shown to be limited by discreteness noise and galaxy clustering behind the lens. Galaxy clustering sets a lower limit to the error along the redshift direction, but a clustering independent lensing signature may be obtained from the magnitude distribution at fixed redshift. Statistical techniques are developed for estimating the surface mass density of the cluster. We extend these methods to account for any obscuration by cluster halo dust, which may be mapped independently of the dark matter. We apply the method to a series of numerical simulations and show the feasibility of the approach. We consider approximate redshift information, and show how the mass estimates are degraded.


title: Variable Smoothing Lengths and Energy Conservation in Smoothed Particle Hydrodynamics

abstract: We present a new formulation of the equations of motion used in smoothed particle hydrodynamics (SPH). The spatial resolution in SPH is determined by the smoothing length, $h$, and it has become common practice for each particle to be given its own adaptive smoothing length, $h_i$. Consequently, the dynamic range that may be spatially resolved is greatly increased, but additional ($\nabla h$) terms, which account for the variability of the smoothing lengths, should appear in the equations of motion in order to satisy conservation requirements. Previous implementations of SPH have neglected these additional terms, whereas we have included them. This is achieved by defining a functional form for the $h_i$s, that depends on inter-particle distances only, and then deriving the equations of motion using a Hamiltonian formalism. A number of test calculations are presented. We find that including the $\nabla h$ terms has no detrimental effect on the ability of SPH to model known problems with reasonable accuracy. For problems where the energy conservation was rather poor, including the $\nabla h$ terms results in a dramatic improvement. In particluar non-conservation of energy during a collision between two polytropes can occur at the level of 10% when the $\nabla h$ terms are neglected. When the $\nabla h$ terms are included, then this error reduces to 0.8%.


title: Circumstellar Interaction in SN 1993J

abstract: The radio and X-ray observations of SN 1993J during the first year can be consistently explained as a result of interaction of the expanding ejecta with a circumstellar medium. The density of the circumstellar gas can be deduced from the free-free absorption of the radio emission and from the X-ray luminosity. During the first two weeks, both sets of observations indicate a mass loss rate of $\sim 4\EE{-5} \Msun ~\rm yr^{-1}$ for a wind velocity of $10 \kms$. The subsequent radio and X-ray observations indicate a density gradient $\propto r^{-s}$, with $1.5 \lsim s \lsim 1.7$, as opposed to the $r^{-2}$ gradient expected for a steady, spherically symmetric wind. This may either be caused by a variation of the mass loss rate from the progenitor system, or by a non-spherically symmetric geometry. To explain the properties of the X-ray emission, a steep density gradient in the ejecta is needed. During the first months most of the observed X-ray emission originates from the circumstellar shock, which is adiabatic, while the reverse shock is radiative. To avoid excessive Comptonization in the X-ray range collisionless heating must be ineffective. The soft X-rays observed at 220 days probably originate from the reverse shock. The ionization and temperature structures of the circumstellar gas are calculated; we find that the temperature is in excess of $10^5 \rm~K$ and the medium is nearly completely ionized by the shock radiation after the formation of the shocks. Preacceleration of the circumstellar gas by the radiation from the outbreak can explain


title: Relativistic Radiative Transfer for Spherical Flows

abstract: We present a new complete set of Lagrangian relativistic hydrodynamical equations describing the transfer of energy and momentum between a standard fluid and a radiation fluid in a general non-stationary spherical flow. The new set of equations has been derived for a particular application to the study of the cosmological Quark--Hadron transition but can also be used in other contexts.


title: Phenomenological Aspects of the Pre-Big-Bang Scenario in String Cosmology

abstract: I review various aspects of the pre-big-bang scenario and of its main open problems, with emphasis on the role played by the dilaton. Since the dilaton is a compelling consequence of string theory, tests of this scenario are direct tests of string theory and also, more generally, of Planck scale physics.


title: Ubc/Laval 2.7-Meter Liquid Mirror Telescope

abstract: We report the successful development and operation of a large astronomical liquid-mirror telescope. Employing a rotating 2.7-meter diameter mirror with a surface of liquid mercury, the telescope images a half-degree diameter field centered at the zenith. Located near Vancouver, British Columbia, it is equipped with a low-noise 2048x2048-pixel CCD detector, operating in TDI mode, which produces continuous imaging of a 20 arcminute-widestrip of sky with 2 minute integration time. Images with FWHM of 2 arcseconds or less are regularly obtained. This image quality is limited only by atmospheric seeing and star-trail curvature. The telescope is equipped with a series of narrow-band filters, designed to produce 40-point spectral energy distributions from 0.4 to 1.0 microns of all detected objects. These will allow classification and redshift estimation of approximately 25,000 galaxies and 1,000 quasars to a limiting magnitude of R~21.


title: Dwarf Spheroidals and the Evolution of Not-Too-Old Population II Stars

abstract: Previous evolutionary computations concerning stars with masses in the range of the ``Red Giants Transition Phase'' are extended to lower metallicities. Clusters isochrones for two values of the metallicity Z (Z=10E-4, 4E-4) and for ages down to 1 billion years are presented, discussing selected evolutionary features characterizing the structure of the stars during both H or He burning phases. One finds that clusters with Z about equal to 1E-4 an age around 1-2 billion years could be characterized by the occurrence of anomalous He burning variables at a luminosity of the order of LogL= 2.0-2.2 (in solar units).


title: Analytic Error Estimates

abstract: I present an analytic method for estimating the errors in fitting a distribution. A well-known theorem from statistics gives the minimum variance bound (MVB) for the uncertainty in estimating a set of parameters $\l_i$, when a distribution function $F(z;\l_1 ... \l_m)$ is fit to $N$ observations of the quantity(ies) $z$. For example, a power-law distribution (of two parameters $A$ and $\gaml$) is $F(z;A,\gaml) = A z^{-\gaml}$. I present the MVB in a form which is suitable for estimating uncertainties in problems of astrophysical interest. For many distributions, such as a power-law distribution or an exponential distribution in the presence of a constant background, the MVB can be evaluated in closed form. I give analytic estimates for the variances in several astrophysical problems including the gallium solar-neutrino experiments and the measurement of the polarization induced by a weak gravitational lens. I show that it is possible to make significant improvements in the accuracy of these experiments by making simple adjustments in how they are carried out or analyzed. The actual variance may be above the MVB because of the form of the distribution function and/or the number of observations. I present simple methods for recognizing when this occurs and for obtaining a more accurate estimate of the variance than the MVB when it does.


title: Quoting Experimental Information

abstract: There was a question raised at this meeting as to the best way to present the results of experiments measuring anisotropies in the CMB. Here we will make some simple comments about the 3 main competing methods and some suggestions. In particular we will give an over-simplified but hopefully useful method of comparing GACF numbers with other power spectra, and a simple rule-of-thumb for converting between the two main conventions.


title: Gravitational Wave Signatures of Superdense Objects in a Nonsingular Solution of Gravity Theory

abstract: A massive gravitationally bound object with a radius $r \leq 2GM/c^2$, which occurs in the non-singular version of the nonsymmetric gravitational theory (NGT), replaces the black hole in Einstein gravity theory. This object is kept stable by the attractive and repulsive forces generated by NGT, as well as standard matter pressures, and is called a superdense object (SDO). The luminosity of gravitational waves emitted by a SDO with a red-shift of order unity is calculated and it is found that it could be a strong source of gravitational radiation at low frequencies. The active galactic nucleus of M87 is identified with a SDO and the recent observational results obtained by the Hubble Space Telesope are used to estimate the amplitude of gravitational radiation.


title: Blue Stragglers in Galactic Globular Clusters: Playing with Specific Quantities

abstract: We present preliminary results obtained from the comparison of the specific frequencies of Blue Straggler Stars (BSS) detected so far in a sample of 26 Galactic globular clusters. The number of BSS seems to increase almost linearly with increasing the amount of sampled light in loose clusters, while it drops abruptly for clusters having intermediate-high central densities. In particular, a simple interpretative scenario where the BSS in loose clusters are produced from primordial binaries and those detected in high density globulars are due to star interactions leading to binary formation, merging. etc. seems compatible with these early results. The possibility that this observational evidence could be ascribed to systematic biases mostly related to the increasing difficulty to detect BSS candidates with increasing cluster concentration, is also discussed.


title: Carbon Abundances in F and G Dwarfs

abstract: We have determined carbon abundances or upper limits from the forbidden [C I] line at 8727.13 A for 85 F and early G type main sequence stars, with metallicities ranging between [Fe/H] = -1.0 and [Fe/H] = +0.25. The [C I] line has not been previously used for a study of this size. We find that the C/Fe abundance ratio is slowly decreasing with time and increasing metallicity in the disk, but with considerable observational scatter. A discussion of results in the literature supports this result. Our data do not reveal any gradients with galactocentric birth distance for our stars, but more accurate observations are desirable since possible gradients should be important to the question of the synthesis of C and Fe in the disk. (Accepted for publication in A&A. Postscript text and figures available via anonymous ftp to "chiron.astro.uu.se", in directory "pub/articles/atmos/P89")


title: Relativistic Jets and the Fanaroff-Riley Classification of Radio Galaxies

abstract: Recently, Owen(1993) and Owen and Ledlow (1994) have shown that the dividing line between Fanaroff-Riley class~I and class~II Radio Sources is very sharp when the sources are represented as points in the radio-optical luminosity plane. It is shown that if one accepts the propositions that the sources in the vicinity of this dividing line are characterized by (1) transonic Mach numbers and (2) mildly relativistic velocities (consistent with deceleration of an initially moderately relativistic or ultrarelativistic jet), then the slope of the dividing line is readily obtained using simple physical arguments and established empirical relationships between the X-Ray luminosities, core radii, velocity dispersions and absolute magnitudes of elliptical galaxies. The intercept of the dividing line depends upon parameters which are known perhaps to within factors of order unity and agrees with the data to within an order of magnitude. ROSAT observations of elliptical galaxies will be important in constraining the central pressures and X-Ray core radii of radio ellipticals. Knowledge of these two parameters will assist in a more detailed assessment of the physics which is proposed here as being relevant. High resolution observations of jets within a kpc of the radio core will also be useful for determining the spreading rates of jets in this region.


title: Capanic: A Parallel Tree N-Body Code for Inhomogeneous Clusters of Processors

abstract: We have implemented a parallel version of the Barnes-Hut 3-D N-body tree algorithm under PVM 3.2.5, adopting an SPMD paradigm. We parallelize the problem by decomposing the physical domain by means of the {\bf Orthogonal Recursive Bisection} oct-tree scheme suggested by Salmon (1991), but we modify the original hypercube communication pattern into an incomplete hypercube, which is more suitable for a generic inhomogenous cluster architecture.\\ We address dynamical load balancing by assigning different "weights" to the spawned tasks according to the dynamically changing workloads of each task. The weights are determined by monitoring the local platforms where the tasks are running and estimating the performance of each task. The monitoring scheme is flexible and allows us to address at the same time cluster and intrinsic sources of load imbalance. We then show measurements of the performance of our code on a test case of astrophysical interest in order to test the performance of our implementation.


title: The Dynamical Survival of Small-Scale Substructure in Relaxed Galaxy Clusters

abstract: We consider the dynamical evolution of small-scale substructure in clusters within two extreme alternate scenarios for their possible origin: 1) the accretion of groups (or small clusters) on quasi-radial orbits, and 2) the merger of clusters of similar masses, followed by the decoupling of their dense cores. Using simple analytical arguments and checking with numerical computations, we show that objects are destroyed by the tidal field of the global cluster potential if their mean density is small compared to the mean cluster density within the radius of closest approach of the group or detached core. Accreted groups and small clusters are thus tidally disrupted in one cluster crossing. Since the cores of clusters are much denser than groups, they are considerably more robust to tides, but the least massive are destroyed or severely stripped by tides, while the others are brought to the cluster center by dynamical friction (and subsequently merge) in less than one orbit. The longest lived substructures are detached cores, roughly ten times less massive than the cluster, starting in near-circular orbits beyond $1 \, h^{-1} \, \rm Mpc$ from the cluster center.


title: Beyond the Small-Angle Approximation For MBR Anisotropy from Seeds

abstract: In this paper we give a general expression for the energy shift of massless particles travelling through the gravitational field of an arbitrary matter distribution as calculated in the weak field limit in an asymptotically flat space-time. It is {\it not} assumed that matter is non-relativistic. We demonstrate the surprising result that if the matter is illuminated by a uniform brightness background that the brightness pattern observed at a given point in space-time (modulo a term dependent on the oberver's velocity) depends only on the matter distribution on the observer's past light-cone. These results apply directly to the cosmological MBR anisotropy pattern generated in the immediate vicinity of of an object like a cosmic string or global texture. We apply these results to cosmic strings, finding a correction to previously published results for in the small-angle approximation. We also derive the full-sky anisotropy pattern of a collapsing texture knot.


title: Dynamics of Silent Universes

abstract: We investigate the local non--linear dynamics of irrotational dust with vanishing magnetic part of the Weyl tensor, $H_{ab}$. Once coded in the initial conditions, this dynamical restriction is respected by the relativistic evolution equations. Thus, the outcome of the latter are {\it exact solutions} for special initial conditions with $H_{ab}=0$, but with no symmetries: they describe inhomogeneous triaxial dynamics generalizing that of a fluid element in a Tolman--Bondi, Kantowski--Sachs or Szekeres geometry. A subset of these solutions may be seen as (special) perturbations of Friedmann models, in the sense that there are trajectories in phase--space that pass arbitrarily close to the isotropic ones. We find that the final fate of ever--expanding configurations is a spherical void, locally corresponding to a Milne universe. For collapsing configurations we find a whole family of triaxial attractors, with vanishing local density parameter $\Omega$. These attractors locally correspond to Kasner vacuum solutions: there is a single physical configuration collapsing to a degenerate {\it pancake}, while the generic configuration collapses to a triaxial {\it spindle} singularity. These {\it silent universe} models may provide a fair representation of the universe on super horizon scales. Moreover, one might conjecture that the non--local information carried by $H_{ab}$ becomes negligible during the late highly non--linear stages of collapse, so that the attractors we find may give all of the relevant expansion or collapse configurations of irrotational dust.


title: CBR Anisotropy in an Open Inflation, CDM Cosmogony

abstract: We compute the cosmic background radiation anisotropy, produced by energy-density fluctuations generated during an early epoch of inflation, in an open cosmological model based on the cold dark matter scenario. At $\Omega_0 \sim$ 0.3 -- 0.4, the COBE normalized open model appears to be consistent with most observations.


title: Reverberation mapping of active galactic nuclei : The SOLA method for time-series inversion

abstract: In this paper a new method is presented to find the transfer function of the broad-line region in active galactic nuclei. The subtractive optimally localized averages (SOLA) method is a modified version of the Backus-Gilbert method and is presented as an alternative to the more often used maximum-entropy method. The SOLA method has been developed for use in helioseismology. It has been applied to the solar oscillation frequency splitting data currently available to deduce the internal rotation rate of the sun. The original SOLA method is reformulated in the present paper to cope with the slightly different problem of inverting time series. We use simulations to test the viability of the method and apply the SOLA method to the real data of the Seyfert-1 galaxy NGC 5548. We investigate the effects of measurement errors and how the resolution of the TF critically depends upon both the sampling rate and the photometric accuracy of the data. A uuencoded compressed postscript file of the paper which includes the figures is available by anonymous ftp at ftp://solaris.astro.uu.se/pub/articles/atmos/frank/PijWan.uue


title: The Nature versus Nurture of Anisotropies

abstract: With the rapidly growing number of cosmic microwave background measurements on various scales, there is real hope that the number of acceptable models for structure formation will be limited to a very few in the near future. Yet any given model can always be saved by introducing and tuning extraneous free parameters. To better understand this question of ``nature versus nurture'' for temperature fluctuations, it is useful to know not only the general features of anisotropy predictions but also their causes. Extracting the physical content of our other works, we present here a {\it simple} account of cosmic microwave background anisotropies on all scales. In particular, we show that analytic approximations can trace the structure of the so-called ``Doppler peaks,'' which arise due to the {\it adiabatic} oscillations in the photon-baryon fluid. We also show how the finite thickness of the last scattering surface and the Silk damping mechanism can be described in a unified way by photon diffusion. In order to present a specific example, we focus on comparing the primordial isocurvature baryon (PIB) model with the standard cold dark matter model (CDM). In particular, we explain why PIB generically predicts larger {\it non}-oscillatory anisotropies from the 1$^\circ$ to 10$^\circ$ scale which may already be in conflict with experiments.


title: Neutrino-Neutrino Scattering and Matter-Enhanced Neutrino Flavor Transformation in Supernovae

abstract: We examine matter-enhanced neutrino flavor transformation ($\nu_{\tau(\mu)}\rightleftharpoons\nu_e$) in the region above the neutrino sphere in Type II supernovae. Our treatment explicitly includes contributions to the neutrino-propagation Hamiltonian from neutrino-neutrino forward scattering. A proper inclusion of these contributions shows that they have a completely negligible effect on the range of $\nu_e$-$\nu_{\tau(\mu)}$ vacuum mass-squared difference, $\delta m^2$, and vacuum mixing angle, $\theta$, or equivalently $\sin^22\theta$, required for enhanced supernova shock re-heating. When neutrino background effects are included, we find that $r$-process nucleosynthesis from neutrino-heated supernova ejecta remains a sensitive probe of the mixing between a light $\nu_e$ and a $\nu_{\tau(\mu)}$ with a cosmologically significant mass. Neutrino-neutrino scattering contributions are found to have a generally small effect on the $(\delta m^2,\ \sin^22\theta)$ parameter region probed by $r$-process nucleosynthesis. We point out that the nonlinear effects of the neutrino background extend the range of sensitivity of $r$-process nucleosynthesis to smaller values of $\delta m^2$.


title: Comment on The Preprint ``Neutrino Flavor Evolution Near A Supernova's Core'' (Astro-Ph 9405008) by J. Pantaleone

abstract: The revised version of the widely circulated preprint ``Neutrino Flavor Evolution Near A Supernova's Core'' by J. Pantaleone (astro-ph 9405008 on the Bulletin Board, Indiana University preprint IUHET-276) is wrong. It contains two errors which lead to incorrect conclusions regarding neutrino flavor transformation in the supernova environment. In this short note we discuss these errors.