Source author record

Ira Wasserman

Ira Wasserman appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

16works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

16 published item(s)

preprint2021arXiv

Non-Axisymmetric Precession of Magnetars and Fast Radio Bursts

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.

preprint2020arXiv

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

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.

preprint2015arXiv

Supergiant Pulses from Extragalactic Neutron Stars

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.

preprint2013arXiv

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

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.

preprint2013arXiv

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

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.

preprint2012arXiv

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

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.

preprint2012arXiv

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

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.

preprint2011arXiv

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

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.

preprint2010arXiv

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

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.

preprint2006arXiv

Bounds on Cosmic Strings from WMAP and SDSS

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.

preprint2006arXiv

Observational Constraints on Cosmic String Production During Brane Inflation

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.

preprint2005arXiv

Precession of the Isolated Neutron Star PSR B1828-11

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.

preprint2005arXiv

Scaling of Multi-Tension Cosmic Superstring Networks

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.

preprint2004arXiv

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

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.

preprint2003arXiv

Saturation of the R-mode Instability

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.

preprint2001arXiv

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

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.