title: Testing the Hypothesis of the Morphological Transformation from Field Spiral to Cluster S0

abstract: Hubble Space Telescope observations of distant clusters have suggested a steep increase in the proportion of S0 galaxies between clusters at high redshifts and similar systems the present-day. It has been proposed that this increase results from the transformation of the morphologies of accreted field galaxies from spirals to S0s. We have simulated the evolution of the morphological mix in clusters based on a simple phenomenological model where the clusters accrete a mix of galaxies from the surrounding field, the spiral galaxies are transformed to S0s (through an unspecified process) and are added to the existing cluster population. We find that in order to reproduce the apparently rapid increase in the ratio of S0 galaxies to ellipticals in the clusters, our model requires that: (1) galaxy accretion rate has to be high (typically, more than half of the present-day cluster population must have been accreted since z~0.5), and (2) most of the accreted spirals, with morphological types as late as Scdm, must have transformed to S0's. Although the latter requirement may be difficult to meet, it is possible that such bulge-weak spirals have already been `pre-processed' into the bulge-strong galaxies prior to entering the cluster core and are eventually transformed into S0's in the cluster environment. On the basis of the evolution of the general morphological mix in clusters our model suggests that the process responsible for the morphological transformation takes a relatively long time (1-3Gyr) after the galaxy has entered the cluster environment.


title: A Deep ROSAT HRI Observation of NGC 1313

abstract: We describe a series of observations of NGC 1313 using the ROSAT HRI with a combined exposure time of 183.5 ksec. The observations span an interval between 1992 and 1998; the purpose of observations since 1994 was to monitor the X-ray flux of SN1978K, one of several luminous sources in the galaxy. No diffuse emission is detected in the galaxy to a level of ~1-2x10^37 ergs/s/arcmin^-2. A total of eight sources are detected in the summed image within the D_25 diameter of the galaxy. The luminosities of five of the eight range from \~6x10^37 to ~6x10^38 erg/s; these sources are most likely accreting X-ray binaries, similar to sources obseved in M31 and M33. The remaining three sources all emit above 10^39 erg/s. We present light curves of the five brightest sources. Variability is detected at the 99.9% level from four of these. We identify one of the sources as an NGC 1313 counterpart of a Galactic X-ray source. The light curve, though crudely sampled, most closely resembles that of a Galactic black hole candidate such as GX339-4, but with considerably higher peak X-ray luminosity. An additional seven sources lie outside of the D_25 diameter and are either foreground stars or background AGN.


title: Detection of a broad iron emission line and sub-millisecond quasiperiodic oscillations from the type I X-ray burster 4U1728-34 in a high state

abstract: We report results from simultaneous RossiXTE and BeppoSAX observations of the neutron-star x-ray binary and type I X-ray burster 4U1728-34. The source was found in a high luminosity state, L_X near 0.1 L_Edd, and quasiperiodic oscillations were detected at 1284 +/- 6 Hz, the highest frequency ever observed from this source. The x-ray spectrum shows a broad, FWHM = 0.8 keV, iron K-alpha fluorescence line. We discuss interpretations of the broad line and the quasiperiodic oscillations and how future simultaneous spectral and timing observations can be used to test these interpretations and, potentially, to estimate the mass of the compact object.


title: What Kind of Obscured AGN Make the Cosmic X-ray Background?

abstract: While it is understood that the origin of the Cosmic X-ray Background (CXRB) can be explained in terms of a combination of obscured and unobscured AGN, the identity of these obscured AGN showing up in large numbers at faint fluxes remains unclear. Deep surveys are showing up increasingly large numbers of AGN at the faintest fluxes which have broad optical lines but flat X-ray spectra. It is commonly assumed that there must a new population of sources making a substantial contribution to the CXRB. We show that we need not look further than our own backyard to find familiar examples of such sources and there is no compelling reason to invoke a mysterious, new population.


title: A search for clusters at high redshift - II. A proto cluster around a radio galaxy at z=2.16

abstract: VLT spectroscopic observations of Ly alpha--excess objects in the field of the clumpy radio galaxy 1138-262 at z=2.16 have led to the discovery of 14 galaxies and one QSO at approximately the same distance as the radio galaxy. All galaxies have redshifts in the range 2.16 \pm 0.02, centered around the redshift of the radio galaxy, and are within a projected physical distance of 1.5 Mpc from it. The velocity distribution suggests that there are two galaxy subgroups having velocity dispersions of ~500 km/s and ~300 km/s and a relative velocity of 1800 km/s. If these are virialized structures, the estimated dynamical masses for the subgroups are ~9 and ~3x10E13 solar masses respectively, implying a total mass for the structure of more than 10E14 solar masses. The new observations, together with previous results, suggest that the structure of galaxies around 1138-262 is likely to be a forming cluster.


title: HST Observations of Black Hole X-ray Transient Outbursts

abstract: The Hubble Space Telescope has now observed four black hole X-ray transients in outburst under an ongoing Target-of-Opportunity program: GRO J0422+32 (1994), GRO J1655-40 (1996), XTE J1859+226 (1999-2000) and XTE J1118+480 (2000). We summarise the highlights of this program to date, including a diversity of spectral energy distributions, correlated X-ray/UV/optical rapid variability and detailed UV spectroscopy.


title: Doppler Tomography of XTE J2123-058 and Other Neutron Star LMXBs

abstract: We describe Doppler tomography obtained in the 1998 outburst of the neutron star low mass X-ray binary (LMXB) XTE J2123-058. This analysis, and other aspects of phase-resolved spectroscopy, indicate similarities to SW Sex systems, except that anomalous emission kinematics are seen in HeII, whilst phase 0.5 absorption is confined to H alpha. This separation of these effects may provide tighter constraints on models in the LMXB case than is possible for SW Sex systems. We will compare results for other LMXBs which appear to show similar kinematics and discuss how models for the SW Sex phenomenon can be adapted to these systems. Finally we will summarise the limited Doppler tomography performed on the class of neutron star LMXBs as a whole, and discuss whether any common patterns can yet be identified.


title: Dust Grain Size Distributions and Extinction in the Milky Way, LMC, and SMC

abstract: We construct size distributions for carbonaceous and silicate grain populations in different regions of the Milky Way, LMC, and SMC. The size distributions include sufficient very small carbonaceous grains (including polycyclic aromatic hydrocarbon molecules) to account for the observed infrared and microwave emission from the diffuse interstellar medium. Our distributions reproduce the observed extinction of starlight, which varies depending upon the interstellar environment through which the light travels. As shown by Cardelli, Clayton & Mathis in 1989, these variations can be roughly parameterized by the ratio of visual extinction to reddening, R_V. We adopt a fairly simple functional form for the size distribution, characterized by several parameters. We tabulate these parameters for various combinations of values for R_V and b_C, the C abundance in very small grains. We also find size distributions for the line of sight to HD 210121, and for sightlines in the LMC and SMC. For several size distributions, we evaluate the albedo and scattering asymmetry parameter, and present model extinction curves extending beyond the Lyman limit.


title: On the magnetism of stars and planets

abstract: A self-consistent statistical approach to the problem of planetary and stellar magnetism is suggested. The mechanism of magnetic field generation in the astronomical objects, where the existence of fields is associated with the axial rotation of objects, is discussed. In the general case different forces produce partial R-separation of the charges. As a result of the system rotation, the magnetic fields of the currents of these charges are not compensated. The influence of various factors on the magnetic field of some object is analysed.


title: A weak turbulence theory for incompressible MHD

abstract: We derive a weak turbulence formalism for incompressible MHD. Three-wave interactions lead to a system of kinetic equations for the spectral densities of energy and helicity. We find energy spectra solution of the kinetic equations. The constants of the spectra are computed exactly and found to depend on the amount of correlation between the velocity and the magnetic field. Comparison with several numerical simulations and models is also made.


title: Theoretical Models of Photodissociation Fronts

abstract: Observations of H2 line emission have revealed higher-than-expected gas temperatures in a number of photodissociation fronts. We discuss the heating and cooling processes in photodissociation regions. Observations of NGC 2023 are compared to a theoretical model in which there is substantial gas at temperatures T = 500-1000K heated by photoelectric emission and collisional deexcitation of H2. In general the model successfully reproduces the observed H2 line emission from a wide range of energy levels. The observed [SiII]34.8um emission appears to indicate substantial depletion of Si in NGC 2023.


title: Interstellar Grains

abstract: Our current understanding of interstellar dust is summarized at an introductory level. Submicron-sized interstellar dust grains absorb and scatter light, and reradiate the absorbed energy in the infrared. The grain population spans a range of sizes, from molecules containing only tens of atoms, to particles containing 10^{10} atoms. Most of the grain mass appears to be due to two types of solid, in approximately equal amounts: amorphous silicate mineral, and carbonaceous material. Approximately 2/3 of the interstellar carbon in diffuse clouds is in solid form. Interstellar grains are important for many reasons: they obscure our view (in a wavelength-dependent manner); the infrared emission from dust provides a valuable probe of dense regions; and because dust grains play a central role in interstellar chemistry, gas dynamics, and heating and cooling of interstellar gas.


title: User Guide for the Discrete Dipole Approximation Code DDSCAT (Version 5a10)

abstract: DDSCAT.5a is a freely available software package which applies the "discrete dipole approximation" (DDA) to calculate scattering and absorption of electromagnetic waves by targets with arbitrary geometries and complex refractive index. The DDA approximates the target by an array of polarizable points. DDSCAT.5a requires that these polarizable points be located on a cubic lattice. DDSCAT.5a10 allows accurate calculations of electromagnetic scattering from targets with "size parameters" 2 pi a/lambda < 15 provided the refractive index m is not large compared to unity (|m-1| < 1). The DDSCAT package is written in Fortran and is highly portable. The program supports calculations for a variety of target geometries (e.g., ellipsoids, regular tetrahedra, rectangular solids, finite cylinders, hexagonal prisms, etc.). Target materials may be both inhomogeneous and anisotropic. It is straightforward for the user to import arbitrary target geometries into the code, and relatively straightforward to add new target generation capability to the package. DDSCAT automatically calculates total cross sections for absorption and scattering and selected elements of the Mueller scattering intensity matrix for specified orientation of the target relative to the incident wave, and for specified scattering directions. This User Guide explains how to use DDSCAT.5a10 to carry out EM scattering calculations. CPU and memory requirements are described.


title: Contemporaneous Observations of Direct and Raman Scattered O VI in Symbiotic Stars

abstract: Symbiotic stars are binary systems consisting of a hot star, typically a white dwarf, and a cool giant companion. The wind from the cool star is ionized by the radiation from the hot star, resulting in the characteristic combination of sharp nebular emission lines and stellar molecular absorption lines in the optical spectrum. Most of the emission lines are readily identifiable with common ions. However, two strong, broad emission lines at $\lambda\lambda$ 6825, 7082 defied identification with known atoms and ions. In 1989 Schmid made the case that these long unidentified emission lines resulted from the Raman scattering of the O VI resonance photons at $\lambda\lambda$ 1032, 1038 by neutral hydrogen. We present contemporaneous far-UV and optical observations of direct and Raman scattered O VI lines for nine symbiotic stars obtained with the Hopkins Ultraviolet Telescope (Astro-2) and various ground-based telescopes. The O VI emission lines are present in every instance in which the $\lambda\lambda$ 6825, 7082 lines are present in support of the Schmid Raman scattering model. We calculate scattering efficiencies and interpret the results in terms of the Raman models. Additionally, we measure the flux of the Fe II fluorescence line at $\lambda$1776, which is excited by the O VI line at $\lambda$1032, and calculate the first estimates of the conversion efficiencies of this process.


title: Gas/Star Offsets in Tidal Tails

abstract: We use numerical simulations to study the development of gas/star offsets in the tidal tails of merging galaxies. These offsets are shown to be a natural consequence of the radially extended HI spatial distribution in disk galaxies, coupled with internal dissipation in the gaseous component driven by the interaction. This mechanism explains the observed gas/star offsets in interacting galaxies without invoking interactions with a hot (unseen) gaseous component.


title: Comptonization of Infrared Radiation from Hot Dust by Relativistic Jets in Quasars

abstract: We demonstrate the importance of near-infrared radiation from hot dust for Compton cooling of electrons/positrons in quasar jets. In our model, we assume that the non-thermal radiation spectra observed in OVV quasars are produced by relativistic electrons/positrons accelerated in thin shells which propagate down the jet with relativistic speeds. We show that the Comptonization of the near-IR flux is likely to dominate the radiative output of OVV quasars in the energy range from tens of keV up to hundreds of MeV, where it exceeds that produced by Comptonization of the UV radiation reprocessed and rescattered in the Broad Emission Line region. The main reason for this lies in the fact that the jet encounters the ambient IR radiation over a relatively large distance as compared to the distance where the energy density of the broad emission line light peaks. In the soft - to mid energy X-ray band, the spectral component resulting from Comptonization of the near-IR radiation joins smoothly with the synchrotron-self-Compton component, which may be responsible for the soft X-ray flux. At the highest observed gamma-ray energies, in the GeV range, Comptonization of broad emission lines dominates over other components.


title: Dynamical Mass Estimates for the Halo of M31 from Keck Spectroscopy

abstract: The last few months have seen the measurements of the radial velocities of all of the dwarf spheroidal companions to the Andromeda galaxy (M31) using the spectrographs (HIRES and LRIS) on the Keck Telescope. This paper summarises the data on the radial velocities and distances for all the companion galaxies and presents new dynamical modelling to estimate the mass of extended halo of M31. The best fit values for the total mass of M31 are between 7 and 10 x 10^{11} solar masses, depending on the details of the modelling. The mass estimate is accompanied by considerable uncertainty caused by the smallness of the dataset; for example, the upper bound on the total mass is roughly 24 x 10^{11} solar masses, while the lower bound is about 3 x 10^{11} solar masses. These values are less than the most recent estimates of the most likely mass of the Milky Way halo. Bearing in mind all the uncertainties, a fair conclusion is that the M31 halo is roughly as massive as that of the Milky Way halo. There is no dynamical evidence for the widely held belief that M31 is more massive -- it may even be less massive.


title: Triaxial Haloes and Particle Dark Matter Detection

abstract: This paper presents the properties of a family of scale-free triaxial haloes. We adduce arguments to suggest that the velocity ellipsoids of such models are aligned in conical coordinates. We provide an algorithm to find the set of conically aligned velocity second moments that support a given density against the gravity field of the halo. The case of the logarithmic ellipsoidal model -- the simplest triaxial generalisation of the familiar isothermal sphere -- is examined in detail. The velocity dispersions required to hold up the self-consistent model are analytic. The velocity distribution of the dark matter can be approximated as a triaxial Gaussian with semiaxes equal to the velocity dispersions. There are roughly twenty experiments worldwide that are searching for evidence of scarce interactions between weakly-interacting massive-particle dark matter (WIMPs) and detector nuclei. The annual modulation signal, caused by the Earth's rotation around the Sun, is a crucial discriminant between WIMP events and the background. The greatest rate is in June, the least in December. We compute the differential detection rate for energy deposited by the rare WIMP-nucleus interactions in our logarithmic ellipsoidal halo models. Triaxiality and velocity anisotropy change the total rate by up to 40 %, and have a substantial effect on the amplitude of the annual modulation signal. The overall rate is greatest, but the amplitude of the modulation is weakest, in our radially anisotropic halo models. Even the sign of the signal can be changed. Restricting attention to low energy events, the models predict that the maximum rate occurs in December, and not in June.


title: The Orbital Structure of Dark Matter Halos with Gas

abstract: With the success of the Chandra and XMM missions and the maturation of gravitational lensing techniques, powerful constraints on the orbital structure of cluster dark matter halos are possible. I show that the X-ray emissivity and mass of a galaxy cluster uniquely specify the anisotropy and velocity dispersion profiles of its dark matter halo. I consider hydrostatic as well as cooling flow scenarios, and apply the formalism to the lensing cluster CL0024+16 and the cooling flow cluster Abell 2199. In both cases, the model predicts a parameter-free velocity dispersion profile that is consistent with independent optical redshift surveys of the clusters.


title: Searching for MACHOs in Galaxy Clusters

abstract: If cluster dark matter is in the form of compact objects it will introduce fluctuations into the light curves of distant sources. Current searches for MACHOs in clusters of galaxies focus on monitoring quasars behind nearby systems. This paper considers the effect of such a compact population on the surface brightness distribution of giant gravitationally lensed arcs. As the microlensing optical depth is significant in these clusters, the expected fluctuations are substantial and are observable. Focusing on the giant arc seen in Abell 370, we demonstrate that several `extreme' events would be visible in a comparison of HST observations at two epochs. Utilizing NGST, long term monitoring should reveal a ubiquitous twinkling of brightness over the surface of the arcs.


title: A study of the shadowing of galactic cosmic rays by the sun in a quiet phase of solar activity with the Tibet air shower array

abstract: We have shown that the Sun's shadow by high energy cosmic rays moves year by year and its behavior is correlated with a time variation of the large-scale structure of the solar and interplanetary magnetic fields. The solar activity was near minimum in the period from 1996 through 1997. Using the data obtained with the Tibet air shower array, we examined the shadowing of cosmic rays by the Sun in this quiet phase of solar cycle, and found that the Sun's shadow was just in the apparent direction of the Sun, though it was observed at the position considerably away from the Sun to the south-west in the period between 1990 and 1993. It is known that the magnetic pole of equivalent solar dipole was reversed during the previous active phase, and near solar minimum the dipole was aligned with the rotating axis, preserving its N-pole on the north pole side of the Sun. This causes the solar magnetic field to shift the Sun's shadow to the east. Thus, the observed results suggest that the shift of the Sun's shadow due to the solar magnetic field was pushed back by the effect of the geomagnetic field, since the geomagnetic field always make the shadow shift to the west. We discuss the Sun's shadow observed during the period near solar minimum in 1996-1997 and compare it with the simulation results.


title: A Photometric Study of the Ages and Metallicities of Early-type Galaxies in A2218

abstract: We present deep optical and near-infrared imaging of the rich cluster A2218 at z=0.17. Our optical imaging comes from new multicolour Hubble Space Telescope WFPC2 observations in the F450W (B), F606W (V) and F814W (I) passbands. These observations are complemented by deep near-infrared, Ks-band, imaging from the new INGRID imager on the 4.2-m William Herschel Telescope. This combination provides unique high-precision multicolour optical-infrared photometry and morphological information for a large sample of galaxies in the core of this rich cluster at a lookback time of ~3Gyrs. We analyse the (B-I), (V-I) and (I-Ks) colours of galaxies spanning a range of a factor of 100 in K-band luminosity in this region and compare these with grids of stellar population models. We find that the locus of the colours of the stellar populations in the luminous (>0.5L*) early-type galaxies, both ellipticals and S0s, traces a sequence of varying metallicity at a single age. At fainter luminosities (<0.1L*), this sequence is extended to lower metallicities by the morphologically-classified ellipticals. However, the faintest S0s exhibit very different behaviour, showing a wide range in colours, including a large fraction (30%) with relatively blue colours which appear to have younger luminosity-weighted ages for their stellar populations, 2-5Gyrs. We show that the proportion of these young S0s in the cluster population is consistent with the observed decrease in the S0 population seen in distant clusters when interpreted within the framework of a two-step spectroscopic and morphological transformation of accreted spiral field galaxies into cluster S0s.


title: Superfluidity in relativistic neutron stars

abstract: The purpose of these notes is to give a brief review of superfluidity in neutron stars. After a short presentation explaining why and how superfluidity is expected in the crust and core of neutron stars, consequences on thermal evolution and rotational dynamics are discussed. The second part summarizes a formalism that has been recently developped to describe the hydrodynamics of superfluids or superconductors in the framework of general relativity. As an application, one can compute the oscillations of a two-component relativistic neutron star.


title: Stromgren Photometry from z=0 to z~1. The Method

abstract: We use rest-frame Stromgren photometry to observe clusters of galaxies in a self-consistent manner from z=0 to z=0.8. Stromgren photometry of galaxies is an efficient compromise between standard broad-band photometry and spectroscopy, in the sense that it is more sensitive to subtle variations in spectral energy distributions than the former, yet much less time-consuming than the latter. Principal Component Analysis (PCA) is used to extract maximum information from the Stromgren data. By calibrating the Principal Components using well-studied galaxies (and stellar population models), we develop a purely empirical method to detect, and subsequently classify, cluster galaxies at all redshifts smaller than 0.8. Interlopers are discarded with unprecedented efficiency (up to 100%). The first Principal Component essentially reproduces the Hubble Sequence, and can thus be used to determine the global star formation history of cluster members. The (PC2, PC3) plane allows us to identify Seyfert galaxies (and distinguish them from starbursts) based on photometric colors alone. In the case of E/S0 galaxies with known redshift, we are able to resolve the age-dust- metallicity degeneracy, albeit at the accuracy limit of our present observations. This technique will allow us to probe galaxy clusters well beyond their cores and to fainter magnitudes than spectroscopy can achieve. We are able to directly compare these data over the entire redshift range without a priori assumptions because our observations do not require k-corrections. The compilation of such data for different cluster types over a wide redshift range is likely to set important constraints on the evolution of galaxies and on the clustering process.


title: The ING Instrumentation Conference Discussion, Options for a Competitive Observatory

abstract: An expert panel initiated discussion on a number of key questions facing the role of 4-m and small telescopes in the new era of 8-m telescopes. The panel and audience agreed that the 4-m class telescope role would necessarily evolve, but would still be important in the coming years. The need for an active development programme of competitive instrumentation for 4-m class telescopes, and in particular the William Herschel Telescope (WHT) was stressed. In conjunction with this, the need to de-commission instrumentation made redundant by 8-m class telescopes was noted. New operational modes, including greater emphasis on survey programmes, and possibly queue scheduling, coupled with changes to the procedures for allocating time were seen as desirable. The panel and audience supported the Isaac Newton Group's emphasis on the development of instrumentation to exploit its imminent deployment of the WHT's facility Adaptive Optics system.


title: The Globular Cluster Omega Centauri and the Oosterhoff Dichotomy

abstract: CCD observations obtained by the OGLE team for 128 RR Lyrae variables in Omega Cen have been analysed. The P-L and P-A plots indicate that, in addition to fundamental (RRab) and first overtone (RRc) pulsators, the Omega Cen RR Lyrae population seems to include second overtone (RRe) and possibly third overtone pulsators. The mean periods derived for the RRab and RRc stars are typical values for an Oosterhoff type II (OoII) cluster. Nevertheless, the P-A plot plot also shows that some of the RR Lyrae variables have `Oosterhoff type I' (OoI) characteristics. Relative luminosities derived for the RRc variables from Fourier coefficients correlate with the observed apparent magnitudes. Masses for the RRc stars have been calculated from Fourier coefficients. A comparison of the derived masses for RRc stars in the four OoII clusters Omega Cen, M15, M55 and M68 indicates that the masses of the RRc stars in M15 and M68 are almost 0.2 solar masses greater than those in the other two. Since M15 and M68 have a high frequency of RRd stars among their first overtone pulsators, while none have been identified in Omega Cen or M55, this suggests that the double-mode pulsation phenomenon may be associated with mass. Among the RRc variables in Omega Cen, the OoII variables have lower derived masses and higher luminosities than the OoI variables. An application of the period-density law to pairs of OoI and OoII RRab stars selected according to their position in the P-A plot also indicates that the OoII variables in general have lower masses and higher luminosities. These findings support the hypothesis that RR Lyrae variables in OoII systems are evolved HB stars that spend their ZAHB phase on the blue side of the instability strip.


title: Global fits for the spectral index of the cosmological curvature perturbation

abstract: Best-fit values of the spectral index of the curvature perturbation are presented, assuming the $\Lambda$CDM cosmology. Apart from the spectral index, the parameters are the Hubble parameter, the total matter density and the baryon density. The data points are intended to represent all measurements which are likely to significantly affect the result. The cosmic microwave anisotropy is represented by the COBE normalization, and heights of the first and second peaks given by the latest Boomerang and Maxima data. The slope of the galaxy correlation function and the matter density contrast on the $8h^{-1}\Mpc$ scale are each represented by a data point, as are the expected values of the Hubble parameter and matter density. The `low-deuterium' nucleosynthesis value of the baryon density provides a final data point, the fit giving a value about one standard deviation higher. The reionization epoch is calculated from the model by assuming that it corresponds to the collapse of a fraction $f\gsim 10^{-4}$ of matter. We consider the case of a scale-independent spectral index, and also the scale-dependent spectral index predicted by running mass models of inflation. In the former case, the result is compared with the prediction of models of inflation based on effective field theory, in which the field value is small on the Planck scale. Detailed comparison is made with other fits, and other approaches to the comparison with theory.


title: Putting The Together A Cyclical Baryonic Universe

abstract: There are multiple examples of gravitational losses in neutron stars and black holes. Protons and neutrons have been found to have enormous repulsive pressures that highly squeezed collapsing matter cannot overcome. The case against singularities follows. Galactic black hole gravitational losses can supply the missing dark energy. With highly squeezed nucleons, the big bang could begin as a hot core and a cold dark matter shell. The 21 cm. radiation data has identified baryon sized particles as cold dark matter. Highly squeezed nucleons will not decompose to produce antimatter. The flatness of the Universe is due to a baryonic bounce. The highly correlated galaxies originated from primordial black holes capturing hot core gasses.There is evidence that galaxies have not grown nor merged significantly since formation.


title: Towards a refined cosmic concordance model: joint 11-parameter constraints from CMB and large-scale structure

abstract: We present a method for calculating large numbers of power spectra C_l and P(k) that accelerates CMBfast by a factor around 10^3 without appreciable loss of accuracy, then apply it to constrain 11 cosmological parameters from current Cosmic Microwave Background (CMB) and Large Scale Structure (LSS) data. While the CMB alone still suffers from several degeneracies, allowing, e.g., closed models with strong tilt and tensor contributions, the shape of the real space power spectrum of galaxies from the IRAS Point Source Catalogue Redshift (PSCz) survey breaks these degeneracies and helps place strong constraints on most parameters. At 95% confidence, the combined CMB and LSS data imply a baryon density 0.020 < omega_b < 0.037, dark matter density 0.10 < omega_d < 0.32 with a neutrino fraction f_nu < 38%, vacuum density Omega_Lambda < 0.76, curvature -0.19 < Omega_k < 0.10, scalar tilt 0.86 < n_s < 1.16 and reionization optical depth tau < 0.44. These joint constraints are quite robust, changing little when we impose priors on the Hubble parameter, tilt, flatness, gravity waves or reionization. Adding nucleosynthesis and neutrino priors on the other hand tightens parameters considerably, requiring Omega_Lambda > 0.49 and a red-tilt. The analysis allows a number of consistency tests to be made, all of which pass. At the 95% level, the flat scalar ``concordance model'' with Omega_Lambda=0.62, omega_d=0.13, omega_b=0.02, f_nu~0, n_s=0.9, tau=0.1, h=0.63 is consistent with external constraints. The inferred PSCz bias b~1.2 agrees with the value estimated independently from redshift space distortions. The inferred cosmological constant value agrees with the one independently derived from SN 1a studies. Cosmology seems to be on the right track.


title: Evidence for Polarized Synchrotron Components in Radio-Optical Aligned Quasars

abstract: A tendency for the axes of double-lobed radio quasars to be aligned with the electric vectors of optical polarization in the active galactic nuclei is well-known. However, the origin of the polarization and reason behind its correlation with radio morphology is not yet established. From accurate spectropolarimetry of 7 quasars which show this alignment effect, we find that the polarization is universally confined to the continuum, and not shared by the line emission or 3000A bump. Over the observed region 4000-8000A, the spectral indices of polarized flux, P x Flambda, are consistent with uniform polarization applied to the light of the big blue bump or with synchrotron emission. However, electron scattering from an optically-thick, geometrically- thin, accretion disk is well-known to give rise to a polarization position angle that is perpendicular to the disk axis. Optical synchrotron emission akin to that shown to exist in the miniblazar 3C273 is a more attractive explanation, and supporting evidence can be found for some of the targets in the form of polarimetric variability over intervals of years. Properties of the most strongly-polarized radio-optical aligned quasars can be explained by misdirected blazar core components that have net polarizations of P~10% and provide ~10% of the total optical light.


title: An accelerated Monte Carlo method to solve two-dimensional radiative transfer and molecular excitation

abstract: We present a numerical method and computer code to calculate the radiative transfer and excitation of molecular lines. Formulating the Monte Carlo method from the viewpoint of cells rather than photons allows us to separate local and external contributions to the radiation field. This separation is critical to accurate and fast performance at high optical depths (tau>100). The random nature of the Monte Carlo method serves to verify the independence of the solution to the angular, spatial, and frequency sampling of the radiation field. These features allow use of our method in a wide variety of astrophysical problems without specific adaptations: in any axially symmetric source model and for all atoms or molecules for which collisional rate coefficients are available. Continuum emission and absorption by dust is explicitly taken into account but scattering is neglected. We illustrate these features in calculations of (i) the HCO+ J=1-0 and 3-2 emission from a flattened protostellar envelope with infall and rotation, (ii) the CO, HCO+, CN and HCN emission from a protoplanetary disk and (iii) HCN emission from a high-mass young stellar object, where infrared pumping is important. The program can be used for optical depths up to 1000-10,000, depending on source model. We expect this program to be an important tool in analysing data from present and future infrared and (sub) millimetre telescopes.


title: Analysis of Energy Spectra with Low Photon Counts via Bayesian Posterior Simulation

abstract: Over the past 10 years Bayesian methods have rapidly grown more popular as several computationally intensive statistical algorithms have become feasible with increased computer power. In this paper, we begin with a general description of the Bayesian paradigm for statistical inference and the various state-of-the-art model fitting techniques that we employ (e.g., Gibbs sampler and Metropolis- Hastings). These algorithms are very flexible and can be used to fit models that account for the highly hierarchical structure inherent in the collection of high-quality spectra and thus can keep pace with the accelerating progress of new space telescope designs. The methods we develop, which will soon be available in the CIAO software package, explicitly model photon arrivals as a Poisson process and, thus, have no difficulty with high resolution low count X-ray and gamma-ray data. We expect these methods to be useful not only for the recently launched Chandra X-ray observatory and XMM but also new generation telescopes such as Constellation X, GLAST, etc. In the context of two examples (Quasar S5 0014+813 and Hybrid-Chromosphere Supergiant Star alpha TrA) we illustrate a new highly structured model and how Bayesian posterior sampling can be used to compute estimates, error bars, and credible intervals for the various model parameters.


title: The Size Distribution of the Kuiper Belt Objects

abstract: The size distribution is determined directly from the current known Kuiper Belt Objects (KBOs). We found that there is a peak for the size distribution around 230 km in diameter. For the objects larger than 230 km, the KBOs populate as $N(s) \propto s^{-4}$ as we see in Jewitt et al. (1998). For the objects smaller than 230 km, we found $N(s) \propto s^{3}$. This result from the current observational data is new and very different from the conventional concept of the KBO size distribution. Though the true size distribution might change from what we have after more KBOs be discovered, we argue that the existence of a peak around 150 km to 250 km is likely.


title: The Evolution and Space Density of Damped Lyman-alpha Galaxies

abstract: The results of a new spectroscopic survey of 66 $z \simgt 4$ quasars for Damped Lyman-alpha absorption systems are presented. The search led to the discovery of 30 new DLA candidates which are analysed in order to compute the comoving mass density of neutral gas in a non-zero lambda Universe. The possible sources of uncertainty are discussed and the implications of our results for the theories of galaxy formation and evolution are emphasized. A subsequent paper will present details of the calculations summarised here and a more extensive explanation of the consequences of our observations for the understanding of the nature of DLAs.


title: Sub-parsec-scale Accleration of the Radio Jet in the Powerful Radio Galaxy NGC 6251

abstract: In order to investigate the genesis of powerful radio jet, we have mapped the central 10 pc region of the nearby radio galaxy NGC 6251 with a 0.2 pc resolution using Very Long Baseline Interferometer (VLBI) at two radio frequencies, 5 GHz and 15 GHz, we have found the sub-parsec-scale counterjet for the first time in this radio galaxy. This discovery allows us to investigate the jet acceleration based on the relativistic beaming model.


title: The Planetary System of Upsilon Andromedae

abstract: The bright F8 V solar-type star upsilon Andromedae has recently been reported to have a system of three planets of Jovian masses. In order to investigate the orbital stability and mutual gravitational interactions among these extrasolar planets, both forward and backward integrations from the latest observed orbital elements for all three planets' orbits have been performed under the coplanar assumption. We reconfirm that the middle and the outer planet have strong interaction leading to large time variations in the eccentricities of these planets, which was shown by the previous studies. However, we discuss the validity of the ignorance of the innermost planet. We argue that this planetary system is likely to be stable and oscillate around current orbital elements since it was formed. We suggest that one possible way to produce these orbital elements: the innermost planet has very low eccentricity but the outermost planet has high eccentricity could be the interaction with the protostellar disc.


title: Feasibility of Reconstructing the Quintessential Potential Using SNIa Data

abstract: We investigate the feasibility of the method for reconstructing the equation of state and the effective potential of the quintessence field from SNIa data. We introduce a useful functional form to fit the luminosity distance with good accuracy (the relative error is less than 0.1%). We assess the ambiguity in reconstructing the equation of state and the effective potential which originates from the uncertainty in $\Omega_M$. We find that the equation of state is sensitive to the assumed $\Omega_M$, while the shape of the effective potential is not. We also demonstrate the actual reconstruction procedure using the data created by Monte-Carlo simulation. Future high precision measurements of distances to thousands of SNIa could reveal the shape of the quintessential potential.


title: RXTE observations of XTE J2012+381 during its 1998 outburst

abstract: We present an analysis of the RXTE observations of X-ray transient source XTE J2012+381 during its 1998 outburst. The spectral and timing properties of the source emission and their evolution during the outburst are very similar to what of the X-ray Novae.


title: Afterglows from jetted GRB remnants

abstract: We have found that the conventional generic dynamical model for gamma-ray bursts (GRBs) cannot reproduce the Sedov solution in the non-relativistic limit. Based on our refined generic dynamical model, we investigate afterglows from jetted GRB remnants numerically. Many new results are reached. For example, we find no obvious break in the optical light curve during the relativistic phase itself. But an obvious break does exist at the transition from the relativistic phase to the non-relativistic phase, which typically occurs at time 10 to 30 days. It is very interesting that the break is affected by many parameters, such as the electron energy fraction, $\xi_{\rm e}$, the magnetic energy fraction, $\xi_{\rm B}^2$, the initial opening angle of the jet,$\theta_0$, and the medium density, $n$. We also find that afterglows from jetted GRB remnants are uniformly characterized by a quick decay during the non-relativistic phase, with power law timing index greater than 2.1 . Afterglows from GRB 970228, 980326, 980519, 990123, 990510 and 991208 can be satisfactorily fitted if the corresponding GRB ejecta are highly collimated.


title: Cosmic Shear with the CFHT

abstract: We present results of our cosmic shear survey in progress at the Canada- France-Hawaii Telescope (CFHT). We analysed 1.7 sq. degrees of data, out of which we measured a significant correlation of galaxy shape-orientation over several arcmin scale. We present measurements of the variance of the shear $<\gamma^2>$ and of the correlation functions $<e_r(\theta)e_r(0)>$, $<e_t(\theta)e_t(0)>$. The signal is consistent with gravitational lensing by large scale structures predictions. The level of residual systematics after PSF correction is small compared to our measured signal. We outline several evolutions of our work by using additional data covering 9 sq. degrees.


title: Cosmic Shear with the VLT

abstract: We report on the detection of cosmic shear on angular scales of 1.3-6.5 arcmin using 45 independent empty fields observed with the Very Large Telescope (VLT). This result confirms previous measurements obtained with the CFH Telescope at the same angular scales. We present the data analysis and the first preliminary results.


title: A Numerical Study of Cosmic Shear Statistics

abstract: We explore the stability of the variance and skewness of the cosmic gravitational convergence field, using two different approaches: first we simulate a whole MEGACAM survey (100 sq. degrees). The reconstructed mass map, obtained from a shear map, shows that the state-of-the-art data analysis methods can measure weak-lensing statistics at angular scales ranging from 2.5' to 25'. We looked also at the influence of a varying signal-to-noise ratio over the shear map (due to local variations of source density) on the mass reconstruction, by means of Monte-Carlo simulation. The effect at small scales can easily be corrected-for in most of the relevant cases. These results enhance the confidence in the capability of future large surveys to measure accurately cosmologically interesting quantities.


title: Hot dust in two hard Chandra X-ray sources

abstract: The two brightest hard X-ray sources discovered serendipitously by Chandra in the field of the lensing cluster A2390 are found to have ISOCAM counterparts at 6.7 and 15 micron. We use this fact, together with their non-detection by SCUBA at 850 micron, as the basis for DUSTY radiative transfer modelling of their infrared spectral energy distributions. For the best-fitting models, we find that the dust which reprocesses the optical--UV light in these Compton-thin AGN is heated to near its sublimation temperature (above 1000K), with an inner radius within a parsec of the nucleus. Some warm dust models with inner temperatures of 200K are also formally acceptable. These findings strongly support the obscured AGN hypothesis for the new hard X-ray Chandra sources, which lack both strong emission lines and 850 micron SCUBA detections.


title: Chandra observations of NGC 253: New insights into the nature of starburst-driven superwinds

abstract: Arcsecond-resolution X-ray imaging of the nucleus of the nearby starburst galaxy NGC 253 with Chandra reveals a well-collimated, strongly limb-brightened, kiloparsec-scale conical outflow from the central starburst region. The outflow is very similar in morphology to the known H-alpha outflow cone, on scales down to <= 20 pc. This provides, for the first time, robust evidence that both X-ray and H-alpha emission come from low volume filling factor regions of interaction between the fast energetic wind of SN-ejecta and the denser ambient interstellar medium (ISM), and not from the wind fluid itself. We provide estimates of the (observationally and theoretically important) filling factor of the X-ray emitting gas, of between 4 and 40 per cent, consistent with an upper limit of ~40 per cent based directly on the observed limb-brightened morphology of the outflow. Only <= 20 per cent of the observed X-ray emission can come from the volume-filling, metal-enriched, wind fluid itself. Spatially-resolved spectroscopy of the soft diffuse thermal X-ray emission reveals that the predominant source of spectral variation along the outflow cones is due to strong variation in the absorption, on scales of < 60 pc, there being little change in the characteristic temperature of the emission. We show that these observations are easily explained by, and fully consistent with, the standard model of a superwind driven by a starburst of NGC 253's observed power. If these results are typical of all starburst-driven winds, then we do not directly see all the energy and gas (in particular the hot metal-enriched gas) transported out of galaxies by superwinds, even in X-ray emission.


title: New theoretical Mass-loss rates of O and B Stars

abstract: We have calculated mass-loss rates for a grid of wind models covering a wide range of stellar parameters and have derived a mass-loss recipe for two ranges of effective temperature at either side of the bi-stability jump around spectral type B1. For a large sample of O stars, it is shown that there is now good agreement between these new theoretical mass-loss rates that take multiple scattering into account and observations. Agreement between the observed and new theoretical wind momenta increases confidence in the possibility to derive distances to luminous stars in distant stellar systems using the Wind momentum Luminosity Relation. For the winds of the B stars there is an inconsistency in the literature between various mass-loss rate determinations from observations by different methods. One group of the determinations for B stars does follow the new theoretical relation, while another group does not. The lack of agreement between the observed mass-loss rates derived by different methods may point to systematic errors in mass-loss determinations from observations for B stars. We show that our theoretical mass-loss recipe is reliable and recommend it be used in evolutionary calculations.


title: The Edinburgh/Durham Southern Galaxy Catalogue - IX. The Galaxy Catalogue

abstract: We announce here the public availability of the Edinburgh/Durham Southern Galaxy Catalogue (EDSGC, http://www.edsgc.org). This objective galaxy catalogue was constructed using the COSMOS micro-densitometer at the Royal Observatory Edinburgh and constitutes one of the largest digitized galaxy surveys currently in existence. The EDSGC contains a total of 1,495,877 galaxies (each with 27 image parameters) covering an contiguous area of 1182 sq deg centered on the South Galactic Pole. The data consists of photographic bj magnitudes calibrated via CCD sequences which provide a plate-to-plate accuracy of bj~0.1. Extensive external checks have demonstrated that the global EDSGC photometry is free of large-scale systematic gradients and is therefore, ideal for studying the distribution of galaxies on large angular scales. Independent spectroscopy of EDSGC galaxies has shown that the accuracy of the star-galaxy separation is consistent with earlier visual checks and that only 12% of EDSGC galaxies are potentially mis-classified. This paper is intended to provide a summary of the essential details of the catalogue's design and construction as well as provide a brief summary of the main scientific achievements using the EDSGC over the last ten years.


title: The Age of the Galaxy from Thorium Cosmochronometry

abstract: The long-lived radioactive element thorium can potentially provide a simple and direct clock to determine the age of our Galaxy. Spectroscopic investigations of thorium in metal-poor stars have yielded some promising initial results. We discuss the major observational and theoretical limitations in thorium cosmochronometry, and point out the ways in which the implied Galactic ages from thorium abundances can be made more accurate.


title: Visualization of Large Multi-Dimensional Datasets

abstract: Visualization techniques are well developed for many problem domains, but these systems break down for datasets which are very large or multidimensional. Techniques for data which is discrete rather than continuous are also less well studied. Astronomy datasets like the Sloan Digital Sky Survey are very much in this category. We propose the extension of information visualization techniques to these very large record-oriented datasets. Specifically, we describe the possible adaptation of the ``Visage'' information visualization tool to terabyte astronomy datasets.


title: Fast Algorithms and Efficient Statistics: Density Estimation in Large Astronomical Datasets

abstract: In this paper, we outline the use of Mixture Models in density estimation of large astronomical databases. This method of density estimation has been known in Statistics for some time but has not been implemented because of the large computational cost. Herein, we detail an implementation of the Mixture Model density estimation based on multi-resolutional KD-trees which makes this statistical technique into a computationally tractable problem. We provide the theoretical and experimental background for using a mixture model of Gaussians based on the Expectation Maximization (EM) Algorithm. Applying these analyses to simulated data sets we show that the EM algorithm - using the AIC penalized likelihood to score the fit - out-performs the best kernel density estimate of the distribution while requiring no ``fine--tuning'' of the input algorithm parameters. We find that EM can accurately recover the underlying density distribution from point processes thus providing an efficient adaptive smoothing method for astronomical source catalogs. To demonstrate the general application of this statistic to astrophysical problems we consider two cases of density estimation: the clustering of galaxies in redshift space and the clustering of stars in color space. From these data we show that EM provides an adaptive smoothing of the distribution of galaxies in redshift space (describing accurately both the small and large-scale features within the data) and a means of identifying outliers in multi-dimensional color-color space (e.g. for the identification of high redshift QSOs). Automated tools such as those based on the EM algorithm will be needed in the analysis of the next generation of astronomical catalogs (2MASS, FIRST, PLANCK, SDSS) and ultimately in in the development of the National Virtual Observatory.


title: BOOMERanG Data Suggest a Purely Baryonic Universe

abstract: The amplitudes of peaks in the angular power spectrum of anisotropies in the microwave background radiation depend on the mass content of the universe. The second peak should be prominent when cold dark matter is dominant, but is depressed when baryons dominate. Recent microwave background data are consistent with a purely baryonic universe with Omega(matter) = Omega(baryon) ~ 0.03 and Omega(Lambda) ~ 1.


title: On a Probable Contribution of the Einstein-Podolsky-Rosen Effect to Solving the Problem of Homogeneity of the Universe

abstract: Evading formation of the domain walls in cosmological phase transitions is one of key problems to be solved for explanation of the observed large-scale homogeneity of the Universe. The previous attempts to get around this obstacle led to imposing severe observational constraints on the parameters of the fields involved. Aim of the present report is to show that yet another way to overcome the above problem is accounting for Einstein-Podolsky-Rosen (EPR) effect. Namely, if the scalar (Higgs) field was presented by a single coherent quantum state at the initial instant of time, then its reduction during a phase transition at some later instant should be correlated even at distances exceeding the local cosmological horizon. Consideration of a simplest one-dimensional cosmological model with Z_2 Higgs field demonstrates that EPR correlations can substantially reduce the probability of spontaneous creation of the domain walls, and such effect seems to be even more promising in the cases of two- and three-dimensional geometry.


title: Influence of Cooled Interstellar Gas on the Fundamental Plane for Elliptical Galaxies

abstract: We explore the possibly important influence of cooled interstellar gas on the fundamental plane of elliptical galaxies. Interstellar cooling is described by a parameterized sink term in the equation of continuity. Parameters that give the best fits to the X-ray observations of NGC 4472 are used as a template for the radial distribution of interstellar cooling in structurally homologous elliptical galaxies of lower mass. Gas that cools within an effective radius can contribute an additional 10 - 30 percent to the mass of the old stellar population. If the cooled gas forms into stars of very low mass, $\ll M_{\odot}$, as is commonly assumed, the cooled mass is optically dark. As a result, the mass to light ratios determined from stellar velocities systematically overestimate that of the old stellar population. Moreover, the total mass and spatial distribution of the optically dark young stellar population does not scale homologously with galactic luminosity or radius and the total stellar mass to light ratio varies with galactic radius. We investigate the non-homologous perturbations of cooled gas on the mass to light ratio for several idealized homologous elliptical galaxies and show that they appear to be incompatible with the observed thinness of the fundamental plane. If optically luminous young stars formed from the cooled gas, the disturbance of the fundamental plane would be lessened.


title: Young massive star clusters in M51

abstract: A search for young massive star clusters (YMCs) in the nearby face-on spiral galaxy M51 (NGC 5194) has been carried out using UBV CCD images from the prime focus camera on the Lick 3 meter Shane telescope. The YMC population is found to be quite rich with a specific U-band luminosity T_L(U) = 1.4, consistent with the high current star formation rate of this galaxy. The brightest clusters have M_V about -12.5, far brighter than any young clusters currently known in the Milky Way and even surpassing the luminosity of the R136 cluster in the 30 Dor complex in the Large Magellanic Cloud. A few of the YMCs are examined on archive HST/WFPC2 images, confirming their cluster nature and providing estimates of their effective radii of 2-3 pc. The number of YMCs in M51 is compatible with extrapolation of a power-law luminosity function with exponent roughly equal to -2 from a Milky Way-like population of open clusters. Both the SFR and T_L(U) of M51 are similar to those of other cluster-rich spiral galaxies like NGC 1313 and M83.


title: The Giant Radio Galaxy MRC B0319-454: Circumnuclear Structure of the Host Galaxy ESO248-G10

abstract: We present optical and near-infrared images and spectra of ESO248-G10, the host of giant radio galaxy MRC B0319-454. From near-infrared colours, the active nucleus is shown to be reddened by hot dust emission or dust extinction. Star forming regions are identified beyond a radius of 5" (8.1 kpc) where hot gas is prevalent. The optical spectrum shows gas rotation at speeds of up to 350km/s out to >=32kpc along the radio axis. A model is proposed in which the giant elliptical is triaxial with the radio axis along the short axis, and the figure rotation is around the long axis. From the model, the angles of the principal axes are psi=34 degrees, theta=65 degrees, and phi=19 degrees with axis ratios of B/A=0.75 and C/A=0.69. A second velocity feature from 5" to 15" to the north-east of the nucleus is proposed to be a merging gas-rich galaxy inducing star formation while settling into an orbit about the figure rotation axis.


title: Linear Coasting in Cosmology and SNe Ia

abstract: A strictly linear evolution of the cosmological expansion scale factor is a characteristic feature in several classes of alternative gravity theories as also in the standard (big-bang) model with specially chosen equations of state of matter. Such an evolution has no free parameters as far as the classical cosmological tests are concerned and should therefore be easily falsifiable. In this article we demonstrate how such models present very good fits to the current supernovae 1a data. We discuss the overall viability of such models.


title: Are X-ray properties of loose groups different from those of compact groups?

abstract: We compare the X-ray properties of loose and compact galaxy groups, using a combined sample of 42 groups. We find that we are unable to separate loose and compact groups on the luminosity-temperature relation, the luminosity-velocity dispersion relation or the velocity dispersion-temperature relation using equally weighted errors. This suggests that the distinction between compact and loose groups is not a fundamental one, and we argue that a more useful distinction is that between X-ray bright and X-ray faint systems. Given their similarity in X-ray properties, we combine the loose and compact subsamples to derive relations based on the full sample. This provides the highest statistical quality results to date on the way in which the correlations in X-ray properties of low mass systems depart from those seen in rich clusters.


title: Evolution of Galaxy morphologies in Clusters

abstract: We have studied the evolution of galaxian morphologies from ground-based, good-seeing images of 9 clusters at z=0.09-0.25. The comparison of our data with those relative to higher redshift clusters (Dressler et al. 1997) allowed us to trace for the first time the evolution of the morphological mix at a look-back time of 2-4 Gyr, finding a dependence of the observed evolutionary trends on the cluster properties.


title: Detection of a Compact X-ray Source in the Supernova Remnant G29.6+0.1: A Variable Anomalous X-ray Pulsar?

abstract: We present follow-up observations of the serendipitously discovered 7-s X-ray pulsar AX J1845-0258, which displays characteristics similar to those observed in the anomalous X-ray pulsars (AXPs). We find a dramatic reduction in its 3-10 keV flux in both new ASCA and RXTE datasets. Within the pulsar's position-error locus, we detect a faint point source, AX J184453-025640, surrounded by an arc of diffuse X-ray emission. This arc is coincident with the South-East quadrant of the radio shell of the newly discovered supernova remnant G29.6+0.1, reported in our companion paper (Gaensler et al. 1999). Lack of sufficient flux from the source prevents us from confirming the 7-s pulsed emission observed in the bright state; hence, at present we cannot definitively resolve whether AX J1845-0258 and AX J184453-025640 are one and the same. If they are the same, then the peak-to-peak luminosity changes recorded for AX J1845-0258 may be larger than seen in other AXPs; closer monitoring of this pulsar might lead to a resolution on the mechanism that drives AXPs.


title: The Current Status of CLASS

abstract: I give a brief overview of the current status of some aspects of the Cosmic Lens All-Sky Survey (CLASS): description of the survey, current list of lens systems, cosmological parameters from lensing statistics, H_0 from time delays, dark lenses, wide-separation lenses.


title: The Lens-Redshift Test Revisited

abstract: Kochanek (1992) suggested that the redshifts of gravitational lens galaxies rule out a large cosmological constant. This result was questioned by Helbig & Kayser (1996), who pointed out that selection effects related to the brightness of the lens can bias the results of this test against a high lambda value; however, we did not claim that the observations favoured a high lambda value, merely that current observational data were not sufficient to say either way, using the test as proposed by Kochanek (1992) but corrected for selection effects. Kochanek (1996) pointed out that additional information (fraction of measured lens redshifts) provides additional information which restores the sensitivity of the test to the cosmological model, at least somewhat. Here, I consider three aspects. First, I examine the accuracy of the correction to the test proposed by Kochanek (1996). Second, I compare the slightly different statistical methods which have been used in connection with this test. Third, I discuss what results can be obtained today now that more and better-defined observations are available.


title: Isolated Neutron Stars in the Galaxy

abstract: In this article we briefly review our recent results on evolution and properties of isolated neutron stars (INSs) in the Galaxy. As the first step we calculate a {\it census} of INSs in our Galaxy. We infer a lower bound for the mean kick velocity of NSs, $ <V> \sim $(200-300) ${\rm km s^{-1}}$. The same conclusion is reached for both a constant magnetic field ($B\sim 10^{12}$ G) and for a magnetic field decaying exponentially with a timescale $\sim 10^9$ yr. These results, moreover, constrain the fraction of low velocity stars, which could have escaped pulsar statistics, to $\sim$few percents. Then we show that the range of minimum value of magnetic moment, $\mu_b$: $\sim 10^{29.5}\ge \mu_b \ge 10^{28} {\rm G} {\rm cm}^3$, and the characteristic decay time, $t_d$: $\sim 10^8\ge t_d \ge 10^7 {\rm yrs}$, can be excluded assuming the standard initial magnetic momentum, $\mu_0=10^{30} {\rm G} {\rm cm}^3$, if accreting INSs are observed. For these parameters an INS would never reach the stage of accretion from the interstellar medium even for a low space velocity of the star and high density of the ambient plasma. The range of excluded parameters increases for lower values of $\mu_0$. It is shown that old accreting INSs become more abundant than young cooling INSs at X-ray fluxes below $\sim 10^{-13}$ erg cm$^{-2}$ s$^{-1}$. We can predict that about one accreting INS per square degree should be observed at the {\it Chandra} and {\it Newton} flux limits of $\sim 10^{-16}$ erg cm$^{-2}$ s$^{-1}.$ The weak {\it ROSAT} sources, associated with INSs, can be young cooling objects, if the NSs birth rate in the solar vicinity during the last $\sim 10^6$ yr was much higher than inferred from radiopulsar observations.


title: Deep Halpha imagery of the Eridanus shells

abstract: A deep \ha image of interlocking filamentary arcs of nebulosity has been obtained with a wide-field ($\approx$ 30\degree diameter) narrow-band filter camera combined with a CCD as a detector. The resultant mosaic of images, extending to a galactic latitude of 65$^{o}$, has been corrected for field distortions and had galactic coordinates superimposed on it to permit accurate correlations with the most recent H{\sc i} (21 cm), X-ray (0.75 kev) and FIR (IRAS 100 $\mu$m) maps. Furthermore, an upper limit of 0.13 arcsec/yr to the expansion proper motion of the primary 25\degree long nebulous arc has been obtained by comparing a recent \ha image obtained with the San Pedro Martir telescope of its filamentary edge with that on a POSS E plate obtained in 1951. It is concluded that these filamentary arcs are the superimposed images of separate shells (driven by supernova explosions and/or stellar winds) rather than the edges of a single `superbubble' stretching from Barnard's Arc (and the Orion Nebula) to these high galactic latitudes. The proper motion measurement argues against the primary \ha emitting arc being associated with the giant radio loop (Loop 2) except in extraordinary circumstances.


title: Galaxy Groups at Intermediate Redshift

abstract: Galaxy groups likely to be virialized are identified within the CNOC2 intermediate redshift galaxy survey using an iterative method. The number-velocity dispersion relation is in agreement with the low-mass extrapolation of the cluster normalized Press-Schechter function. The two-point group-group correlation function has r_0=6.8+/- 0.3 Mpc, which is larger than the correlations of individual galaxies at the level predicted from n-body calibrated halo clustering. The groups are stacked in velocity and position to create a sample large enough for measurement of a density and velocity dispersion profile. The stacked mean galaxy density profile falls nearly as a power law with r^{-2.5} and has no well-defined core. The projected velocity dispersion is examined for a variety of samples with different methods and found to be either flat or slowly rising outwards. The combination of a steeper-than-isothermal density profile and the outward rising velocity dispersion implies that the mass-to-light ratio of groups rises with radius. The M/L can be kept nearly constant if the galaxy orbits are nearly circular, although such strong tangential anisotropy is not supported by other evidence. The segregation of mass and light is not dependent on galaxy luminosity but is far more prominent in the red galaxies than the blue. The M/L gradient could arise from orbital ``sloshing'' of the galaxies in the group halos, dynamical friction acting on the galaxies in a background of ``classical'' collisionless dark matter, or, more speculatively, the dark matter may have a true core.


title: Infrared Helium-Hydrogen Line Ratios as a Measure of Stellar Effective Temperature

abstract: We have observed a large sample of compact planetary nebulae in the near-infrared to determine how the 2^1P-2^1S HeI line at 2.058um varies as a function of stellar effective temperature, Teff. The ratio of this line with HI Br g at 2.166um has often been used as a measure of the highest Teff present in a stellar cluster, and hence on whether there is a cut-off in the stellar initial mass function at high masses. However, recent photoionisation modelling has revealed that the behaviour of this line is more complex than previously anticipated. Our work shows that in most aspects the photoionisation models are correct. In particular, we confirm the weakening of the 2^1P-2^1S as Teff increases beyond 40000K. However, in many cases the model underpredicts the observed ratio when we consider the detailed physical conditions in the individual planetary nebulae. Furthermore, there is evidence that there is still significant 2^1P-2^1S HeI line emission even in the planetary nebulae with very hot (Teff>100000K) central stars. It is clear from our work that this ratio cannot be considered as a reliable measure of effective temperature on its own.


title: The Structure and Emission of the Accretion Shock in T Tauri Stars II: the Ultraviolet Continuum Emission

abstract: We compare accretion shock models with optical and ultraviolet spectra of pre-main sequence stars to (1) make the first determinations of accretion rates in intermediate mass T Tauri stars from continuum emission, and (2) derive improved estimates of accretion rates and extinctions for continuum T Tauri stars. Our method extends the shock models developed by Calvet and Gullbring to enable comparisons with optical and archival International Ultraviolet Explorer ultraviolet spectra. We find good agreement between the observations and the model predictions, supporting the basic model of magnetospheric accretion shocks as well as previous determinations of accretion rates and interstellar reddening for the low mass T Tauri stars. The accretion rates determined for the intermediate-mass T Tauri stars agree well with values obtained through other methods that use near-infrared hydrogen line strengths (Muzerolle et al. 1998).


title: The Disk and Dark Halo Mass of the Barred Galaxy NGC 4123. I. Observations

abstract: The non-circular streaming motions in barred galaxies are sensitive to the mass of the bar and can be used to lift the degeneracy between disk and dark matter halo encountered when fitting axisymmetric rotation curves of disk galaxies. In this paper, we present photometric and kinematic observations of NGC 4123, a barred galaxy of modest size (V_rot = 130 km/sec, L = 0.7 L_*), which reveal strong non-circular motions. The bar has straight dust lanes and an inner Lindblad resonance. The disk of NGC 4123 has no sign of truncation out to 10 scale lengths, and star-forming regions are found well outside R_25. A Fabry-Perot H-alpha velocity field shows velocity jumps of >100 km/sec at the location of the dust lanes within the bar, indicating shocks in the gas flow. VLA observations yield the velocity field of the H I disk. Axisymmetric mass models yield good fits to the rotation curve outside the bar regionfor disk I-band M/L of 2.25 or less, and dark halos with either isothermal or power-law profiles can fit the data well. In a companion paper, we model the full 2-D velocity field, including non-circular motions, to determine the stellar M/L and the mass of the dark halo.


title: The Disk and Dark Halo Mass of the Barred Galaxy NGC 4123. II. Fluid-Dynamical Models

abstract: We report a dynamical determination of the separate contributions of disk and dark halo masses to the rotation curve of a spiral galaxy. We use fluid-dynamical models of gas flow in the barred galaxy NGC 4123 to constrain the dynamical properties of the galaxy: disk M/L, bar pattern speed, and the central density and scale radius of the dark halo. We derive a realistic barred potential directly from the light distribution. For each model we assume a value of the stellar M/L and a bar pattern speed Omega_p and add a dark halo to fit the rotation curve. We then compute the gas flow velocities with a 2-D gas dynamical code, and compare the model flow patterns to a 2-D velocity field derived from Fabry-Perot observations. The strong shocks and non-circular motions in the observed gas flow require a high stellar M/L and a fast-rotating bar. Models with I-band disk M/L of 2.0 -- 2.5 h_75, or 80 -- 100% of the maximum disk value, are highly favored. The corotation radius of the bar must be <= 1.5 times the bar semi-major axis. These results contradict some recent claimed ``universal'' galaxy disk/halo relations, since NGC 4123 is of modest size (rotation curve maximum 145 km/sec, and V_flat = 130 km/sec) yet is quite disk-dominated. The dark halo of NGC 4123 is less concentrated than favored by current models of dark halos based on cosmological simulations. Since some 30% of bright disk galaxies are strongly barred and have dust lanes indicating shock morphology similar to that of NGC 4123, it is likely that they also have high stellar M/L and low density halos. We suggest that luminous matter dominates inside the optical radius R_25 of high surface brightness disk galaxies.


title: A Measurement of the Average Longitudinal Development Profile of Cosmic Ray Air Showers Between 10^{17}eV and 10^{18}eV

abstract: The average extensive air shower longitudinal development profile as a function of shower age in the energy range from 10$^{17}$ to 10$^{18}$eV is measured using data from the hybrid HiRes/MIA experiment. An angular bin signal based Cerenkov light component subtraction method and a shower maximum fitting method using a local parabolic function are used to correct and normalize the data. The Gaisser-Hillas and Greisen functions work equally well for describing the shower profile. The Gaussian function is a poor fit. A simple scale-free function is proposed and fits the data equally well. The best-fit parameters for the above functions are determined.


title: Structured Mixed Phase is Favored in Neutron Stars

abstract: We give a general thermodynamical argument showing that in neutron stars, the Coulomb structured mixed phase is always favored for any first order phase transition involving systems in equilibrium with baryon number and electric charge as the two independent components. This finding is likely to have important consequences for many neutron star properties, e.g., glitch phenomena, transport and superfluid properties, r-mode instabilities, and the braking index.


title: Emergence of power-law correlation in 1-dimensional self-gravitating system

abstract: A new phase of temporal evolution of the one-dimensional self-gravitating system is numerically discovered. Fractal structure is dynamically created from non-fractal initial conditions. Implication to astrophysics and mathematical physics is discussed.


title: Acceleration of Particles by Oblique Shocks and Cosmic Ray Spectra around the Knee Region

abstract: We examine the first order Fermi acceleration on the presumption that supernova remnant shocks cross ambient magnetic fields with various angles. These oblique shocks accelerate particles more efficiently than the parallel shocks and elevate the maximum energies achievable by the particles. The primary cosmic ray spectrum is strongly dependent upon these energies. We also consider the dependence of the injection efficiency and of spectral indices on obliquity. When indices and absolute fluxes at $10^{12}$ eV are given for six nuclear groups from balloon-borne data, each energy spectrum develops a smooth rigidity dependent knee structure. The resultant total spectrum also behaves similarly and fits well with ground-based experimental data up to several $10^{17}$ eV. It is shown as well that the chemical composition changes significantly from lighter to heavier nuclei as the energies of particles exceed the knee region. Other predicted curves are compared with the experimental data which they reproduce rather well.


title: Interpretation of Extreme Scattering Events

abstract: Extreme Scattering Events are sometimes manifest in the light-curves of compact radio-quasars at frequencies of a few GHz. These events are not understood. The model which appears to offer the best explanation requires a new population of AU-sized, neutral gas clouds; these clouds would then make up a large fraction of the Galaxy's dark matter. Independent of the question of which theoretical model is correct, if we extrapolate the observed behaviour to low radio-frequencies, we expect that the sky should be criss-crossed by a network of narrow caustics, at frequencies below about 700 MHz. Consequently at these frequencies sources should typically manifest additional, faint images which are substantially delayed with respect to the primary image. Although some examples of this type of behaviour are already known, it is expected that these are just the tip of the iceberg, with strong selection biases having been imposed by the instrumentation employed to date.


title: The TeV Energy Spectrum of Mkn 501 Measured with the Stereoscopic Telescope System of HEGRA during 1998 and 1999

abstract: During 1997, the BL Lac object Mkn 501 went into an extraordinary state of high X-ray and TeV gamma-ray activity, lasting more than 6 months. In this paper we report on the TeV emission characteristics of the source in the subsequent years of 1998 and 1999 as measured with the Stereoscopic Cherenkov Telescope System of HEGRA (La Palma, Canary Islands). Our observations reveal a 1998-1999 mean emission level at 1 TeV of 1/3 of the flux of the Crab Nebula, a factor of 10 lower than during the year of 1997. A dataset of 122 observations hours with the HEGRA telescope system makes it possible to assess for the first time the Mkn 501 TeV energy spectrum for a mean flux level substantially below that of the Crab Nebula with reasonable statistical accuracy. Excluding the data of a strong flare, we find evidence that the 1998--1999 low-flux spectrum is substantially softer (by 0.44+-0.1(stat) in spectral index) than the 1997 time averaged spectrum. The 500 GeV to 10 TeV energy spectrum can well be described by a power law model with exponential cutoff: dN/dE ~ E^(-alpha) exp(-E/E0) with alpha=2.31+-0.22(stat), and E0=5.1 (-2.3+7.8)(stat) TeV. Within statistical accuracy, also a pure power law model gives an acceptable fit to the data: dN/dE ~ E^(-Gamma) with Gamma=2.76+-0.08(stat). After presenting the 1998-1999 TeV characteristics of the source we discuss the implications of the results.


title: Vortex solitons - Mass, Energy and Angular momentum bunching in relativistic electron-positron plasmas

abstract: It is shown that the interaction of large amplitude electromagnetic waves with a hot electron-positron (e-p) plasma (a principal constituent of the universe in the MeV epoch) leads to a bunching of mass, energy, and angular momentum in stable, long-lived structures. Electromagnetism in the MeV epoch, then, could provide a possible route for seeding the observed large-scale structure of the universe.


title: CO band emission from MWC349: I. First overtone bands from a disk or from a wind?

abstract: We observed the near infrared emission in the wavelength range 2.28-2.50 micron from the peculiar B[e]-star MWC349. The spectra contain besides the strong IR continuum the first overtone CO bands and most of the hydrogen recombination lines of the Pfund series, both in emission. We also modeled the spectra. The Pfund lines have a gaussian profile with a FWHM of ~100 km/s, and it turned out that their emission is in LTE and optically thin. To explain the CO bands, several scenarios were investigated. We found that the CO band heads are formed under LTE and that the gas must have a temperature of 3500 to 4000 K. The width of the 2 -> 0 band head indicates kinematical broadening of 50 to 60 km/s. We can obtain fits to the measured spectra assuming that the CO gas has a column density of 5 10^{20} cm^{-2} and is located either at the inner edge of the rotating circumstellar disk. In this case, the disk must have a bulge which partly blocks the radiation so that the observer sees only a sector on the far side where the radial velocities are small. Or the CO emission originates in a wind with gaussian line profiles. Both fits are of equal quality and satisfactory. In a third alternative where the fit is less convincing, the CO emission is optically thin and comes from an extended Keplerian disk.


title: Particle-Mesh Simulations of the Lyman-Alpha Forest

abstract: Numerical hydrodynamical simulations have proven a successful means of reproducing many of the statistical properties of the Lyman-Alpha forest as measured in high redshift quasar spectra. Pseudo-hydrodynamical methods based only on simulating the dark matter component have been claimed to yield a comparable level of success. We investigate the degree to which two pseudo-methods, with and without allowing for a pseudo-gas pressure, are able to match the predictions of fully hydrodynamical plus dark matter simulations. We also address the requirements for convergence to the statistics of the spectra and the inferred properties of the Lyman-Alpha forest as a function of resolution and box size. Generally we find it is possible to reach agreement with full hydrodynamic simulations at the 10% level in the cumulative distributions of the flux and absorption line parameter statistics for readily achievable particle and grid numbers, but difficult to do much better.


title: Evolution of Buoyant Bubbles in M87

abstract: The morphology of the X-ray and radio emitting features in the central $\sim$ 50 kpc region around the galaxy M87 strongly suggests that buoyant bubbles of cosmic rays (inflated by an earlier nuclear active phase of the galaxy) rise through the cooling gas at roughly half the sound speed. In the absence of strong surface tension, initially spherical bubbles will transform into tori as they rise through an external medium. Such structures can be identified in the radio images of the halo of M87. During their rise, bubbles will uplift relatively cool X-ray emitting gas from the central regions of the cooling flow to larger distances. This gas is colder than the ambient gas and has a higher volume emissivity. As a result, rising ``radio'' bubbles may be trailed by elongated X-ray features as indeed is observed in M87. We performed simple hydrodynamic simulations to qualitatively illustrate the evolution of buoyant bubbles in the M87 environment.


title: Modelling the Extinction Properties of Galaxies

abstract: Recently (Granato, Lacey, Silva et al. 2000, astro-ph/0001308) we have combined our spectrophotometric galaxy evolution code which includes dust reprocessing (GRASIL, Silva et al. 1998) with semi-analytical galaxy formation models (GALFORM, Cole et al. 1999). One of the most characteristic features of the former is that the dust is divided in two main phases: molecular cloud complexes, where stars are assumed to be born, and the diffuse interstellar medium. As a consequence, stellar populations of different ages have different geometrical relationships with the two phases, which is essential in understanding several observed properties of galaxies, in particular those undergoing major episodes of star formation at any redshift. Indeed, our merged GRASIL+GALFORM model reproduces fairly well the SEDs of normal spirals and starbursts from the far-UV to the sub-mm and their internal extinction properties. In particular in the model the observed starburst attenuation law (Calzetti 1999) is accounted for as an effect of geometry of stars and dust, and has nothing to do with the optical properties of dust grains.


title: Self-similar spherical collapse with non-radial motions

abstract: We derive the asymptotic mass profile near the collapse center of an initial spherical density perturbation, $\delta \propto M^{-\epsilon}$, of collision-less particles with non-radial motions. We show that angular momenta introduced at the initial time do not affect the mass profile. Alternatively, we consider a scheme in which a particle moves on a radial orbit until it reaches its turnaround radius, r_*. At turnaround the particle acquires an angular momentum $L={\cal L} \sqrt{GM_* r_*}$ per unit mass, where M_* is the mass interior to r_*. In this scheme, the mass profile is $M\propto r^{3/(1+3\epsilon)}$ for all $\epsilon >0$, in the region $r/r_t\ll {\cal L}$, where r_t is the current turnaround radius. If ${\cal L} \ll 1$ then the profile in the region ${\cal L} \ll r/r_t \ll$ is $M\propto r$ for $\epsilon <2/3$. The derivation relies on a general property of non-radial orbits which is that ratio of the pericenter to apocenter is constant in a force field k(t) r^{n} with k(t) varying adiabatically.


title: Stellar Populations and Surface Brightness Fluctuations: New Observations and Models

abstract: We examine the use of surface brightness fluctuations (SBF) for both stellar population and distance studies. New V-band SBF data are reported for five Fornax cluster galaxies and combined with literature data to define a new V-band SBF distance indicator. We use new stellar population models, based on the latest Padua isochrones transformed empirically to the observational plane, to predict SBF magnitudes and integrated colours for a wide range of population ages and metallicities. We examine the sensitivity of the predictions to changes in the isochrones, transformations, and IMF. The new models reproduce the SBF data for globular clusters fairly well, especially if higher metallicity globulars are younger. The models also give a good match to the "fluctuation colors" of elliptical galaxies. In order to obtain theoretical calibrations of the SBF distance indicators, we combine our single-burst models into composite population models. These models reproduce the observed behavior of the SBF magnitudes as a function of stellar population parameters, including the steep colour dependence found for HST/WFPC2 F814W SBF data. Because the theoretical SBF calibrations are fairly sensitive to uncertain details of stellar evolution, the empirical calibrations are more secure. However, the sensitivity of SBF to these finer details potentially makes it a powerful constraint for stellar evolution and population synthesis. [abbridged]


title: RXTE observations of the Seyfert galaxy NGC 5506: evidence for reflection from disk and torus

abstract: We report on observations of the narrow line Seyfert galaxy NGC 5506 with the Rossi X-ray TIming Explorer. The observations cover a time interval of ~1000 days during which the source showed strong flux and spectral variability. The spectrum clearly shows iron fluorescence emission at 6.4 keV and significant reflection features. Both the equivalent width of the iron line and the relative strength of the reflected continuum are higher during low flux states. The variability of the reflection features can be explained by the presence of two reflected components: one that responds rapidly to flux changes of the primary continuum and a second, slowly variable, component originating from a distant reflector, e.g. a molecular torus.


title: On the large-scale structure of galactic disks

abstract: The disk positions for galaxies of various morphological and nuclear-activity types (normal galaxies, QSO, Sy, E/S0, low-surface brightness galaxies, etc.) on the "central surface brightness - exponential scalelength" plane are considered. The stellar disks are shown to form a single sequence on this plane over a wide range of surface brightnesses (\mu_0(I)=12-25) and sizes (h=10 pc - 100 kpc). The existence of this observed sequence can probably be explaned by a combination of three factors: a disk-stability requirement, a limited total disk luminosity, and observational selection. The model by Mo et al. (1998) for disk formation in the CDM hierarchical-clustering scenario is shown to satisfactorily reproduce the salient features of the galaxy disk distribution on the "\mu_0 - h" plane.


title: Physical parameters of erupting Luminous Blue Variables: NGC 2363-V1 caught in the act

abstract: A quantitative study of the Luminous Blue Variable NGC 2363-V1 in the Magellanic galaxy NGC 2366 (D = 3.44 Mpc) is presented, based on ultraviolet and optical HST/STIS spectroscopy. Contemporary WFPC2 and WHT imaging reveals a modest V-band brightness increase of ~ 0.2 mag per year between 1996 January and 1997 November, reaching V=17.4 mag, corresponding to Mv=-10.4 mag. Subsequently, V1 underwent a similar decrease in V-band brightness, together with a UV brightening of 0.35 mag from 1997 November to 1999 November. The optical spectrum of V1 is dominated by H emission lines, with Fe II, He I and Na I also detected. In the ultraviolet, a forest of Fe absorption features and numerous absorption lines typical of mid-B supergiants are observed. From a spectral analysis with the non-LTE, line-blanketed code of Hillier & Miller (1998), we derive stellar parameters of T*=11kK, R*=420Ro, log(L/Lo)=6.35 during 1997 November, and T*=13kK, R*=315Ro, log(L/Lo)=6.4 for 1999 July. The wind properties of V1 are also exceptional, with Mdot ~ 4.4 x 10e-4 Mo/yr and $v_{\infty} ~ $300 km/s, allowing for a clumped wind (filling factor = 0.3), and assuming H/He ~ 4 by number. The presence of Fe lines in the UV and optical spectrum of V1 permits an estimate of the heavy elemental abundance of NGC 2363 from our spectral synthesis. Although some deficiencies remain, allowance for charge exchange reactions in our calculations supports a SMC-like metallicity, that has previously been determined for NGC 2363 from nebular oxygen diagnostics. Considering a variety of possible progenitor stars, V1 has definitely undergone a giant eruption, with a substantial increase in stellar luminosity, radius, and almost certainly mass-loss rate.


title: Relativistic Viscous Fluid Description of Microscopic Black Hole Wind

abstract: Microscopic black holes explode with their temperature varying inversely as their mass. Such explosions would lead to the highest temperatures in the present universe, all the way to the Planck energy. Whether or not a quasi-stationary shell of matter undergoing radial hydrodynamic expansion surrounds such black holes is been controversial. In this paper relativistic viscous fluid equations are applied to the problem. It is shown that a self-consistent picture emerges of a fluid just marginally kept in local thermal equilibrium; viscosity is a crucial element of the dynamics.


title: Mining the Metal-Rich Stars for Planets

abstract: We examine the correlation between stellar metallicity and short period planets. It appears that approximately 1% of dwarf stars in the solar neighborhood harbor short-period planets characterized by near-circular orbits and orbital periods P<20 days. However, among the most metal-rich stars (defined as having [Fe/H]>0.2 dex), it appears that the fraction increases to 10%. Using the Hipparcos database and the Hauck & Mermilliod (1998) compilation of Stromgren uvby photometry, we identify a sample of 206 metal-rich stars of spectral type K, G, and F which have an enhanced probability of harboring short-period planets. Many of these stars would be excellent candidates for addition to radial velocity surveys. We have searched the Hipparcos epoch photometry for transiting planets within our 206 star catalog. We find that the quality of the Hipparcos data is not high enough to permit unambiguous transit detections. It is, however, possible to identify candidate transit periods. We then discuss various ramifications of the stellar metallicity - planet connection. First, we show that there is preliminary evidence for increasing metallicity with increasing stellar mass among known planet-bearing stars. This trend can be explained by a scenario in which planet-bearing stars accrete an average of 30 Earth Masses of rocky material after the gaseous protoplanetary disk phase has ended. We present dynamical calculation which suggest that a survey of metallicities of spectroscopic binary stars can be used to understand the root cause of the stellar metallicity - planet connection.


title: Cosmic Statistics of Statistics: N-point Correlations

abstract: The fully general calculation of the cosmic error on N-point correlation functions and related quantities is presented. More precisely, the variance caused by the finite volume, discreteness, and edge effects is determined for {\em any} estimator which is based on a general function of N-tuples, such as multi-point correlation functions and multi-spectra. The results are printed explicitly for the two-point correlation function (or power-spectrum), and for the three-point correlation (or bispectrum). These are the most popular statistics in the study of large scale structure, yet, the a general calculation of their variance has not been performed until now.


title: A New Underground Laboratory in the USA for a Neutrino Factory Detector and Other Scientific Projects

abstract: A neutrino factory storage ring can provide beams to various locations. We discuss the ICANOE detector (at LNGS) at one such site. We then describe the prospects for the use of the underground location at Carlsbad, NM for a neutrino factory detector. A brief discussion is given of a simple magnetized Fe detector of 10 50 kT for this site.


title: Stellar Collisions and Pulsar Planets

abstract: I describe models for the formation of planetary systems surrounding the remnants of stellar mergers and collisions. I focus primarily on models for the viscous evolution of disks suitable for the formation of the planets surrounding the pulsar B1257+12. I show that the adaptation of models for traditional protoplanetary disks which invoke quiescent or `dead' zones are quite successful in producing disks appropriate for the formation of the pulsar planets. I also briefly describe some even more exotic possibilities that may arise from compact object mergers.


title: Fast Convolution on the Sphere

abstract: We propose fast, exact and efficient algorithms for the convolution of two arbitrary functions on the sphere which speed up computations by a factor \order{\sqrt{N}} compared to present methods where $N$ is the number of pixels. No simplifying assumptions are made other than bandlimitation. This reduces typical computation times for convolving the full sky with the asymmetric beam pattern of a megapixel Cosmic Microwave Background (CMB) mission from months to minutes. Our methods enable realistic simulation and careful analysis of data from such missions, taking into account the effects of asymmetric ``point spread functions'' and far side lobes of the physical beam. While motivated by CMB studies, our methods are general and hence applicable to the convolution or filtering of any scalar field on the sphere with an arbitrary, asymmetric kernel. We show in an appendix that the same ideas can be applied to the inverse problems of map-making and beam reconstruction by similarly accelerating the {\em transpose convolution} which is needed for the iterative solution of the normal equations.


title: All-sky convolution for polarimetry experiments

abstract: We discuss all-sky convolution of the instrument beam with the sky signal in polarimetry experiments, such as the Planck mission which will map the temperature anisotropy and polarization of the cosmic microwave background (CMB). To account properly for stray light (from e.g. the galaxy, sun, and planets) in the far side-lobes of such an experiment, it is necessary to perform the beam convolution over the full sky. We discuss this process in multipole space for an arbitrary beam response, fully including the effects of beam asymmetry and cross-polarization. The form of the convolution in multipole space is such that the Wandelt-Gorski fast technique for all-sky convolution of scalar signals (e.g. temperature) can be applied with little modification. We further show that for the special case of a pure co-polarized, axisymmetric beam the effect of the convolution can be described by spin-weighted window functions. In the limits of a small angle beam and large Legendre multipoles, the spin-weight 2 window function for the linear polarization reduces to the usual scalar window function used in previous analyses of beam effects in CMB polarimetry experiments. While we focus on the example of polarimetry experiments in the context of CMB studies, we emphasise that the formalism we develop is applicable to anisotropic filtering of arbitrary tensor fields on the sphere.


title: [Ne II] 12.8 Micron Images of Four Galactic Ultracompact H II Regions: Ionized Neon Abundance as a Tracer of the Ionizing Stars

abstract: We present the results of ground-based imaging spectroscopy of the [Ne II] 12.8 micron line emitted from the ultracompact (UC) H II regions; W51d, G45.12+0.13, G35.20-1.74 and Monoceros R2, with 2arcsec spatial resolution. We found that the overall distribution of the [Ne II] emission is generally in good agreement with the radio (5 or 15 GHz) VLA distribution for each source. The Ne+ abundance ([Ne+/H+]) distributions are also derived from the [Ne II] and the radio maps. As for G45.12+0.13 and W51d, the Ne+ abundance decreases steeply from the outer part of the map toward the radio peak. On the other hand, the Ne+ abundance distributions of G35.20-1.74 and Mon R2 appear rather uniform. These results can be interpreted by the variation of ionizing structures of neon, which is primarily determined by the spectral type of the ionizing stars. We have evaluated the effective temperature of the ionizing star by comparing the Ne+ abundance averaged over the whole observed region with that calculated by H II region models based on recent non-LTE stellar atmosphere models: 39,100 (+1100, -500) K (O7.5V-O8V) for W51d, 37,200 (+1000, -700) K (O8V-O8.5V) for G45.12+0.13, 35,000-37,600 (+1500, -600) K (O8V-O9V) for G35.20-1.74, and < 34,000 K (< B0V) for Mon R2. These effective temperatures are consistent with those inferred from the observed Ne+ abundance distributions.


title: The Dynamics of the Formation of Degenerate Heavy Neutrino Stars

abstract: Through reformulating the cold, self-gravitating fermion gas as a Bose condensate by identifying their mutual Thomas-Fermi limits, the dissipationless formation of a heavy neutrino star in gravitational collapse is numerically demonstrated. Such stars offer an alternative to supermassive black holes for the compact dark objects at the centers of galaxies.


title: Is the deficit of z > 1 field ellipticals real ?

abstract: The results of a 700 arcmin^2 survey for EROs are presented. The sky distribution of EROs is very inhomogeneous and their two point correlation function shows a very high amplitude, a factor of 10 larger than that of the field galaxy population. Such a clustering can explain the strong variance found in previous works, and provides evidence that ERO are mostly made of high-z ellipticals. The surface density of EROs found in our large survey is in good agreement with that expected for passively evolving ellipticals formed at high redshift (z\simgt2.5).