Showing 1–50 of 70 results for author: Wasserman, I

Magnetohydrodynamic stability of magnetars in the ultrastrong field regime II: The crust

Authors: Peter B. Rau, Ira Wasserman

Abstract: We study the stability of Hall MHD with strong magnetic fields in which Landau quantization of electrons is important. We find that the strong-field Hall modes can be destabilized by the dependence of the differential magnetic susceptibility on magnetic field strength. This instability is studied using linear perturbation theory, and is found to have typical growth time of order $\lesssim 10^3$ yr… ▽ More We study the stability of Hall MHD with strong magnetic fields in which Landau quantization of electrons is important. We find that the strong-field Hall modes can be destabilized by the dependence of the differential magnetic susceptibility on magnetic field strength. This instability is studied using linear perturbation theory, and is found to have typical growth time of order $\lesssim 10^3$ yrs, with the growth time decreasing as a function of wavelength of the perturbation. The instability is self-limiting, turning off following a period of local field growth by a few percent of the initial value. Finite temperature is also shown to limit the instability, with sufficiently high temperatures eliminating it altogether. Alfvén waves can show similar unstable behaviour on shorter timescales. We find that Ohmic heating due to the large fields developed via the instability and magnetic domain formation is not large enough to account for observed magnetar surface temperatures. However, Ohmic heating is enhanced by the oscillatory differential magnetic susceptibility of Landau-quantized electrons, which could be important to magneto-thermal simulations of neutron star crusts. △ Less

Submitted 19 February, 2023; v1 submitted 11 October, 2022; originally announced October 2022.

Comments: 20 pages, 18 figures

Report number: INT preprint: INT-PUB-22-026

Journal ref: Monthly Notices of the Royal Astronomical Society 520, 1173 (2023)

Non-Axisymmetric Precession of Magnetars and Fast Radio Bursts

Authors: Ira Wasserman, James M. Cordes, Shami Chatterjee, Gauri Batra

Abstract: The repeating FRBs 180916.J0158 and 121102 are visible during periodically-occuring windows in time. We consider the constraints on internal magnetic fields and geometry if the cyclical behavior observed for FRB~180916.J0158 and FRB 121102 is due to precession of magnetars. In order to frustrate vortex line pinning we argue that internal magnetic fields must be stronger than about $10^{16}$ Gauss,… ▽ More The repeating FRBs 180916.J0158 and 121102 are visible during periodically-occuring windows in time. We consider the constraints on internal magnetic fields and geometry if the cyclical behavior observed for FRB~180916.J0158 and FRB 121102 is due to precession of magnetars. In order to frustrate vortex line pinning we argue that internal magnetic fields must be stronger than about $10^{16}$ Gauss, which is large enough to prevent superconductivity in the core and destroy the crustal lattice structure. We conjecture that the magnetic field inside precessing magnetars has three components, (1) a dipole component with characteristic strength $\sim 10^{14}$ Gauss; (2) a toroidal component with characteristic strength $\sim 10^{15}-10^{16}$ Gauss which only occupies a modest fraction of the stellar volume; and (3) a disordered field with characteristic strength $\sim 10^{16}$ Gauss. The disordered field is primarily responsible for permitting precession, which stops once this field component decays away, which we conjecture happens after $\sim 1000$ years. Conceivably, as the disordered component damps bursting activity diminishes and eventually ceases. We model the quadrupolar magnetic distortion of the star, which is due to its ordered components primarily, as triaxial and very likely prolate. We address the question of whether or not the spin frequency ought to be detectable for precessing, bursting magnetars by constructing a specific model in which bursts happen randomly in time with random directions distributed in or between cones relative to a single symmetry axis. Within the context of these specific models, we find that there are precession geometries for which detecting the spin frequency is very unlikely. △ Less

Submitted 27 July, 2021; originally announced July 2021.

Comments: 29 pages, 6 figures

Empirical Assessment of Aperiodic and Periodic Radio Bursts from Young Precessing Magnetars

Authors: J. M. Cordes, I. Wasserman, Shami Chatterjee, Gauri Batra

Abstract: We analyze the slow periodicities identified in burst sequences from FRB 121102 and FRB 180916 with periods of about 16 and 160 d, respectively, while also addressing the absence of any fast periodicity that might be associated with the spin of an underlying compact object. Both phenomena can be accounted for by a young, highly magnetized, precessing neutron star that emits beamed radiation with s… ▽ More We analyze the slow periodicities identified in burst sequences from FRB 121102 and FRB 180916 with periods of about 16 and 160 d, respectively, while also addressing the absence of any fast periodicity that might be associated with the spin of an underlying compact object. Both phenomena can be accounted for by a young, highly magnetized, precessing neutron star that emits beamed radiation with significant imposed phase jitter. Sporadic narrow-beam emission into an overall wide solid angle can account for the necessary phase jitter, but the slow periodicities with 25 to 55 % duty cycles constrain beam traversals to be significantly smaller. Instead, phase jitter may result from variable emission altitudes that yield large retardation and aberration delays. A detailed arrival-time analysis for triaxial precession includes wobble of the radio beam and the likely larger, cyclical torque resulting from the changes in the spin-magnetic moment angle. These effects will confound identification of the fast periodicity in sparse data sets longer than about a quarter of a precession cycle unless fitted for and removed as with orbital fitting. Stochastic spin noise, likely to be much larger than in radio pulsars, may hinder detection of any fast-periodicity in data spans longer than a few days. These decoherence effects will dissipate as FRB sources age, so they may evolve into objects with properties similar to Galactic magnetars. △ Less

Submitted 20 April, 2022; v1 submitted 27 July, 2021; originally announced July 2021.

Comments: 24 pages, 12 figures

Journal ref: J. M. Cordes, I. Wasserman, S. Chatterjee, G. Batra 2022, Astrophysical Journal, 929, 97

Magnetohydrodynamic stability of magnetars in the ultrastrong field regime I: The core

Authors: Peter B. Rau, Ira Wasserman

Abstract: We study magnetohydrodynamic stability of neutron star core matter composed of neutrons, protons and leptons threaded by a magnetar-strength magnetic field $10^{14}$--$10^{17}$ G, where quantum electrodynamical effects and Landau quantization of fermions are important. Stability is determined using the Friedman--Schutz formalism for the canonical energy of fluid perturbations, which we calculate f… ▽ More We study magnetohydrodynamic stability of neutron star core matter composed of neutrons, protons and leptons threaded by a magnetar-strength magnetic field $10^{14}$--$10^{17}$ G, where quantum electrodynamical effects and Landau quantization of fermions are important. Stability is determined using the Friedman--Schutz formalism for the canonical energy of fluid perturbations, which we calculate for a magnetizable fluid with $H\neq B$. Using this and the Euler--Heisenberg--Fermi--Dirac Lagrangian for a strongly magnetized fluid of Landau-quantized charged fermions, we calculate the local stability criteria for a fluid slab as a stand-in for a segment of a neutron star core, accounting for magnetic and composition gradient buoyancy. The slab is threaded by a field orthogonal to the gravitational field, the Cartesian analogy to a toroidal field. We find that, for sufficiently strong fields $B\gtrsim10^{15}$ G, the magnetized fluid is unstable to a magnetosonic-type instability with growth times of order $10^{-3}$ s. The instability is triggered by sharp changes in the second-order field derivative of the Euler--Heisenberg--Fermi--Dirac Lagrangian which occur where additional Landau levels start being populated. These sharp changes are divergent at zero temperature, but are finite for nonzero temperature, so realistic neutron star core temperatures $5\times10^7$ K$<T<5\times10^8$ K are used. We conjecture that this mechanism could promote the formation of magnetic domains as predicted by Blandford and Hernquist (1982) and Suh and Mathews (2010). △ Less

Submitted 24 August, 2021; v1 submitted 17 April, 2021; originally announced April 2021.

Comments: 28 pages. Matches published version, including supplemental material

Journal ref: MNRAS 506, 4632 (2021)

Relativistic finite temperature multifluid hydrodynamics in a neutron star from a variational principle

Authors: Peter B. Rau, Ira Wasserman

Abstract: 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 hy… ▽ More 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. △ Less

Submitted 10 September, 2020; v1 submitted 16 April, 2020; originally announced April 2020.

Comments: 27 pages, version accepted by Physical Review D

Journal ref: Phys. Rev. D 102, 063011 (2020)

Compressional modes in two-superfluid neutron stars with leptonic buoyancy

Authors: P. B. Rau, I. Wasserman

Abstract: 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 ad… ▽ More 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äisälä 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. △ Less

Submitted 12 October, 2018; v1 submitted 23 February, 2018; originally announced February 2018.

Comments: 19 pages, 12 figures

Journal ref: MNRAS 481, 4427-4444 (2018)

Lensing of Fast Radio Bursts by Plasma Structures in Host Galaxies

Authors: J. M. Cordes, I. Wasserman, J. W. T. Hessels, T. J. W. Lazio, S. Chatterjee, R. S. Wharton

Abstract: 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… ▽ More 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~μ$s to tens of ms and they will show different apparent dispersion measures, $δ{\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. △ Less

Submitted 19 March, 2017; originally announced March 2017.

Comments: 11 pages, 7 figures, submitted to the Astrophysical Journal

Radio Wave Propagation and the Provenance of Fast Radio Bursts

Authors: J. M. Cordes, R. S. Wharton, L. G. Spitler, S. Chatterjee, I. Wasserman

Abstract: 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 $τ$-DM relation for Galactic pulsars and use it as a benchmark for discussing FRB scattering. Though scattering times $τ$ for FRBs are large in the majority… ▽ More 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 $τ$-DM relation for Galactic pulsars and use it as a benchmark for discussing FRB scattering. Though scattering times $τ$ 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 $τ$ and DM for FRBs suggests that the intergalactic medium (IGM) cannot account for both $τ$ and DM. We therefore consider mixed models involving the IGM and host galaxies. If the IGM contributes significantly to DM while host galaxies dominate $τ$, the scattering deficit with respect to the mean Galactic trend can be explained with a $τ$-DM relation in the host that matches that for the Milky Way. However, it is possible that hosts dominate both $τ$ 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. △ Less

Submitted 19 May, 2016; originally announced May 2016.

Comments: 15 pages, 10 figures, submitted to ApJ

Magnetic Field Amplification via Protostellar Disc Dynamos

Authors: Sergei Dyda, Richard V. E. Lovelace, Galina V. Ustyugova, Alexander V. Koldoba, Ira Wasserman

Abstract: We numerically investigate the generation of a magnetic field in a protostellar disc via an $αΩ$-dynamo and the resulting magnetohydrodynamic (MHD) driven outflows. We find that for small values of the dimensionless dynamo parameter $α_d$ the poloidal field grows exponentially at a rate $σ\propto Ω_K \sqrt{α_d}$, before saturating to a value $\propto \sqrt{α_d}$. The dynamo excites dipole and octu… ▽ More We numerically investigate the generation of a magnetic field in a protostellar disc via an $αΩ$-dynamo and the resulting magnetohydrodynamic (MHD) driven outflows. We find that for small values of the dimensionless dynamo parameter $α_d$ the poloidal field grows exponentially at a rate $σ\propto Ω_K \sqrt{α_d}$, before saturating to a value $\propto \sqrt{α_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 $αΩ$-dynamos may be responsible for generating magnetic fields strong enough to launch observed outflows. △ Less

Submitted 2 March, 2018; v1 submitted 30 July, 2015; originally announced July 2015.

Comments: 13 pages, 8 figures

Supergiant Pulses from Extragalactic Neutron Stars

Authors: J. M. Cordes, Ira Wasserman

Abstract: 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… ▽ More 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. △ Less

Submitted 12 November, 2015; v1 submitted 4 January, 2015; originally announced January 2015.

Comments: Five figures. New version resubmitted to MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society 2016 457 (1): 232-257

Nonlinear Development of the R Mode Instability and the Maximum Rotation Rate of Neutron Stars

Authors: Ruxandra Bondarescu, Ira Wasserman

Abstract: 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… ▽ More 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. △ Less

Submitted 5 September, 2013; v1 submitted 10 May, 2013; originally announced May 2013.

Comments: 14 pages, 3 figures. Accepted by ApJ

Journal ref: ApJ 778, 9, 2013

Poloidal Magnetic Fields In Superconducting Neutron Stars

Authors: K. T. Henriksson, I. Wasserman

Abstract: 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 tha… ▽ More 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ε_b/3\simeq 3\times 10^{12}$ G, where $H_b$ is the magnetic field at the outer edge of the core and $ε_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}$. △ Less

Submitted 26 September, 2014; v1 submitted 23 December, 2012; originally announced December 2012.

Comments: 17 pages, 4 figures; v4 revised version accepted for publication

Journal ref: MNRAS (June 1, 2013) 431 (4): 2986-3002

Luminosity distance in Swiss cheese cosmology with randomized voids and galaxy halos

Authors: Éanna É. Flanagan, Naresh Kumar, Ira Wasserman

Abstract: 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 she… ▽ More 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. △ Less

Submitted 16 July, 2012; originally announced July 2012.

Comments: 16 pages, 13 figures

Multilevel Bayesian framework for modeling the production, propagation and detection of ultra-high energy cosmic rays

Authors: Kunlaya Soiaporn, David Chernoff, Thomas Loredo, David Ruppert, Ira Wasserman

Abstract: 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… ▽ More 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. △ Less

Submitted 28 November, 2013; v1 submitted 20 June, 2012; originally announced June 2012.

Comments: Published in at http://dx.doi.org/10.1214/13-AOAS654 the Annals of Applied Statistics (http://www.imstat.org/aoas/) by the Institute of Mathematical Statistics (http://www.imstat.org)

Report number: IMS-AOAS-AOAS654

Journal ref: Annals of Applied Statistics 2013, Vol. 7, No. 3, 1249-1285

Guilt by Association: Finding Cosmic Ray Sources Using Hierarchical Bayesian Clustering

Authors: Kunlaya Soiaporn, David Chernoff, Thomas Loredo, David Ruppert, Ira Wasserman

Abstract: 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 phenome… ▽ More 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. △ Less

Submitted 15 June, 2012; originally announced June 2012.

Luminosity distance in Swiss cheese cosmology with randomized voids. II. Magnification probability distributions

Authors: Éanna É. Flanagan, Naresh Kumar, Ira Wasserman, R. Ali Vanderveld

Abstract: 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 Λ-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 densit… ▽ More 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 Λ-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. △ Less

Submitted 14 November, 2011; v1 submitted 8 September, 2011; originally announced September 2011.

Journal ref: Phys.Rev.D85:023510,2012

Lemaitre-Tolman-Bondi cosmological models, smoothness, and positivity of the central deceleration parameter

Authors: R. Ali Vanderveld, Eanna E. Flanagan, Ira Wasserman

Abstract: 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… ▽ More 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. △ Less

Submitted 7 May, 2009; v1 submitted 28 April, 2009; originally announced April 2009.

Comments: 4 pages, no figures, revtex4

Modification to the Luminosity Distance Redshift Relation in Modified Gravity Theories

Authors: Éanna É. Flanagan, Eran Rosenthal, Ira M. Wasserman

Abstract: 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 c… ▽ More 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. △ Less

Submitted 11 February, 2009; v1 submitted 2 October, 2008; originally announced October 2008.

Comments: 15 pages, we discuss in more detail how the luminosity distance expression can be used to differentiate among various theories of gravity

Journal ref: Phys.Rev.D79:044032,2009

Spinning down newborn neutron stars: nonlinear development of the r-mode instability

Authors: Ruxandra Bondarescu, Saul A. Teukolsky, Ira Wasserman

Abstract: 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 i… ▽ More 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. △ Less

Submitted 13 March, 2009; v1 submitted 19 September, 2008; originally announced September 2008.

Comments: 23 pages, 13 figures, 1 table. Submitted to Phys. Rev. D. Detectability discussion expanded. Includes referee input

Journal ref: Phys.Rev.D79:104003,2009

Luminosity distance in "Swiss cheese" cosmology with randomized voids: I. Single void size

Authors: R. Ali Vanderveld, Eanna E. Flanagan, Ira Wasserman

Abstract: 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 inho… ▽ More 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. △ Less

Submitted 30 September, 2008; v1 submitted 7 August, 2008; originally announced August 2008.

Comments: 6 pages, 2 figures, revtex4. v2: final version accepted by PRD, minor changes that do not qualitatively affect results

Journal ref: Phys.Rev.D78:083511,2008

Cosmic Strings as the Source of Small-Scale Microwave Background Anisotropy

Authors: Levon Pogosian, S. -H. Henry Tye, Ira Wasserman, Mark Wyman

Abstract: 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… ▽ More 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. △ Less

Submitted 23 March, 2010; v1 submitted 7 April, 2008; originally announced April 2008.

Comments: 5 pages, 2 figures; improved fitting formula at 1000<l<3500, typos fixed

Journal ref: JCAP 0902:013,2009

Systematic corrections to the measured cosmological constant as a result of local inhomogeneity

Authors: R. Ali Vanderveld, Eanna E. Flanagan, Ira Wasserman

Abstract: 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… ▽ More 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. △ Less

Submitted 19 June, 2007; v1 submitted 13 June, 2007; originally announced June 2007.

Comments: 28 pages, 3 figures, revtex4; v2: corrected comments and the section on previous work; v3: clarified wording. References added

Journal ref: Phys.Rev.D76:083504,2007

Toroidal Magnetic Fields in Type II Superconducting Neutron Stars

Authors: Taner Akgun, Ira Wasserman

Abstract: 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 p… ▽ More 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. △ Less

Submitted 1 November, 2007; v1 submitted 15 May, 2007; originally announced May 2007.

Comments: 32 pages, 6 figures; derivations shortened, comments and references added; accepted for publication in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.383:1551-1580,2008

Spin Evolution of Accreting Neutron Stars: Nonlinear Development of the R-mode Instability

Authors: Ruxandra Bondarescu, Saul A. Teukolsky, Ira Wasserman

Abstract: 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… ▽ More 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 $ν_{max} \sim 800 Hz [S_{ns}/(M_{1.4} R_6)]^{4/11} T_8^{-2/11}$. We find that for $ν> 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. △ Less

Submitted 21 September, 2007; v1 submitted 5 April, 2007; originally announced April 2007.

Comments: 18 pages, 9 figures. Error Corrected. Includes referee input. Corrected References

Journal ref: Phys.Rev.D76:064019,2007

On vector mode contribution to CMB temperature and polarization from local strings

Authors: Levon Pogosian, Ira Wasserman, Mark Wyman

Abstract: 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… ▽ More 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μ< 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. △ Less

Submitted 6 April, 2006; originally announced April 2006.

Comments: 6 pages, 2 figures. Our code is available at http://physics.syr.edu/~lepogosi/cmbact.html

Mimicking Dark Energy with Lemaitre-Tolman-Bondi Models: Weak Central Singularities and Critical Points

Authors: R. Ali Vanderveld, Eanna E. Flanagan, Ira Wasserman

Abstract: 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… ▽ More 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. △ Less

Submitted 14 July, 2006; v1 submitted 21 February, 2006; originally announced February 2006.

Comments: 18 pages, 4 figures, revtex4; very minor revisions

Journal ref: Phys.Rev. D74 (2006) 023506

Rotational Broadening of Atomic Spectral Features from Neutron Stars

Authors: P. Chang, S. Morsink, L. Bildsten, I. Wasserman

Abstract: 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$α$ line formed above the photosphere of a bursting NS is intrinsically broad. We now include rotational broadening withi… ▽ More 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$α$ 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. △ Less

Submitted 8 November, 2005; originally announced November 2005.

Comments: 4 pages, 4 figures, to appear in ApJ Letters

Journal ref: Astrophys.J.636:L117-L120,2006

Precession of the Isolated Neutron Star PSR B1828-11

Authors: Taner Akgun, Bennett Link, Ira Wasserman

Abstract: 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 ri… ▽ More 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. △ Less

Submitted 30 September, 2005; v1 submitted 24 June, 2005; originally announced June 2005.

Comments: 22 pages, 14 figures, submitted to MNRAS; added references, corrected typos

Journal ref: Mon.Not.Roy.Astron.Soc.365:653-672,2006

Spherical Gravitational Collapse of Annihilating Dark Matter and the Minimum Mass of CDM Black Holes

Authors: R. Ali Vanderveld, Ira Wasserman

Abstract: Spherical gravitational collapse of a cold gas of annihilating particles involves a competition between the free-fall rate $\propto\sqrtρ$ and the (s-wave) annihilation rate $\proptoρ$. 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… ▽ More Spherical gravitational collapse of a cold gas of annihilating particles involves a competition between the free-fall rate $\propto\sqrtρ$ and the (s-wave) annihilation rate $\proptoρ$. 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πG^3M^2\lesssim\rhoann$, or $M\gtrsim\Mann\equiv (3/32π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 $ν{\barν}$ 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. △ Less

Submitted 25 May, 2005; originally announced May 2005.

Journal ref: JCAP 0511 (2005) 004

Formation of Resonant Atomic Lines during Thermonuclear Flashes on Neutron Stars

Authors: P. Chang, L. Bildsten, I. Wasserman

Abstract: Motivated by the measurement of redshifted Fe H$α$ 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$α$, Ly… ▽ More Motivated by the measurement of redshifted Fe H$α$ 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$α$, Ly$α$ and P$α$ lines. From the observed equivalent width of the Fe H$α$ 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$α$ EW $\approx$ 15-20 eV (redshifted 11-15 eV) and the P$α$ EW $\approx$ 4-7 eV (redshifted 3-5 eV) when the H$α$ EW is 10 eV (redshifted 8 eV). The Ly$α$ 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. △ Less

Submitted 3 May, 2005; originally announced May 2005.

Comments: 10 pages, 8 figures, to be published in ApJ

Journal ref: Astrophys.J. 629 (2005) 998-1007

Scaling of Multi-Tension Cosmic Superstring Networks

Authors: S. -H. Henry Tye, Ira Wasserman, Mark Wyman

Abstract: 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… ▽ More 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. △ Less

Submitted 17 June, 2005; v1 submitted 23 March, 2005; originally announced March 2005.

Comments: v2: 30 pages, 6 figures, references added to match version submitted to Phys. Rev. D v3: typos corrected on pp. 7-8

Journal ref: Phys.Rev.D71:103508,2005; Erratum-ibid.D71:129906,2005

Bounds on Cosmic Strings from WMAP and SDSS

Authors: Mark Wyman, Levon Pogosian, Ira Wasserman

Abstract: 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 ($μ$), we also varied the amount of small-scale structure on the strings. Our results indicate that cosmic strings can account for up to 7 (1… ▽ More 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 ($μ$), 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μ< 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. △ Less

Submitted 7 April, 2006; v1 submitted 16 March, 2005; originally announced March 2005.

Comments: 11 pages, 7 figures. A few errors in the computer code used to calculated CMB anistotropy from strings are fixed, resulting in a somewhat tighter bound on Gμand an enhanced B-mode polarization. Details of the corrected errors and their implications can be found in astro-ph/0604141. Go to http://physics.syr.edu/~lepogosi/cmbact.html for our now-corrected cosmic string CMB and LSS code

Journal ref: Phys.Rev.D72:023513,2005; Erratum-ibid.D73:089905,2006

A Nonlinear Coupling Network to Simulate the Development of the r-mode Instablility in Neutron Stars II. Dynamics

Authors: Jeandrew Brink, Saul A. Teukolsky, Ira Wasserman

Abstract: 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 numeric… ▽ More 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. △ Less

Submitted 15 October, 2004; originally announced October 2004.

Comments: 18 Pages 14 Figures

Journal ref: Phys.Rev. D71 (2005) 064029

A Nonlinear Coupling Network to Simulate the Development of the r-mode Instablility in Neutron Stars I. Construction

Authors: Jeandrew Brink, Saul A. Teukolsky, Ira Wasserman

Abstract: 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… ▽ More 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. △ Less

Submitted 10 September, 2004; originally announced September 2004.

Comments: 15 pages 3 figures

Journal ref: Phys.Rev. D70 (2004) 124017

Bunching instability of rotating relativistic electron layers and coherent synchrotron radiation

Authors: Bjoern S. Schmekel, Richard V. E. Lovelace, Ira M. Wasserman

Abstract: 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 pulsar… ▽ More 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. △ Less

Submitted 18 March, 2005; v1 submitted 27 September, 2004; originally announced September 2004.

Comments: 15 pages, 9 figures, LaTeX; presented at the April Meeting in Denver, Colorado 2004; numerous typos corrected, one approximation removed, one new proof added, accepted for publication in Phys. Rev. E

Journal ref: Phys.Rev. E71 (2005) 046502

Nonlinear Couplings of R-modes: Energy Transfer and Saturation Amplitudes at Realistic Timescales

Authors: Jeandrew Brink, Saul A Teukolsky, Ira Wasserman

Abstract: 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 time… ▽ More 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. △ Less

Submitted 23 June, 2004; v1 submitted 22 June, 2004; originally announced June 2004.

Comments: 3 figures 4 pages Submitted to PRD

Journal ref: Phys.Rev. D70 (2004) 121501

Observational constraints on cosmic strings: Bayesian analysis in a three dimensional parameter space

Authors: Levon Pogosian, Mark Wyman, Ira Wasserman

Abstract: 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 pres… ▽ More 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. △ Less

Submitted 11 January, 2005; v1 submitted 10 March, 2004; originally announced March 2004.

Comments: v3: 4 pages, 5 figures, slight modifications to match published version

Journal ref: JCAP 0409:008,2004

Explosive Disruption of Polytropes: a One Dimensional Hydrodynamic Calculation

Authors: Mark Wyman, David Chernoff, Ira Wasserman

Abstract: 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 inaccu… ▽ More 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. △ Less

Submitted 27 July, 2004; v1 submitted 1 May, 2003; originally announced May 2003.

Comments: v2: 12 pages, 10 figures, MNRAS LaTeX Macros. Paper contents significantly extended, with one new section and several new figures v3: 13 pages. Version to appear in MNRAS. Text expanded, particularly sections 5 and 6; corrected typos; added two new references

Journal ref: Mon.Not.Roy.Astron.Soc. 354 (2004) 1053-1064

Observational Constraints on Cosmic String Production During Brane Inflation

Authors: Levon Pogosian, S. -H. Henry Tye, Ira Wasserman, Mark Wyman

Abstract: 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 contr… ▽ More 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μ\lesssim 3.5\times 10^{-7}(λ/0.25)\sqrt{B/0.1}$, where $λ$ 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. △ Less

Submitted 6 April, 2006; v1 submitted 22 April, 2003; originally announced April 2003.

Comments: 16 pages, 4 figures. A few errors in the computer code used to calculated CMB anistotropy from strings are fixed, resulting in a somewhat tighter bound on Gμand an enhanced B-mode polarization. Details of the corrected errors and their implications can be found in astro-ph/0604141

Journal ref: Phys.Rev.D68:023506,2003; Erratum-ibid.D73:089904,2006

Spin-perpendicular kicks from evanescent binaries formed in the aftermath of rotational core-collapse and the nature of the observed bimodal distribution of pulsar peculiar velocities

Authors: Monica Colpi, Ira Wasserman

Abstract: 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.… ▽ More 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. △ Less

Submitted 17 February, 2003; originally announced February 2003.

Comments: To appear in the proceedings of the International Workshop on " Pulsars, AXPs, and SGRs observed with BeppoSAX and other Observatories" (Marsala, Sicily, 23-25 September 2002); 6 pages, 1 figure

The Evolution Of LMC X-4 Flares: Evidence For Super-Eddington Radiation Oozing Through Inhomogeneous Polar Cap Accretion Flows ?

Authors: Dae-Sik Moon, Stephen S. Eikenberry, Ira M. Wasserman

Abstract: 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 du… ▽ More 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. △ Less

Submitted 9 December, 2002; originally announced December 2002.

Comments: 39 pages (including figures and tables), accepted for publication in ApJ

Journal ref: Astrophys.J. 586 (2003) 1280-1296

SMC X-1 As An Intermediate-Stage Flaring X-ray Pulsar

Authors: Dae-Sik Moon, Stephen S. Eikenberry, Ira M. Wasserman

Abstract: 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. W… ▽ More 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. △ Less

Submitted 19 September, 2002; originally announced September 2002.

Comments: 14 pages (5 figures), To appear in ApJ Letters

Journal ref: Astrophys.J. 582 (2003) L91-L94

Precession of Isolated Neutron Stars II: Magnetic Fields and Type II Superconductivity

Authors: Ira Wasserman

Abstract: 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 t… ▽ More 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. △ Less

Submitted 23 August, 2002; v1 submitted 20 August, 2002; originally announced August 2002.

Comments: 40 pages, 1 figure. Discussion added on vortex pinning and compatibility with glitch models. References added and corrected. Typo corrected (Eq. 58)

Journal ref: Mon.Not.Roy.Astron.Soc. 341 (2003) 1020

Formation of an evanescent proto-neutron star binary and the origin of pulsar kicks

Authors: Monica Colpi, Ira Wasserman

Abstract: 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 explod… ▽ More 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. △ Less

Submitted 23 August, 2002; v1 submitted 16 July, 2002; originally announced July 2002.

Comments: 11 pages, 3 Figures, To appear in The Astrophysical Journal

Journal ref: Astrophys.J. 581 (2002) 1271-1279

On the Signature of Short Distance Scale in the Cosmic Microwave Background

Authors: Gary Shiu, Ira Wasserman

Abstract: 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φ^2/M^4$ where $φ$ is the {\it classical background} of the inflaton field. Therefore, the imprints of short distance p… ▽ More 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φ^2/M^4$ where $φ$ 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. △ Less

Submitted 12 March, 2002; originally announced March 2002.

Journal ref: Phys.Lett. B536 (2002) 1-8

On the Degeneracy Inherent in Observational Determination of the Dark Energy Equation of State

Authors: Ira Wasserman

Abstract: 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, $Ω_M$, somehow can be found independently. A phenomenological model is… ▽ More 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, $Ω_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. △ Less

Submitted 12 March, 2002; v1 submitted 8 March, 2002; originally announced March 2002.

Comments: 12 pages, 3 figures. Several references added

Journal ref: Phys.Rev. D66 (2002) 123511

Saturation of the R-mode Instability

Authors: Phil Arras, Eanna E. Flanagan, Sharon M. Morsink, A. Katrin Schenk, Saul A. Teukolsky, Ira Wasserman

Abstract: 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 am… ▽ More 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. △ Less

Submitted 10 February, 2003; v1 submitted 18 February, 2002; originally announced February 2002.

Comments: 34 pages. 6 figures. Accepted at ApJ

Journal ref: Astrophys.J.591:1129-1151,2003

Search for high-frequency periodicities in time-tagged event data from gamma ray bursts and soft gamma repeaters

Authors: Adam T. Kruger, Thomas J. Loredo, Ira Wasserman

Abstract: 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… ▽ More 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. △ Less

Submitted 7 December, 2001; originally announced December 2001.

Comments: 22 pages, 7 figures, submitted to ApJ

Nonlinear mode coupling in rotating stars and the r-mode instability in neutron stars

Authors: A. Katrin Schenk, Phil Arras, Eanna E. Flanagan, Saul A. Teukolsky, Ira Wasserman

Abstract: 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… ▽ More 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. △ Less

Submitted 11 July, 2001; v1 submitted 23 January, 2001; originally announced January 2001.

Comments: 62 pages, no figures. Corrected error in computation of coupling coefficients, added new selection rule and an appendix on energy and angular momentum of modes

Journal ref: Phys.Rev. D65 (2002) 024001

The tensor virial method and its applications to self-gravitating superfluids

Authors: A. Sedrakian, I. Wasserman

Abstract: 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 t… ▽ More 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. △ Less

Submitted 30 January, 2001; originally announced January 2001.

Comments: 30 pages, 3 figures; Chapter contributed to "Physics of Neutron Star Interiors", Lecture Notes in Physics (Springer 2001), eds. D. Blaschke, N. K. Glendenning, A. Sedrakian