We develop a relativistic multifluid dynamics appropriate for describing neutron star cores at finite temperatures based on Carter's convective variational procedure. The model includes seven fluids, accounting for both normal and superfluid/superconducting neutrons and protons, leptons (electrons and muons) and entropy. The formulation is compared to the non-variational relativistic multifluid hydrodynamics of Gusakov and collaborators and shown to be equivalent. Vortex lines and flux tubes, mutual friction, vortex pinning, heat conduction and viscosity are incorporated into the model in steps after the basic hydrodynamics is described. The multifluid system is then considered at the mesoscopic scale where the currents around individual vortex lines and flux tubes are important, and this mesoscopic theory is averaged to determine the detailed vortex line/flux tube contributions to the macroscopic ``effective'' theory. This matching procedure is partially successful, though obtaining full agreement between the averaged mesoscopic and macroscopic effective theory requires discarding subdominant terms. The matching procedure allow us to interpret the magnetic $H$-field inside a neutron star in a way that is consistent with condensed matter physics literature, and to clarify the difference between this interpretation and that in previous astrophysical works.

We investigate the compressional modes of cold neutron stars with cores consisting of superfluid neutrons, superconducting protons and normal fluid electrons and muons, and crusts that contain superfluid neutrons plus a normal fluid of (spherical) nuclei and electrons. We develop a two-fluid formalism for the core that accounts for leptonic buoyancy, and an analogous treatment for the crust. We adopt the Cowling approximation, neglecting gravitational perturbations, but include all effects of the background space-time. We introduce a phenomenological, easily-modified nuclear equation of state which contains all of the thermodynamic information required to compute the coupled fluid oscillations, with parameters that are constrained by nuclear physics and the requirement that the maximum mass of a neutron star is $\geq 2M_{\odot}$. Using four parametrizations of this equation of state with nuclear compressibilities $K=230$-$280$ MeV, we calculate the Brunt-V\"{a}is\"{a}l\"{a} frequency due to leptonic buoyancy, and find the corresponding $g$-mode frequencies and eigenfunctions. We find that the WKB approximation reproduces $g$-mode frequencies closely. We examine the dependence of $g$-mode frequencies on stellar mass, nuclear compressibility and strength of neutron-proton entrainment, and compare to previous calculations of $g$-mode frequencies due to leptonic buoyancy. We also compute the $p$-mode spectra, confirming previous findings that the two fluids behave as if uncoupled except the case of large entrainment, and show the existence of nearly resonant mode pairs which could lead to nonlinear $p$-$g$ instabilities even at zero temperature.

We numerically investigate the generation of a magnetic field in a protostellar disc via an $\alpha \Omega$-dynamo and the resulting magnetohydrodynamic (MHD) driven outflows. We find that for small values of the dimensionless dynamo parameter $\alpha_d$ the poloidal field grows exponentially at a rate $\sigma \propto \Omega_K \sqrt{\alpha_d}$, before saturating to a value $\propto \sqrt{\alpha_d}$. The dynamo excites dipole and octupole modes, but quadrupole modes are suppressed, because of the symmetries of the seed field. Initial seed fields too weak to launch MHD outflows are found to grow sufficiently to launch winds with observationally relevant mass fluxes of order $10^{-9} M_{\odot}/\rm{yr}$ for T Tauri stars. This suggests $\alpha \Omega$-dynamos may be responsible for generating magnetic fields strong enough to launch observed outflows.

Plasma lenses in the host galaxies of fast radio bursts (FRBs) can strongly modulate FRB amplitudes for a wide range of distances, including the $\sim $ Gpc distance of the repeater FRB121102. To produce caustics, the lens' dispersion-measure depth (${\rm DM}_{\ell}$), scale size ($a$), and distance from the source ($d_{\rm sl}$) must satisfy ${\rm DM}_{\ell} d_{\rm sl} / a^2 \gtrsim 0.65~ {\rm pc^2 \ AU^{-2} \ cm^{-3}}$. Caustics produce strong magnifications ($\lesssim 10^2$) on short time scales ($\sim$ hours to days and perhaps shorter) along with narrow, epoch dependent spectral peaks (0.1 to 1~GHz). However, strong suppression also occurs in long-duration ($\sim$ months) troughs. For geometries that produce multiple images, the resulting burst components will arrive differentially by $< 1~\mu$s to tens of ms and they will show different apparent dispersion measures, $\delta{\rm DM}_{\rm apparent} \sim 1$ pc cm$^{-3}$. Arrival time perturbations may mask any underlying periodicity with period $\lesssim 1$ s. When arrival times differ by less than the burst width, interference effects in dynamic spectra are expected. Strong lensing requires source sizes smaller than $({\rm Fresnel~scale)^2} / a$, which can be satisfied by compact objects such as neutron star magnetospheres but not by AGNs. Much of the phenomenology of the repeating fast radio burst source FRB121102 is similar to lensing effects. The overall picture can be tested by obtaining wideband spectra of bursts (from $<1$ to 10 GHz and possibly higher), which can also be used to characterize the plasma environment near FRB sources. A rich variety of phenomena is expected from an ensemble of lenses near the FRB source. We discuss constraints on densities, magnetic fields, and locations of plasma lenses related to requirements for lensing to occur.

We analyze plasma dispersion and scattering of fast radio bursts (FRBs) to identify the dominant locations of free electrons along their lines of sight and thus constrain the distances of the burst sources themselves. We establish the average $\tau$-DM relation for Galactic pulsars and use it as a benchmark for discussing FRB scattering. Though scattering times $\tau$ for FRBs are large in the majority of the 17 events we analyze, they are systematically smaller than those of Galactic pulsars that have similar dispersion measures (DMs). The lack of any correlation between $\tau$ and DM for FRBs suggests that the intergalactic medium (IGM) cannot account for both $\tau$ and DM. We therefore consider mixed models involving the IGM and host galaxies. If the IGM contributes significantly to DM while host galaxies dominate $\tau$, the scattering deficit with respect to the mean Galactic trend can be explained with a $\tau$-DM relation in the host that matches that for the Milky Way. However, it is possible that hosts dominate both $\tau$ and DM, in which case the observed scattering deficits require free electrons in the host to be less turbulent than in the Galaxy, such as if they are in hot rather than warm ionized regions. Our results imply that distances or redshifts of FRB sources can be significantly overestimated if they are based on the assumption that the extragalactic portion of DM is dominated by the IGM.

We consider radio bursts that originate from extragalactic neutron stars (NSs) by addressing three questions about source distances. What are the physical limitations on coherent radiation at GHz frequencies? Do they permit detection at cosmological distances? How many bursts per NS are needed to produce the inferred burst rate $\sim 10^3$-$10^4 $sky$^{-1}$ day$^{-1}$? The burst rate is comparable to the NS formation rate in a Hubble volume, requiring only one per NS if they are bright enough. However, radiation physics causes us to favor a closer population. More bursts per NS are then required but repeats in 10 to 100 yr could still be negligible. Bursts are modeled as sub-ns, coherent shot pulses superposed incoherently to produce ms-duration $\sim 1$ Jy amplitudes; each shot-pulse can be much weaker than the burst amplitude, placing less restrictive requirements on the emission process. Nonetheless, single shot pulses are similar to the extreme, unresolved ($< 0.4$ ns) MJy shot pulse seen from the Crab pulsar, which is consistent with coherent curvature radiation emitted near the light cylinder by an almost neutral clump with net charge $\sim \pm 10^{21}e$ and total energy $\gtrsim 10^{23}$ ergs. Bursts from Gpc distances require incoherent superposition of $\sim 10^{12}d_{\rm Gpc}^2$ shot pulses or a total energy $\gtrsim 10^{35} d_{\rm Gpc}^2$ erg. The energy reservoir near the light cylinder limits the detection distance to $\lesssim {\rm few} \times 100$ Mpc for a fluence $\sim 1$ Jy ms unless conditions are more extreme than for the Crab pulsar. Similarly, extreme single pulses from ordinary pulsars and magnetars could be detectable from throughout the Local Group and perhaps farther. Contributions to dispersion measures from galaxy clusters will be significant for some of the bursts. We discuss tests for the signatures of bursts associated with extragalactic NSs.

We develop the formalism for computing the magnetic field within an axisymmetric neutron star with a strong Type II superconductor core surrounded by a normal conductor. The formalism takes full account of the constraints imposed by hydrostatic equilibrium with a barotropic equation of state. We specialize to purely poloidal magnetic fields and develop the "most dipolar case" for which we find that the surface field strength is $\simeq H_b\epsilon_b/3\simeq 3\times 10^{12}$ G, where $H_b$ is the magnetic field at the outer edge of the core and $\epsilon_b R$ is the thickness of the normal shell. This value only depends on the equation of state of nuclear matter. We also find that the quadrupolar distortion of the star is $\sim 10^{-9}$.

Ultra-high energy cosmic rays (UHECRs) are atomic nuclei with energies over ten million times energies accessible to human-made particle accelerators. Evidence suggests that they originate from relatively nearby extragalactic sources, but the nature of the sources is unknown. We develop a multilevel Bayesian framework for assessing association of UHECRs and candidate source populations, and Markov chain Monte Carlo algorithms for estimating model parameters and comparing models by computing, via Chib's method, marginal likelihoods and Bayes factors. We demonstrate the framework by analyzing measurements of 69 UHECRs observed by the Pierre Auger Observatory (PAO) from 2004-2009, using a volume-complete catalog of 17 local active galactic nuclei (AGN) out to 15 megaparsecs as candidate sources. An early portion of the data ("period 1," with 14 events) was used by PAO to set an energy cut maximizing the anisotropy in period 1; the 69 measurements include this "tuned" subset, and subsequent "untuned" events with energies above the same cutoff. Also, measurement errors are approximately summarized. These factors are problematic for independent analyses of PAO data. Within the context of "standard candle" source models (i.e., with a common isotropic emission rate), and considering only the 55 untuned events, there is no significant evidence favoring association of UHECRs with local AGN vs. an isotropic background. The highest-probability associations are with the two nearest, adjacent AGN, Centaurus A and NGC 4945. If the association model is adopted, the fraction of UHECRs that may be associated is likely nonzero but is well below 50%. Our framework enables estimation of the angular scale for deflection of cosmic rays by cosmic magnetic fields; relatively modest scales of $\approx\!3^{\circ}$ to $30^{\circ}$ are favored. Models that assign a large fraction of UHECRs to a single nearby source (e.g., Centaurus A) are ruled out unless very large deflection scales are specified a priori, and even then they are disfavored. However, including the period 1 data alters the conclusions significantly, and a simulation study supports the idea that the period 1 data are anomalous, presumably due to the tuning. Accurate and optimal analysis of future data will likely require more complete disclosure of the data.

We describe how the nonlinear development of the R mode instability of neutron stars influences spin up to millisecond periods via accretion. Our arguments are based on nearly-resonant interactions of the R mode with pairs of "daughter modes". The amplitude of the R mode saturates at the lowest value for which parametric instability leads to significant excitation of a particular pair of daughters. The lower bound on this limiting amplitude is proportional to the damping rate of the daughter modes that are excited parametrically. Based on this picture, we show that if modes damp because of dissipation in a very thin boundary layer at the crust-core boundary then spin up to frequencies larger than about 300 Hz does not occur. Within this conventional scenario the R mode saturates at an amplitude that is too large for angular momentum gain from accretion to overcome gravitational loss to gravitational radiation. We conclude that lower dissipation is required for spin up to frequencies much higher than 300 Hz. We conjecture that if the transition from the fluid core to the crystalline crust occurs over a distance much longer than 1 cm then a sharp viscous boundary layer fails to form. In this case, damping is due to shear viscosity dissipation integrated over the entire star; the rate is slower than if a viscous boundary layer forms. We use statistical arguments and scaling relations to estimate the lowest parametric instability threshold from first principles. The resulting saturation amplitudes are low enough to permit spin up to higher frequencies. Further, we show that the requirement that the lowest parametric instability amplitude be small enough to allow continued spin up imposes an upper bound to the frequencies that may be attained via accretion that may plausibly be about 750 Hz. Within this framework, the R mode is unstable for all millisecond pulsars, whether accreting or not.

We study the fluctuations in luminosity distance due to gravitational lensing produced both by galaxy halos and large scale voids. Voids are represented via a "Swiss cheese" model consisting of a \LambdaCDM Friedman-Robertson-Walker background in which a number of randomly distributed, spherical regions of comoving radius 35 Mpc are removed. A fraction of the removed mass is then placed on the shells of the spheres, in the form of randomly located halos, modeled with Navarro-Frenk-White profiles. The remaining mass is placed in the interior of the spheres, either smoothly distributed, or as randomly located halos. We compute the distribution of magnitude shifts using a variant of the method of Holz & Wald (1998), which includes the effect of lensing shear. In the two models we consider, the standard deviation of this distribution is 0.065 and 0.072 magnitudes and the mean is -0.0010 and -0.0013 magnitudes, for voids of radius 35 Mpc, sources at redshift 1.5, with the voids chosen so that 90% of the mass is on the shell today. The standard deviation due to voids and halos is a factor ~ 3 larger than that due to 35 Mpc voids alone with a 1 Mpc shell thickness which we studied in our previous work. To a good approximation, the variance of the distribution depends only on the mean column depth and concentration of halos and on the fraction of the mass density that is in the form of halos (as opposed to smoothly distributed): it is independent of how the halos are distributed in space. We derive an approximate analytic formula for the variance that agrees with our numerical results to \lesssim 20% out to z\simeq 1.5.

The Earth is continuously showered by charged cosmic ray particles, naturally produced atomic nuclei moving with velocity close to the speed of light. Among these are ultra high energy cosmic ray particles with energy exceeding 5x10^19 eV, which is ten million times more energetic than the most energetic particles produced at the Large Hadron Collider. Astrophysical questions include: what phenomenon accelerates particles to such high energies, and what sort of nuclei are energized? Also, the magnetic deflection of the trajectories of the cosmic rays makes them potential probes of galactic and intergalactic magnetic fields. We develop a Bayesian hierarchical model that can be used to compare different association models between the cosmic rays and source population, using Bayes factors. A measurement model with directional uncertainties and accounting for non-uniform sky exposure is incoporated into the model. The methodology allows us to learn about astrophysical parameters, such as those governing the source luminosity function and the cosmic magnetic field.

We study the fluctuations in luminosity distances due to gravitational lensing by large scale (> 35 Mpc) structures, specifically voids and sheets. We use a simplified "Swiss cheese" model consisting of a \Lambda -CDM Friedman-Robertson-Walker background in which a number of randomly distributed non-overlapping spherical regions are replaced by mass compensating comoving voids, each with a uniform density interior and a thin shell of matter on the surface. We compute the distribution of magnitude shifts using a variant of the method of Holz & Wald (1998), which includes the effect of lensing shear. The standard deviation of this distribution is ~ 0.027 magnitudes and the mean is ~ 0.003 magnitudes for voids of radius 35 Mpc, sources at redshift z_s=1.0, with the voids chosen so that 90% of the mass is on the shell today. The standard deviation varies from 0.005 to 0.06 magnitudes as we vary the void size, source redshift, and fraction of mass on the shells today. If the shell walls are given a finite thickness of ~ 1 Mpc, the standard deviation is reduced to ~ 0.013 magnitudes. This standard deviation due to voids is a factor ~ 3 smaller than that due to galaxy scale structures. We summarize our results in terms of a fitting formula that is accurate to ~ 20%, and also build a simplified analytic model that reproduces our results to within ~ 30%. Our model also allows us to explore the domain of validity of weak lensing theory for voids. We find that for 35 Mpc voids, corrections to the dispersion due to lens-lens coupling are of order ~ 4%, and corrections to due shear are ~ 3%. Finally, we estimate the bias due to source-lens clustering in our model to be negligible.

Cosmic string networks generate cosmological perturbations actively throughout the history of the universe. Thus, the string sourced anisotropy of the cosmic microwave background is not affected by Silk damping as much as the anisotropy seeded by inflation. The spectrum of perturbations generated by strings does not match the observed CMB spectrum on large angular scales (l<1000) and is bounded to contribute no more than 10% of the total power on those scales. However, when this bound is marginally saturated, the anisotropy created by cosmic strings on small angular scales l>2000 will dominate over that created by the primary inflationary perturbations. This range of angular scales in the CMB is presently being measured by a number of experiments; their results will test this prediction of cosmic string networks soon.

We argued in a previous paper [R. A. Vanderveld et al. 2006, arXiv:astro-ph/0602476] that negative deceleration parameters at the center of symmetry in Lemaitre-Tolman-Bondi cosmological models can only occur if the model is not smooth at the origin. Here we demonstrate explicitly the connection between non-smoothness and the failure of positivity theorems for deceleration. We also address some confusion that has arisen in the literature and respond to some recent criticisms of our arguments.

We model the nonlinear saturation of the r-mode instability via three-mode couplings and the effects of the instability on the spin evolution of young neutron stars. We include one mode triplet consisting of the r-mode and two near resonant inertial modes that couple to it. We find that the spectrum of evolutions is more diverse than previously thought. The evolution of the star is dynamic and initially dominated by fast neutrino cooling. Nonlinear effects become important when the r-mode amplitude grows above its first parametric instability threshold. The balance between neutrino cooling and viscous heating plays an important role in the evolution. Depending on the initial r-mode amplitude, and on the strength of the viscosity and of the cooling this balance can occur at different temperatures. If thermal equilibrium occurs on the r-mode stability curve, where gravitational driving equals viscous damping, the evolution may be adequately described by a one-mode model. Otherwise, nonlinear effects are important and lead to various more complicated scenarios. Once thermal balance occurs, the star spins-down oscillating between thermal equilibrium states until the instability is no longer active. For lower viscosity we observe runaway behavior in which the r-mode amplitude passes several parametric instability thresholds. In this case more modes need to be included to model the evolution accurately. In the most optimistic case, we find that gravitational radiation from the r-mode instability in a very young, fast spinning neutron star within about 1 Mpc of Earth may be detectable by advanced LIGO for years, and perhaps decades, after formation. Details regarding the amplitude and duration of the emission depend on the internal dissipation of the modes of the star, which would be probed by such detections.

We derive an expression for the luminosity distance as a function of redshift for a flat Robertson-Walker spacetime perturbed by arbitrary scalar perturbations possibly produced by a modified gravity theory with two different scalar perturbation potentials. Measurements of the luminosity distance as function of redshift provide a constraint on a combination of the scalar potentials and so they can complement weak lensing and other measurements in trying to distinguish among the various alternative theories of gravity.

Recently there have been suggestions that the Type Ia supernova data can be explained using only general relativity and cold dark matter with no dark energy. In "Swiss cheese" models of the Universe, the standard Friedmann-Robertson-Walker picture is modified by the introduction of mass compensating spherical inhomogeneities, typically described by the Lemaitre-Tolman-Bondi metric. If these inhomogeneities correspond to underdense cores surrounded by mass-compensating overdense shells, then they can modify the luminosity distance-redshift relation in a way that can mimic accelerated expansion. It has been argued that this effect could be large enough to explain the supernova data without introducing dark energy or modified gravity. We show that the large apparent acceleration seen in some models can be explained in terms of standard weak field gravitational lensing together with insufficient randomization of void locations. The underdense regions focus the light less than the homogeneous background, thus dimming supernovae in a way that can mimic the effects of acceleration. With insufficient randomization of the spatial location of the voids and of the lines of sight, coherent defocusing can lead to anomalously large demagnification effects. We show that a proper randomization of the voids and lines of sight reduces the effect to the point that it can no longer explain the supernova data.

We determine constraints on the form of axisymmetric toroidal magnetic fields dictated by hydrostatic balance in a type II superconducting neutron star with a barotropic equation of state. Using Lagrangian perturbation theory, we find the quadrupolar distortions due to such fields for various models of neutron stars with type II superconducting and normal regions. We find that the star becomes prolate and can be sufficiently distorted to display precession with a period of the order of years. We also study the stability of such fields using an energy principle, which allows us to extend the stability criteria established by R. J. Tayler for normal conductors to more general media with magnetic free energy that depends on density and magnetic induction, such as type II superconductors. We also derive the growth rate and instability conditions for a specific instability of type II superconductors, first discussed by P. Muzikar, C. J. Pethick and P. H. Roberts, using a local analysis based on perturbations around a uniform background.

The nonlinear saturation of the r-mode instability and its effects on the spin evolution of Low Mass X-ray Binaries (LMXBs) are modeled using the triplet of modes at the lowest parametric instability threshold. We solve numerically the coupled equations for the three mode amplitudes in conjunction with the spin and temperature evolution equations. We observe that very quickly the mode amplitudes settle into quasi-stationary states. Once these states are reached, the mode amplitudes can be found algebraically and the system of equations is reduced from eight to two equations: spin and temperature evolution. Eventually, the system may reach thermal equilibrium and either (1) undergo a cyclic evolution with a frequency change of at most 10%, (2) evolve toward a full equilibrium state in which the accretion torque balances the gravitational radiation emission, or (3) enter a thermogravitational runaway on a very long timescale of about $10^6$ years. Alternatively, a faster thermal runaway (timescale of about 100 years) may occur. The sources of damping considered are shear viscosity, hyperon bulk viscosity and boundary layer viscosity. We vary proprieties of the star such as the hyperon superfluid transition temperature T_c, the fraction of the star that is above the threshold for direct URCA reactions, and slippage factor, and map the different scenarios we obtain to ranges of these parameters. For all our bound evolutions the r-mode amplitude remains small $\sim 10^{-5}$. The spin frequency is limited by boundary layer viscosity to $\nu_{max} \sim 800 Hz [S_{ns}/(M_{1.4} R_6)]^{4/11} T_8^{-2/11}$. We find that for $\nu > 700$ Hz the r-mode instability would be active for about 1 in 1000 LMXBs and that only the gravitational waves from LMXBs in the local group of galaxies could be detected by advanced LIGO interferometers.

We calculate the systematic inhomogeneity-induced correction to the cosmological constant that one would infer from an analysis of the luminosities and redshifts of Type Ia supernovae, assuming a homogeneous universe. The calculation entails a post-Newtonian expansion within the framework of second order perturbation theory, wherein we consider the effects of subhorizon density perturbations in a flat, dust dominated universe. Within this formalism, we calculate luminosity distances and redshifts along the past light cone of an observer. The resulting luminosity distance-redshift relation is fit to that of a homogeneous model in order to deduce the best-fit cosmological constant density Omega_Lambda. We find that the luminosity distance-redshift relation is indeed modified, by a small fraction of order 10^{-5}. When fitting this perturbed relation to that of a homogeneous universe, we find that the inferred cosmological constant can be surprisingly large, depending on the range of redshifts sampled. For a sample of supernovae extending from z=0.02 out to z=0.15, we find that Omega_Lambda=0.004. The value of Omega_Lambda has a large variance, and its magnitude tends to get larger for smaller redshifts, implying that precision measurements from nearby supernova data will require taking this effect into account. However, we find that this effect is likely too small to explain the observed value of Omega_Lambda=0.7. There have been previous claims of much larger backreaction effects. By contrast to those calculations, our work is directly related to how observers deduce cosmological parameters from astronomical data.

There has been much debate over whether or not one could explain the observed acceleration of the Universe with inhomogeneous cosmological models, such as the spherically-symmetric Lemaitre-Tolman-Bondi (LTB) models. It has been claimed that the central observer in these models can observe a local acceleration, which would contradict general theorems. We resolve the contradiction by noting that many of the models that have been explored contain a weak singularity at the location of the observer which makes them unphysical. In the absence of this singularity, we show that LTB models must have a positive central deceleration parameter $q_{0}$, in agreement with the general theorems. We also show that it is possible to achieve a negative apparent deceleration parameter at nonzero redshifts in LTB models that do not contain this singularity. However, we find other singularities that tend to arise in LTB models when attempting to match luminosity distance data, and these generally limit the range of redshifts for which these models can mimic observations of an accelerating Universe. Exceptional models do exist that can extend to arbitrarily large redshift without encountering these pathologies, and we show how these may be constructed. These special models exhibit regions with negative effective equation of state parameter, which may fall below negative one, but we have failed to find any singularity-free models that agree with observations. Moreover, models based on dust-filled LTB metrics probably fail to reproduce observed properties of large scale structure.

We find the constraints from WMAP and SDSS data on the fraction of cosmological fluctuations sourced by local cosmic strings using a Markov Chain Monte Carlo (MCMC) analysis. In addition to varying the usual 6 cosmological parameters and the string tension ($\mu$), we also varied the amount of small-scale structure on the strings. Our results indicate that cosmic strings can account for up to 7 (14)% of the total power of the microwave anisotropy at 68 (95)% confidence level. The corresponding bound on the string mass per unit length, within our string model, is $G\mu < 1.8 (2.7) \times 10^{-7}$ at 68 (95)% c.l., where this constraint has been altered from what appears below following the correction of errors in our cosmic string code outlined in a recent erratum, astro-ph/0604141. We also calculate the B-type polarization spectra sourced by cosmic strings and discuss the prospects of their detection.

Overall, brane inflation is compatible with the recent analysis of the WMAP data. Here we explore the constraints of WMAP and 2dFGRS data on the various brane inflationary scenarios. Brane inflation naturally ends with the production of cosmic strings, which may provide a way to distinguish these models observationally. We argue that currently available data cannot exclude a non-negligible contribution from cosmic strings definitively. We perform a partial statistical analysis of mixed models that include a sub-dominant contribution from cosmic strings. Although the data favor models without cosmic strings, we conclude that they cannot definitively rule out a cosmic-string-induced contribution of $\sim 10 %$ to the observed temperature, polarization and galaxy density fluctuations. These results imply that $G\mu \lesssim 3.5\times 10^{-7}(\lambda/0.25)\sqrt{B/0.1}$, where $\lambda$ is a dimensionless parameter related to the interstring distance, and $B$ measures the importance of perturbations induced by cosmic strings. We argue that, conservatively, the data available currently still permit $B\lesssim 0.1$. Precision measurements sensitive to the B-mode polarization produced by vector density perturbation modes driven by the string network could provide evidence for these models. Accurate determinations of $n_s(k)$, the scalar fluctuation index, could also distinguish among various brane inflation models.

In a recent publication, we used the data from WMAP and SDSS to constrain the primordial perturbations and to predict the B-mode polarization sourced by cosmic string networks. We have been alerted by A. Slosar to the existence of errors in the code we used to calculate the Cosmic Microwave Background anisotropies from cosmic strings. Correcting the errors leads to a significant increase in the vector mode contribution to the CMB temperature and polarization anisotropies as well as an overall renormalization of the various string spectra. In these notes we explain the nature of the errors and discuss their implications for previously published constraints on cosmic strings based on this code. The chief change in our results is that our derived limit for the cosmic string tension is strengthened: G\mu < 1.8 (2.7) * 10^{-7} at 68 (95)% confidence. We also note that the newly-enhanced vector mode contribution produces a greatly-increased amplitude for B-mode polarization in the CMB which could exceed the B-mode power produced by the lensing of primordial E-mode polarization into B-mode polarization.

The discovery of the first gravitationally redshifted spectral line from a neutron star (NS) by Cottam, Paerels and Mendez has triggered theoretical studies of the physics of atomic line formation in NS atmospheres. Chang, Bildsten and Wasserman showed that the hydrogenic Fe H$\alpha$ line formed above the photosphere of a bursting NS is intrinsically broad. We now include rotational broadening within general relativity and compare the resulting profile to that observed during Type I bursts from EXO 0748-676. We show that the fine structure splitting of the line precludes a meaningful constraint on the radius. Our fitting of the data show that the line forming Fe column is ${\rm log}_{10} (N_{\rm Fe, n=2}/{\rm cm^{-2}})=17.9_{-0.42}^{+0.27}$ and gravitational redshift $1+z =1.345_{-0.008}^{+0.005}$ with 95% confidence. We calculate the detectability of this spectral feature for a large range of spins and inclinations assuming that the emission comes from the entire surface. We find that at 300 (600) Hz only 10-20% (5-10%) of NSs would have spectral features as deep as that seen in EXO 0748-676.

Stairs, Lyne & Shemar have found that arrival time residuals from PSR B1828-11 vary periodically with a period of 500 days. This behavior can be accounted for by precession of the radiopulsar, an interpretation that is reinforced by the detection of variations in its pulse profile on the same timescale. Here, we model the period residuals from PSR B1828-11 in terms of precession of a triaxial rigid body. We include two contributions to the residuals: (i) the geometric effect, which arises because the times at which the pulsar emission beam points toward the observer varies with precession phase; (ii) the spindown contribution, which arises from any dependence of the spindown torque acting on the pulsar on the angle between its spin and magnetic axes. We use the data to probe numerous properties of the pulsar, most notably its shape, and the dependence of its spindown torque on the angle between its spin and magnetic axes, for which we assume a sum of a spin-aligned component (with a weight 1-a) and a dipolar component perpendicular to the magnetic beam axis (weight a), rather than the vacuum dipole torque (a=1). We find that a variety of shapes are consistent with the residuals, with a slight statistical preference for a prolate star. Moreover, a range of torque possibilities fit the data equally well, with no strong preference for the vacuum model. In the case of a prolate star we find evidence for an angle-dependent spindown torque. Our results show that the combination of geometrical and spin-down effects associated with precession can account for the principal features of PSR B1828-11's timing behavior, without fine tuning of the parameters.

Brane inflation in superstring theory ends when branes collide, initiating the hot big bang. Cosmic superstrings are produced during the brane collision. The cosmic superstrings produced in a D3-brane-antibrane inflationary scenario have a spectrum: $(p,q)$ bound states of $p$ fundamental (F) strings and $q$ D-strings, where $p$ and $q$ are coprime. By extending the velocity-dependent one-scale network evolution equations for abelian Higgs cosmic strings to allow a spectrum of string tensions, we construct a coupled (infinite) set of equations for strings that interact through binding and self-interactions. We apply this model to a network of $(p,q)$ superstrings. Our numerical solutions show that $(p,q)$ networks rapidly approach a stable scaling solution. We also extract the relative densities of each string type from our solutions. Typically, only a small number of the lowest tension states are populated substantially once scaling is reached. The model we study also has an interesting new feature: the energy released in $(p,q)$ string binding is by itself adequate to allow the network to reach scaling. This result suggests that the scaling solution is robust. To demonstrate that this result is not trivial, we show that choosing a different form for string interactions can lead to network frustration.

Spherical gravitational collapse of a cold gas of annihilating particles involves a competition between the free-fall rate $\propto\sqrt{\rho}$ and the (s-wave) annihilation rate $\propto\rho$. Thus, there is a critical density $\rhoann$ above which annihilation proceeds faster than free fall. Gravitational collapse of a cloud of (initial) mass $M$ to a black hole is only possible if $3/32\pi G^3M^2\lesssim\rhoann$, or $M\gtrsim\Mann\equiv (3/32\pi G^3\rhoann)^{1/2}$. For a particle mass $m$ and freeze-out temperature $T_f=m/x_f$, the minimum black hole mass is $\Mann\approx 10^{10}\msun \times(x_f\sqrt{g_\star}/100\omcdm g_{\star S}m({\rm Gev}))$, where $g_{\star S}$ and $g_\star$ are degeneracy factors. The formation of a black hole of initial mass $M_{BH}$ is accompanied by the annihilation of about $M_{ann}$ released in a burst lasting a time $\sim GM_{BH}$ that could reach a total annihilation luminosity $\sim 10^{59} {\rm erg s^{-1}}$. The absence of astronomical observations of such spectacular events suggests either: (i) the branching ratio for CDM annihilation to $e^{+}e^{-}$ pairs or quarks $\lesssim 10^{-10}$, while the branching ratio to $\nu{\bar{\nu}}$ is $\lesssim 10^{-5}$; or (ii) CDM is not made of annihilating particles, but may be in some non-annihilating form, such as axions; or (iii) CDM black holes never form.

Motivated by the measurement of redshifted Fe H$\alpha$ lines during type I X-ray bursts on EXO 0748-676 (Cottam, Paerels & Mendez), we study the formation of atomic Fe lines above the photosphere of a bursting neutron star ($k_BT_{\rm eff} \approx 1-2 {\rm keV}$). We discuss the effects of Stark broadening, resonant scattering and NLTE (level population) on the formation of hydrogenic Fe H$\alpha$, Ly$\alpha$ and P$\alpha$ lines. From the observed equivalent width of the Fe H$\alpha$ line, we find an implied Fe column of $1-3 \times 10^{20} {\rm cm}^{-2}$, which is 3-10 times larger than the Fe column calculated from the accretion/spallation model of Bildsten, Chang & Paerels. We also estimate that the implied Fe column is about a factor of 2-3 larger than a uniform solar metallicity atmosphere. We discuss the effects of rotational broadening and find that the rotation rate of \EXO must be slow, as confirmed by the recent measurement of a 45 Hz burst oscillation by Villarreal & Strohmayer. We also show that the Fe Ly$\alpha$ EW $\approx$ 15-20 eV (redshifted 11-15 eV) and the P$\alpha$ EW $\approx$ 4-7 eV (redshifted 3-5 eV) when the H$\alpha$ EW is 10 eV (redshifted 8 eV). The Ly$\alpha$ line is rotationally broadened to a depth of $\approx 10%$, making it difficult to observe with {\it Chandra}. We also show that radiative levitation can likely support the Fe column needed to explain the line.

We study the stability of a collisionless, relativistic, finite-strength, cylindrical layer of charged particles in free space by solving the linearized Vlasov-Maxwell equations and compute the power of the emitted electromagnetic waves. The layer is rotating in an external magnetic field parallel to the layer. This system is of interest to understanding the high brightness temperature of pulsars which cannot be explained by an incoherent radiation mechanism. Coherent synchrotron radiation has also been observed recently in bunch compressors used in particle accelerators. We consider equilibrium layers with a `thermal' energy spread and therefore a non-zero radial thickness. The particles interact with their retarded electromagnetic self-fields. The effect of the betatron oscillations is retained. A short azimuthal wavelength instability is found which causes a modulation of the charge and current densities. The growth rate is found to be an increasing function of the azimuthal wavenumber, a decreasing function of the Lorentz factor, and proportional to the square root of the total number of electrons. We argue that the growth of the unstable perturbation saturates when the trapping frequency of electrons in the wave becomes comparable to the growth rate. Owing to this saturation we can predict the radiation spectrum for a given set of parameters. Our predicted brightness temperatures are proportional to the square of the number of particles and scale by the inverse five-third power of the azimuthal wavenumber which is in rough accord with the observed spectra of radio pulsars.

Current data exclude cosmic strings as the primary source of primordial density fluctuations. However, in a wide class of inflationary models, strings can form at later stages of inflation and have potentially detectable observational signatures. We study the constraints from WMAP and SDSS data on the fraction of primordial fluctuations sourced by local cosmic strings. The Bayesian analysis presented in this brief report is restricted to the minimal number of parameters. Yet it is useful for two reasons. It confirms the results of Pogosian et al (2003) using an alternative statistical method. Secondly, it justifies the more costly multi-parameter analysis. Already, varying only three parameters -- the spectral index and the amplitudes of the adiabatic and string contributions -- we find that the upper bound on the cosmic string contribution is of order 10%. We expect that the full multi-parameter study, currently underway, will likely loosen this bound.

Two mechanisms for nonlinear mode saturation of the r-mode in neutron stars have been suggested: the parametric instability mechanism involving a small number of modes and the formation of a nearly continuous Kolmogorov-type cascade. Using a network of oscillators constructed from the eigenmodes of a perfect fluid incompressible star, we investigate the transition between the two regimes numerically. Our network includes the 4995 inertial modes up to n<= 30 with 146,998 direct couplings to the r-mode and 1,306,999 couplings with detuning< 0.002 (out of a total of approximately 10^9 possible couplings). The lowest parametric instability thresholds for a range of temperatures are calculated and it is found that the r-mode becomes unstable to modes with 13<n<15. In the undriven, undamped, Hamiltonian version of the network the rate to achieve equipartition is found to be amplitude dependent, reminiscent of the Fermi-Pasta-Ulam problem. More realistic models driven unstable by gravitational radiation and damped by shear viscosity are explored next. A range of damping rates, corresponding to temperatures 10^6K to 10^9K, is considered. Exponential growth of the r-mode is found to cease at small amplitudes, approximately 10^-4. For strongly damped, low temperature models, a few modes dominate the dynamics. The behavior of the r-mode is complicated, but its amplitude is still no larger than about 10^-4 on average. For high temperature, weakly damped models the r-mode feeds energy into a sea of oscillators that achieve approximate equipartition. In this case the r-mode amplitude settles to a value for which the rate to achieve equipartition is approximately the linear instability growth rate.

R-modes of a rotating neutron star are unstable because of the emission of gravitational radiation. We explore the saturation amplitudes of these modes determined by nonlinear mode-mode coupling. Modelling the star as incompressible allows the analytic computation of the coupling coefficients. All couplings up to n=30 are obtained, and analytic values for the shear damping and mode normalization are presented. In a subsequent paper we perform numerical simulations of a large set of coupled modes.

We study explosions of stars using a one-dimensional Lagrangian hydrodynamics code. We calculate how much mass is liberated as a function of the energy of explosion for a variety of pre-explosion stellar structures and for equations of state with a range of radiation-to-gas pressure ratios. The results show that simple assumptions about the amount of mass lost in an explosion can be quite inaccurate, and that even one-dimensional stars exhibit a rich phenomenology. The mass loss fraction rises from about 50 to 100 per cent as a function of the explosion energy in an approximately discontinuous manner. Combining our results with those of other, more realistic models, we suggest that Nova Scorpii (J1655-40) may have experienced significant mass fallback because the explosion energy was less than the critical value. We infer that the original progenitor was less than twice the mass of today's remnant.

Non-linear interactions among the inertial modes of a rotating fluid can be described by a network of coupled oscillators. We use such a description for an incompressible fluid to study the development of the r-mode instability of rotating neutron stars. A previous hydrodynamical simulation of the r-mode reported the catastrophic decay of large amplitude r-modes. We explain the dynamics and timescale of this decay analytically by means of a single three mode coupling. We argue that at realistic driving and damping rates such large amplitudes will never actually be reached. By numerically integrating a network of nearly 5000 coupled modes, we find that the linear growth of the r-mode ceases before it reaches an amplitude of around 10^(-4). The lowest parametric instability thresholds for the r-mode are calculated and it is found that the r-mode becomes unstable to modes with 13<n<15 if modes up to n=30 are included. Using the network of coupled oscillators, integration times of 10^6 rotational periods are attainable for realistic values of driving and damping rates. Complicated dynamics of the modal amplitudes are observed. The initial development is governed by the three mode coupling with the lowest parametric instability. Subsequently a large number of modes are excited, which greatly decreases the linear growth rate of the r-mode.

If rotating core collapse leads to the formation of a proto-neutron star binary in super-close orbit, then the lighter star, propelled toward the minimum stable mass, explodes. The neutron star (or black hole) that remains acquires a spin-perpendicular kick of very large amplitude. Kicks of the type are required to explain the geodesic precession in double neutron star binaries such as B1913+16. On the contrary, spin-kick alignment has been claimed for the Vela and Crab pulsars whose kicks are relatively small. We propose that the larger kick component, when present in a pulsar, results from the disruption of an evanescent proto-neutron star binary and is spin perpendicular; the smaller kick component is associated to some other mechanism that leads to less vigorous kicks, predominantly parallel to the spin, because of phase averaging. This gives rise to a "bimodal distribution" in the peculiar velocities of neutron stars, as observed in the pulsar sample. The kick velocity is found to scale with the mass of the star that remains as M^(-2/3). This scenario can explain the run-away black hole GRO J1655-40, the first to show evidence for a natal kick.

Rossby waves (r-modes) in rapidly rotating neutron stars are unstable because of the emission of gravitational radiation. We study saturation of this instability by nonlinear transfer of energy to stellar "inertial" oscillation modes. We present detailed calculations of stellar inertial modes in the WKB limit, their linear damping by bulk and shear viscosity, and the nonlinear coupling forces among these modes. The saturation amplitude is derived in the extreme limits of strong or weak driving by radiation reaction, as compared to the damping rate of low order inertial modes. We find the saturation energy is {\it extremely small}, at least four orders of magnitude smaller than that found by previous investigators. We discuss the consequences of this result for spin evolution of young neutron stars, and neutron stars being spun up by accretion in Low Mass X-ray Binaries.We also discuss the detection of these gravitational waves by LIGO.

We present the results of two extensive Rossi X-ray Timing Explorer observations of large X-ray flaring episodes from the high-mass X-ray binary pulsar LMC X-4. Light curves during the flaring episodes comprise bright peaks embedded in relatively fainter regions, with complex patterns of recurrence and clustering of flares. We identify precursors preceding the flaring activity. Pulse profiles during the flares appear to be simple sinusoids, and pulsed fractions are proportional to the flare intensities. We fit Gaussian functions to flare peaks to estimate the mean full-width-half-maximum to be $\sim$68 s. Significant rapid aperiodic variability exists up to a few hertz during the flares, which is related to the appearance of narrow, spiky peaks in the light curves. While spectral fits and softness ratios show overall spectral softening as the flare intensity increases, the narrow, spiky peaks do not follow this trend. The mean fluence of the flare peaks is (3.1 $\pm$ 2.9) $\times$ 10$^{40}$ ergs in the 2.5--25 keV energy range, with its maximum at $\sim$1.9 $\times$ 10$^{41}$ ergs. The flare peak luminosity reaches up to (2.1 $\pm$ 0.2) $\times$ 10$^{39}$ ergs s$^{-1}$, far above the Eddington luminosity of a neutron star. We discuss possible origins of the flares, and we also propose that inhomogeneous accretion columns onto the neutron star polar caps are responsible for the observed properties.

We present Rossi X-ray Timing Explorer observations of the X-ray pulsar SMC X-1. The source is highly variable on short time scales (< 1 h), exhibiting apparent X-ray flares occupying a significant fraction (~3 %) of the total observing time, with a recurrence time of ~100 s. The flares seem to occur over all binary orbital phases, and correlate with the overall variability in the light curve. We find a total of 323 discrete flares which have a mean full width half maximum of ~18 s. The detailed properties of SMC X-1 do not vary significantly between the flares and the normal state, suggesting that the flare may be an extension of the normal state persistent emission with increased accretion rates. The flares resemble Type II X-ray bursts from GRO J1744--28. We discuss the origin of the SMC X-1 flares in terms of a viscous instability near the inner edge of the accretion disk around a weakly magnetized X-ray pulsar, and find this is consistent with the interpretation that SMC X-1 is in fact an intermediate-stage source like GRO J1744--28.

We consider the physics of free precession of a rotating neutron star with an oblique magnetic field. We show that if the magnetic stresses are large enough, then there is no possibility of steady rotation, and precession is inevitable. Even if the magnetic stresses are not strong enough to prevent steady rotation, we show that the minimum energy state is one in which the star precesses. Since the moment of inertia tensor is inherently triaxial in a magnetic star, the precession is periodic but not sinusoidal in time, in agreement with observations of PSR 1828-11. However, the problem we consider is {\it not} just precession of a triaxial body. If magnetic stresses dominate, the amplitude of the precession is not set just by the angle between the rotational angular velocity and any principal axis, which allows it to be small without suppressing oscillations of timing residuals at harmonics of the precession frequency. We argue that magnetic distortions can lead to oscillations of timing residuals of the amplitude, period, and relative strength of harmonics observed in PSR 1828-11 if magnetic stresses in its core are about 200 times larger than the classical Maxwell value for its dipole field, and the stellar distortion induced by these enhanced magnetic stresses is about 100-1000 times larger than the deformation of the neutron star's crust. Magnetic stresses this large can arise if the core is a Type II superconductor, or from toroidal fields $\sim 10^{14}$ G if the core is a normal conductor. The observations of PSR 1828-11 appear to require that the neutron star is slightly prolate.

If core collapse leads to the formation of a rapidly rotating bar-unstable proto-neutron star surrounded by fall-back material, then we might expect it to cool and fragment to form a double (proto)-neutron star binary into a super-close orbit. The lighter star should survive for awhile, until tidal mass loss propels it toward the minimum stable mass of a (proto)-neutron star, whereupon it explodes. Imshennik and Popov have shown that the explosion of the unstable, cold star can result in a large recoil velocity of the remaining neutron star. Here, we consider several factors that mitigate the effect and broaden the range of final recoil speeds, in particular the finite velocity and gravitational deflection of the ejecta, a range of original masses for the low mass companion and its cooling history, rotational phase averaging of the momentum impulse from non-instantaneous mass loss, and the possibility of a common envelope phase. In spite of these mitigating factors, we argue that this mechanism can still lead to substantial neutron star recoil speeds, close to, or even above, 1000 km/s.

Using a specific model for the expansion rate of the Universe as a function of scale factor, it is demonstrated that the equation of state of the dark energy cannot be determined uniquely from observations at redshifts $z\lesssim{\rm a few}$ unless the fraction of the mass density of the Universe in nonrelativistic particles, $\Omega_M$, somehow can be found independently. A phenomenological model is employed to discuss the utility of additional constraints from the formation of large scale structure and the positions of CMB peaks in breaking the degeneracy among models for the dark energy.

We discuss the signature of the scale of short distance physics in the Cosmic Microwave Background. In addition to effects which depend on the ratio of Hubble scale H during inflation to the energy scale M of the short distance physics, there can be effects which depend on $\dot{\phi}^2/M^4$ where $\phi$ is the {\it classical background} of the inflaton field. Therefore, the imprints of short distance physics on the spectrum of Cosmic Microwave Background anisotropies generically involve a {\it double expansion}. We present some examples of a single scalar field with higher order kinetic terms coupled to Einstein gravity, and illustrate that the effects of short distance physics on the Cosmic Microwave Background can be substantial even for H << M, and generically involve corrections that are not simply powers of H/M. The size of such effects can depend on the short distance scale non-analytically even though the action is local.

We analyze the Time-Tagged Event (TTE) data from observations of gamma ray bursts (GRBs) and soft gamma repeaters (SGRs) by the Burst and Transient Source Experiment (BATSE). These data provide the best available time resolution for GRBs and SGRs. We have performed an extensive search for weak periodic signals in the frequency range 400 Hz to 2500 Hz using the burst records for 2203 GRBs and 152 SGR flares. The study employs the Rayleigh power as a test statistic to evaluate the evidence for periodic emissions. We find no evidence of periodic emissions from these events at these frequencies. In all but a very few cases the maximum power values obtained are consistent with what would be expected by chance from a non-periodic signal. In those few instances where there is marginal evidence for periodicity there are problems with the data that cast doubt on the reality of the signal. For classical GRBs, the largest Rayleigh power occurs in bursts whose TTE data appear to be corrupted. For SGRs, our largest Rayleigh power, with a significance of 1%, occurs in one record for SGR 1900+14 (at 2497 Hz), and in no other outbursts associated with this source; we thus consider it unlikely to represent detection of a real periodicity. From simulations, we deduce that the Rayleigh test would have detected significant oscillations with relative amplitude ~10% about half the time. Thus, we conclude that high frequency oscillations, if present, must have small relative amplitudes.

We develop the formalism required to study the nonlinear interaction of modes in rotating Newtonian stars in the weakly nonlinear regime. The formalism simplifies and extends previous treatments. At linear order, we elucidate and extend slightly a formalism due to Schutz, show how to decompose a general motion of a rotating star into a sum over modes, and obtain uncoupled equations of motion for the mode amplitudes under the influence of an external force. Nonlinear effects are added perturbatively via three-mode couplings. We describe a new, efficient way to compute the coupling coefficients, to zeroth order in the stellar rotation rate, using spin-weighted spherical harmonics. We apply this formalism to derive some properties of the coupling coefficients relevant to the nonlinear interactions of unstable r-modes in neutron stars, postponing numerical integrations of the coupled equations of motion to a later paper. From an astrophysical viewpoint, the most interesting result of this paper is that many couplings of r-modes to other rotational modes (modes with zero frequencies in the non-rotating limit) are small: either they vanish altogether because of various selection rules, or they vanish to lowest order in the angular velocity. In zero-buoyancy stars, the coupling of three r-modes is forbidden entirely and the coupling of two r-modes to one hybrid rotational mode vanishes to zeroth order in rotation frequency. In incompressible stars, the coupling of any three rotational modes vanishes to zeroth order in rotation frequency.

This review starts with a discussion of the hierarchy of scales, relevant to the description of superfluids in neutron stars, which motivates a subsequent elementary exposition of the Newtonian superfluid hydrodynamics. Starting from the Euler equations for a superfluid and a normal fluid we apply the tensor virial method to obtain the virial equations of the first, second, and third order and to compute their Eulerian perturbations. Special emphasis is put on the computation of perturbations of the new terms due to mutual gravitational attraction and mutual friction between the two fluids. The oscillation modes of superfluid Maclaurin spheroids are derived from the first and second order perturbed virial equations. We discuss two generic classes of oscillation modes which correspond to the co-moving and relative oscillations of two fluids. These modes decouple if the normal fluid is inviscid. We also discuss the mixing of these modes (when the normal fluid is viscous) and its effect on the dynamical and secular instabilities of the co-moving modes and their damping.

We elaborate on the recently proposed static brane world scenario, where the effective 4-D cosmological constant is exponentially small when parallel 3-branes are far apart. We extend this result to a compactified model with two positive tension branes. Besides an exponentially small effective 4-D cosmological constant, this model incorporates a Randall-Sundrum-like solution to the hierarchy problem. Furthermore, the exponential factors for the hierarchy problem and the cosmological constant problem obey an inequality that is satisfied in nature. This inequality implies that the cosmological constant problem can be explained if the hierarchy problem is understood. The basic idea generalizes to the multibrane world scenario. We discuss models with piecewise adjustable bulk cosmological constants (to be determined by the 5-dimensional Einstein equation), a key element of the scenario. We also discuss the global structure of this scenario and clarify the physical properties of the particle (Rindler) horizons that are present. Finally, we derive a 4-D effective theory in which all observers on all branes not separated by particle horizons measure the same Newton's constant and 4-D cosmological constant.

We study the perturbation modes of rotating superfluid ellipsoidal figures of equilibrium in the framework of the two-fluid superfluid hydrodynamics and Newtonian gravity. Our calculations focus on linear perturbations of background equilibria in which the two fluids move together, the total density is uniform, and the densities of the two components are proportional to one another, with ratios that are independent of position. The motions of the two fluids are coupled by mutual friction, as formulated by Khalatnikov. We show that there are two general classes of modes for small perturbations: one class in which the two fluids move together and the other in which there is relative motion between them. The former are identical to the modes found for a single fluid, except that the rate of viscous dissipation, when computed in the secular (or ``low Reynolds number'') approximation under the assumption of a constant kinematic viscosity, is diminished by a factor $f_N$, the fraction of the total mass in the normal fluid. The relative modes are completely new, and are studied in detail for a range of values for the phenomenological mutual friction coefficients, relative densities of the superfluid and normal components, and, for Roche ellipsoids, binary mass ratios. We find that there are no new secular instabilities connected with the relative motions of the two fluid components. Moreover, although the new modes are subject to viscous dissipation (a consequence of viscosity of the normal matter), they do not emit gravitational radiation at all.

We consider a model with two parallel (positive tension) 3-branes separated by a distance $L$ in 5-dimensional spacetime. If the interbrane space is anti-deSitter, or is not precisely anti-deSitter but contains no event horizons, the effective 4-dimensional cosmological constant seen by observers on one of the branes (chosen to be the visible brane) becomes exponentially small as $L$ grows large.

We establish constraints on the mass and abundance of black holes in the Galactic halo by determining their impact on globular clusters which are conventionally considered to be little evolved. Using detailed Monte Carlo simulations and simple analytic estimates, we conclude that, at Galactocentric radius R~8 kpc, black holes with masses M_bh >~(1-3) x 10^6 M_sun can comprise no more than a fraction f_bh ~ 0.025-0.05 of the total halo density. This constraint significantly improves those based on disk heating and dynamical friction arguments as well as current lensing results. At smaller radius, the constraint on f_bh strengthens, while, at larger radius, an increased fraction of black holes is allowed.

We develop a possible cosmology for a Universe with n additional spatial dimensions of variable scale, and an associated scalar field, the radion, which is distinct from the field responsible for inflation, the inflaton. Based on a particular ansatz for the effective potential for the inflaton and radion (which may emerge in string theory), we show that the early expansion of the Universe may proceed in three stages. First, the radion becomes trapped at a value much smaller than the size of the extra dimensions today. Second, the Universe expands exponentially, but with a Planck mass smaller than its present value. Because the Planck mass during inflation is small, density fluctuations in agreement with observations can arise naturally. Third, when inflation ends, the Universe reheats, and the radion becomes free to expand once more. During the third phase the Universe is ``radiation-dominated'' and tends toward a fixed-point evolutionary model in which the radion grows but the temperature remains unchanged. Ultimately, the radion becomes trapped once again at its present value, and a short period of exponential expansion, which we identify with the electroweak phase transition, ensues. Once this epoch is over, the Universe reheats to a temperature of order the electroweak scale, and the mature Universe evolves as in standard cosmological models. We show that the present day energy density in radions can be smaller than the closure density if the second inflationary epoch lasts about 8 e-foldings or more; the present-day radion mass turns out to be small (less than or of the order of an eV, depending on parameters). We argue that although our model envisages considerable time evolution in the Planck mass, substantial spatial fluctuations in Newton's constant are not produced.

We explore the possible role of evolution in the analysis of data on SNe Ia at cosmological distances. First, using a variety of simple sleuthing techniques, we find evidence that the properties of the high and low redshift SNe Ia observed so far differ from one another. Next, we examine the effects of including simple phenomenological models for evolution in the analysis. The result is that cosmological models and evolution are highly degenerate with one another, so that the incorporation of even very simple models for evolution makes it virtually impossible to pin down the values of $\Omega_M$ and $\Omega_\Lambda$, the density parameters for nonrelativistic matter and for the cosmological constant, respectively. Moreover, we show that if SNe Ia evolve with time, but evolution is neglected in analyzing data, then, given enough SNe Ia, the analysis hones in on values of $\Omega_M$ and $\Omega_\Lambda$ which are incorrect. Using Bayesian methods, we show that the probability that the cosmological constant is nonzero (rather than zero) is unchanged by the SNe Ia data when one accounts for the possibility of evolution, provided that we do not discriminate among open, closed and flat cosmologies a priori. The case for nonzero cosmological constant is stronger if the Universe is presumed to be flat, but still depends sensitively on the degree to which the peak luminosities of SNe Ia evolve as a function of redshift. The estimated value of $H_0$, however, is only negligibly affected by accounting for possible evolution.

We demonstrate that models of resonant cyclotron radiation transfer in a strong field (i.e. cyclotron scattering) can account for spectral lines seen at two epochs, denoted S1 and S2, in the Ginga data for GRB870303. Using a generalized version of the Monte Carlo code of Wang et al. (1988,1989b), we model line formation by injecting continuum photons into a static plane-parallel slab of electrons threaded by a strong neutron star magnetic field (~ 10^12 G) which may be oriented at an arbitrary angle relative to the slab normal. We examine two source geometries, which we denote "1-0" and "1-1," with the numbers representing the relative electron column densities above and below the continuum photon source plane. We compare azimuthally symmetric models, i.e. models in which the magnetic field is parallel to the slab normal, with models having more general magnetic field orientations. If the bursting source has a simple dipole field, these two model classes represent line formation at the magnetic pole, or elsewhere on the stellar surface. We find that the data of S1 and S2, considered individually, are consistent with both geometries, and with all magnetic field orientations, with the exception that the S1 data clearly favor line formation away from a polar cap in the 1-1 geometry, with the best-fit model placing the line-forming region at the magnetic equator. Within both geometries, fits to the combined (S1+S2) data marginally favor models which feature equatorial line formation, and in which the observer's orientation with respect to the slab changes between the two epochs. We interpret this change as being due to neutron star rotation, and we place limits on the rotation period.

It has recently been suggested that the maximum observed quasi-periodic oscillation (QPO) frequencies, $\nu_{max}$, for several low-mass X-ray binaries, particularly 4U 1820-30, correspond to the orbital frequency at the inner-most stable orbit of the accretion disk. This would imply that the neutron stars in these systems have masses $\go 2~M_\odot$, considerably larger than any well-measured neutron star mass. We suggest that the levelling off of $\nu_{QPO}$ may be also understood in terms of a steepening magnetic field which, although possibly dipolar at the stellar surface, is altered substantially by disk accretion, and presents a ``wall'' to the accretion flow that may be outside the innermost stable orbit. General relativistic effects add to the flattening of the $\nu_{QPO}-\mdot$ relation at frequencies below the Kepler frequency at the innermost stable orbit. We offer two other possible ways to reconcile the low value of $\nu_{max}$ ($\approx 1060$ Hz for 4U 1820-30) with a moderate neutron star mass, $\approx 1.4\msun$: at sufficiently large $\mdot$, either (i) the disk terminates in a very thin boundary layer near the neutron star surface, or (ii) $\nu_{QPO}$ is not the orbital frequency right at the inner edge of the disk, but rather at a somewhat larger radius, where the emissivity of the disk peaks.

We explore the possibility that massive black holes comprise a significant fraction of the dark matter of our galaxy by studying the dissolution of galactic globular clusters bombarded by them. In our simulations, we evolve the clusters along a sequence of King models determined by changes of state resulting from collisions with the black holes. The results divide naturally into regimes of `small' and `large' black hole mass. `Small' black holes do not destroy clusters in single collisions; their effect is primarily cumulative, leading to a relation between $\mbh$ and $\fhalo$, the fraction of the halo in black holes of mass $\mbh$, which is $\fhalo\mbh < $ constant (up to logarithmic corrections). For $\fhalo=1$, we find $\mbh \simless 10^{3} \msun$ by requiring survival of the same clusters studied by Moore (1993), who neglected cluster evolution, mass loss, and stochasticity of energy inputs in his estimates, but reached a similar conclusion. `Large' black holes may not penetrate a cluster without disrupting it; their effect is mainly catastrophic (close collisions), but also partly cumulative (distant collisions). In the large $\mbh$ limit, $\fhalo$ (but not $\mbh$) can be constrained by computing the probability that a cluster survives a combination of close, destructive encounters and distant, nondestructive encounters. We find that it is unlikely that $\fhalo \simgreat 0.3$ by requiring 50 per cent survival probability for Moore's clusters over $10^{10}$ years.

From considerations of spin evolution and kinematics in the galactic potential, we argue that the pulsars B1913+16, B1534+12, and B2127+11C may be younger than previously assumed, and we find that a lower bound on the formation and merger rate of close double neutron star binaries is $10^{-6.7} f_b^{-1} f_d^{-1}$ yr$^{-1}$, where $f_b$ is the beaming fraction and $f_d$ accounts for the possibility that the known NS-NS binaries are atypical of the underlying population (e.g., if most such binaries are born with shorter orbital periods). If we assume no prior knowledge of the detectable lifetimes of such systems, the rate could be as large as $\simeq 10^{-5.0} f_b^{-1} f_d^{-1}$ yr$^{-1}$. From both plausible bounds on $f_b$ and $f_d$, and a revision of the independently derived limit proposed by Bailes (1996), we argue that a firm upper bound is $10^{-4}$ yr$^{-1}$. We also present a unifying empirical overview of the spin-up of massive binary pulsars based on their distribution in spin period $P$ and spin-down rate $\dot P$, finding evidence for two distinct spin-up processes, one dominated by disk accretion, the other by wind accretion. We argue that the positions of binary pulsars in the $P$-$\dot P$ diagram can be understood if (1) there exists a Galactic population of pulsars in double neutron star systems with combinations of spin and orbital periods that have prevented their detection in surveys to date; (2) recycled pulsars in wide-orbit binaries are not born near the canonical spin-up line in the $P$-$\dot P$ diagram because they were predominantly spun up through wind accretion; and (3) there exists a disfavored evolutionary endpoint for radio pulsars with spin periods 5--30 ms and $\dot P > 10^{-19}$ s-s$^{-1}$.

We consider the precession of isolated neutron stars in which superfluid is not pinned to the stellar crust perfectly. In the case of perfect pinning, Shaham (1977) showed that there are no slowly oscillatory, long-lived modes. When the assumption of perfect pinning is relaxed, new modes are found that can be long-lived, but are expected to be damped rather than oscillatory, unless the drag force on moving superfluid vortex lines has a substantial component perpendicular to the direction of relative motion. The response of a neutron star to external torques, such as the spindown torque, is also treated. We find that when computing the response of a star to perturbations, assuming perfect coupling of superfluid to normal matter from the start can miss some effects.

We consider an ultra-relativistic wind consisting of electron-positron pairs and photons with the principal goal of finding the asymptotic Lorentz factor $\gamma_{\infty}$ for zero baryon number. The wind is assumed to originate at radius $r_i$ where it has a Lorentz factor $\gamma_i$ and a temperature $T_i$ sufficiently high to maintain pair equilibrium. As $r$ increases, $T$ decreases and becomes less than the temperature corresponding to the electron mass $m_e$, after which non-equilibrium effects become important. Further out in the flow the optical depth $\tau$ drops below one, but the pairs may still be accelerated by the photons until $\tau$ falls below $\sim 2\times10^{-5} \gamma_{i}^{3/4}$. Radiative transfer calculations show that only at this point do the radiation flux and pressure start to deviate significantly from their blackbody values. The acceleration of the pairs increases $\gamma$ by a factor $\sim 45$ as compared to its value at the photosphere; it is shown to approach $\gamma_{\infty} \sim 1.4\times 10^3 (r_i/10^6\mbox{cm})^{1/4} \gamma_{i}^{3/4} T_i/m_e$.

We present the first fully relativistic investigation of matter fallback in a supernova. We investigate spherically symmetric supernova fallback using a relativistic radiation hydrodynamics Lagrangian code that handles radiation transport in all regimes. Our goal is to answer the fundamental question: did SN1987A produce a black hole and, if so, when will the hole become detectable ? We compute the light curve, assuming that a black hole has been formed during the explosion, and compare it with the observations. Our preliminary calculations lack radioactive energy input and adopt a very simple chemical composition (pure hydrogen). As a result, our computed models cannot fit the observed data of SN1987A in detail. Nevertheless, we can show that, during the first hours, the accretion flow is self--regulated and the accretion luminosity stays very close to the Eddington limit. The light curve is completely dominated, during the first few weeks, by the emission of the stellar envelope thermal energy, and resembles that obtained in ``standard'' supernova theory. Only long after hydrogen recombination takes place is there even a chance to actually detect radiation from the accreting black hole above the emission of the expanding envelope. The presence of a black hole is thus not inconsistent with observations to date. Because of the exponential decay of the $^{44}$Ti radioactive heating rate, the date of the emergence of the black hole is not very sensitive to the actual parameters of the models and turns out to be about 1000 years. The bulk of the emission then is expected to be in the visible band, but will be unobservable with present instrumentation. We discuss the implications of our results in connection with the possible emergence of a black hole in other supernovae.

We present results from a fully relativistic investigation of matter fallback in a supernova, assuming that the central compact star is a black hole (BH). Our goal is to answer the following questions: can we infer from the light curve whether a stellar black hole (BH) has formed in the aftermath of a supernova explosion ? Did SN1897A produce a BH and, if so, when will the hole become detectable ? We show that it should not be possible to infer the presence of a BH in the remnant for another several hundred to a thousand years. In the case of SN1987A, we estimate the BH luminosity to be $L\simeq 5\times 10^{34}$ erg s$^{-1}$, which is well below the present day bolometric luminosity of the remnant ($\sim 10^{36}erg/s$; Suntzeff 1997).

A recent study of the rotation of the plane of polarization of light from 160 cosmological sources claims to find significant evidence for cosmological anisotropy. We point out methodological weaknesses of that study, and reanalyze the same data using Bayesian methods that overcome these problems. We find that the data always favor isotropic models for the distribution of observed polarizations over counterparts that have a cosmological anisotropy of the type advocated in the earlier study. Although anisotropic models are not completely ruled out, the data put strong lower limits on the length scale $\lambda$ (in units of the Hubble length) associated with the anisotropy; the lower limits of 95% credible regions for $\lambda$ lie between 0.43 and 0.62 in all anisotropic models we studied, values several times larger than the best-fit value of $\lambda \approx 0.1$ found in the earlier study. The length scale is not constrained from above. The vast majority of sources in the data are at distances closer than 0.4 Hubble lengths (corresponding to a redshift of $\approx$0.8); the results are thus consistent with there being no significant anisotropy on the length scale probed by these data.

We use Bayesian methods to analyze the distribution of gamma ray burst intensities reported in the Third BATSE Catalog (3B catalog) of gamma ray bursts, presuming the distribution of burst sources (``bursters'') is isotropic. We study both phenomenological and cosmological source distribution models, using Bayes's theorem both to infer unknown parameters in the models, and to compare rival models. We analyze the distribution of the time-averaged peak photon number flux, F, measured on both 64 ms and 1024 ms time scales, performing the analysis of data based on each time scale independently. Several of our findings differ from those of previous analyses that modeled burst detection less completely. In particular, we find that the width of the intrinsic luminosity function for bursters is unconstrained, and the luminosity function of the actually observed bursts can be extremely broad, in contrast to the findings of all previous studies. Useful constraints probably require observation of bursts significantly fainter than those visible to BATSE. We also find that the 3B peak flux data do not usefully constrain the redshifts of burst sources; useful constraints require the analysis of data beyond that in the 3B catalog (such as burst time histories), or data from brighter bursts than have been seen by BATSE (such as those observed by the Pioneer Venus Orbiter). In addition, we find that an accurate understanding of the peak flux distributions reported in the 3B almost certainly requires consideration of data on the temporal and spectral properties of bursts beyond that reported in the 3B catalog, and more sophisticated modeling than has so far been attempted.

We use Bayesian methods to study anisotropic models for the distribution of gamma ray burst intensities and directions reported in the Third BATSE Catalog (3B catalog) of gamma ray bursts. We analyze data obtained using both the 64 ms and 1024 ms measuring timescales. We study both purely local models in which burst sources (``bursters'') are presumed to be distributed in extended halos about the Galaxy and M31, and mixed models consisting of a cosmological population of standard candle bursters and a local population distributed throughout a standard Bahcall-Soneira dark matter halo with a 2 kpc core. We find that the purely local models we have studied can account for the 3B data as successfully as cosmological models, provided one considers halos with core sizes significantly larger than those used to model the distribution of dark matter. We infer core sizes for the halo distribution that are smaller than one might expect based on popular semiquantitative arguments, and show why such arguments can lead to unwarranted conclusions. We also find that the 3B data do not constrain the width of power-law luminosity functions for burst sources. Our analysis of mixed models finds two families of models that can successfully account for the data: models with up to 20% of observed bursts in a bright local population visible to ~ 50 kpc; and models with up to 50% of observed bursts in a dim local population visible only nearby (to less than a disk scale height). These models fit as well or better than purely cosmological models. They indicate that a surprisingly large local, anisotropic component could be present whose size is comparable to the sizes of hypothetical classes of bursts inferred from analyses of temporal and spectral characteristics.

This article describes a new, fully adaptive Particle-Multiple-Mesh numerical simulation code developed primarily for simulations of small regions (such as a group of galaxies) in a cosmological context. It integrates the equations of motion of a set of particles subject to their mutual gravitational interaction and to an arbitrary external field. The interactions are computed using a hierarchy of nested grids constructed anew at each integration step to enhance the spatial resolution in high-density regions of interest. Significant effort has gone into supporting isolated boundary conditions at the top grid level. This makes our method also applicable to non-cosmological problems, at the cost of some complications which we discuss. We point out the implications of some differences between our approach and those of other authors of similar codes, in particular with respect to the handling of the interface between regions of different spatial resolution. We present a selection of tests performed to verify the correctness and performance of our implementation. The conclusion suggests possible further improvements in the areas of independent time steps and particle softening lengths.

Clumping of elementary dark matter in the Galaxy halo may be inevitable. If so, the nondetection of certain dark matter candidates could simply mean that the local halo density is low. Conversely, indirect detection of annihilation products could be facilitated, perhaps to an embarrassing degree.

Estimates indicate there may be $\simgreat 10^3$ low luminosity X-ray pulsars ($L \simless 10^{34} \rom{erg\;s^{-1}}$) in the Galaxy undergoing ``low-state'' wind accretion in Be/X-ray binary systems, and $\sim 10^8-10^9$ isolated neutron stars which may be accreting directly from the interstellar medium. Despite their low effective temperatures ($kT_e \simless 300 \rom{eV}$), low luminosity accreting neutron stars with magnetic fields $B \sim (0.7-7) \times 10^{12}$ G could emit a substantial fraction ($0.5 - 5\%$) of their total luminosity in a moderately broadened ($E/\Delta E \sim 2-4$) cyclotron {\it emission} line which peaks in the energy range $ \sim 5-20\rom{keV}$. The bulk of the thermal emission from these stars will be in the extreme ultraviolet/soft X-ray regime, which is subject to strong interstellar absorption and would be difficult to distinguish from spectra of other types of objects. In sharp contrast, the {\it nonthermal} cyclotron component predicted here will not be strongly absorbed, and consequently it may be the only distinguishing signature for the bulk of these low luminosity sources. We propose a search for this cyclotron emission feature in long pointed observations of the newly discovered candidate isolated neutron star MS0317.7-6477, and the Be/X-ray transient