Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userSimulating pulsar glitches: an $N$-body solver for superfluid vortex motion in two dimensions
289
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
A rotating superfluid forms an array of quantized vortex lines which determine its angular velocity. The spasmodic evolution of the array under the influence of deceleration, dissipation, and pinning forces is thought to be responsible for the phenomenon of pulsar glitches, sudden jumps in the spin frequency of rotating neutron stars. We describe and implement an $N$-body method for simulating the motion of up to 5000 vortices in two dimensions and present the results of numerical experiments validating the method, including stability of a vortex ring and dissipative formation of an Abrikosov array. Vortex avalanches occur routinely in the simulations, when chains of unpinning events are triggered collectively by vortex-vortex repulsion, consistent with previous, smaller-scale studies using the Gross-Pitaevskii equation. The probability density functions of the avalanche sizes and waiting times are consistent with both exponential and log-normal distributions. We find weak correlations between glitch sizes and waiting times, consistent with astronomical data and meta-models of pulsar glitch activity as a state-dependent Poisson process or a Brownian stress-accumulation process, and inconsistent with a threshold-triggered stress-release model with a single, global stress reservoir. The spatial distribution of the effective stress within the simulation volume is analysed before and after a glitch.
rate research
Read More
The current-quadrupole gravitational-wave signal emitted during the spin-up phase of a pulsar glitch is calculated from first principles by modeling the vortex dynamics observed in recent Gross-Pitaevskii simulations of pinned, decelerating quantum condensates. Homogeneous and inhomogeneous unpinning geometries, representing creep- and avalanche-like glitches, provide lower and upper bounds on the gravitational wave signal strength respectively. The signal arising from homogeneous glitches is found to scale with the square root of glitch size, whereas the signal from inhomogeneous glitches scales proportional to glitch size. The signal is also computed as a function of vortex travel distance and stellar angular velocity. Convenient amplitude scalings are derived as functions of these parameters. For the typical astrophysical situation, where the glitch duration (in units of the spin period) is large compared to the vortex travel distance (in units of the stellar radius), an individual glitch from an object $1,rm{kpc}$ from Earth generates a wave strain of $10^{-24} [(Deltaomega/omega) / 10^{-7}] (omega/10^2 rm{rad s}^{-1})^3 (Delta r / 10^{-2} rm{m})^{-1}$, where $Delta r$ is the average distance travelled by a vortex during a glitch, $Deltaomega/omega$ is the fractional glitch size, and $omega$ is the pulsar angular velocity. The non-detection of a signal from the 2006 Vela glitch in data from the fifth science run conducted by the Laser Interferometer Gravitational-Wave Observatory implies that the glitch duration exceeds $sim 10^{-4},rm{ms}$. This represents the first observational lower bound on glitch duration to be obtained.
The scale-invariant glitch statistics observed in individual pulsars (exponential waiting-time and power-law size distributions) are consistent with a critical self-organization process, wherein superfluid vortices pin metastably in macroscopic domains and unpin collectively via nearest-neighbor avalanches. Macroscopic inhomogeneity emerges naturally if pinning occurs at crustal faults. If, instead, pinning occurs at lattice sites and defects, which are macroscopically homogeneous, we show that an alternative, noncritical self-organization process operates, termed coherent noise, wherein the global Magnus force acts uniformly on vortices trapped in a range of pinning potentials and undergoing thermal creep. It is found that vortices again unpin collectively, but not via nearest-neighbor avalanches, and that, counterintuitively, the resulting glitch sizes are scale invariant, in accord with observational data. A mean-field analytic theory of the coherent noise process, supported by Monte-Carlo simulations, yields a power-law size distribution, between the smallest and largest glitch, with exponent $a$ in the range $-2leq a leq 0$. When the theory is fitted to data from the nine most active pulsars, including the two quasiperiodic glitchers PSR J0537$-$6910 and PSR J0835$-$4510, it directly constrains the distribution of pinning potentials in the star, leading to two conclusions: (i) the potentials are broadly distributed, with the mean comparable to the standard deviation; and (ii) the mean potential decreases with characteristic age. An observational test is proposed to discriminate between nearest-neighbor avalanches and coherent noise.
Pulsars are known for their superb timing precision, although glitches can interrupt the regular timing behavior when the stars are young. These glitches are thought to be caused by interactions between normal and superfluid matter in the crust of the star. However, glitching pulsars such as Vela have been shown to require a superfluid reservoir that greatly exceeds that available in the crust. We examine a model in which glitches tap the superfluid in the core. We test a variety of theoretical superfluid models against the most recent glitch data and find that only one model can successfully explain up to 45 years of observational data. We develop a new technique for combining radio and X-ray data to measure pulsar masses, thereby demonstrating how current and future telescopes can probe fundamental physics such as superfluidity near nuclear saturation.
Pulsar-like compact stars provide us a unique laboratory to explore properties of dense matter at supra-nuclear densities. One of the models for pulsar-like stars is that they are totally composed of strangeons, and in this paper we studied the pulsar glitches in a strangeon star model. Strangeon stars would be solidified during cooling, and the solid stars would be natural to have glitches as the result of starquakes. Based on the starquake model established before, we proposed that when the starquake occurs, the inner motion of the star which changes the moment of inertia and has impact on the glitch sizes, is divided into plastic flow and elastic motion. The plastic flow which is induced in the fractured part of the outer layer, would move tangentially to redistribute the matter of the star and would be hard to recover. The elastic motion, on the other hand, changes its shape and would recover significantly. Under this scenario, we could understand the behaviors of glitches without significant energy releasing, including the Crab and the Vela pulsars, in an uniform model. We derive the recovery coefficient as a function of glitch size, as well as the time interval between two successive glitches as the function of the released stress. Our results show consistency with observational data under reasonable ranges of parameters. The implications on the oblateness of the Crab and the Vela pulsars are discussed.
Few statistically compelling correlations are found in pulsar timing data between the size of a rotational glitch and the time to the preceding glitch (backward waiting time) or the succeeding glitch (forward waiting time), except for a strong correlation between sizes and forward waiting times in PSR J0537-6910. This situation is counterintuitive, if glitches are threshold-triggered events, as in standard theories (e.g. starquakes, superfluid vortex avalanches). Here it is shown that the lack of correlation emerges naturally, when a threshold trigger is combined with secular stellar braking slower than a critical, calculable rate. The Pearson and Spearman correlation coefficients are computed and interpreted within the framework of a state-dependent Poisson process. Specific, falsifiable predictions are made regarding what objects currently targeted by long-term timing campaigns should develop strong size-waiting-time correlations, as more data are collected in the future.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
