title: SN~1991T: Reflections of Past Glory

abstract: We have obtained photometry and spectra of SN~1991T which extend more than 1000 days past maximum light, by far the longest a SN~Ia has been followed. Although SN~1991T exhibited nearly normal photometric behavior in the first 400 days following maximum, by 600 days its decline had slowed, and by 950~days the supernova brightness was consistent with a constant apparent magnitude of $m_B=21.30$. Spectra near maximum showed minor variations on the SN~Ia theme which grew less conspicuous during the exponential decline. At 270 days the nebular spectrum was composed of Fe and Co lines common to SNe~Ia. However, by 750 days past maximum light, these lines had shifted in wavelength, and were superimposed on a strong blue continuum. The luminosity of SN~1991T at 950 days is more than $9.0\times10^{38}(D/13~{\rm Mpc})^2$~ergs~s$^{-1}$ with a rate of decline of less than $0.04$ mags/100~days. We show that this emission is likely to be light that was emitted by SN~1991T near maximum light which has reflected from foreground dust, much like the light echos observed around SN~1987A.


title: The Distances to Five Type~II Supernovae Using the Expanding Photosphere Method and the Value of $H_0$

abstract: We have used observations gathered at CTIO to measure distances by the Expanding Photosphere Method (EPM) to 5 Type II supernovae. These supernovae lie at redshifts from cz = 1100 km/s to z = 5500 km/s, and increase to 18 the number of distances measured using EPM. We compare distances derived to 11 Type II supernovae with distances to their host galaxies measured using the Tully-Fisher method. We find that the Tully-Fisher distances average 11%+-7% smaller. The comparison shows no significant evidence of any large distance-dependent bias in the Tully-Fisher distances. We employ the sample of EPM distances from 4.5~Mpc to 180~Mpc to derive a value for the Hubble Constant. We find H_0 = 73 +- 6 (statistical) +- 7 (systematic}) km/s/Mpc.


title: Implications of a High Population II B/Be Ratio

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


title: Estimation of Peculiar Velocity from the Inverse Tully-Fisher Relation

abstract: We present a method for deriving a smoothed estimate of the peculiar velocity field of a set of galaxies with measured circular velocities $\eta\equiv {\rm log} \Delta v$ and apparent magnitudes $m$. The method is based on minimizing the scatter of a linear inverse Tully-Fisher relation $\eta= \eta(M)$ where the absolute magnitude of each galaxy is inferred from its redshift $z$, corrected by a peculiar velocity field, $M \propto m - 5\log(z-u)$. We describe the radial peculiar velocity field $u({\bf z})$ in terms of a set of orthogonal functions which can be derived from any convenient basis set; as an example we take them to be linear combinations of low order spherical harmonic and spherical Bessel functions. The model parameters are then found by maximizing the likelihood function for measuring a set of observed $\eta$. The predicted peculiar velocities are free of Malmquist bias in the absence of multi-streaming, provided no selection criteria are imposed on the measurement of circular velocities. This procedure can be considered as a generalized smoothing algorithm of the peculiar velocity field, and is particularly useful for comparison to the smoothed gravity field derived from full-sky galaxy redshift catalogs such as the IRAS surveys. We demonstrate the technique using a catalog of ``galaxies" derived from an N-body simulation. Increasing the resolution of the velocity smoothing beyond a certain level degrades the correlation of fitted velocities against the velocities calculated from linear theory methods, which have finite resolution,


title: Polarization-Temperature Correlations in the Microwave Background

abstract: The polarization of the CMB offers a potentially valuable probe of the physics of the early universe, specifically of the velocity of the photon-baryon fluid on the last scattering surface. Here we discuss qualitatively the correlations that arise between the temperature anisotropies and polarization of the CMB and how they might be observed.


title: Dark Matter in the Light of COBE

abstract: The observations of all three COBE instruments are examined for the effects of dark matter. The anisotropy measured by the DMR, and especially the degree-scale ground- and balloon-based experiments, is only compatible with large-scale structure formation by gravity if the Universe is dominated by non-baryonic dark matter. The FIRAS instrument measures the total power radiated by cold dust, and thus places tight limits on the absorption of starlight by very cold gas and dust in the outer Milky Way. The DIRBE instrument measures the infrared background, and will place tight limits on the emission by low mass stars in the Galactic halo.


title: Comparison of Spectral Index Determinations

abstract: The index $n$ of a power law power spectrum of primordial density fluctuations, $P(k) \propto k^n$, has been estimated using many different techniques. The most precise compare the COBE DMR large angular scale $\Delta T$ to the amplitude of the large scale structure, but these are also the most model-dependent. The COBE DMR $\Delta T$ has also been compared to the degree-scale $\Delta T$ from several experiments. And finally, a relatively model-independent value of $n$ can derived from the COBE data alone, but the small range of angular scales covered by COBE limits the precision of these methods.


title: Self-Lensing By A Stellar Disk

abstract: I derive a general expression for the optical depth $\tau$ for gravitational lensing of stars in a disk by Massive Compact Objects (Machos) in the same disk. For the more restricted case where the disk is self-gravitating and the stars and Machos have the same distribution function, I find $\tau = 2\VEV{v^2}/c^2\sec^2 i$ where $\VEV{v^2}$ is the mass-weighted vertical velocity dispersion, and $i$ is the angle of inclination. This result does not depend on any assumptions about the velocity distribution. As an example, if stars within the bar of the Large Magellanic Cloud (LMC) account for the observed optical depth $\tau\sim 8\times 10^{-8}$ as has recently been suggested, then $v\gsim 60\,\kms$. This is substantially larger than the measured dispersions of known LMC populations.


title: The Difference between Radio-Loud and Radio-Quiet Active Galaxies

abstract: The recent development of unified theories of active galactic nuclei (AGN) has indicated that there are two physically distinct classes of these objects - radio-loud and radio-quiet. The primary observational distinctions between the two types are: (1) The radio-loud objects produce large scale radio jets and lobes, with the kinetic power of the jets being a significant fraction of the total bolometric luminosity. On the other hand, the weak radio ejecta of the radio-quiet objects are energetically insignificant. (2) The radio-loud objects are associated with elliptical galaxies which have undergone recent mergers, while the radio-quiets prefer spiral hosts. (3) The space density of the radio-louds at a given optical luminosity is $\approx$ 10 times lower than that of the radio-quiets. Despite these differences, the (probably) thermal emissions from the AGN (continua and lines from X-ray to infrared wavelengths) are quite similar in the two classes of object. We argue that this last result suggests that the black hole masses and mass accretion rates in the two classes are not greatly different, and that the difference between the classes is associated with the spin of the black hole. We assume that the normal process of accretion through a disk does not lead to rapidly spinning holes, and propose instead that galaxies (e.g. spirals) which have not suffered a recent major merger event contain non-rotating or only slowly rotating black holes. When two such galaxies merge, the two black holes are known to form a binary and we assume that they eventually coalesce. In the small fraction of mergers in which the two ``parent'' galaxies contain very massive holes of roughly equal mass, a rapidly spinning, very massive hole results. It is proposed that such mergers are the progenitors of powerful radio


title: Response to Comments on "Neutrino Flavor Evolution Near a Supernova's Core"

abstract: In a recent paper I concluded that the neutrino background is important for calculating neutrino flavor evolution in supernovae, and including it weakens the connection between r-process nucleosynthesis and cosmologically relevant neutrino masses. Comments on this paper have been posted to this bulletin board; one by J. Cline and the other by Y.Z. Qian and G.M. Fuller. Here I briefly discuss these Comments and demonstrate that they are inaccurate and ill-considered. No corrections to my paper are necessary.


title: Light Element Abundance Inhomogeneities in Globular Clusters: Probing Star Formation and Evolution in the Early Milky Way

abstract: Abundance patterns of the elements C, N, and O are sensitive probes of stellar nucleosynthesis processes and, in addition, O abundances are an important input for stellar age determinations. Understanding the nature of the observed distribution of these elements is key to constraining protogalactic star formation history. Patterns deduced from low-resolution spectroscopy of the CN, CH, NH, and CO molecules for low-mass stars in their core-hydrogen or first shell-hydrogen burning phases in the oldest ensembles known, the Galactic globular star clusters, are reviewed. New results for faint stars in NGC 104 (47 Tuc, C0021-723) reveal that the bimodal, anticorrelated pattern of CN and CH strengths found among luminous evolved stars is also present in stars nearing the end of their main-sequence lifetimes. In the absence of known mechanisms to mix newly synthesized elements from the interior to the observable surface layers of such unevolved stars, those particular inhomogeneities imply that the original material from which the stars formed some 15 billion years ago was chemically inhomogeneous in the C and N elements. However, in other clusters, observations of abundance ratios and C isotope ratios suggest that alterations to surface chemical compositions are produced as stars evolve from the main sequence through the red giant branch. Thus, the current observed distributions of C, N, and O among the brightest stars (those also observed most often) may not reflect the true distribution from which the protocluster cloud formed. The picture which is emerging of the C, N and O abundance patterns within globular clusters may be one which


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 difference between them in the way the duration distribution of faint bursts deviates from that of bright ones, assuming the luminosity function of gamma-ray bursts is independent of distance. 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: The Width of the Gamma-ray Burst Luminosity Function

abstract: We examine the width of the gamma-ray burst (GRB) luminosity function through the distribution of GRB peak count rates, C$_{\rm{peak}}$, as detected by BATSE (\cite{batse:93}). In the context of galactic corona spatial distribution models, we attempt to place constraints on the characteristic width of the luminosity function by comparing the observed intensity distribution with those produced by a range of density and luminosity functions. We find that the intrinsic width of the luminosity function cannot be very well restricted. However, the distribution of intrinsic luminosities of {\it detected bursts} can be limited: we find that most observed bursts have luminosities that are in a range of one to two decades, but a significant population of undetected less luminous bursts cannot be excluded. These findings demonstrate that the assumption that GRB are standard candles is sufficient but not necessary to explain the observed intensity distribution. We show that the main reason for the relatively poor constraints is the fact that the bright-end part of the GRB flux distribution is not yet sampled by BATSE, and better sampling in the future may lead to significantly stronger constraints on the width of the luminosity function.


title: The Reliability of X-ray Constraints of Intrinsic Cluster Shapes

abstract: Using the simulation of Katz & White we have tested the viability of X-ray analysis for constraining the intrinsic shapes of clusters of galaxies considering the effects of both substructure and steep temperature gradients. We restrict analysis to aggregate cluster shapes on scales of $r\sim 1-2$ Mpc to reduce our sensitivity to core subclustering. For low redshifts $(z\lesssim 0.25)$ the X-ray method accurately measures the true ellipticity of the three-dimensional cluster dark matter; assuming the gas is isothermal adds only small errors. At higher redshifts the X-ray method yields unreliable results since the gas does not trace the cluster gravitational potential. We proffer necessary conditions for the reliability of X-ray methods characterized by the amount of substructure in the X-ray surface brightness images and the isophotal shapes. We conclude that measurements of aggregate cluster shapes on scales $r\sim 1-2$ Mpc are insensitive to core substructure on scales of a few hundred kpc. Therefore our results suggest that the X-ray measurements of aggregate cluster shapes by Fabricant, Rybicki, & Gorenstein and Buote & Canizares are valid provided that they do not suffer from serious projection effects. A substantial number of Abell clusters observed with the ROSAT PSPC will be amenable to X-ray shape analysis.


title: The Near-Infrared Structure and Spectra of the Bipolar Nebulae M 2--9 and Afgl 2688: The Role of UV-Pumping and Shocks in Molecular Hydrogen Excitation

abstract: High-resolution near-infrared images and moderate resolution spectra were obtained of the bipolar nebulae M~2--9 and AFGL 2688. The ability to spatially and spectrally resolve the various components of the nebulae has proved to be important in determining their physical structure and characteristics. In M~2--9, the lobes are found to have a double-shell structure. Analysis of \h2\ line ratios indicates that the \h2\ emission is radiatively excited. A well-resolved photodissociation region is observed in the lobes. The spectrum of the central source is dominated by H recombination lines and a strong continuum rising towards longer wavelengths consistent with a $T = 795$ K blackbody. In AFGL 2688, the emission from the bright lobes is mainly continuum reflected from the central star. Several molecular features from C$_2$ and CN are present. In the extreme end of the N lobe and in the E equatorial region, the emission is dominated by lines of \h2 in the 2--2.5 \microns region. The observed \h2 line ratios indicate that the emission is collisionally excited, with an excitation temperature $T_{ex} \approx 1600\pm 100$ K.


title: Cosmic Ray Hydrodynamics at Shock Fronts

abstract: When a magnetised, thermal plasma containing a shock wave efficiently accelerates Cosmic Rays (CR) a reaction is exerted on the fluid structure. The simplest description of this process involves three conservation equations for the background fluid, with the CR pressure gradient included as an extra force, plus a hydrodynamic equation for the CR energy density. Since the CR spectrum is not determined in this approach two quantities are introduced, the ratio of CR pressure to energy density and a mean diffusion coefficient, that must be specified to close the system of equations. Physical arguments are presented for how these quantities should evolve and then compared with time dependent numerical solutions of acceleration at a plane parallel, piston driven shock.


title: Gamma-Rays from the Funnels of Thick Accretion Disks Around Active Galactic Nuclei

abstract: We analyse the propagation of particles in the narrow funnel of a thick accretion disk. It is assumed that: (1) the funnel walls emit black body radiation with temperature decreasing outwards; (2) the magnetic and electric fields are longitudinal in the funnel. Such a scenario is consistent with models in which a large potential drop is induced by a rotating massive black hole threaded by a magnetic field (Blandford & Znajek 1977, Macdonald & Thorne 1982). The interaction of relativistic protons with thermal photons from the funnel results in direct production of pairs and/or pions. We discuss the relative importance of these processes for different conditions in the funnel (temperature and electric field profiles). Injected pairs interact with thermal photons in the Thomson or Klein-Nishina regime. Under some conditions, a pair avalanche results, which we assume saturates in the production of stable bunches, containing almost equal numbers of electrons and positrons. As a consequence, highly collimated gamma-ray photons are produced. We obtain gamma-ray spectra from our model in order to test its applicability to the gamma-ray emitting AGNs recently reported by EGRET.


title: Evidences of Secular Dynamical Evolution in Disk Galaxies

abstract: After a recall of fundamental concepts used in galactic dynamics, we review observational facts as well as results of orbit theory and numerical simulations which suggest long-term evolution of galaxies. Dynamical interactions between galaxy constituents (bulges, disks, bars, haloes, dark matter) are discussed and effects of external perturbations on internal structures are examined. We report recent developments on the connection between dynamical interaction processes and efficiency of star formation in galaxies of various morphological types. The Hubble sequence as an evolution sequence is revisited.


title: Formalising the Slow-Roll Approximation in Inflation

abstract: The meaning of the inflationary slow-roll approximation is formalised. Comparisons are made between an approach based on the Hamilton-Jacobi equations, governing the evolution of the Hubble parameter, and the usual scenario based on the evolution of the potential energy density. The vital role of the inflationary attractor solution is emphasised, and some of its properties described. We propose a new measure of inflation, based upon contraction of the comoving Hubble length as opposed to the usual e-foldings of physical expansion, and derive relevant formulae. We introduce an infinite hierarchy of slow-roll parameters, and show that only a finite number of them are required to produce results to a given order. The extension of the slow-roll approximation into an analytic slow-roll expansion, converging on the exact solution, is provided. Its role in calculations of inflationary dynamics is discussed. We explore rational-approximants as a method of extending the range of convergence of the slow-roll expansion up to, and beyond, the end of inflation.


title: The Andante Regime of Scalar Field Dynamics

abstract: The andante regime of scalar field dynamics in the chaotic inflationary Universe is defined as the epoch when the field is rolling moderately slowly down its interaction potential, but at such a rate that first-order corrections to the slow-roll approximation become important. These conditions should apply towards the end of inflation as the field approaches the global minimum of the potential. Solutions to the Einstein-scalar field equations for the class of power law potentials $V(\phi) \propto \phi^{2n}$ are found in this regime in terms of the inverse error function.


title: Evolution of the Potential in Cosmological Gravitational Clustering

abstract: In recent years there has been a developing realization that the interesting large--scale structure of voids, ``pancakes", and filaments in the Universe is a consequence of the efficacy of an approximation scheme for cosmological gravitation clustering proposal by Zel'dovich in 1970. However, this scheme was only supposed to apply to smoothed initial conditions, and its efficacy for hierarchical clustering has not been justified. We show that this is explained by the fact that the gravitationally evolved potential from N--body simulations closely resembles the smoothed potential of the initial conditions. The resulting ``hierarchical pancaking" picture combines features of the former Soviet and Western theoretical pictures for galaxy and large--scale structure formation.


title: Velocity Fields in Non--Gaussian Cold Dark Matter Models

abstract: We analyse the large--scale velocity field obtained by N--body simulations of cold dark matter (CDM) models with non--Gaussian primordial density fluctuations, considering models with both positive and negative primordial skewness in the density fluctuation distribution. We study the velocity probability distribution and calculate the dependence of the bulk flow, one--point velocity dispersion and Cosmic Mach Number on the filtering size. We find that the sign of the primordial skewness of the density field provides poor discriminatory power on the evolved velocity field. All non--Gaussian models here considered tend to have lower velocity dispersion and bulk flow than the standard Gaussian CDM model, while the Cosmic Mach Number turns out to be a poor statistic in characterizing the models. Next, we compare the large--scale velocity field of a composite sample of optically selected galaxies as described by the Local Group properties, bulk flow, velocity correlation function and Cosmic Mach Number, with the velocity field of mock catalogues extracted from the N--body simulations. The comparison does not clearly permit to single out a best model: the standard Gaussian model is however marginally preferred by the maximum likelihood analysis.


title: Hierarchical Pancaking: Why the Zel'dovich Approximation Describes Coherent Large-Scale Structure in N-Body Simulations of Gravitational Clustering

abstract: To explain the rich structure of voids, clusters, sheets, and filaments apparent in the Universe, we present evidence for the convergence of the two classic approaches to gravitational clustering, the ``pancake'' and ``hierarchical'' pictures. We compare these two models by looking at agreement between individual structures -- the ``pancakes'' which are characteristic of the Zel'dovich Approximation (ZA) and also appear in hierarchical N-body simulations. We find that we can predict the orientation and position of N-body simulation objects rather well, with decreasing accuracy for increasing large-$k$ (small scale) power in the initial conditions. We examined an N-body simulation with initial power spectrum $P(k) \propto k^3$, and found that a modified version of ZA based on the smoothed initial potential worked well in this extreme hierarchical case, implying that even here very low-amplitude long waves dominate over local clumps (although we can see the beginning of the breakdown expected for $k^4$). In this case the correlation length of the initial potential is extremely small initially, but grows considerably as the simulation evolves. We show that the nonlinear gravitational potential strongly resembles the smoothed initial potential. This explains why ZA with smoothed initial conditions reproduces large-scale structure so well, and probably why our Universe has a coherent large-scale structure.


title: Are BL Lacertae Objects Beamed QSO Remnants?

abstract: This paper considers the hypothesis that BL Lacertae objects (BLLs) are the beamed remnants of Quasi Stellar Objects. The hypothesis explains why BLLs do not undergo the strong evolution seen in other active galactic nuclei since it naturally predicts that the space density of BLLs should increase with cosmic time, as shown by recent observations. Numerical models reproduce, with reasonable parameters, the known redshift and magnitude counts of BL Lac objects. It is assumed that radio-quiet as well as radio-loud quasars are capable of generating jets but that jets are snuffed in young radio-quiet objects and only emerge in aged ones. I argue that the observations allow this assumption.


title: BIMA observations of possible microwave background sources

abstract: We present sensitive upper limits on the 90 GHz flux of known radio and infrared sources in regions associated with possible cosmic microwave background fluctuations at 0.5--1 degree scales. Specifically we look at the MAX GUM region and the region of strongest fluctuation in the MSAM scan. None of the known sources can account for the levels of anisotropy seen.


title: The Orientation of Spin Vectors of Galaxies in the Ursa Major Filament

abstract: Spin vectors of 60 galaxies in the ``Ursa Major Filament'' are obtained from CCD images and spectroscopic determinations of the directions of rotation. These data are used to remove the four-fold degeneracy introduced by the project ed images of galaxies, making possible the complete reconstruction of the galaxy spin vectors. Several possible organizations of spin vectors are investigated by means of statistical analyses of the three-dimensional spin vector maps. The results indicate that there exists no significant alignment of galaxy spin vectors with respect to the supergalactic plane, the supergalactic axis, or the cylindrical axis of the filament itself. No tendency for linear alignment of spin vectors towards any specific direction was detected. We also obtain a statistically negative result for azimuthal, radial, and circul ar-helical orientation of spin vectors.


title: Big Bang Nucleosynthesis Constraints and Light Element Abundance Estimates

abstract: To elucidate the significance of the effect of systematic uncertainties in light element abundance estimates on cosmological bounds derivable from Big Bang Nucleosynthesis (BBN) we present tables giving bounds on $\Omega_{baryon}$ and $N_{\nu}$ as one changes the limits on primordial \he and $^7$Li. This allows us to derive new relations between these estimates and constraints on $\Omega_{baryon}$ and $N_{\nu}$. For example, only if the helium mass fraction, $Y_p \ge .245$ does $^7$Li (or D) presently play a role in placing an upper limit on the baryon density, and only if $Y_p \ge .250$ does $^4$He cease to play a role in bounding $\etaten$. All the elements combine together tend to give a stringent upper bound of $0.16$ on $\Omega_{baryon}$. We also find that $Y_p$ must exceed .239 for consistency between theory and observation if D+$^3$He/H is less than $ 10^{-4}$. Updated nuclear reaction rates, an updated neutron half life, Monte Carlo techniques, and correlations between the predicted abundances are incorporated in our analysis. We also discuss the handling of systematic uncertainties in the context of statistical analyses of BBN predictions.


title: Symplectic Integration without Roundoff Error

abstract: Most numerical integration algorithms are not designed specifically for Hamiltonian systems and do not respect their characteristic properties, which include the preservation of phase space volume with time. This can lead to spurious damping or excitation. Methods that do preserve all the Hamiltonian properties, i.e., for which the time-forward map is symplectic, are called symplectic integration algorithms or SIAs. Although such integrators are symplectic in theory, they are not symplectic if implemented using finite-precision arithmetic. This paper explains how to eliminate this problem by using ``Lattice SIAs'' and shows that these methods yield significant advantages when the computational error is dominated by roundoff. Using a lattice SIA is equivalent to evolving the exact solution of a problem with a Hamiltonian that is slightly different from the original. Lattice methods are useful for studies of the long-term evolution of Hamiltonian dynamical systems.

categories: astro-ph chao-dyn comp-gas nlin.CD nlin.CG nlin.SI solv-int


title: Finite Disks with Power-Law Potentials

abstract: We describe a family of circular, and elliptical, finite disks with a disk potential that is a power of the radius. These are all flattened ellipsoids, obtained by squashing finite spheres with a power-law density distribution, and cutoff at some radius Ro. First we discuss circular disks whose circular rotation speed v is proportional to r^alpha, with any alpha> -1/2. The surface-density of the disks is expressed in terms of hypergeometric functions of 1-(Ro/r)^2. We give closed expressions for the full 3-D potentials in terms of hypergeometric functions of two variables. We express the potential and acceleration in the plane at r>Ro, and along the rotation axis, in terms of simple hypergeometric functions. All the multipoles of the disk are given. We then generalize to non-axisymmetric disks. The potential in the midplane is given in terms of the hypergeometric function of two variables. For integer values of 2 alpha the above quantities are given in more elementary terms. All these results follow straightforwardly from formulae we derive for the general, cutoff, power-law, triaxial ellipsoid.


title: The Optical Gravitational Lensing Experiment. the Early Warning System

abstract: The discoveries of 17 microlensing event candidates have been reported over the last year by three teams conducting unprecedented mass photometric searches in the direction of the Galactic bulge and the Magellanic Clouds. These include 10 events found by the OGLE collaboration, 5 by the MACHO team and 2 by the EROS team. All searches have the main goal to detect dark matter in our Galaxy. The detection of 17 event candidates proves that the microlensing is a powerful tool in the search for dark matter, and it may be used for reliable mass determination when the geometry of the event is known. Here we present the first microlensing event, OGLE~\#11, discovered in real time, using the newly implemented "Early Warning System". We describe our system which makes it possible to monitor and study in great details any very rare phenomena, not only lensing events, with a broad array of instruments almost immediately after they have changed their brightness.


title: Tomography of Collisionless Stellar Systems

abstract: In this paper the concept of tomography of a collisionless stellar system of general shape is introduced, and a generalization of the Projected Virial Theorem is obtained. Applying the tomographic procedure we then derive a new family of virial equations which coincides with the already known ones for spherically symmetric systems. This result is obtained without any use of explicit expressions for the line-of-sight velocity dispersion, or spherical coordinate system.


title: Cold Dark Matter I: The Formation of Dark Halos

abstract: We use numerical simulations of critically-closed cold dark matter (CDM) models to study the effects of numerical resolution on observable quantities. We study simulations with up to $256^3$ particles using the particle-mesh (PM) method and with up to $144^3$ particles using the adaptive particle-particle --particle-mesh (P$^3$M) method. Comparisons of galaxy halo distributions are made among the various simulations. We also compare distributions with observations and we explore methods for identifying halos, including a new algorithm that finds all particles within closed contours of the smoothed density field surrounding a peak. The simulated halos show more substructure than predicted by the Press-Schechter theory. We are able to rule out all $\Omega=1$ CDM models for linear amplitude $\sigma_8\gsim 0.5$ because the simulations produce too many massive halos compared with the observations. The simulations also produce too many low mass halos. The distribution of halos characterized by their circular velocities for the P$^3$M simulations is in reasonable agreement with the observations for $150\kms\lsim V_{\rm circ} \lsim 350\kms$.}}


title: Cold Dark Matter II: Spatial and Velocity Statistics

abstract: We examine high-resolution gravitational N-body simulations of the $\Omega=1$ cold dark matter (CDM) model in order to determine whether there is any normalization of the initial density fluctuation spectrum that yields acceptable results for galaxy clustering and velocities. Dense dark matter halos in the evolved mass distribution are identified with luminous galaxies; the most massive halos are also considered as sites for galaxy groups, with a range of possibilities explored for the group mass to light ratios. We verify the earlier conclusions of White et al. (1987) for the low amplitude (high bias) CDM model --- the galaxy correlation function is marginally acceptable but that there are too many galaxies. We also show that the peak biasing method does not accurately reproduce the results obtained using dense halos identified in the simulations themselves. The COBE anisotropy implies a higher normalization, resulting in problems with excessive pairwise galaxy velocity dispersion unless a strong velocity bias is present. Although we confirm the strong velocity bias of halos reported by Couchman \& Carlberg (1992), we show that the galaxy motions are still too large on small scales. We find no amplitude for which the CDM model can reconcile simultaneously the galaxy correlation function, the low pairwise velocity dispersion, and the richness distribution of groups and clusters. With the normalization implied by COBE, the CDM spectrum has too much power on small scales if $\Omega=1$.


title: Detecting Bimodality in Astronomical Datasets

abstract: We discuss statistical techniques for detecting and quantifying bimodality in astronomical datasets. We concentrate on the KMM algorithm, which estimates the statistical significance of bimodality in such datasets and objectively partitions data into sub-populations. By simulating bimodal distributions with a range of properties we investigate the sensitivity of KMM to datasets with varying characteristics. Our results facilitate the planning of optimal observing strategies for systems where bimodality is suspected. Mixture-modeling algorithms similar to the KMM algorithm have been used in previous studies to partition the stellar population of the Milky Way into subsystems. We illustrate the broad applicability of KMM by analysing published data on globular cluster metallicity distributions, velocity distributions of galaxies in clusters, and burst durations of gamma-ray sources. PostScript versions of the tables and figures, as well as FORTRAN code for KMM and instructions for its use, are available by anonymous ftp from kula.phsx.ukans.edu.


title: A Survey of Galaxy Redshifts from Low-Resolution Slitless Spectra: Suggestion of Large-Scale Structures at Z ~0.1

abstract: We have conducted a galaxy survey based on low-resolution slitless spectra taken from the automated CFHT-Laval survey. We present redshift distributions for 522 galaxies distributed in 4 distinct regions of the sky. Redshifts are determined from the shifted positions of the 4000 A. stellar break using an automated break-finding algorithm. The redshifts so obtained have a precision of < 3000 km s-1, good enough to trace the large-scale distribution of galaxies. The most striking feature of the survey is an apparent excess of galaxies at z ~ 0.1 in 8 of the 10 fields probed. Although the statistical significance of the peaks appears marginal (~ 2.5s per field) after taken in account the small-scale galaxy clustering, the peaks may nevertheless reveal the existence of large-scale structures at z ~ 0.1, indicating that the known large structures in the local universe are not unusual.


title: Macho Parallaxes From A Single Satellite

abstract: Massive Compact Objects (Machos) are currently being discovered at substantially higher rates than would be expected from standard models of known stellar populations. To determine whether they are due to non-standard distri- butions of known populations or to a heretofore unknown (`dark') population, one must acquire more information about the individual events. Space-based parallaxes are potentially the best tool for extracting additional information. To leading order, parallax measurements by a single satellite result in a four- fold degeneracy: two possible values of $\tilde v$ and two possible signs for the component of motion perpendicular to the projected satellite-Earth vector. It had been believed that a second satellite would be required to break this degeneracy. I show that the velocity difference between the satellite and the Earth allows one to partially or totally break the degeneracy using a single satellite. For most Macho events it is possible to measure $\tilde v$ and $\tilde r_e$. For some it is also possible to measure $\Phi$. The proposed Space Infrared Telescope could measure $\sim 100$ parallaxes per year by applying $\sim 400\,$hr of telescope time.


title: Anticorrelated CN and CH Variations on the 47 Tucanae Main-Sequence Turnoff

abstract: Observations of CN and CH band strengths among a random sample of main-sequence turn-off stars (+3.9 < MV < +4.6) in the globular cluster 47 Tuc (NGC 104, C0021-723) were made with the CTIO Argus fiber spectrograph for the purpose of determining the ratio of CN-strong to CN-weak stars and investigating the behavior of CH relative to CN. Of the 20 turn-off stars, 12 were found to be CN-strong while 8 appear to be CN-weak. This ratio of CN-strong to CN-weak stars is similar to the ratios found among the more luminous 47 Tuc stars and implies little change in the overall distribution of CN with evolutionary state, although the present sample size is small. A general anticorrelation between CN and CH is also observed in that the CN-weak main-sequence stars all (with one possible exception) exhibit strong CH bands - a trend similar to that found among the brighter stars. That these variations occur among such relatively un-evolved stars and that the overall CN distribution appears to be independent of evolutionary state presents serious challenges to internal or mixing theories of their origin. We therefore suggest that at least some component of the C and N abundance inhomogeneities observed among the brighter (more evolved) stars of this cluster appears to have been established prior to the commencement of evolution up the red giant branch.


title: Cosmic Strings and Cosmic Variance

abstract: By using a simple analytical model based on counting random multiple impulses inflicted on photons by a network of cosmic strings we show how to construct the general q-point temperature correlation function of the Cosmic Microwave Background radiation. Our analysis is sensible specially for large angular scales where the Kaiser-Stebbins effect is dominant. Then we concentrate our study on the four-point function and in particular on its zero-lag limit, namely, the excess kurtosis parameter, for which we obtain a predicted value of $\sim 10^{-2}$. In addition, we estimate the cosmic variance for the kurtosis due to a Gaussian fluctuation field, showing its dependence on the primordial spectral index of density fluctuations $n$ and finding agreement with previous published results for the particular case of a flat Harrison-Zel'dovich spectrum. Our value for the kurtosis compares well with previous analyses but falls below the threshold imposed by the cosmic variance when commonly accepted parameters from string simulations are considered. In particular the non-Gaussian signal is found to be inversely proportional to the scaling number of defects, as could be expected by the central limit theorem.


title: Current Status of Axion Cosmology

abstract: We review the current status of axion cosmology in the light of recent work. We conclude that in the standard scenario, in which a global string network forms at the Peccei-Quinn phase transition, the largest contribution to the axion density is from loops, even when uncertainties in the nature of the QCD phase transition are taken into account. We also briefly discuss the implications of non-standard scenarios, such as an inflationary epoch and entropy production.


title: Observations of the Cosmic Microwave Background

abstract: The properties of the cosmic microwave background radiation provide unique constraints on cosmological models, i.e. on the content, history, and evolution of the Universe. I discuss the latest measurements of the spectral and spatial properties of the cosmic microwave background radiation. Recent measurements from NASA's Cosmic Background Explorer (COBE) satellite, and from balloon-borne and ground-based platforms, are summarized and their cosmological implications are discussed.


title: CMB Maps at 0.5 degree Resolution I: Full-Sky Simulations and Basic Results

abstract: We have simulated full-sky maps of the Cosmic Microwave Background (CMB) anisotropy expected from Cold Dark Matter (CDM) models at 0.5 and 1.0 degree angular resolution. Statistical properties of the maps are presented as a function of sky coverage, angular resolution, and instrument noise, and the implications of these results for observability of the Doppler peak are discussed. The rms fluctuations in a map are not a particularly robust probe of the existence of a Doppler peak, however, a full correlation analysis can provide reasonable sensitivity. We find that sensitivity to the Doppler peak depends primarily on the fraction of sky covered, and only secondarily on the angular resolution and noise level. Color plates and one-dimensional scans of the maps are presented to visually illustrate the anisotropies.


title: CMB Maps at 0.5 Degree Resolution II: Unresolved Features

abstract: High-contrast peaks in the cosmic microwave background (CMB) anisotropy can appear as unresolved sources to observers. We fit simulated CMB maps generated with a cold dark matter model to a set of unresolved features at instrumental resolution 0.5 degrees to 1.5 degrees and present the integral density per steradian n(>|T|) of unresolved features brighter than threshold temperature |T|. A typical medium-scale experiment observing 0.001 sr at 0.5 degree resolution would expect to observe one unresolved feature brighter than 85 \muK after convolution with the beam profile, with less than 5% probability to observe a source brighter than 150 \muK. Increasing the power-law index of primordial density perturbations n from 1 to 1.5 raises these temperature limits |T| by a factor of 2. The results are only weakly dependent on the assumed values of the Hubble constant and baryon density.


title: BD+30 3639: The Infrared Spectrum During Post-AGB Stellar Evolution

abstract: We present a radiative-transfer calculation which reproduces the infrared spectrum of the planetary nebula BD~+30$^{\circ}$3639. We calculate the transfer process through absorption and scattering in a spherical-symmetric multi-grain dust shell. The emission of transiently heated particles is taken into account, as well as polycyclic aromatic hydrocarbons. We obtain an acceptable fit to most of the spectrum, including the PAH infrared bands. At submillimetre wavelengths the observed emission is larger than the model predicts, indicating that large dust conglomerates (``f{}luffy grains'') may be needed as an additional constituent. The fit favours a distance of $\ge 2 \,$kpc, which implies that BD~+30$^\circ$3639 has evolved from a massive progenitor of several solar masses. A low dust-to-gas mass ratio is found in the ionised region. The calculations yield an original mass-loss rate of $2\times10^{-4} \msolar \peryr$ on the Asymptotic Giant Branch. Using this mass-loss rate, we calculate how the infrared spectrum has evolved during the post-AGB evolution. We show in particular the evolution of the IRAS colours during the preceding post-AGB evolution.


title: Semi-analytic Continuum Spectra of Type II Supernovae and the Expanding Photosphere Method of Distance Determination

abstract: We extend the approximate radiative transfer analysis of Hershkowitz, Linder, and Wagoner (1986) to a more general class of supernova model atmospheres, using a simple fit to the effective continuum opacity produced by lines (Wagoner, Perez, and Vasu 1991). At the low densities considered, the populations of the excited states of hydrogen are governed mainly by photoionization and recombination, and scattering dominates absorptive opacity. We match the asymptotic expressions for the spectral energy density $\Jnu$ at the photosphere, whose location at each frequency is determined by a first-order calculation of the deviation of $\Jnu$ from the Planck function $\Bnu$. The emergent spectral luminosity then assumes the form $ L_\nu = 4\pi^2 \r*^2 \zeta^2 \Bnu(T_p) $, where $T_p(\nu)$ is the photospheric temperature, $\zeta$ is the dilution factor, and $\r*$ is a fiducial radius [ultimately taken to be the photospheric radius $r_p (\nu)$]. The atmosphere is characterized by an effective temperature $T_e$ ($\propto L^{1/4} r_*^{-1/2}$) and hydrogen density $n_H=n_* (r_*/r)^\alpha$ ; and less strongly by the heavy element abundance and velocity gradient. Our major result is the dependence of $\zeta$ on frequency $\nu$ and the parameters $T_p$, $r_p$, and $\alpha$.


title: Magnetic Merging in Colliding Flux Tubes

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


title: Detection of Compact Ultraviolet Nuclear Emission in LINER Galaxies

abstract: Low-ionization nuclear emission-line regions (LINERs), which exist in a large fraction of galaxies, may be the least luminous manifestation of quasar activity. The nature of LINERs has, however, remained controversial because an AGN-like nonstellar continuum source has not been directly observed in them. We report the detection of bright, unresolved (FWHM $\ltorder 0.1''$) point sources of UV ($\sim 2300$ \AA) emission in the nuclei of nine nearby galaxies from a complete sample of 110 nearby galaxies imaged with {\it HST}. Ground-based optical spectroscopy reveals that five of the nuclei are LINERs, three are starburst nuclei, and one is a Seyfert nucleus. The observed UV flux in each of the five LINERs implies an ionizing flux that is sufficient to account for the observed emission lines through photoionization. The detection of a strong UV continuum in the LINERs argues against shock excitation as the source of the observed emission lines, and supports the idea that photoionization excites the lines in at least some objects of this class. Among the Northern-hemisphere galaxies in the sample, 26 are LINERs, of which only the above five LINERs have a detected nuclear UV source. There are no obvious differences in the optical line intensity ratios between the UV-bright and UV-dark LINERs. If all LINERs are photoionized, then the continuum source is unobscured along our line of sight in $5/26\approx 20\%$ of LINERs. Alternatively, spectrally-similar LINERs may be produced by various excitation mechanisms, with photoionization responsible in only about 20\% of the cases. The high angular resolution allows us to set upper limits, typically several parsecs, on the physical size of the compact star-cluster or AGN-type continuum source that is emitting the UV light in these objects.


title: Velocity Bias in MS1224+20

abstract: Velocity bias offers a single parameter ($b_v$) description of the likely segregation of mass and light in a galaxy cluster. The relation between the projected mass profile and the light profile, normalized by the virial $M/L$, is presented as a function of $b_v\equiv\sigma_\ell/\sigma_\rho$, the ratio of the RMS velocity dispersions of the galaxies and the underlying dark matter. The light distribution does not trace the mass at any radius within the cluster. At the apparent harmonic radius the gravitational mass can exceed the virial mass calculated from the light by about a factor of 2 (or more). The data for the MS1224+20 cluster, covering nearly a decade in radius, are well described by a velocity bias value of $b_v=0.85\pm0.05$, comparable to the n-body estimate of $b_v=0.8\pm0.1$. This particular value will remain tentative until more data are available.


title: Primordial Spectral Indices from Generalized Einstein Theories

abstract: Primordial spectral indices are calculated to second order in slow-roll parameters for three closely-related models of inflation, all of which contain a scalar field non-minimally coupled to the Ricci curvature scalar. In most cases, $n_s$ may be written as a function of the non-minimal curvature coupling strength $\xi$ alone, with $n_s (\xi) \leq 1$, although the constraints on $\xi$ differ greatly between \lq new inflation' and \lq chaotic inflation' initial conditions. Under \lq new inflation' initial conditions, there are discrepancies between the values of $n_s$ as calculated in the Einstein frame and the Jordan frame. The sources for these discrepancies are addressed, and shown to have negligible effects on the numerical predictions for $n_s$. No such discrepancies affect the calculations under \lq chaotic inflation' initial conditions.


title: Bias Minimisation in POTENT

abstract: Much current interest in cosmology has revolved around the peculiar velocity field of galaxies. Many methods have been proposed to derive this, probably the most prominent being POTENT (Dekel et al. 1990). However, this is prone to many forms of systematic error or bias. Here, we examine some of the methods currently used to minimise these biases, particularly Malmquist corrections, and highlight some of their pitfalls. As an alternative we propose an iterative scheme based on Monte Carlo simulations of POTENT that shows noticeable improvement over the other methods considered.


title: Likelihood Analysis of Large-Scale Flows

abstract: We apply a likelihood analysis to the data of \markcite{Lauer \& Postman 1994} With $P(k)$ parametrized by $(\sigma_8, \Gamma)$, the likelihood function peaks at $\sigma_8\simeq0.3$, $\Gamma\lesssim0.025$, indicating at face value very strong large-scale power, though at a level incompatible with COBE\@. There is, however, a ridge of likelihood such that more conventional power spectra do not seem strongly disfavored. The likelihood calculated using as data only the components of the bulk flow solution peaks at higher $\sigma_8$, in agreement with other analyses, but is rather broad. The likelihood incorporating both bulk flow and shear gives a different picture. The components of the shear are all low, and this pulls the peak to lower amplitudes as a compromise. The Lauer \& Postman velocity data alone are therefore {\em consistent}\/ with models with very strong large scale power which generates a large bulk flow, but the small shear (which also probes fairly large scales) requires that the power would have to be at {\em very}\/ large scales, which is strongly disfavored by COBE\@. The velocity data also seem compatible with more conventional $P(k)$ with $0.2\lesssim\Gamma\lesssim0.5$, and the likelihood is peaked around $\sigma_8\sim1$, in which case the bulk flow is a moderate, but not extreme, statistical fluctuation. Applying the same techniques to the data of \markcite{Riess, Press, \& Kirshner 1995}, the results are quite different. The flow is not inconsistent with the microwave dipole and we derive only an upper limit to the amplitude of the power spectrum: $\sigma_8\lesssim1.5$ at roughly 99\%.


title: Microlensing Events of the LMC are Better Explained by Stars within the LMC Than by Machos

abstract: The recently reported microlensing events of the LMC have caused much excitement, and have been interpreted as due to `dark objects' (MACHOs) in the halo of our Galaxy. It is shown here that the stars within the LMC play a dominant role as gravitational lenses and can indeed account for the observed events. The MACHO event and one of the EROS events lie within the bar for which the probability of microlensing is consistent with being caused by an object within the LMC. If the stars within the LMC play a dominant role as lenses, the events should be concentrated towards the center of the LMC. If MACHOs play a dominant role as lenses then, for a given number of monitored stars, the events should be uniformly distributed over the whole of the LMC. Thus the galactic and the LMC lenses can be statistically distinguished in most cases. It is further shown that, under certain conditions, the light curve of an event caused by a star within the LMC would be different from the one caused by a MACHO. This can also be a distinguishing signature, and seems to have been observed in case of the MACHO event. The fit discrepancy near the peak which the authors say ``is not yet understood" would be a natural consequence if the event is caused by an object within the LMC, which further indicates that the lensing is due to a low mass star within the LMC itself.


title: X-ray Variability and Period Determinations in the Eclipsing Polar DP Leo

abstract: An analysis of ROSAT observations for the eclipsing magnetic cataclysmic binary DP Leo provides constraints on the origin, size, temperature, variability and structure of the soft X-ray emission region on the surface of the white dwarf. These data, when combined with prior observations, show a progression of approximately 2 degrees per year in the impact position of the accretion stream onto the white dwarf. One explanation for the observed drift in stream position is that a magnetic activity cycle on the secondary produces orbital period oscillations. These oscillations result in an orbital period which cycles above and below the rotational period of the nearly synchronous white dwarf. The accretion stream and X-ray emission regions are modeled to fit the observational data. A distance to the system is also calculated. [An erroneous value for the cyclotron luminosity, included in an earlier paper version of the preprint, is corrected here.]


title: Gross Spectral Differences between Bright and Dim Gamma-Ray Bursts

abstract: We find that dim gamma-ray bursts (GRBs) are softer than bright GRBs, as indicated on average by data from the Burst and Transient Source Experiment (BATSE) on board the Compton Gamma Ray Observatory. We show that this correlation is statistically significant with respect to variations due to random differences between GRBs. This effect is discernable using a variety of methods and data sets, including public domain data. We analyze several types of systematic errors and selection effects in the BATSE data and conclude that the observed effect is not dominated by any of them. We therefore assert that this dim/soft effect is a real property of GRBs. It is possible that this correlation is a consequence of the time dilation detected by Norris et al. (1994) and that this is additional evidence that burst sources are located at cosmological distances.


title: Steps towards Nonlinear Cluster Inversion Through Gravitational Distortions. II. Generalization of the Kaiser & Squires Method

abstract: The weak distortions of high-redshift galaxies caused by gravitational light deflection near clusters of galaxies can be used to reconstruct the projected (two-dimensional) surface mass density of intermediate redshift clusters. This technique, pioneered by Tyson, Valdes & Wenk, and Kaiser & Squires, is reconsidered in the present paper, where we generalize the inversion equation found by Kaiser & Squires (KS) in several respect. Adopting a different smoothing procedure for the discreetly sampled data (individual galaxy images), we effectively reduce the shot noise in the KS procedure. In particular we show that the best density reconstructions are obtained if the smoothing scale is adopted to the `strength of the signal', which yields a better resolution near the center of the cluster where the distortions are strongest. Furthermore, we point out the importance of boundary effects and demonstrate their disastrous impact on rectangular data fields (CCDs) with large side ratio. Most important, however, is the generalization of the KS method to critical clusters, i.e., to such clusters which are capable of producing multiple images and giant luminous arcs. The corresponding modifications of the inversion procedure are severe; in particular, the resulting inversion equation is much more difficult to solve. As we pointed out in a previous paper (Schneider & Seitz), there exists a local degeneracy if the cluster is critical. We have developed an iteration procedure to solve the inversion equation, which we demonstrate to yield a very accurate


title: Super-Eddington Radiation Transfer in Soft Gamma Repeaters

abstract: Bursts from soft gamma repeaters have been shown to be super-Eddington by a factor of 1000 and have been persuasively associated with compact objects. Here, a model of super-Eddington radiation transfer on the surface of a strongly magnetic ($\geq 10^{13}$ gauss) neutron star is studied and related to the observational constraints on soft gamma repeaters. In strong magnetic fields, the cross-section to electron scattering is strongly suppressed in one polarization state, so super-Eddington fluxes can be radiated while the plasma remains in hydrostatic equilibrium. The model offers a somewhat natural explanation for the observation of similarity between spectra from bursts of varying intensity. The radiation produced in the model is found to be linearly polarized to about 1 part in 1000 in a direction determined by the local magnetic field, and the large intensity variations between bursts are understood as a change in the radiating area on the source. Therefore, the polarization may vary as a function of burst intensity, since the complex structure of the magnetic field may be more apparent for larger radiating areas. It is shown that for radiation transfer calculations in this limit of super-strong magnetic fields it is sufficient to solve the radiation transfer equations for the low opacity state rather than the coupled equations for both. With this approximation, standard stellar atmosphere techniques are utilized to calculate the model energy spectrum.


title: A New Method for Finding Clusters of Galaxies at z>1

abstract: At large redshifts, a cluster or group may be too distant for the galaxies within the cluster to be detected individually. However, the light from these ``undetected'' galaxies still modulates the surface brightness of the background sky. Clusters can appear as 10"-1.5' sized fluctuations in the surface brightness of the EBL. The fluctuations have central surface brightnesses between roughly 26 and 28 mag/sqr-arcsec (in V) for clusters between z=1 and z=2, and are brighter than the fluctuations produced by background field galaxies. While such low surface brightnesses are difficult to achieve with direct high-resolution imaging, we demonstrate that they are easily reached in short exposures through smoothing the sky in very flat CCD images. For a reasonable extrapolation of the properties and space densities of clusters and groups, we find that for a wide range of cosmological assumptions there should be tens of clusters per square degree visible in the extragalactic background. The detection rate can range between 0.1 and 100 clusters per square degree, for extreme assumptions about the rate of cluster and galaxy evolution. Drift scans provide sufficient accuracy of flat-fielding to make the method discussed in this paper an efficient techinque for finding candidate high redshift clusters.


title: Gravitational Radiation from Rapidly Rotating Nascent Neutron Stars

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


title: Hydrodynamics of Coalescing Binary Neutron Stars: Ellipsoidal Treatment

abstract: We employ an approximate treatment of dissipative hydrodynamics in three dimensions to study the coalescence of binary neutron stars driven by the emission of gravitational waves. The stars are modeled as compressible ellipsoids obeying a polytropic equation of state; all internal fluid velocities are assumed to be linear functions of the coordinates. The hydrodynamic equations then reduce to a set of coupled ordinary differential equations for the evolution of the principal axes of the ellipsoids, the internal velocity parameters and the binary orbital parameters. Gravitational radiation reaction and viscous dissipation are both incorporated. We set up exact initial binary equilibrium configurations and follow the transition from the quasi-static, secular decay of the orbit at large separation to the rapid dynamical evolution of the configurations just prior to contact. A hydrodynamical instability resulting from tidal interactions significantly accelerates the coalescence at small separation, leading to appreciable radial infall velocity and tidal lag angles near contact. This behavior is reflected in the gravitational waveforms and may be observable by gravitational wave detectors under construction.


title: Resolving the Question of Time for Semiclassical Gravity

abstract: Hamilton-Jacobi theory provides a natural starting point for a covariant description of the gravitational field. Using a spatial gradient expansion, one may solve for the phase of the wavefunction by using a line-integral in superspace. Each contour of integration corresponds to a particular choice of the time-hypersurface, and each yields the same answer. In this way, one can describe all time choices simultaneously. As a demonstration of the formalism, I will compute large-angle microwave background anisotropies and the galaxy-galaxy correlation function associated with scalar and tensor fluctuations of power-law inflation.


title: Effect of Adiabatic Deceleration on the Focused Transport of Solar Cosmic Rays

abstract: In the framework of focused transport theory, adiabatic deceleration arises from adiabatic focusing in the solar wind frame and from differential solar wind convection. An explicit formula is given for the deceleration of individual particles as a function of the pitch angle. Deceleration and other first-order effects of the solar wind, including convection, are incorporated into a numerical code for simulating the transport of energetic particles along the interplanetary magnetic field. We use this code to model the transport of solar flare protons. We find that including deceleration can increase the decay rate of the near-Earth intensity by 75\% more than would be expected based on advection from higher momenta, due to an interplay with diffusive processes. Improved response functions are derived for the impulsive injection of particles near the Sun, and it is found that neglecting deceleration leads to incorrect estimates of the scattering mean free path based on the intensity decay alone, especially for lower-energy particles.


title: A comparison of the evolution of the density field in perturbation theory and numerical simulations - II Counts in cells analysis

abstract: We present a detailed comparison of the predictions of perturbation theory for the averaged J-point correlation functions, $\xibar_J$, with the results of numerical simulations of gravitational clustering. We have carried out a systematic analysis of this method using ensembles of simulations with different numbers of particles, different box sizes and using different particle arrangements and clustering amplitudes in the initial conditions. We estimate $\xibar_J$, for $J=2-10$, from moments of counts-in-cells. We find significant non-linear effects in the variance, $J=2$, even at scales as large as $R \sim 30 \Mpc$. Perturbation theory gives remarkable agreement at large scales, where $\xibar_2 \simlt 1$, with the measured hierarchical amplitudes $S_J=\xibar_J/\xibar_2^{J-1}$. We have followed the evolution of $\xibar_J$ in time and find that for a change in the effective redshift of at least $\Delta z \simeq 2$ the amplitudes $S_J$ remain unchanged, despite the fact that the $\xibar_J$ have evolved by large factors, $\simeq 10^{J-1}$. We illustrate how these results can be applied to interpret the clustering in galaxy surveys and conclude that the observed hierarchical pattern in the APM is compatible with gravitational evolution in unbiased, initially Gaussian, models.


title: Stellar Models with Microscopic Diffusion and Rotational Mixing I: Application to the Sun

abstract: The Yale stellar evolution code has been modified to include the combined effects of diffusion and rotational mixing on $^1$H, $^4$He and the trace elements $^3$He, $^6$Li, $^7$Li, and $^9$Be. The interaction between rotational mixing and diffusion is studied by calculating a number of calibrated solar models. The rotational mixing inhibits the diffusion in the outer parts of the models, leading to a decrease in the envelope diffusion by 25 -- 50\%. Conversely, diffusion leads to gradients in mean molecular weight which can inhibit the rotational mixing. The degree to which gradients in mean molecular weight inhibits the rotational mixing is somewhat uncertain. A comparison to the observed solar oblateness suggests that gradients in the mean molecular weight play a smaller role in inhibiting the rotational mixing previously believed. This is reinforced by the fact that the model with the standard value for the inhibiting effect of mean molecular weight on the rotational mixing depletes no Li on the main sequence. This is in clear contrast to the observations. A reduction in the inhibiting effect of mean molecular weight gradients by a factor of ten leads to noticeable main sequence Li depletion.


title: Stellar Models with Microscopic Diffusion and Rotational Mixing II: Application to Open Clusters

abstract: Stellar models with masses ranging from 0.5 to $1.3~M_\odot$ were constructed in order to compare to young cluster observations of Li and of rotation velocities. The amount of Li depletion in cool stars is sensitive to the amount of overshoot at the base of the surface convection zone, and the exact metallicity of the models. Even when this uncertainty is taken into account, the Li observations are a severe constraint for the models and rule out standard models, and pure diffusion models. Stellar models which include diffusion and rotational mixing in the radiative regions of stars are able to simultaneously match the Li abundances observed in the Pleiades, UMaG, Hyades, Praesepe, NGC 752 and M67. They also match the observed rotation periods in the Hyades. However, these models are unable to simultaneously explain the presence of the rapidly rotating late G and K stars in the Pleiades and the absence of rapidly rotating late F and early G stars.


title: Microlensing Events: End of the Dark Halo?

abstract: I obtain an upper limit for the optical depth to microlensing toward Baade's Window of $3\times 10^{-6}$ by assuming that all of the mass of the Galaxy interior to the Sun (and not in the bulge) is in a disk. The exponential scale height of the disk is left as a completely arbitrary function of radius and is varied to maximize the optical depth. I take account of the relatively small corrections induced by the fact that the bulge is not axisymmetric. If initial estimates by the OGLE collaboration of an observed optical depth $\tau\sim 3.3\times 10^{-6}$ are confirmed, then essentially all of the dark matter interior to the Sun must be in a disk with a scale height of a few hundred parsecs.


title: A Merger Model and Globular Cluster Formation

abstract: We propose a self-consistent model for globular cluster formation in, but not limited to, our Galaxy, based on the merger model of Mathews and Schramm (1993). Stars and star clusters form in bursts at the merging interfaces as protogalactic clouds collide. We describe the formation of those star clusters with a simple schematic ansatz which takes into account the thermal and Kelvin-Helmholtz instabilities. It is shown that this model is consistent with many observational properties such as the age and metallicity distributions of globular clusters, the overall number of globular clusters, and the near constancy of the number of globular clusters in different size host galaxies. Most of the features of this merger model are insensitive to choices of parameters. However, the model does not exhibit two distinct populations of globular clusters, i.e. halo clusters and disk clusters. Possible explanations for this are presented.


title: The Deuterium Abundance and Nucleocosmochronology

abstract: We examine galactic chemical evolution models which reproduce the present-day and pre-solar values of deuterium starting with a primordial value which is consistent with a baryon-to-photon ratio of $3 \times 10^{-10}$. We consider various galactic chemical evolution models to determine the viability of significant deuterium destruction and which provide a consistent age of the galaxy at the time of the formation of the solar system and consequently its present day age from nuclear chronometers. These models generally require some amount of infall which we take with rates proportional to the gas mass as well as exponentially decreasing rates and some initial disk enrichment which we limit to the range of 0\% to 30\%. We present those models which give the observed pre-solar value and present-day value of D/H and which lead to a present-day gas fraction of $\sigma = .05\quad - \quad .2$. These models result in a braod range for the age of galaxy between $9.8 - 21.6$ Gyrs.


title: Duration Distributions of Bright and Dim BATSE Gamma-Ray Bursts

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


title: The complex theory of gravitational lensing: Beltrami equation and cluster lensing

abstract: We introduce a new complex formalism to describe arclet fields in clusters of galaxies and derive the appropriate inversion techniques to find the related mass distribution of the lensing cluster. Applying the complex formalism to the statistics of lensed background sources shows that the mean of the axial ratios of locally lensed isotropically distributed elliptical sources is {\em not} influenced by the intrinsic ellipticity distribution, if the cluster is subcritical in that region.


title: Galaxy Formation in a Variety of Hierarchical Models

abstract: We predict the observable properties of the galaxy population in popular hierarchical models of galaxy formation. We employ a detailed semianalytic procedure which incorporates the formation and merging of dark matter halos, the shock heating and radiative cooling of gas, self-regulated star formation, the merging of galaxies within dark matter halos, and the spectral evolution of the stellar populations. We contrast the standard CDM cosmogony with variants of the CDM model having either a low value of H_0, or a low value of Omega with or without a cosmological constant. In addition, we compare galaxy formation in these CDM universes with a CHDM model. We find that although the models have some success in remedying the shortcomings of the standard CDM cosmogony, none of these new models produce broad agreement with the whole range of observations. Although the low-Omega and Omega+Lambda=1 CDM models reduce the discrepancy between the predicted and observed Tully-Fisher relations (the main weakness of galaxy formation in standard CDM), these models predict an inverted colour-magnitude relation and do not produce an exponential cut-off at the bright end of the galaxy luminosity function. All of our models predict recent star formation and exhibit galaxy colours bluer than observed, but this problem is far more severe in the CHDM model which produces colours about two magnitudes too blue in B-K. Unlike in the variants of the CDM model in the CHDM case this result is not dependent on our model of stellar feedback, but is instead directly caused by the late epoch of structure formation in this model.


title: Interpreting Large Scale Structure Observations

abstract: The standard model of large scale structure is considered, in which the structure originates as a Gaussian adiabatic density perturbation with a nearly scale invariant spectrum. The basic theoretical tool of cosmological perturbation theory is described, as well as the possible origin of the density perturbation as a vacuum fluctuation during inflation. Then, after normalising the spectrum to fit the cosmic microwave background anisotropy measured by COBE, some versions of the standard model are compared with a variety of data coming from observations of galaxies and galaxy clusters. The recent COBE analysis of G\'{o}rski and collaborators is used, which gives a significantly higher normalization than earlier ones. The comparison with galaxy and cluster data is done using linear theory, supplemented by the Press-Schechter formula when discussing object abundances of rich clusters and of damped Lyman alpha systems. By focussing on the smoothed density contrast as a function of scale, the observational data can be conveniently illustrated on a single figure, facilitating easy comparison with theory. The spectral index is constrained to $0.6<n<1.1$, and in particle physics motivated models that predict significant gravitational waves the lower limit is tightened to $0.8$. [To appear, Proceedings of Journee Cosmologie, Paris, June 1994]


title: The Assembly of Galaxies in a Hierarchically Clustering Universe

abstract: We study the formation of galaxies by using $N$-body/hydrodynamics simulations to investigate how baryons collect at the centre of dark matter halos. We treat the dark matter as a collisionless fluid and the baryons as an ideal gas. We include the effects of gravity, pressure gradients, hydrodynamical shocks, and radiative energy losses, but we neglect star formation. Our initial conditions assume a flat universe dominated by cold dark matter with a mean baryon abundance of 10\% by mass. Typical halos form through the merging of a few smaller systems which had themselves formed in a similar manner at higher redshift. The gas collects at the bottom of dark matter potential wells as soon as these are properly resolved by our simulations. There it settles into cold, tightly bound disks, and it remains cold during subsequent evolution. As their halos coalesce, these disks merge on a timescale that is consistent with dynamical friction estimates based on their {\it total} (gas + surrounding dark matter) mass. Both the merger rates of the disks and their mass spectrum are in remarkably good agreement with recent analytic models that describe the evolution of {\it dark halos} in a hierarchical universe. This very simple model of galaxy formation suffers from serious shortcomings. It predicts that most baryons should be locked up in galaxies, whereas in the real universe most baryons are thought to lie outside visible galaxies. In addition, it predicts the specific angular momentum of a disk to be only about 20\% that of its surrounding halo, corresponding to a radius smaller than that of observed spiral galaxy disks.


title: The Damping and Excitation of Galactic Warps by Dynamical Friction

abstract: We investigate the dynamical interaction of galactic warps with the surrounding dark matter halo, using analytic perturbation theory. A precessing warp induces a density wake in the collisionless dark matter, which acts back on the original warp, transferring energy and angular momentum between the warp and halo (dynamical friction). In most cases dynamical friction damps the warp, but in unusual circumstances (such as a halo that rotates in the same direction as the precession of the warp, or a warp in the equatorial plane of an axisymmetric prolate halo) friction can excite a warp. The damping/excitation time is usually short compared to the Hubble time for realistic systems. Thus most warps cannot be primordial; they must be maintained by some ongoing excitation mechanism.


title: Simulations of X-ray Clusters

abstract: We present simulations of the formation and evolution of galaxy clusters in the Cold Dark Matter cosmogony. Clusters with a wide range of mass were selected from previous N-body models, and were resimulated at higher resolution using a combined N-body/Smooth Particle Hydrodynamics code. The effects of radiative cooling on the gas are neglected. While many present-day clusters are predicted to be undergoing mergers, the density profiles of those that are approximately in equilibrium are all very similar, both for the gas and for the dark matter. These profiles show no sign of a uniform density core and steepen gradually from the centre outwards. The standard $\beta$-model is a reasonable fit over most of the radius range observable in real clusters. However, the value obtained for the slope parameter $\beta_f$ increases with the outermost radius of the fit. Temperature profiles of different simulated clusters are also similar. Typically the temperature is almost uniform in the regions which emit most of the X-ray flux but drops at larger radii. The gas temperature and dark matter velocity dispersion in equilibrium clusters give values of $\beta_T\equiv \mu m_p\sigma_{DM}^2/kT$ which are consistent with unity provided an X-ray emission-weighted temperature is used. Larger values of $\beta_T$ are found in merging objects where there is a transient boost in the velocity dispersion of the system. Thus $\beta_T >1$ may be an observational indicator of merging in real clusters. The similar structure of clusters of differing mass results in scaling relations between the X-ray and dynamical properties of clusters identified at


title: Gaussian Statistics of the Cosmic Microwave Background: Correlation of Temperature Extrema In the COBE DMR Two-Year Sky Maps

abstract: We use the two-point correlation function of the extrema points (peaks and valleys) in the COBE Differential Microwave Radiometers (DMR) 2-year sky maps as a test for non-Gaussian temperature distribution in the cosmic microwave background anisotropy. A maximum likelihood analysis compares the DMR data to n=1 toy models whose random-phase spherical harmonic components $a_{\ell m}$ are drawn from either Gaussian, $\chi^2$, or log-normal parent populations. The likelihood of the 53 GHz (A+B)/2 data is greatest for the exact Gaussian model. All non-Gaussian models tested are ruled out at 90\% confidence, limited by type II errors in the statistical inference. The extrema correlation function is a stronger test for this class of non-Gaussian models than topological statistics such as the genus.


title: The Clustering Properties of Faint Galaxies

abstract: The two-point angular correlation function of galaxies, \wte, has been computed from a new survey of faint galaxies covering a 2 deg$^2$ area near the North Galactic Pole. This survey, which is complete to limiting magnitudes \jmag=24 and \fmag=23, samples angular scales as large as $1\degpoint5$. Faint galaxies are found to be more weakly clustered (by a factor of at least two) compared to galaxies observed locally. Clustering amplitudes are closer to model predictions in the red than in the blue. The weak clustering of faint galaxies cannot be explained by any plausible model of clustering evolution with redshift. However, one possible explanation of the clustering properties of intermediate redshift galaxies is that they resemble those of starburst galaxies and H II region galaxies, which are observed locally to possess weak clustering amplitudes. Our clustering amplitudes are also similar to those of nearby late-type galaxies, which are observed to be more weakly clustered than early-type galaxies A simple, self-consistent model is presented that predicts the fraction of galaxies in the ``excess'' population at intermediate redshifts, and correctly matches observed color distributions. The available data on the clustering properties of faint galaxies are consistent with this model if the ``excess'' population of faint blue galaxies is also the weakly clustered population. Evidence is presented that the power-law slope of the angular correlation function becomes shallower at fainter magnitudes. A similar effect is seen locally both for dwarf galaxies and for galaxies with late morphological type;


title: The Medium Redshift Clusters CL0017-20 and CL0500-24

abstract: We present magnitudes and redshifts of galaxies in the cluster of galaxies CL0017-20 and in the double cluster CL0500-24. Luminosity functions and velocity distributions have been derived. CL0017-20 possesses a virial mass of $2.9 \times 10^{14}$ M$_\odot$ and a \vmag band luminosity of $2.3 \pm 0.5 \times 10^{12}$ L$_\odot$ ($M/L_{\rm V} = 127$), much of it in a very compact core of size 113 kpc. Our new velocities and photometry of CL0500-24 confirm its double nature. Our group algorithm clearly detects two subconcentrations with masses 1.9 and $2.2 \times 10^{14}$ M$_{\odot}$. Its total \vmag band luminosity amounts to $2.7 \pm 0.8 \times 10^{12}$ L$_{\odot}$ ($M/L_{\rm V}$ of 157). Since the crossing and collapse times in both clusters are smaller than the age of the Universe, they appear bound and virialized. It is surprising that although $M/L$ is small in these clusters, giant arcs have been detected in both.


title: Counts and Colours of Faint Stars in 5 Fields Near the North Galactic Pole

abstract: Faint star number counts in the photographic \jmag band and \bvcol colour distributions are presented for a 1.08 deg$^2$ field near the North Galactic Pole. Due to the excellent star/galaxy discrimination we count stars as faint as \jmag = 23. We compare the number counts and colour distributions in 5 adjacent fields near SA57. The number counts and colour distributions are in very good agreement with previous data. However, we find that the large field-to-field scatter in the colour distributions, which we argue is real, might prevent us setting strong limits on Galactic structure. A simple two component standard model, Bahcall and Soneira (1984), fits the number counts reasonably well at the bright \jmag $< 21$, but fails notably at the faint end, even if a third component is added, as in Reid and Majewski (1993). The standard models are in good agreement with both the number counts and colour distribution at $20 < $\vmag$< 21$. Although the standard models bimodal shape of the colour distribution compares well with the data at $21 < $\vmag$< 21.5$, the number counts of stars are underestimated.


title: The Asymmetric Drift of the Thick Disk Population

abstract: We investigate the kinematics of the intermediate population using the photometric and astrometric sample surveys towards two opposite directions (galactic centre and anticentre) at intermediate latitude. A multivariate discriminant analysis (MDA) is used to distinguish the `Thick Disk' from other populations with the help of the Besan\c con model of population synthesis. The data constrain the asymmetric drift of the intermediate population, which is found to be 40$\pm$10 km sec$^{-1}$ and does not depend on the galactic radius.


title: Astrometric Plate Reduction with Orthogonal Functions and Milli-Arcseconds Accuracy in Deep Proper Motion Surveys

abstract: We have been doing a sample survey in UBV photometry and proper motions as part of an investigation of galactic structure and evolution. The programme is based on Schmidt plates (from OCA, Tautenburg, Palomar and ESO Schmidt telescopes) digitized with the MAMA machine (CAI, Paris). The high astrometric quality of the MAMA gives access to micronic accuracy (leading to $<$ 2 mas per year) on proper motions with a 35 years time base. The high proper motion accuracy for faint star probes in wide-areas give access to the properties of star samples out of the solar neighborhood. We have analyzed components of the UVW galactic space motions resulting from the accurate proper motion surveys. The kinematical distribution of F and G--type stars have been probed to distances up to 2.5 kpc above the galactic plane. We have derived the constrain on the structural parameters of the thin and thick disk components of the Galaxy.


title: Probing the Era of Galaxy Formation via TeV Gamma Ray Absorption by the Near Infrared Extragalactic Background

abstract: We present models of the extragalactic background light (EBL) based on several scenarios of galaxy formation and evolution. We have treated galaxy formation with the Press-Schecter approximation for both cold dark matter (CDM) and cold+hot dark matter (CHDM) models. Galaxy evolution has been treated by considering a variety of stellar types, different initial mass functions and star formation histories, and with an accounting of dust absorption and emission. We find that the dominant factor influencing the EBL is the epoch of galaxy formation. A recently proposed method for observing the EBL utilizing the absorption of $\sim 0.1$ to 10 TeV gamma-rays from active galactic nuclei (AGN) is shown to be capable of discriminating between different galaxy formation epochs. Observations of TeV gamma-ray bursts (GRB) would constrain or eliminate models in which the GRB sources lie at cosmological distances.


title: Self-Consistent Particle Acceleration in Active Galactic Nuclei

abstract: Adopting the hypothesis that the nonthermal emission of Active Galactic Nuclei (AGN) is primarily due to the acceleration of protons, we construct a simple model in which the interplay of acceleration and losses can be studied together with the formation of the emitted spectrum. The acceleration process is assumed to be of the first order Fermi type, and the proton distribution as well as the injected electrons and photons in the central region of the AGN are described by spatially averaged kinetic equations. The various relevant processes which dominate the three species are incorporated into the equations. The technique used to solve these is presented and several tests of the numerical implementation are presented. We also present results of a sample time-dependent AGN model in which photons appear suddenly as a result of a feedback instability and the system evolves to a steady state, in which the acceleration process is saturated self-consistently by the photons it produces. This example combines an X-ray power law index of about -1.7, together with a break at an energy between 50 and 500~keV.


title: String Radiative Backreaction

abstract: We discuss radiative backreaction for global strings described by the Kalb-Ramond action with an analogous derivation to that for the point electron in classical electrodynamics. We show how local corrections to the equations of motion allow one to separate the self-field of the string from that of the radiation-field. Modifications to this 'local backreaction approximation' circumvent the runaway solutions familiar from electrodynamics, allowing these corrections to be used to evolve string trajectories numerically. Comparisons are made with analytic and numerical radiation calculations from previous work and the merits and limitations of this approach are discussed. Finally, we comment on the relevance of this work to simulations of a cosmological network of global or axion strings. These methods can also be applied to describe graviational, electromagnetic and other forms of string radiative backreaction.


title: Paramagnetic alignment of fractal grains

abstract: Paramagnetic alignment of fractal suprathermally rotating grains is discussed. It is shown that if the concentration of H$_{2}$formation sites is low and resurfacing is active, fractal structure of grains enhances their alignment. Studying the influence of grain surface physics and chemistry on the alignment we found that there exist two critical values of grain sizes, and the alignment of smaller grains is expected to decrease as compared to the predictions of the Purcell's theory (1979). One of the critical sizes is due to intensified poisoning of active sites, while the other is a result of a finite value of the imaginary part of magnetic susceptibility in the limit of high angular velocities. On the contrary, if active sites completely cover grain surface, suprathermal rotation, and therefore efficient alignment, is possible only for a limited range of grain sizes.


title: The Cloud-in-Cloud Problem in the Press-Schechter Formalism of Hierarchical Structure Formation

abstract: The formalism by Press and Schechter (PS) is often used to infer number densities of virialized objects of mass M (e.g. quasars, galaxies, clusters of galaxies, etc.) from a count of initially overdense regions in a Gaussian density perturbation field. We reanalyze the PS-formalism by explicitly counting underdense regions which are embedded within overdense regions, so called cloud-in-clouds. In contrast to the original PS-formalism, our revised analysis automatically accounts for all the cosmic material. We find that mass distribution functions for virialized objects are altered by the proper solution of the cloud-in-cloud problem. These altered distribution functions agree much better with distribution functions inferred form N-body simulations than the original PS-distribution functions.


title: Baryon Number Transport in a Cosmic QCD-Phase Transition

abstract: We investigate the transport of baryon number across phase boundaries in a putative first order QCD-phase transition. Two independent phenomenological models are employed to estimate the baryon penetrability at the phase boundary: chromoelectric flux tube models; and an analogy to baryon-baryon coalescence in nuclear physics. Our analysis indicates that baryon transport across phase boundaries may be order of magnitude more efficient than other work has suggested. We discuss the substantial uncertainties involved in estimating baryon penetrability at phase boundaries.


title: Arc Statistics with Realistic Cluster Potentials. II. Influence of Cluster Asymmetry and Substructure

abstract: We construct a sample of numerical models for clusters of galaxies and employ these to investigate their capability of imaging background sources into long arcs. Emphasis is laid on the statistics of these arcs. We study cross sections for arc length and length-to-width ratio and optical depths for these arc properties, we examine the distribution of arc widths and curvature radii among long arcs, and we compare these results to predictions based on simplified (radially symmetric) cluster models. We find that the capability of the numerically modeled clusters to produce long arcs is larger by about two orders of magnitude than that of spherically symmetric cluster models with the same observable parameters (core radii and velocity dispersions), and that they are similarly efficient as singular isothermal spheres with the same velocity-dispersion distribution. The influence of source ellipticity is also investigated; we find that the optical depth for arcs with a length-to-width ratio $\ga10$ is significantly larger for elliptical than for circular sources. Given these results, we conclude that spherically symmetric lens models for galaxy clusters, adapted to the observable parameters of these clusters, grossly underestimate the frequency of long arcs. We attribute this difference between numerically constructed and simplified analytical lens models to the abundance and the extent of intrinsic asymmetry and to substructure in galaxy clusters.


title: Deep Imaging of the Double Quasar 0957+561: New Constraints on H_0

abstract: We present new results from extremely deep, high-resolution images of the field around the double quasar QSO 0957+561. A possible gravitational arc system near the double quasar has recently been reported, which, if real, would set strong constraints on determinations of the Hubble constant from the time delay in the double quasar. We find that both the morphology and the colors of the claimed arc systems suggest that they are chance alignments of three and two different objects, and not gravitationally lensed arcs. Hence, the constraints on $H_0$-determinations from the arcs are not valid. Also, a small group of galaxies at $z=0.5$ near the line-of-sight which was required to have a very large mass in the physically interesting arc models, is most likely insignificant. From our deep images we are able to use weak lensing of faint background galaxies in the field to map the gravitational potential in the main cluster. This sets new constraints on determinations of $H_0$. We find that the Hubble constant is constrained to be less than 70km/(s Mpc), if the time delay between the two images of the QSO is equal to or larger than 1.1 years.


title: Parameter degeneracy in models of the quadruple lens system Q2237+0305

abstract: The geometry of the quadruple lens system Q2237+0305 is modeled with a simple astigmatic lens: a power-law mass distribution: m \propto r^\beta with an external shear \gamma. The image positions can be reproduced with an accuracy better than 0.01 arcsec for any 0.0 < \beta < 1.85 and the corresponding value of \gamma = 0.1385 - 0.0689 \beta. This is a factor of about 4 more precise than what can be achieved by the best constant M/L lens models. The image intensity ratios and the time delay ratios are almost constant along our one parameter family of models, but the total magnification varies from 8 to > 1000, and the maximum time delay (between leading image B and trailing image C) for H_0 of 75 km/sec/Mpc varies from more than 20 hours to about 1.5 hours, while \beta increases from 0.0 to 1.85.


title: Sky Subtraction with Fiber-Spectrographs

abstract: The sky subtraction performances of multi-fiber spectrographs are discussed, analysing in detail the case of the OPTOPUS system at the 3.6 meter ESO telescope at La Silla. A standard technique, based on flat-fields obtained with a uniformly illuminated screen on the dome, provides poor results. A new method has been developed, using the [OI] emission line at 5577 \AA\ as a calibrator of the fiber transmittance, taking into account the diffuse light and the influence of each fiber on the adjacent ones, and correcting for the effects of the image distortions on the sky sampling. In this way the accuracy of the sky subtraction improves from $2-8 \%$ to $1.3-1.6 \%$.