We propose that the periodic fast radio bursts of FRB 180916.J0158+65 are sourced by axion emission (mass $m_{a} \sim 10^{-14}$ eV) from cosmic superstrings. Some of the emitted axions are converted to photons by magnetic fields as they travel along the line of sight to Earth. An impulsive burst of axion emission generates a photon signal typically lasting for milliseconds and varying with frequency in the observed manner. We find a range of parameters in our cosmic string network model consistent with the properties of FRB 180916.J0158+65. We suggest followup gravitational wave observations to test our model.

The existence of cosmic superstrings may be probed by astronomical time domain surveys. When crossing the line of sight to point-like sources, strings produce a distinctive microlensing signature. We consider two avenues to hunt for a relic population of superstring loops: frequent monitoring of (1) stars in Andromeda, lensed by loops in the haloes of the Milky-Way and Andromeda and (2) supernovae at cosmological distances, lensed by loops in the intergalactic medium. We assess the potential of such experiments to detect and/or constrain strings with a range of tensions, $10^{-15} \lesssim G \mu/c^2 \lesssim 10^{-6}$. The practical sensitivity is tied to cadence of observations which we explore in detail. We forecast that high-cadence monitoring of $\sim 10^5$ stars on the far side of Andromeda over a year-long period will detect microlensing events if $G\mu/c^2 \sim 10^{-13}$, while $\sim 10^6$ stars will detect events if $10^{-13.5} <G\mu/c^2 < 10^{-11.5}$; the upper and lower bounds of the accessible tension range continue to expand as the number of stars rises. We also analyse the ability to reject models in the absence of fluctuations. While challenging, these studies are within reach of forthcoming time-domain surveys. Supernova observations can hypothetically constrain models with $10^{-12} < G\mu/c^2 < 10^{-6}$ without any optimisation of the survey cadence. However, the event rate forecast suggests it will be difficult to reject models of interest. As a demonstration, we use observations from the Pantheon Type Ia supernova cosmology data-set to place modest constraints on the number density of cosmic superstrings in a poorly tested region of the parameter space.

We describe the simulated data sample for the "Photometric LSST Astronomical Time Series Classification Challenge" (PLAsTiCC), a publicly available challenge to classify transient and variable events that will be observed by the Large Synoptic Survey Telescope (LSST), a new facility expected to start in the early 2020s. The challenge was hosted by Kaggle, ran from 2018 September 28 to 2018 December 17, and included 1,094 teams competing for prizes. Here we provide details of the 18 transient and variable source models, which were not revealed until after the challenge, and release the model libraries at https://doi.org/10.5281/zenodo.2612896. We describe the LSST Operations Simulator used to predict realistic observing conditions, and we describe the publicly available SNANA simulation code used to transform the models into observed fluxes and uncertainties in the LSST passbands (ugrizy). Although PLAsTiCC has finished, the publicly available models and simulation tools are being used within the astronomy community to further improve classification, and to study contamination in photometrically identified samples of type Ia supernova used to measure properties of dark energy. Our simulation framework will continue serving as a platform to improve the PLAsTiCC models, and to develop new models.

We describe work in progress by a collaboration of astronomers and statisticians developing a suite of Bayesian data analysis tools for extrasolar planet (exoplanet) detection, planetary orbit estimation, and adaptive scheduling of observations. Our work addresses analysis of stellar reflex motion data, where a planet is detected by observing the "wobble" of its host star as it responds to the gravitational tug of the orbiting planet. Newtonian mechanics specifies an analytical model for the resulting time series, but it is strongly nonlinear, yielding complex, multimodal likelihood functions; it is even more complex when multiple planets are present. The parameter spaces range in size from few-dimensional to dozens of dimensions, depending on the number of planets in the system, and the type of motion measured (line-of-sight velocity, or position on the sky). Since orbits are periodic, Bayesian generalizations of periodogram methods facilitate the analysis. This relies on the model being linearly separable, enabling partial analytical marginalization, reducing the dimension of the parameter space. Subsequent analysis uses adaptive Markov chain Monte Carlo methods and adaptive importance sampling to perform the integrals required for both inference (planet detection and orbit measurement), and information-maximizing sequential design (for adaptive scheduling of observations). We present an overview of our current techniques and highlight directions being explored by ongoing research.

Cosmic superstrings of string theory differ from conventional cosmic strings of field theory. We review how the physical and cosmological properties of the macroscopic string loops influence experimental searches for these relics from the epoch of inflation. The universe's average density of cosmic superstrings can easily exceed that of conventional cosmic strings having the same tension by two or more orders of magnitude. The cosmological behavior of the remnant superstring loops is qualitatively distinct because the string tension is exponentially smaller than the string scale in flux compactifications in string theory. Low tension superstring loops live longer, experience less recoil (rocket effect from the emission of gravitational radiation) and tend to cluster like dark matter in galaxies. Clustering enhances the string loop density with respect to the cosmological average in collapsed structures in the universe. The enhancement at the Sun's position is $\sim 10^5$. We develop a model encapsulating the leading order string theory effects, the current understanding of the string network loop production and the influence of cosmological structure formation suitable for forecasting the detection of superstring loops via optical microlensing, gravitational wave bursts and fast radio bursts. We evaluate the detection rate of bursts from cusps and kinks by LIGO- and LISA-like experiments. Clustering dominates rates for $G \mu < 10^{-11.9}$ (LIGO cusp), $G \mu<10^{-11.2}$ (LISA cusp), $G \mu < 10^{-10.6}$ (LISA kink); we forecast experimentally accessible gravitational wave bursts for $G \mu>10^{-14.2}$ (LIGO cusp), $G \mu>10^{-15}$ (LISA cusp) and $G \mu>10^{- 14.1}$ (LISA kink).

Pulsar timing arrays (PTAs) are used to search for long-wavelength gravitational waves (GWs) by monitoring a set of spin-stable millisecond pulsars. Most theoretical analyses assume that the relevant GW sources are much more distant from Earth than the pulsars comprising the array. Unlike ground- or solar system-based GW detectors, PTAs might well contain embedded GW sources. We derive the PTA response from sources at any distance, with a specific focus on GW bursts with memory (BWMs). We consider supernovae and compact binary mergers as potential Galactic BWM sources and evaluate the signature for an array with pulsars in globular clusters or in the Galactic center. Understanding the response of PTAs to nearby sources of BWM is a step towards investigating other more complex Galactic sources.

We explore the collision of two cylindrical bubbles in classical general relativity with a scalar field stress-energy tensor. Inside each bubble the field rests at a local minimum of the potential with non-negative energy density. Outside the field rests at zero potential, the global minimum. The calculation resolves the connection from the inner de-Sitter region to the asymptotically flat Minkowski spacetime. We choose initial conditions such that the two bubbles collide and study the full nonlinear evolution by means of a two-dimensional numerical simulation of Einstein's equations. The collision generates a strongly interacting region with spatially varying fields and potentials. These circumstances promote dynamical exploration of the potential's landscape. No horizon is present and the scalar curvature invariants eventually diverge. We speculate that Schwarzschild-like horizons will encompass only part of the complicated, interesting regions of spacetime in the analogous case of colliding spherical bubbles.

At its very beginning, the universe is believed to have grown exponentially in size via the mechanism of inflation. The almost scale-invariant density perturbation spectrum predicted by inflation is strongly supported by cosmological observations, in particular the cosmic microwave background radiation. However, the universe's precise inflationary scenario remains a profound problem for cosmology and for fundamental physics. String theory, the most-studied theory as the final physical theory of nature, should provide an answer to this question. Some of the proposals on how inflation is realized in string theory are reviewed. Since everything is made of strings, some string loops of cosmological sizes are likely to survive in the hot big bang that followed inflation. They appear as cosmic strings, which can have intricate properties. Because of the warped geometry in flux compactification of the extra spatial dimensions in string theory, some of the cosmic strings may have tensions substantially below the Planck or string scale. Such strings cluster in a manner similar to dark matter leading to hugely enhanced densities. As a result, numerous fossil remnants of the low tension cosmic strings may exist within the galaxy. They can be revealed through the optical lensing of background stars in the near future and studied in detail through gravitational wave emission. We anticipate that these cosmic strings will permit us to address central questions about the properties of string theory as well as the birth of our universe.

The discovery of a binary comprising a black hole (BH) and a millisecond pulsar (MSP) would yield insights into stellar evolution and facilitate exquisitely sensitive tests of general relativity. Globular clusters (GCs) are known to harbor large MSP populations and recent studies suggest that GCs may also retain a substantial population of stellar mass BHs. We modeled the formation of BH+MSP binaries in GCs through exchange interactions between binary and single stars. We found that in dense, massive clusters most of the dynamically formed BH+MSP binaries will have orbital periods of 2 to 10 days, regardless of the mass of the BH, the number of BHs retained by the cluster, and the nature of the GC's binary population. The size of the BH+MSP population is sensitive to several uncertain parameters, including the BH mass function, the BH retention fraction, and the binary fraction in GCs. Based on our models, we estimate that there are $0.6\pm0.2$ dynamically formed BH+MSP binaries in the Milky Way GC system, and place an upper limit on the size of this population of $\sim 10$. Interestingly, we find that BH+MSP binaries will be rare even if GCs retain large BH populations.

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.

Cosmic superstring loops within the galaxy microlens background point sources lying close to the observer-string line of sight. For suitable alignments, multiple paths coexist and the (achromatic) flux enhancement is a factor of two. We explore this unique type of lensing by numerically solving for geodesics that extend from source to observer as they pass near an oscillating string. We characterize the duration of the flux doubling and the scale of the image splitting. We probe and confirm the existence of a variety of fundamental effects predicted from previous analyses of the static infinite straight string: the deficit angle, the Kaiser-Stebbins effect, and the scale of the impact parameter required to produce microlensing. Our quantitative results for dynamical loops vary by O(1) factors with respect to estimates based on infinite straight strings for a given impact parameter. A number of new features are identified in the computed microlensing solutions. Our results suggest that optical microlensing can offer a new and potentially powerful methodology for searches for superstring loop relics of the inflationary era.

We present a new method to calculate formation of cosmological structure in the Newtonian limit. The method is based on Lagrangian perturbation theory plus two key theoretical extensions. One advance involves identifying and fixing a previously ignored gauge-like degree of freedom relating quantities calculated in LPT to those measured by a preferred Friedmann-Robertson-Walker (FRW) observer. Handling this connection between calculational and observer frames is physically essential and ensures a momentum conserving description. The second extension is to systematically re-expand the equations of motion to increase LPT's radius of convergence to the maximum future time prior to orbit crossing. The paper implements a complete algorithm and performs extensive "proof of principle" tests of the new method, including direct comparison to known solutions, evaluation of conserved quantities and formal convergence studies. All are satisfactory. We show convergence is exponential in grid size and Lagrangian order and polynomial in step size. There are three {\it powerful advantages} afforded by the new technique: (1) it employs a smooth representation of all fields and the results are not limited by particle induced shot-noise errors, (2) it permits the numerical error to be controlled by changing Lagrangian order and/or number of steps allowing, in principle, arbitrarily small errors to be achieved prior to orbit crossing and (3) it handles generic cold initial data (any periodic density and velocity fields, including those with initial rotational components). Together, these properties make the new technique well-suited to handle quasi-linear scales where analytic methods and/or numerical simulations fail to provide suitably accurate answers.

We model the formation of black hole-neutron star (BH-NS) binaries via dynamical interactions in globular clusters. We find that in dense, massive clusters, 16-61% of the BH-NS binaries formed by interactions with existing BH binaries will undergo mergers driven by the emission of gravitational radiation. If the BHs are retained by the cluster after merging with a NS, the BHs acquire subsequent NS companions and undergo several mergers. Thus, the merger rate depends critically upon whether or not the BH is retained by the cluster after the merger. Results from numerical relativity suggest that kick imparted to a ~7 M_sun BH after it merges with a NS will greatly exceed the cluster's escape velocity. In this case, the models suggest that the majority of BH-NS mergers in globular clusters occur within 4 Gyrs of the cluster's formation and would be unobservable by Advanced LIGO. For more massive BHs, on the other hand, the post merger kick is suppressed and the BH is retained. Models with 35 M_sun BHs predict Advanced LIGO detection rates in the range 0.04 - 0.7 per year. On the pessimistic end of this range, BH-NS mergers resulting from binary-single star interactions in globular clusters could account for an interesting fraction of all BH-NS mergers. On the optimistic end, this channel may dominate the rate of detectable BH-NS mergers.

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 investigate convergence of Lagrangian Perturbation Theory (LPT) by analyzing the model problem of a spherical homogeneous top-hat in an Einstein-deSitter background cosmology. We derive the formal structure of the LPT series expansion, working to arbitrary order in the initial perturbation amplitude. The factors that regulate LPT convergence are identified by studying the exact, analytic solution expanded according to this formal structure. The key methodology is to complexify the exact solution, demonstrate that it is analytic and apply well-known convergence criteria for power series expansions of analytic functions.This analysis fully explains the previously reported observation that LPT fails to predict the evolution of an underdense, open region beyond a certain time. It also implies the existence of other examples, including overdense, closed regions, for which LPT predictions should also fail. We show that this is indeed the case by numerically computing the LPT expansion in these problematic cases. The formal limitations to the validity of LPT expansion are considerably more complicated than simply the first occurrence of orbit crossings as is often assumed. Evolution to a future time generically requires re-expanding the solution in overlapping domains that ultimately link the initial and final times, each domain subject to its own convergence criterion. We demonstrate that it is possible to handle all the problematic cases by taking multiple steps (LPT re-expansion). We characterize how the leading order numerical error for a solution generated by LPT re- expansion varies with the choice of Lagrangian order and of time step size. Convergence occurs when the Lagrangian order increases and/or the time step size decreases in a simple, well-defined manner. We develop a recipe for time step control for LPT re-expansion based on these results.

We study the numerical solution of the non-relativistic Schr\"{o}dinger equation for two-electron atoms in ground and excited S-states using pseudospectral (PS) methods of calculation. The calculation achieves convergence rates for the energy, Cauchy error in the wavefunction, and variance in local energy that are exponentially fast for all practical purposes. The method requires three separate subdomains to handle the wavefunction's cusp-like behavior near the two-particle coalescences. The use of three subdomains is essential to maintaining exponential convergence. A comparison of several different treatments of the cusps and the semi-infinite domain suggest that the simplest prescription is sufficient. For many purposes it proves unnecessary to handle the logarithmic behavior near the three-particle coalescence in a special way. The PS method has many virtues: no explicit assumptions need be made about the asymptotic behavior of the wavefunction near cusps or at large distances, the local energy is exactly equal to the calculated global energy at all collocation points, local errors go down everywhere with increasing resolution, the effective basis using Chebyshev polynomials is complete and simple, and the method is easily extensible to other bound states. This study serves as a proof-of-principle of the method for more general two- and possibly three-electron applications.

The fragmentation of horizon crossing superstrings creates loops that accrete and bind to collapsing matter perturbations if the string tension is less than a critical value and if the loops are large-scale. For typical model parameters (velocity dispersion of newly formed loops, loop length distribution, etc.) the galactic loop density is greatly enhanced as compared to the universe's mean loop density. The limiting enhancement of loop (energy) density is 0.25-0.4 times the enhancement of cold dark matter (ratio of galaxy to universe's dark matter densities) and is fully achieved for dimensionless string tensions less than 10^-13 at galactic scales less than 100 kpc. The enhancement is substantial for any tension less than the critical value which is of order 10^-9. The fact that the string loop enhancement roughly tracks that of cold dark matter is a robust result for small tension and large-scale fragmentation. Experiments sensitive to the local Galactic population of loops, especially microlensing, should enjoy characteristic detection rates greatly boosted over homogeneous estimates.

Cosmic superstrings are produced towards the end of the brane inflation. If the string tension is low enough, loops tend to be relatively long-lived. The resultant string network is expected to contain many loops which are smaller than typical Galactic scales. Cosmic expansion damps the center of mass motion of the loops which then cluster like cold dark matter. Loops will lens stars within the galaxy and local group. We explore microlensing of stars as a tool to detect and to characterize some of the fundamental string and string network properties, including the dimensionless string tension $G \mu/c^2$ and the density of string loops within the Galaxy. As $G \mu \to 0$ the intrinsic microlensing rate diverges as $1/\sqrt{G \mu}$ but experimental detection will be limited by shortness of the lensing timescale and/or smallness of the bending angle which each vary $\propto G \mu$. We find that detection is feasible for a range of tensions. As an illustration, the planned optical astrometric survey mission, Gaia, should be able to detect numerous microlensing events for string networks with tensions $10^{-10} \simless G \mu \simless 10^{-6}$. A null result for optical microlensing implies $G \mu \simless 10^{-10}$. If lensing of a given source is observed it will repeat because the internal motions of the loop are relativistic but the center of mass motion may be much smaller, of order the halo velocity. This distinctive hallmark $\sim 1000$ repetitions, suggests a useful method for confirmation of a potential lensing detection.

We discuss the signature of a peculiar constellation of very massive stars at a projected distance of 2-3 pc around R136a. We discuss various scenarios for its possible origin, such as independent clusters, triggered star formation, supernovae, and ejections via dynamical interactions. If the latter scenario were the correct one this would have significant implications on the way to probe the conditions in dense stellar cores, and on the evolution of massive clusters in general.

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.

(Abridged) We infer the velocity distribution of radio pulsars by modelling their birth, evolution, and detection in large-scale 0.4 GHz pulsar surveys, and by comparing model distributions of measurable pulsar properties with survey data using a likelihood function. We test models that characterize a population's birth rate, luminosity, shutoff of radio emission, birth locations, and birth velocities. We infer that the radio beam luminosity (i) is comparable to the energy flux of relativistic particles in models for spin-driven magnetospheres, signifying that radio emission losses reach nearly 100% for the oldest pulsars; and (ii) scales approximately as sqrt(Edot) which, in magnetosphere models, is proportional to the voltage drop available for acceleration of particles. We find that a two-component velocity distribution with characteristic velocities of 90 km/s and 500 km/s is greatly preferred to any one-component distribution. We explore some consequences of the preferred birth velocity distribution: (i) roughly 50% of pulsars in the solar neighborhood will escape the Galaxy, while ~15% have velocities greater than 1000 km/s; (ii) an important low-velocity population exists that increases the fraction of neutron stars retained by globular clusters and is consistent with the number of old objects that accrete from the interstellar medium; (iii) under standard assumptions for supernova remnant expansion and pulsar spin-down, ~10% of pulsars younger than 20 kyr will appear to lie outside of their host remnants. Finally, we comment on the ramifications of our birth velocity distribution for binary survival and the population of inspiraling binary neutron stars relevant to some GRB models and potential sources for LIGO.

We study implications for the apparent alignment of the spin axes, proper-motions, and polarization vectors of the Crab and Vela pulsars. The spin axes are deduced from recent Chandra X-ray Observatory images that reveal jets and nebular structure having definite symmetry axes. The alignments indicate these pulsars were born either in isolation or with negligible velocity contributions from binary motions. We examine the effects of rotation and the conditions under which spin-kick alignment is produced for various models of neutron star kicks. If the kick is generated when the neutron star first forms by asymmetric mass ejection or/and neutrino emission, then the alignment requires that the protoneutron star possesses an original spin with period $P_s$ much less than the kick timescale, thus spin-averaging the kick forces. The kick timescale ranges from 100 ms to 10 s depending on whether the kick is hydrodynamically driven or neutrino-magnetic field driven. For hydrodynamical models, spin-kick alignment further requires the rotation period of an asymmetry pattern at the radius near shock breakout (>100 km) to be much less than ~100 ms; this is difficult to satisfy unless rotation plays a dynamically important role in the core collapse and explosion ($P_s\lo 1$ ms). Aligned kick and spin vectors are inherent to the slow process of asymmetric electromagnetic radiation from an off-centered magnetic dipole. We reassess the viability of this effect, correcting a factor of 4 error in Harrison and Tademaru's calculation that increases the size of the effect. To produce a kick velocity of order a few hundred km/s requires that the neutron star be born with an initial spin close to 1 ms and that spindown due to r-mode driven gravitational radiation be inefficient compared to standard magnetic braking.

We discuss high-resolution simulations of one-dimensional, plane-parallel shock waves with mean speeds between 150 and 240 km/s propagating into gas with Alfven velocities up to 40 km/s and outline the conditions under which these radiative shocks experience an oscillatory instability in the cooling length, shock velocity, and position of the shock front. We investigate two forms of postshock cooling: a truncated single power law and a more realistic piecewise power law. The degree of nonlinearity of the instability depends strongly on the cooling power law and the Alfven Mach number: for power-law indices \alpha < 0 typical magnetic field strengths may be insufficient either to stabilize the fundamental oscillatory mode or to prevent the oscillations from reaching nonlinear amplitudes.

We use astrometric, distance and spindown data on pulsars to: (1) estimate three-dimensional velocity components, birth distances from the galactic plane, and ages of individual objects; (2) determine the distribution of space velocities and the scale height of pulsar progenitors; (3) test spindown laws for pulsars; (4) test for correlations between space velocities and other pulsar parameters; and (5) place empirical requirements on mechanisms than can produce high velocity neutron stars. Our approach incorporates measurement errors, uncertainties in distances, deceleration in the Galactic potential, and differential galactic rotation. We find that the scale height of the progenitors is approximately 0.13 kpc, that the 3D velocities are distributed in two components with characteristic speeds of 175(+20,-30) km/s and 700(+200,-150) km/s representing 83% and 17% of the population respectively. These results are insensitive to the explicit relation of chronological and spindown ages. We infer that the most probable chronological ages are typically smaller than conventional spindown ages by factors as large as two. We assess mechanisms for producing high-velocity neutron stars in view of the derived velocity distribution function.

We study molecular hydrogen formation in and on solids. We construct a model with surface sites and bulk sites capable of describing (1) the motion and exchange of H and H_2 between surface and bulk, (2) the recombination of H and dissociation of H_2 in and on the solid, and (3) the injection of H from the gas phase and the loss of H and H_2 from the solid. The basic physical processes include thermally activated reactions, collisionally induced reactions, and tunneling reactions. Our main application is to the astronomical problem of H_2 formation on grains in space but the model has more general applicability. We investigate the steady-state H and H_2 concentrations in and on the solid when gas phase atoms or ions stick to the surface or penetrate the body of the grain. The model identifies ranges of physical parameters for which the solid becomes saturated (surface and/or bulk) with particles (H and/or H_2) and facilitates the calculation of the efficiency of molecule formation (the fraction of the gas phase atoms that leave as molecules). These solutions are highly degenerate in the sense that they depend only on a small number of dimensionless parameters. We find that a variety of recombination pathways operate under a broad range of conditions. As an example we study H_2 formation in and on carbon grains.

We study the field millisecond pulsar (MSP) population to infer its intrinsic distribution in spin period and luminosity and to determine its spatial distribution within the Galaxy. Our likelihood analysis on data from extant surveys (22 pulsars with periods <20 ms) accounts for important selection effects. We infer a minimum period cutoff P(min) > 0.65 ms (99% confidence), a period distribution proportional to P^{-2.0 +- 0.33} and a pseudo-luminosity distribution proportional to L_p^{-2.0 +- 0.2} (where L_p = flux density times distance^2, for L_p >= 1.1 mJy kpc^2). We find a vertical scale height 0.65{+0.16,-0.12} kpc. We use our results to estimate the total number and birthrate of MSPs in the disk of the Galaxy. We limit the density contribution of a diffuse halo-like component to <1% of the midplane value. The MSP velocity dispersion is smaller that that of young, long-period pulsars by about a factor of 5. Our best estimate of the 1D velocity kick that is unique to MSP evolution is approximately 40 km s^-1. We discuss the evolutionary relationship of MSPs and low-mass X-ray binaries and prospects for future searches for MSPs.

High-resolution simulations of the cooling regions of spherically symmetric supernova remnants demonstrate a strong radiative instability. This instability, whose presence is dependent on the shock velocity, causes large-amplitude fluctuations in the shock velocity. The fluctuations begin almost immediately after the radiative phase begins (upon shell formation) if the shock velocity lies in the unstable range; they last until the shock slows to speeds less than approximately 130 km/s. We find that shock-velocity fluctuations from the reverberations of waves within the remnant are small compared to those due to the instability. Further, we find (in plane-parallel simulations) that advected inhomogeneities from the external medium do not interfere with the qualitative nature of the instability-driven fluctuations. Large-amplitude inhomogeneities may alter the phases of shock-velocity fluctuations, but do not substantially reduce their amplitudes.

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.

We investigate a population of primordial binaries during the initial stage of evolution of a star cluster. For our calculations we assume that equal mass stars form rapidly in a tidally truncated gas cloud, that $\sim10$\% of the stars are in binaries and that the resulting star cluster undergoes an epoch of violent relaxation. We study the collisional interaction of the binaries and single stars, in particular, the ionization of the binaries and the energy exchange between binaries and single stars. We find that for large $N$ systems ($N > 10^3$) even the most violent beginning leaves the binary distribution function largely intact. Hence, the binding energy originally tied up in the cloud's protostellar pairs is preserved during the relaxation process and the binaries are available to interact at later times within the virialized cluster.

We investigate asymptotic convergence in the~$\Delta x \!\rightarrow\! 0$ limit as a tool for determining whether numerical computations involving shocks are accurate. We use one-dimensional operator-split finite-difference schemes for hydrodynamics with a von Neumann artificial viscosity. An internal-energy scheme converges to demonstrably wrong solutions. We associate this failure with the presence of discontinuities in the limiting solution. Our extension of the Lax-Wendroff theorem guarantees that certain conservative, operator-split schemes converge to the correct continuum solution. For such a total-energy scheme applied to the formation of a single shock, convergence of a Cauchy error approaches the expected rate slowly. We relate this slowness to the effect of varying diffusion, due to varying linear artificial-viscous length, on small-amplitude waves. In an appendix we discuss the scaling of shock-transition regions with viscous lengths, and exhibit several difficulties for attempts to make extrapolations.

Some bright cD galaxies in cluster cooling flows have Thomson optical depths exceeding 0.01. A few percent of their luminosity is scattered and appears as diffuse polarized emission. We calculate the scattering process for different geometric combinations of luminosity sources and scattering media. We apply our results to galaxies, with and without active nuclei, immersed in cooling flows. We model observations of NGC 1275 and M87 (without active nuclei) in the presence of sky and galactic background fluxes which hinder the measurement of the scattered light at optical wavelengths. Current instruments are unable to detect the scattered light from such objects. However, when a galaxy has an active nucleus of roughly the same luminosity as the remainder of the galaxy in V, both the total and polarized scattered intensity should observable on large scales (5--30kpc), meaning intensity levels greater than 1% of the background level. For typical AGN and galaxy spectral distributions, the scattering is most easily detected at short (U) wavelengths. We point out that a number of such cases will occur. We show that the radiation pattern from the central nuclear region can be mapped using the scattering. We also show that the scattered light can be used to measure inhomogeneities in the cooling flow.

We show how to measure cosmological parameters using observations of inspiraling binary neutron star or black hole systems in one or more gravitational wave detectors. To illustrate, we focus on the case of fixed mass binary systems observed in a single Laser Interferometer Gravitational-wave Observatory (LIGO)-like detector. Using realistic detector noise estimates, we characterize the rate of detections as a function of a threshold signal-to-noise ratio $\rho_0$, the Hubble constant $H_0$, and the binary ``chirp'' mass. For $\rho_0 = 8$, $H_0 = 100$ km/s/Mpc, and $1.4 \msun$ neutron star binaries, the sample has a median redshift of $0.22$. Under the same assumptions but independent of $H_0$, a conservative rate density of coalescing binaries ($8\times10^{-8}\,{\rm yr}^{-1}\,{\rm Mpc}^{-3}$) implies LIGO will observe $\sim 50\,{\rm yr}^{-1}$ binary inspiral events. The precision with which $H_0$ and the deceleration parameter $q_0$ may be determined depends on the number of observed inspirals. For fixed mass binary systems, $\sim 100$ observations with $\rho_0 = 10$ in the LIGO detector will give $H_0$ to 10\% in an Einstein-DeSitter cosmology, and 3000 will give $q_0$ to 20\%. For the conservative rate density of coalescing binaries, 100 detections with $\rho_0 = 10$ will require about 4~yrs.

We investigate the sensitivity of individual LIGO/VIRGO-like interferometers and the precision with which they can determine the characteristics of an inspiralling binary system. Since the two interferometers of the LIGO detector share nearly the same orientation, their joint sensitivity is similar to that of a single, more sensitive interferometer. We express our results for a single interferometer of both initial and advanced LIGO design, and also for the LIGO detector in the limit that its two interferometers share exactly the same orientation. We approximate the evolution of a binary system as driven exclusively by leading order quadrupole gravitational radiation. To assess the sensitivity, we calculate the rate at which sources are expected to be observed, the range to which they are observable, and the precision with which characteristic quantities describing the observed binary system can be determined. Assuming a conservative rate density for coalescing neutron star binary systems we expect that the advanced LIGO detector will observe approximately 69~yr${}^{-1}$ with an amplitude SNR greater than 8. Of these, approximately 7~yr${}^{-1}$ will be from binaries at distances greater than 950~Mpc. We explore the sensitivity of these results to a tunable parameter in the interferometer design (the recycling frequency). The optimum choice of the parameter is dependent on the goal of the observations, e.g., maximizing the rate of detections or maximizing the precision of measurement. We determine the optimum parameter values for these two cases.