Lutz-Kelker bias in pulsar parallax measurements

Lutz & Kelker showed that parallax measurements are systematically overestimated because they do not properly account for the larger volume of space that is sampled at smaller parallax values. We apply their analysis to neutron stars, incorporating the bias introduced by the intrinsic radio luminosity function and a realistic Galactic population model for neutron stars. We estimate the bias for all published neutron star parallax measurements and find that measurements with less than ~95% certainty, are likely to be significantly biased. Through inspection of historic parallax measurements, we confirm the described effects in optical and radio measurements, as well as in distance estimates based on interstellar dispersion measures. The potential impact on future tests of relativistic gravity through pulsar timing and on X-ray--based estimates of neutron star radii is briefly discussed.Comment: 9 pages, 3 tables, 1 figure. Accepted for publication in MNRA

Monte Carlo Neutrino Oscillations

We demonstrate that the effects of matter upon neutrino propagation may be recast as the scattering of the initial neutrino wavefunction. Exchanging the differential, Schrodinger equation for an integral equation for the scattering matrix S permits a Monte Carlo method for the computation of S that removes many of the numerical difficulties associated with direct integration techniques

THE EMERGING FOOD STORE DELICATESSEN: SOME PRELIMINARY ECONOMICS

Agribusiness,

Evidence in Virgo for the Universal Dark Matter Halo

A model is constructed for the mass and dynamics of M87 and the Virgo Cluster. Existing surface photometry of the galaxy, mass estimates from X-ray observations of the hot intracluster gas, and the velocity dispersions of early-type Virgo galaxies, all are used to constrain the run of dark matter density over radii to 2 Mpc in the cluster. The ``universal'' halo advocated by Navarro, Frenk, & White provides an excellent description of the combined data, as does a Hernquist profile. These models are favored over isothermal spheres, and their central structure is preferred to density cusps either much stronger or much weaker than r^{-1}. The galaxies and gas in the cluster trace its total mass distribution, the galaxies' velocity ellipsoid is close to isotropic, and the gas temperature follows the virial temperature profile of the dark halo. The virial radius and mass and the intracluster gas fraction of Virgo are evaluated.Comment: ApJ Letters in pres

Structural parameters for globular clusters in M31 and generalizations for the fundamental plane

The structures of globular clusters (GCs) reflect their dynamical states and past histories. High-resolution imaging allows the exploration of morphologies of clusters in other galaxies. Surface brightness profiles from new Hubble Space Telescope observations of 34 globular clusters in M31 are presented, together with fits of several different structural models to each cluster. M31 clusters appear to be adequately fit by standard King models, and do not obviously require alternate descriptions with relatively stronger halos, such as are needed to fit many GCs in other nearby galaxies. The derived structural parameters are combined with corrected versions of those measured in an earlier survey to construct a comprehensive catalog of structural and dynamical parameters for M31 GCs with a sample size similar to that for the Milky Way. Clusters in M31, the Milky Way, Magellanic Clouds, Fornax dwarf spheroidal and NGC 5128 define a very tight fundamental plane with identical slopes. The combined evidence for these widely different galaxies strongly reinforces the view that old globular clusters have near-universal structural properties regardless of host environment.Comment: AJ in press; 59 pages including 16 figure

Fatigue of notched fiber composite laminates. Part 1: Analytical model

A description is given of a semi-empirical, deterministic analysis for prediction and correlation of fatigue crack growth, residual strength, and fatigue lifetime for fiber composite laminates containing notches (holes). The failure model used for the analysis is based upon composite heterogeneous behavior and experimentally observed failure modes under both static and fatigue loading. The analysis is consistent with the wearout philosophy. Axial cracking and transverse cracking failure modes are treated together in the analysis. Cracking off-axis is handled by making a modification to the axial cracking analysis. The analysis predicts notched laminate failure from unidirectional material fatique properties using constant strain laminate analysis techniques. For multidirectional laminates, it is necessary to know lamina fatique behavior under axial normal stress, transverse normal stress and axial shear stress. Examples of the analysis method are given

Explicit Integration of the Full Symmetric Toda Hierarchy and the Sorting Property

We give an explicit formula for the solution to the initial value problem of the full symmetric Toda hierarchy. The formula is obtained by the orthogonalization procedure of Szeg\"{o}, and is also interpreted as a consequence of the QR factorization method of Symes \cite{symes}. The sorting property of the dynamics is also proved for the case of a generic symmetric matrix in the sense described in the text, and generalizations of tridiagonal formulae are given for the case of matrices with $2M+1$ nonzero diagonals.Comment: 13 pages, Latex

The Gamma Ray Pulsar Population

We apply a likelihood analysis to pulsar detections, pulsar upper limits, and diffuse background measurements from the OSSE and EGRET instruments on the Compton Gamma Ray Observatory to constrain the luminosity law for gamma-ray pulsars and some properties of the gamma-ray pulsar population. We find that the dependence of luminosity on spin period and dipole magnetic field is much steeper at OSSE than at EGRET energies (50-200 keV and >100 MeV, respectively), suggesting that different emission mechanisms are responsible for low- and high-energy gamma-ray emission. Incorporating a spin-down model and assuming a pulsar spatial distribution, we estimate the fraction of the Galactic gamma-ray background due to unidentified pulsars and find that pulsars may be an important component of the OSSE diffuse flux, but are most likely not important at EGRET energies. Using measurements of the diffuse background flux from these instruments, we are able to place constraints on the braking index, initial spin period, and magnetic field of the Galactic pulsar population. We are also able to constrain the pulsar birthrate to be between 1/(25 yr) and 1/(500 yr). Our results are based on a large gamma-ray beam, but they do not scale in a simple way with beam size. With our assumed beam size, the implied gamma-ray efficiency for the EGRET detections is no more than 20%. We estimate that about 20 of the 169 unidentified EGRET sources are probably gamma-ray pulsars. We use our model to predict the pulsar population that will be seen by future gamma-ray instruments and estimate that GLAST will detect roughly 750 gamma-ray pulsars as steady sources, only 120 of which are currently known radio pulsars.Comment: 32 pages, including figures. submitted to Ap

The impact of a stochastic gravitational-wave background on pulsar timing parameters

Gravitational waves are predicted by Einstein's theory of general relativity as well as other theories of gravity. The rotational stability of the fastest pulsars means that timing of an array of these objects can be used to detect and investigate gravitational waves. Simultaneously, however, pulsar timing is used to estimate spin period, period derivative, astrometric, and binary parameters. Here we calculate the effects that a stochastic background of gravitational waves has on pulsar timing parameters through the use of simulations and data from the millisecond pulsars PSR J0437--4715 and PSR J1713+0747. We show that the reported timing uncertainties become underestimated with increasing background amplitude by up to a factor of $\sim10$ for a stochastic gravitational-wave background amplitude of $A=5\times 10^{-15}$, where $A$ is the amplitude of the characteristic strain spectrum at one-year gravitational wave periods. We find evidence for prominent low-frequency spectral leakage in simulated data sets including a stochastic gravitational-wave background. We use these simulations along with independent Very Long Baseline Interferometry (VLBI) measurements of parallax to set a 2--sigma upper limit of $A\le9.1\times 10^{-14}$. We find that different supermassive black hole assembly scenarios do not have a significant effect on the calculated upper limits. We also test the effects that ultralow--frequency (10$^{-12}$--10$^{-9}$ Hz) gravitational waves have on binary pulsar parameter measurements and find that the corruption of these parameters is less than those due to $10^{-9}$--$10^{-7}$ Hz gravitational waves.Comment: 16 pages, 7 figures, accepted to MNRA