-
Measuring glitch recoveries and braking indices with Bayesian model selection
Authors:
Yang Liu,
Michael J. Keith,
Danai Antonopoulou,
Patrick Weltevrede,
Benjamin Shaw,
Benjamin W. Stappers,
Andrew G. Lyne,
Mitchell B. Mickaliger,
Avishek Basu
Abstract:
For a selection of 35 pulsars with large spin-up glitches ($Δν/ν\geq10^{-6}$), which are monitored by the Jodrell Bank Observatory, we analyse 157 glitches and their recoveries. All parameters are measured consistently and we choose the best model to describe the post-glitch recovery based on Bayesian evidence. We present updated glitch epochs, sizes, changes of spin down rate, exponentially recov…
▽ More
For a selection of 35 pulsars with large spin-up glitches ($Δν/ν\geq10^{-6}$), which are monitored by the Jodrell Bank Observatory, we analyse 157 glitches and their recoveries. All parameters are measured consistently and we choose the best model to describe the post-glitch recovery based on Bayesian evidence. We present updated glitch epochs, sizes, changes of spin down rate, exponentially recovering components (amplitude and corresponding timescale) when present, as well as pulsars' second frequency derivatives and their glitch associated changes if detected. We discuss the different observed styles of post-glitch recovery as well as some particularly interesting sources. Several correlations are revealed between glitch parameters and pulsar spin parameters, including a very strong correlation between a pulsar's interglitch $|\ddotν|$ and $\dotν$, as well as between the glitch-induced spin-down rate change $Δ\dotν_{\rm p}$ that does not relax exponentially and $\dotν$. We find that the ratio $\left|Δ\dotν_{\mathrm{p}}/\ddotν\right|$ can be used as an estimate of glitch recurrence times, especially for those pulsars for which there are indications of a characteristic glitch size and interglitch waiting time. We calculate the interglitch braking index $n$ and find that pulsars with large glitches typically have $n$ greater than $3$, suggesting that internal torques dominate the rotational evolution between glitches. The external torque, e.g. from electromagnetic dipole radiation, could dominate the observed $\ddotν$ for the youngest pulsars ($\lesssim10^{4}\;\mathrm{yr}$), which may be expected to display $n\sim3$.
△ Less
Submitted 13 June, 2024;
originally announced June 2024.
-
A new small glitch in Vela discovered with a hidden Markov model
Authors:
L. Dunn,
A. Melatos,
C. M. Espinoza,
D. Antonopoulou,
R. Dodson
Abstract:
A striking feature of the Vela pulsar (PSR J0835$-$4510) is that it undergoes sudden increases in its spin frequency, known as glitches, with a fractional amplitude on the order of $10^{-6}$ approximately every 900 days. Glitches of smaller magnitudes are also known to occur in Vela. Their distribution in both time and amplitude is less well constrained but equally important for understanding the…
▽ More
A striking feature of the Vela pulsar (PSR J0835$-$4510) is that it undergoes sudden increases in its spin frequency, known as glitches, with a fractional amplitude on the order of $10^{-6}$ approximately every 900 days. Glitches of smaller magnitudes are also known to occur in Vela. Their distribution in both time and amplitude is less well constrained but equally important for understanding the physical process underpinning these events. In order to better understand these small glitches in Vela, an analysis of high-cadence observations from the Mount Pleasant Observatory is presented. A hidden Markov model (HMM) is used to search for small, previously undetected glitches across 24 years of observations covering MJD 44929 to MJD 53647. One previously unknown glitch is detected around MJD 48636 (Jan 15 1992), with fractional frequency jump $Δf/f = (8.19 \pm 0.04) \times 10^{-10}$ and frequency derivative jump $Δ\dot{f}/\dot{f} = (2.98 \pm 0.01) \times 10^{-4}$. Two previously reported small glitches are also confidently re-detected, and independent estimates of their parameters are reported. Excluding these events, 90% confidence frequentist upper limits on the sizes of missed glitches are also set, with a median upper limit of $Δf^{90\%}/f = 1.35 \times 10^{-9}$. Upper limits of this kind are enabled by the semi-automated and computationally efficient nature of the HMM, and are crucial to informing studies which are sensitive to the lower end of the glitch size distribution.
△ Less
Submitted 26 April, 2023;
originally announced April 2023.
-
Existence of maximal solutions for the financial stochastic Stefan problem of a volatile asset with spread
Authors:
D. C. Antonopoulou,
D. Farazakis,
G. Karali
Abstract:
In this work, we consider the outer Stefan problem for the short-time prediction of the spread of a volatile asset traded in a financial market. The stochastic equation for the evolution of the density of sell and buy orders is the Heat Equation with a non-smooth noise in the sense of Walsh, posed in a moving boundary domain with velocity given by the Stefan condition. This condition determines th…
▽ More
In this work, we consider the outer Stefan problem for the short-time prediction of the spread of a volatile asset traded in a financial market. The stochastic equation for the evolution of the density of sell and buy orders is the Heat Equation with a non-smooth noise in the sense of Walsh, posed in a moving boundary domain with velocity given by the Stefan condition. This condition determines the dynamics of the spread, and the solid phase $[s^-(t),s^+(t)]$ defines the bid-ask spread area wherein the transactions vanish. We introduce a reflection measure and prove existence and uniqueness of maximal solutions up to stopping times in which the spread $s^+(t)-s^-(t)$ stays a.s. non-negative and bounded. For this, we use a Picard approximation scheme and some of the estimates of \cite{BH} for the Green's function and the associated to the reflection measure obstacle problem. Analogous results are obtained for the equation without reflection corresponding to a signed density. Additionally, we apply some formal asymptotics when the noise depends only on time to derive that the spread is given by the integral of the solution of a linear diffusion stochastic equation.
△ Less
Submitted 20 February, 2023;
originally announced February 2023.
-
Timing six energetic rotation-powered X-ray pulsars, including the fast-spinning young PSR J0058-7218 and Big Glitcher PSR J0537-6910
Authors:
Wynn C. G. Ho,
Lucien Kuiper,
Cristobal M. Espinoza,
Sebastien Guillot,
Paul S. Ray,
D. A. Smith,
Slavko Bogdanov,
Danai Antonopoulou,
Zaven Arzoumanian,
Michal Bejger,
Teruaki Enoto,
Paolo Esposito,
Alice K. Harding,
Brynmor Haskell,
Natalia Lewandowska,
Chandreyee Maitra,
Georgios Vasilopoulos
Abstract:
Measuring a pulsar's rotational evolution is crucial to understanding the nature of the pulsar. Here we provide updated timing models for the rotational evolution of six pulsars, five of which are rotation phase-connected using primarily NICER X-ray data. For the newly-discovered fast energetic young pulsar, PSR J0058-7218, we increase the baseline of its timing model from 1.4 days to 8 months and…
▽ More
Measuring a pulsar's rotational evolution is crucial to understanding the nature of the pulsar. Here we provide updated timing models for the rotational evolution of six pulsars, five of which are rotation phase-connected using primarily NICER X-ray data. For the newly-discovered fast energetic young pulsar, PSR J0058-7218, we increase the baseline of its timing model from 1.4 days to 8 months and not only measure more precisely its spin-down rate nudot = (-6.2324+/-0.0001)x10^-11 Hz s^-1 but also for the first time the second time derivative of spin rate nuddot = (4.2+/-0.2)x10^-21 Hz s^-2. For the fastest and most energetic young pulsar, PSR J0537-6910 (with 16 ms spin period), we detect 4 more glitches, for a total of 15 glitches over 4.5 years of NICER monitoring, and show that its spin-down behavior continues to set this pulsar apart from all others, including a long-term braking index n = -1.234+/-0.009 and interglitch braking indices that asymptote to <~ 7 for long times after a glitch. For PSR J1101-6101, we measure a much more accurate spin-down rate that agrees with a previous value measured without phase-connection. For PSR J1412+7922 (also known as Calvera), we extend the baseline of its timing model from our previous 1-year model to 4.4 years, and for PSR J1849-0001, we extend the baseline from 1.5 years to 4.7 years. We also present a long-term timing model of the energetic pulsar, PSR J1813-1749, by fitting previous radio and X-ray spin frequencies from 2009-2019 and new ones measured here using 2018 NuSTAR and 2021 Chandra data.
△ Less
Submitted 29 July, 2022; v1 submitted 5 May, 2022;
originally announced May 2022.
-
Narrowband searches for continuous and long-duration transient gravitational waves from known pulsars in the LIGO-Virgo third observing run
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
the KAGRA Collaboration,
R. Abbott,
T. D. Abbott,
F. Acernese,
K. Ackley,
C. Adams,
N. Adhikari,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
D. Agarwal,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
L. Aiello,
A. Ain,
P. Ajith,
T. Akutsu,
S. Albanesi,
A. Allocca,
P. A. Altin,
A. Amato
, et al. (1636 additional authors not shown)
Abstract:
Isolated neutron stars that are asymmetric with respect to their spin axis are possible sources of detectable continuous gravitational waves. This paper presents a fully-coherent search for such signals from eighteen pulsars in data from LIGO and Virgo's third observing run (O3). For known pulsars, efficient and sensitive matched-filter searches can be carried out if one assumes the gravitational…
▽ More
Isolated neutron stars that are asymmetric with respect to their spin axis are possible sources of detectable continuous gravitational waves. This paper presents a fully-coherent search for such signals from eighteen pulsars in data from LIGO and Virgo's third observing run (O3). For known pulsars, efficient and sensitive matched-filter searches can be carried out if one assumes the gravitational radiation is phase-locked to the electromagnetic emission. In the search presented here, we relax this assumption and allow the frequency and frequency time-derivative of the gravitational waves to vary in a small range around those inferred from electromagnetic observations. We find no evidence for continuous gravitational waves, and set upper limits on the strain amplitude for each target. These limits are more constraining for seven of the targets than the spin-down limit defined by ascribing all rotational energy loss to gravitational radiation. In an additional search we look in O3 data for long-duration (hours-months) transient gravitational waves in the aftermath of pulsar glitches for six targets with a total of nine glitches. We report two marginal outliers from this search, but find no clear evidence for such emission either. The resulting duration-dependent strain upper limits do not surpass indirect energy constraints for any of these targets.
△ Less
Submitted 27 June, 2022; v1 submitted 21 December, 2021;
originally announced December 2021.
-
Searches for Gravitational Waves from Known Pulsars at Two Harmonics in the Second and Third LIGO-Virgo Observing Runs
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
the KAGRA Collaboration,
R. Abbott,
H. Abe,
F. Acernese,
K. Ackley,
N. Adhikari,
R. X. Adhikari,
V. K. Adkins,
V. B. Adya,
C. Affeldt,
D. Agarwal,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
L. Aiello,
A. Ain,
P. Ajith,
T. Akutsu,
S. Albanesi,
R. A. Alfaidi,
A. Allocca,
P. A. Altin
, et al. (1672 additional authors not shown)
Abstract:
We present a targeted search for continuous gravitational waves (GWs) from 236 pulsars using data from the third observing run of LIGO and Virgo (O3) combined with data from the second observing run (O2). Searches were for emission from the $l=m=2$ mass quadrupole mode with a frequency at only twice the pulsar rotation frequency (single harmonic) and the $l=2, m=1,2$ modes with a frequency of both…
▽ More
We present a targeted search for continuous gravitational waves (GWs) from 236 pulsars using data from the third observing run of LIGO and Virgo (O3) combined with data from the second observing run (O2). Searches were for emission from the $l=m=2$ mass quadrupole mode with a frequency at only twice the pulsar rotation frequency (single harmonic) and the $l=2, m=1,2$ modes with a frequency of both once and twice the rotation frequency (dual harmonic). No evidence of GWs was found so we present 95\% credible upper limits on the strain amplitudes $h_0$ for the single harmonic search along with limits on the pulsars' mass quadrupole moments $Q_{22}$ and ellipticities $\varepsilon$. Of the pulsars studied, 23 have strain amplitudes that are lower than the limits calculated from their electromagnetically measured spin-down rates. These pulsars include the millisecond pulsars J0437\textminus4715 and J0711\textminus6830 which have spin-down ratios of 0.87 and 0.57 respectively. For nine pulsars, their spin-down limits have been surpassed for the first time. For the Crab and Vela pulsars our limits are factors of $\sim 100$ and $\sim 20$ more constraining than their spin-down limits, respectively. For the dual harmonic searches, new limits are placed on the strain amplitudes $C_{21}$ and $C_{22}$. For 23 pulsars we also present limits on the emission amplitude assuming dipole radiation as predicted by Brans-Dicke theory.
△ Less
Submitted 20 July, 2022; v1 submitted 25 November, 2021;
originally announced November 2021.
-
The Jodrell Bank Glitch Catalogue: 106 new rotational glitches in 70 pulsars
Authors:
Avishek Basu,
Benjamin Shaw,
Danai Antonopoulou,
Michael J. Keith,
Andrew G. Lyne,
Mitchell B. Mickaliger,
Benjamin W. Stappers,
Patrick Weltevrede,
Christine A. Jordan
Abstract:
Pulsar glitches are rapid spin-up events that occur in the rotation of neutron stars, providing a valuable probe into the physics of the interiors of these objects. Long-term monitoring of a large number of pulsars facilitates the detection of glitches and the robust measurements of their parameters. The Jodrell Bank pulsar timing programme regularly monitors more than 800 radio pulsars and has ac…
▽ More
Pulsar glitches are rapid spin-up events that occur in the rotation of neutron stars, providing a valuable probe into the physics of the interiors of these objects. Long-term monitoring of a large number of pulsars facilitates the detection of glitches and the robust measurements of their parameters. The Jodrell Bank pulsar timing programme regularly monitors more than 800 radio pulsars and has accrued, in some cases, over 50 years of timing history on individual objects. In this paper we present 106 new glitches in 70 radio pulsars as observed up to the end of 2018. For 70% of these pulsars, the event we report is its only known glitch. For each new glitch we provide measurements of its epoch, amplitude and any detected changes to the spin-down rate of the star. Combining these new glitches with those listed in the Jodrell Bank glitch catalogue we analyse a total sample of 543 glitches in 178 pulsars. We model the distribution of glitch amplitudes and spin-down rate changes using a mixture of two Gaussian components. We corroborate the known dependence of glitch rate and activity on pulsar spin-down rates and characteristic ages, and show that younger pulsars tend to exhibit larger glitches. Pulsars whose spin-down rates between $10^{-14}$ Hz s$^{-1}$ and $10^{-10.5}$ Hz s$^{-1}$ show a mean reversal of 1.8% of their spin-down as a consequence of glitches. Our results are qualitatively consistent with the superfluid vortex unpinning models of pulsar glitches.
△ Less
Submitted 12 November, 2021;
originally announced November 2021.
-
Layer dynamics for the one dimensional $\varepsilon$-dependent Cahn-Hilliard / Allen-Cahn equation
Authors:
Dimitra Antonopoulou,
Georgia Karali,
Konstantinos Tzirakis
Abstract:
We study the dynamics of the one-dimensional $\varepsilon$-dependent Cahn-Hilliard / Allen-Cahn equation within a neighborhood of an equilibrium of $N$ transition layers, that in general does not conserve mass. Two different settings are considered which differ in that, for the second, we impose a mass-conservation constraint in place of one of the zero-mass flux boundary conditions at $x=1$. Moti…
▽ More
We study the dynamics of the one-dimensional $\varepsilon$-dependent Cahn-Hilliard / Allen-Cahn equation within a neighborhood of an equilibrium of $N$ transition layers, that in general does not conserve mass. Two different settings are considered which differ in that, for the second, we impose a mass-conservation constraint in place of one of the zero-mass flux boundary conditions at $x=1$. Motivated by the study of Carr and Pego on the layered metastable patterns of Allen-Cahn in [10], and by this of Bates and Xun in [5] for the Cahn-Hilliard equation, we implement an $N$-dimensional, and a mass-conservative $ N-1$-dimensional manifold respectively; therein, a metastable state with $N$ transition layers is approximated. We then determine, for both cases, the essential dynamics of the layers (ode systems with the equations of motion), expressed in terms of local coordinates relative to the manifold used. In particular, we estimate the spectrum of the linearized Cahn-Hilliard / Allen-Cahn operator, and specify wide families of $\varepsilon$-dependent weights $δ(\varepsilon)$, $μ(\varepsilon)$, acting at each part of the operator, for which the dynamics are stable and rest exponentially small in $\varepsilon$. Our analysis enlightens the role of mass conservation in the classification of the general mixed problem into two main categories where the solution has a profile close to Allen-Cahn, or, when the mass is conserved, close to the Cahn-Hilliard solution.
△ Less
Submitted 23 May, 2021;
originally announced May 2021.
-
Constraints from LIGO O3 data on gravitational-wave emission due to r-modes in the glitching pulsar PSR J0537-6910
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
the KAGRA Collaboration,
R. Abbott,
T. D. Abbott,
S. Abraham,
F. Acernese,
K. Ackley,
A. Adams,
C. Adams,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
D. Agarwal,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
L. Aiello,
A. Ain,
P. Ajith,
T. Akutsu,
K. M. Aleman,
G. Allen,
A. Allocca
, et al. (1574 additional authors not shown)
Abstract:
We present a search for continuous gravitational-wave emission due to r-modes in the pulsar PSR J0537-6910 using data from the LIGO-Virgo Collaboration observing run O3. PSR J0537-6910 is a young energetic X-ray pulsar and is the most frequent glitcher known. The inter-glitch braking index of the pulsar suggests that gravitational-wave emission due to r-mode oscillations may play an important role…
▽ More
We present a search for continuous gravitational-wave emission due to r-modes in the pulsar PSR J0537-6910 using data from the LIGO-Virgo Collaboration observing run O3. PSR J0537-6910 is a young energetic X-ray pulsar and is the most frequent glitcher known. The inter-glitch braking index of the pulsar suggests that gravitational-wave emission due to r-mode oscillations may play an important role in the spin evolution of this pulsar. Theoretical models confirm this possibility and predict emission at a level that can be probed by ground-based detectors. In order to explore this scenario, we search for r-mode emission in the epochs between glitches by using a contemporaneous timing ephemeris obtained from NICER data. We do not detect any signals in the theoretically expected band of 86-97 Hz, and report upper limits on the amplitude of the gravitational waves. Our results improve on previous amplitude upper limits from r-modes in J0537-6910 by a factor of up to 3 and place stringent constraints on theoretical models for r-mode driven spin-down in PSR J0537-6910, especially for higher frequencies at which our results reach below the spin-down limit defined by energy conservation.
△ Less
Submitted 7 January, 2022; v1 submitted 29 April, 2021;
originally announced April 2021.
-
The multi-dimensional Stochastic Stefan Financial Model for a portfolio of assets
Authors:
Dimitra C. Antonopoulou,
Marina Bitsaki,
Georgia Karali
Abstract:
The financial model proposed involves the liquidation process of a portfolio of $n$ assets through sell or (and) buy orders with volatility. We present the rigorous mathematical formulation of this model in a financial setting resulting to an $n$-dimensional outer parabolic Stefan problem with noise. In particular, our aim is to estimate for a short time period the areas of zero trading, and their…
▽ More
The financial model proposed involves the liquidation process of a portfolio of $n$ assets through sell or (and) buy orders with volatility. We present the rigorous mathematical formulation of this model in a financial setting resulting to an $n$-dimensional outer parabolic Stefan problem with noise. In particular, our aim is to estimate for a short time period the areas of zero trading, and their diameter which approximates the minimum of the $n$ spreads of the portfolio assets for orders from the $n$ limit order books of each asset respectively. In dimensions $n=3$, and for zero volatility, this problem stands as a mean field model for Ostwald ripening, and has been proposed and analyzed by Niethammer. Therein, when the initial moving boundary consists of well separated spheres, a first order approximation system of odes had been rigorously derived for the dynamics of the interfaces and the asymptotic profile of the solution. In our financial case, we propose a spherical moving boundaries approach where the zero trading area consists of a union of spherical domains centered at portfolios various prices, while each sphere may correspond to a different market; the relevant radii represent the half of the minimum spread. We apply Itô calculus and provide second order formal asymptotics for the stochastic version dynamics, written as a system of stochastic differential equations for the radii evolution in time. A second order approximation seems to disconnect the financial model from the large diffusion assumption for the trading density. Moreover, we solve the approximating systems numerically.
△ Less
Submitted 24 December, 2020;
originally announced December 2020.
-
Diving below the spin-down limit: Constraints on gravitational waves from the energetic young pulsar PSR J0537-6910
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
the KAGRA Collaboration,
R. Abbott,
T. D. Abbott,
S. Abraham,
F. Acernese,
K. Ackley,
A. Adams,
C. Adams,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
D. Agarwal,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
L. Aiello,
A. Ain,
P. Ajith,
T. Akutsu,
K. M. Aleman,
G. Allen,
A. Allocca
, et al. (1568 additional authors not shown)
Abstract:
We present a search for continuous gravitational-wave signals from the young, energetic X-ray pulsar PSR J0537-6910 using data from the second and third observing runs of LIGO and Virgo. The search is enabled by a contemporaneous timing ephemeris obtained using NICER data. The NICER ephemeris has also been extended through 2020 October and includes three new glitches. PSR J0537-6910 has the larges…
▽ More
We present a search for continuous gravitational-wave signals from the young, energetic X-ray pulsar PSR J0537-6910 using data from the second and third observing runs of LIGO and Virgo. The search is enabled by a contemporaneous timing ephemeris obtained using NICER data. The NICER ephemeris has also been extended through 2020 October and includes three new glitches. PSR J0537-6910 has the largest spin-down luminosity of any pulsar and is highly active with regards to glitches. Analyses of its long-term and inter-glitch braking indices provided intriguing evidence that its spin-down energy budget may include gravitational-wave emission from a time-varying mass quadrupole moment. Its 62 Hz rotation frequency also puts its possible gravitational-wave emission in the most sensitive band of LIGO/Virgo detectors. Motivated by these considerations, we search for gravitational-wave emission at both once and twice the rotation frequency. We find no signal, however, and report our upper limits. Assuming a rigidly rotating triaxial star, our constraints reach below the gravitational-wave spin-down limit for this star for the first time by more than a factor of two and limit gravitational waves from the $l=m=2$ mode to account for less than 14% of the spin-down energy budget. The fiducial equatorial ellipticity is limited to less than about 3e-5, which is the third best constraint for any young pulsar.
△ Less
Submitted 10 June, 2021; v1 submitted 23 December, 2020;
originally announced December 2020.
-
Return of the Big Glitcher: NICER timing and glitches of PSR J0537-6910
Authors:
Wynn C. G. Ho,
Cristobal M. Espinoza,
Zaven Arzoumanian,
Teruaki Enoto,
Tsubasa Tamba,
Danai Antonopoulou,
Michal Bejger,
Sebastien Guillot,
Brynmor Haskell,
Paul S. Ray
Abstract:
PSR J0537-6910, also known as the Big Glitcher, is the most prolific glitching pulsar known, and its spin-induced pulsations are only detectable in X-ray. We present results from analysis of 2.7 years of NICER timing observations, from 2017 August to 2020 April. We obtain a rotation phase-connected timing model for the entire timespan, which overlaps with the third observing run of LIGO/Virgo, thu…
▽ More
PSR J0537-6910, also known as the Big Glitcher, is the most prolific glitching pulsar known, and its spin-induced pulsations are only detectable in X-ray. We present results from analysis of 2.7 years of NICER timing observations, from 2017 August to 2020 April. We obtain a rotation phase-connected timing model for the entire timespan, which overlaps with the third observing run of LIGO/Virgo, thus enabling the most sensitive gravitational wave searches of this potentially strong gravitational wave-emitting pulsar. We find that the short-term braking index between glitches decreases towards a value of 7 or lower at longer times since the preceding glitch. By combining NICER and RXTE data, we measure a long-term braking index n=-1.25+/-0.01. Our analysis reveals 8 new glitches, the first detected since 2011, near the end of RXTE, with a total NICER and RXTE glitch activity of 8.88x10^-7 yr^-1. The new glitches follow the seemingly unique time-to-next-glitch---glitch-size correlation established previously using RXTE data, with a slope of 5 d microHz^-1. For one glitch around which NICER observes two days on either side, we search for but do not see clear evidence of spectral nor pulse profile changes that may be associated with the glitch.
△ Less
Submitted 31 August, 2020;
originally announced September 2020.
-
Small glitches and other rotational irregularities of the Vela pulsar
Authors:
C. M. Espinoza,
D. Antonopoulou,
R. Dodson,
M. Stepanova,
A. Scherer
Abstract:
Glitches are sudden increases in the rotation rate $ν$ of neutron stars, which are thought to be driven by the neutron superfluid inside the star. The Vela pulsar presents a comparatively high rate of glitches, with 21 events reported since observations began in 1968. These are amongst the largest known glitches (17 of them have sizes $Δν/ν\geq10^{-6}$) and exhibit very similar characteristics. Th…
▽ More
Glitches are sudden increases in the rotation rate $ν$ of neutron stars, which are thought to be driven by the neutron superfluid inside the star. The Vela pulsar presents a comparatively high rate of glitches, with 21 events reported since observations began in 1968. These are amongst the largest known glitches (17 of them have sizes $Δν/ν\geq10^{-6}$) and exhibit very similar characteristics. This similarity, combined with the regularity with which large glitches occur, has turned Vela into an archetype of this type of glitching behaviour. The properties of its smallest glitches, on the other hand, are not clearly established. High-cadence observations of the Vela pulsar were taken between 1981 and 2005 at the Mount Pleasant Radio Observatory. An automated systematic search was carried out that investigated whether a significant change of spin frequency $ν$ and/or the spin-down rate $\dotν$ takes place at any given time. We find two new glitches, with respective sizes $Δν/ν$ of $(5.55\pm0.03)\times10^{-9}$ and $(38\pm4)\times10^{-9}$. In addition to these two glitch events, our study reveals numerous events of all possible signatures (i.e. combinations of $Δν$ and $Δ\dotν$ signs), all of them small with $|Δν|/ν<10^{-9}$, which contribute to the Vela timing noise. The Vela pulsar presents an under-abundance of small glitches compared to many other glitching pulsars, which appears genuine and not a result of observational biases. In addition to typical glitches, the smooth spin-down of the pulsar is also affected by an almost continuous activity that can be partially characterised by small step-like changes in $ν$, $\dot{ν,}$ or both. Simulations indicate that a continuous wandering of the rotational phase, following a red spectrum, could mimic such step-like changes in the timing residuals.
△ Less
Submitted 12 April, 2021; v1 submitted 6 July, 2020;
originally announced July 2020.
-
Turbulent, pinned superfluids in neutron stars and pulsar glitch recoveries
Authors:
Brynmor Haskell,
Danai Antonopoulou,
Carlo Barenghi
Abstract:
Pulsar glitches offer an insight into the dynamics of superfluids in the high density interior of a neutron star. To model these phenomena, however, one needs to have an understanding of the dynamics of a turbulent array of superfluid vortices moving through a pinning lattice. In this paper we develop a theoretical approach to describe vortex mediated mutual friction in a pinned, turbulent and rot…
▽ More
Pulsar glitches offer an insight into the dynamics of superfluids in the high density interior of a neutron star. To model these phenomena, however, one needs to have an understanding of the dynamics of a turbulent array of superfluid vortices moving through a pinning lattice. In this paper we develop a theoretical approach to describe vortex mediated mutual friction in a pinned, turbulent and rotating superfluid. Our model is then applied to the study of the post glitch rotational evolution in the Vela pulsar and in PSR J0537-6910. We show that in both cases a turbulent model fits the evolution of the spin frequency derivative better than a laminar one. We also predict that the second derivative of the frequency after a glitch should be correlated with the waiting time since the previous glitch, which we find to be consistent with observational data for these pulsars. The main conclusion of this paper is that in the post-glitch rotational evolution of these two pulsars we are most likely observing the response to the glitch of a pinned turbulent region of the star (possibly the crust) and not the laminar response of a regular straight vortex array.
△ Less
Submitted 6 July, 2020;
originally announced July 2020.
-
Numerical approximation of the Stochastic Cahn-Hilliard Equation near the Sharp Interface Limit
Authors:
Dimitra Antonopoulou,
Lubomir Banas,
Robert Nürnberg,
Andreas Prohl
Abstract:
We consider the stochastic Cahn-Hilliard equation with additive noise term $\varepsilon^γg\, \dot{W}$ ($γ>0$) that scales with the interfacial width parameter $\varepsilon$. We verify strong error estimates for a gradient flow structure-inheriting time-implicit discretization, where $\varepsilon^{-1}$ only enters polynomially; the proof is based on higher-moment estimates for iterates, and a (disc…
▽ More
We consider the stochastic Cahn-Hilliard equation with additive noise term $\varepsilon^γg\, \dot{W}$ ($γ>0$) that scales with the interfacial width parameter $\varepsilon$. We verify strong error estimates for a gradient flow structure-inheriting time-implicit discretization, where $\varepsilon^{-1}$ only enters polynomially; the proof is based on higher-moment estimates for iterates, and a (discrete) spectral estimate for its deterministic counterpart. For $γ$ sufficiently large, convergence in probability of iterates towards the deterministic Hele-Shaw/Mullins-Sekerka problem in the sharp-interface limit $\varepsilon \rightarrow 0$ is shown. These convergence results are partly generalized to a fully discrete finite element based discretization.
We complement the theoretical results by computational studies to provide practical evidence concerning the effect of noise (depending on its 'strength' $γ$) on the geometric evolution in the sharp-interface limit. For this purpose we compare the simulations with those from a fully discrete finite element numerical scheme for the (stochastic) Mullins-Sekerka problem. The computational results indicate that the limit for $γ\geq 1$ is the deterministic problem, and for $γ=0$ we obtain agreement with a (new) stochastic version of the Mullins-Sekerka problem.
△ Less
Submitted 7 July, 2020; v1 submitted 27 May, 2019;
originally announced May 2019.
-
Generation of fine transition layers and their dynamics for the stochastic Allen--Cahn equation
Authors:
Matthieu Alfaro,
Dimitra Antonopoulou,
Georgia Karali,
Hiroshi Matano
Abstract:
We study an $\ep$-dependent stochastic Allen--Cahn equation with a mild random noise on a bounded domain in $\mathbb{R}^n$, $n\geq 2$. Here $\ep$ is a small positive parameter that represents formally the thickness of the solution interface, while the mild noise $ξ^\ep(t)$ is a smooth random function of $t$ of order $\mathcal O(\ep^{-γ})$ with $0<γ<1/3$ that converges to white noise as…
▽ More
We study an $\ep$-dependent stochastic Allen--Cahn equation with a mild random noise on a bounded domain in $\mathbb{R}^n$, $n\geq 2$. Here $\ep$ is a small positive parameter that represents formally the thickness of the solution interface, while the mild noise $ξ^\ep(t)$ is a smooth random function of $t$ of order $\mathcal O(\ep^{-γ})$ with $0<γ<1/3$ that converges to white noise as $\ep\rightarrow 0^+$. We consider initial data that are independent of $\ep$ satisfying some non-degeneracy conditions, and prove that steep transition layers---or interfaces---develop within a very short time of order $\ep^2|\ln\ep|$, which we call the "generation of interface". Next we study the motion of those transition layers and derive a stochastic motion law for the sharp interface limit as $\ep\rightarrow 0^+$. Furthermore, we prove that the thickness of the interface for $\ep$ small is indeed of order $\mathcal O(\ep)$ and that the solution profile near the interface remains close to that of a (squeezed) travelling wave, this means that the presence of the noise does not destroy the solution profile near the interface as long as the noise is spatially uniform. Our results on the motion of interface improve the earlier results of Funaki (1999) and Weber (2010) by considerably weakening the requirements for the initial data and establishing the robustness of the solution profile near the interface that has not been known before.
△ Less
Submitted 19 December, 2018; v1 submitted 10 December, 2018;
originally announced December 2018.
-
Core and crust contributions in pulsar glitches: constraints from the slow rise of the largest glitch observed in the Crab pulsar
Authors:
Brynmor Haskell,
Vadym Khomenko,
Marco Antonelli,
Danai Antonopoulou
Abstract:
Pulsar glitches are attributed to the sudden re-coupling of very weakly coupled large scale superfluid components in the neutron star interior. This process leads to rapid exchange of angular momentum and an increase in spin frequency. The transfer of angular momentum is regulated by a dissipative mutual friction, whose strength defines the spin-up timescale of a glitch. Hence, observations of gli…
▽ More
Pulsar glitches are attributed to the sudden re-coupling of very weakly coupled large scale superfluid components in the neutron star interior. This process leads to rapid exchange of angular momentum and an increase in spin frequency. The transfer of angular momentum is regulated by a dissipative mutual friction, whose strength defines the spin-up timescale of a glitch. Hence, observations of glitch rises can be used to shed light on the dominant microphysical interactions at work in the high density interior of the star. We present a simple analytical model, complemented with more detailed numerical simulations, which produces a fast spin-up followed by a more gradual rise. Such features are observed in some large glitches of the Crab pulsar, including the largest recent glitch of 2017. We also use observations to constrain the mutual friction coefficient of the glitch-driving region for two possible locations: the inner crust and outer core of the star. We find that the features of Crab glitches require smaller values of the mutual friction coefficient than those needed to explain the much faster Vela spin-ups. This suggests a crustal origin for the former but an outer core contribution for the latter.
△ Less
Submitted 26 June, 2018;
originally announced June 2018.
-
Malliavin calculus for the stochastic Cahn-Hilliard / Allen Cahn equation with unbounded noise diffusion
Authors:
D. C. Antonopoulou,
D. Farazakis,
G. D. Karali
Abstract:
The stochastic partial differential equation analyzed in this work, is motivated by a simplified mesoscopic physical model for phase separation. It describes pattern formation due to adsorption and desorption mechanisms involved in surface processes, in the presence of a stochastic driving force. This equation is a combination of Cahn-Hilliard and Allen-Cahn type operators with a multiplicative, w…
▽ More
The stochastic partial differential equation analyzed in this work, is motivated by a simplified mesoscopic physical model for phase separation. It describes pattern formation due to adsorption and desorption mechanisms involved in surface processes, in the presence of a stochastic driving force. This equation is a combination of Cahn-Hilliard and Allen-Cahn type operators with a multiplicative, white, space-time noise of unbounded diffusion. We apply Malliavin calculus, in order to investigate the existence of a density for the stochastic solution $u$. In dimension one, according to the regularity result in \cite{AKM}, $u$ admits continuous paths a.s. Using this property, and inspired by a method proposed in \cite{CW1}, we construct a modified approximating sequence for $u$, which properly treats the new second order Allen-Cahn operator. Under a localization argument, we prove that the Malliavin derivative of $u$ exists locally, and that the law of $u$ is absolutely continuous, establishing thus that a density exists.
△ Less
Submitted 18 February, 2018;
originally announced February 2018.
-
The enigmatic spin evolution of PSR J0537-6910: r-modes, gravitational waves and the case for continued timing
Authors:
N. Andersson,
D. Antonopoulou,
C. M. Espinoza,
B. Haskell,
W. C. G. Ho
Abstract:
We discuss the unique spin evolution of the young X-ray pulsar PSR J0537-6910, a system in which the regular spin down is interrupted by glitches every few months. Drawing on the complete timing data from the Rossi X-ray Timing Explorer (RXTE, from 1999-2011), we argue that a trend in the inter-glitch behaviour points to an effective braking index close to $n=7$, much larger than expected. This va…
▽ More
We discuss the unique spin evolution of the young X-ray pulsar PSR J0537-6910, a system in which the regular spin down is interrupted by glitches every few months. Drawing on the complete timing data from the Rossi X-ray Timing Explorer (RXTE, from 1999-2011), we argue that a trend in the inter-glitch behaviour points to an effective braking index close to $n=7$, much larger than expected. This value is interesting because it would accord with the neutron star spinning down due to gravitational waves from an unstable r-mode. We discuss to what extent this, admittedly speculative, scenario may be consistent and if the associated gravitational-wave signal would be within reach of ground based detectors. Our estimates suggest that one may, indeed, be able to use future observations to test the idea. Further precision timing would help enhance the achievable sensitivity and we advocate a joint observing campaign between the Neutron Star Interior Composition ExploreR (NICER) and the LIGO-Virgo network.
△ Less
Submitted 15 November, 2017;
originally announced November 2017.
-
Pulsar spin-down: the glitch-dominated rotation of PSR J0537-6910
Authors:
D. Antonopoulou,
C. M. Espinoza,
L. Kuiper,
N. Andersson
Abstract:
The young, fast-spinning, X-ray pulsar J0537-6910 displays an extreme glitch activity, with large spin-ups interrupting its decelerating rotation every ~100 days. We present nearly 13 years of timing data from this pulsar, obtained with the {\it Rossi X-ray Timing Explorer}. We discovered 22 new glitches and performed a consistent analysis of all 45 glitches detected in the complete data span. Our…
▽ More
The young, fast-spinning, X-ray pulsar J0537-6910 displays an extreme glitch activity, with large spin-ups interrupting its decelerating rotation every ~100 days. We present nearly 13 years of timing data from this pulsar, obtained with the {\it Rossi X-ray Timing Explorer}. We discovered 22 new glitches and performed a consistent analysis of all 45 glitches detected in the complete data span. Our results corroborate the previously reported strong correlation between glitch spin-up size and the time to the next glitch, a relation that has not been observed so far in any other pulsar. The spin evolution is dominated by the glitches, which occur at a rate ~3.5 per year, and the post-glitch recoveries, which prevail the entire inter-glitch intervals. This distinctive behaviour provides invaluable insights into the physics of glitches. The observations can be explained with a multi-component model which accounts for the dynamics of the neutron superfluid present in the crust and core of neutron stars. We place limits on the moment of inertia of the component responsible for the spin-up and, ignoring differential rotation, the velocity difference it can sustain with the crust. Contrary to its rapid decrease between glitches, the spin-down rate increased over the 13 years, and we find the long-term braking index $n_{\rm l}=-1.22(4)$, the only negative braking index seen in a young pulsar. We briefly discuss the plausible interpretations of this result, which is in stark contrast to the predictions of standard models of pulsar spin-down.
△ Less
Submitted 18 September, 2017; v1 submitted 30 August, 2017;
originally announced August 2017.
-
Pinning down the superfluid and nuclear equation of state and measuring neutron star mass using pulsar glitches
Authors:
Wynn C. G. Ho,
Cristobal M. Espinoza,
Danai Antonopoulou,
Nils Andersson
Abstract:
Pulsars are rotating neutron stars that are renowned for their timing precision, although glitches can interrupt the regular timing behavior when these stars are young. Glitches are thought to be caused by interactions between normal and superfluid matter in the star. We update our recent work on a new technique using pulsar glitch data to constrain superfluid and nuclear equation of state models,…
▽ More
Pulsars are rotating neutron stars that are renowned for their timing precision, although glitches can interrupt the regular timing behavior when these stars are young. Glitches are thought to be caused by interactions between normal and superfluid matter in the star. We update our recent work on a new technique using pulsar glitch data to constrain superfluid and nuclear equation of state models, demonstrating how current and future astronomy telescopes can probe fundamental physics such as superfluidity near nuclear saturation and matter at supranuclear densities. Unlike traditional methods of measuring a star's mass by its gravitational effect on another object, our technique relies on nuclear physics knowledge and therefore allows measurement of the mass of pulsars which are in isolation.
△ Less
Submitted 2 March, 2017;
originally announced March 2017.
-
On the Dirichlet to Neumann Problem for the 1-dimensional Cubic NLS equation on the Half-Line; Zero Initial Data
Authors:
Dimitra C. Antonopoulou,
Spyridon Kamvissis
Abstract:
Initial-boundary value problems for 1-dimensional `completely integrable' equations can be solved via an extension of the inverse scattering method, which is due to Fokas and his collaborators. A crucial feature of this method is that it requires the values of more boundary data than given for a well-posed problem. In the case of cubic NLS, knowledge of the Dirichet data suffices to make the probl…
▽ More
Initial-boundary value problems for 1-dimensional `completely integrable' equations can be solved via an extension of the inverse scattering method, which is due to Fokas and his collaborators. A crucial feature of this method is that it requires the values of more boundary data than given for a well-posed problem. In the case of cubic NLS, knowledge of the Dirichet data suffices to make the problem well-posed but the Fokas method also requires knowledge of the values of Neumann data. The study of the Dirichlet to Neumann map is thus necessary before the application of the `Fokas transform'. In this paper, we provide a rigorous study of this map for a large class of decaying Dirichlet data. We show that the Neumann data are also sufficiently decaying and hence that the Fokas method can be applied.
For simplicity we considered here the case of zero initial data. An addendum will follow, discussing the case of non-zero initial data.
△ Less
Submitted 1 September, 2016; v1 submitted 21 July, 2016;
originally announced July 2016.
-
On the Dirichlet to Neumann Problem for the 1-dimensional Cubic NLS equation on the Half-Line; Non-zero Initial Data
Authors:
Dimitra C. Antonopoulou,
Spyridon Kamvissis
Abstract:
Initial-boundary value problems for 1-dimensional `completely integrable' equations can be solved via an extension of the inverse scattering method, which is due to Fokas and his collaborators. A crucial feature of this method is that it requires the values of more boundary data than given for a well-posed problem. In the case of cubic NLS, knowledge of the Dirichet data suffices to make the probl…
▽ More
Initial-boundary value problems for 1-dimensional `completely integrable' equations can be solved via an extension of the inverse scattering method, which is due to Fokas and his collaborators. A crucial feature of this method is that it requires the values of more boundary data than given for a well-posed problem. In the case of cubic NLS, knowledge of the Dirichet data suffices to make the problem well-posed but the Fokas method also requires knowledge of the values of Neumann data. The study of the Dirichlet to Neumann map is thus necessary before the application of the `Fokas transform'. In this paper, we provide a rigorous study of this map for a large class of decaying Dirichlet data. We show that the Neumann data are also sufficiently decaying and hence that the Fokas method can be applied.
This is an addendum to a previous paper. We consider the case of non-zero initial data, thus completing the discussion of the previous paper.
△ Less
Submitted 6 September, 2016; v1 submitted 21 July, 2016;
originally announced July 2016.
-
Pinning down the superfluid and measuring masses using pulsar glitches
Authors:
Wynn C. G. Ho,
Cristobal M. Espinoza,
Danai Antonopoulou,
Nils Andersson
Abstract:
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…
▽ More
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.
△ Less
Submitted 1 October, 2015;
originally announced October 2015.
-
The Sharp interface limit for the stochastic Cahn-Hilliard equation
Authors:
Dimitra C. Antonopoulou,
Dirk Blömker,
Georgia D. Karali
Abstract:
We study the two and three dimensional stochastic Cahn-Hilliard equation in the sharp interface limit, where the positive parameter $ε$ tends to zero, which measures the width of transition layers generated during phase separation. We also couple the noise strength to this parameter. Using formal asymptotic expansions, we identify the limit. In the right scaling we indicate that the solutions of s…
▽ More
We study the two and three dimensional stochastic Cahn-Hilliard equation in the sharp interface limit, where the positive parameter $ε$ tends to zero, which measures the width of transition layers generated during phase separation. We also couple the noise strength to this parameter. Using formal asymptotic expansions, we identify the limit. In the right scaling we indicate that the solutions of stochastic Cahn-Hilliard converge to a solution of a Hele-Shaw problem with stochastic forcing. In the case when the noise is sufficiently small, we rigorously prove that the limit is a deterministic Hele-Shaw problem. Finally, we discuss which estimates are necessary in order to extend the rigorous result to larger noise strength.
△ Less
Submitted 22 September, 2016; v1 submitted 27 July, 2015;
originally announced July 2015.
-
Motion of a droplet for the mass-conserving stochastic Allen-Cahn equation
Authors:
Dimitra C. Antonopoulou,
Peter W. Bates,
Dirk Blömker,
Georgia D. Karali
Abstract:
We study the stochastic mass-conserving Allen-Cahn equation posed on a bounded two-dimensional domain with additive spatially smooth space-time noise. This equation associated with a small positive parameter describes the stochastic motion of a small almost semicircular droplet attached to domain's boundary and moving towards a point of locally maximum curvature. We apply Itô calculus to derive th…
▽ More
We study the stochastic mass-conserving Allen-Cahn equation posed on a bounded two-dimensional domain with additive spatially smooth space-time noise. This equation associated with a small positive parameter describes the stochastic motion of a small almost semicircular droplet attached to domain's boundary and moving towards a point of locally maximum curvature. We apply Itô calculus to derive the stochastic dynamics of the droplet by utilizing the approximately invariant manifold introduced by Alikakos, Chen and Fusco for the deterministic problem. In the stochastic case depending on the scaling, the motion is driven by the change in the curvature of the boundary and the stochastic forcing. Moreover, under the assumption of a sufficiently small noise strength, we establish stochastic stability of a neighborhood of the manifold of droplets in $L^2$ and $H^1$, which means that with overwhelming probability the solution stays close to the manifold for very long time-scales.
△ Less
Submitted 21 January, 2015;
originally announced January 2015.
-
Probing the neutron star interior and the Equation of State of cold dense matter with the SKA
Authors:
Anna Watts,
Renxin Xu,
Cristobal Espinoza,
Nils Andersson,
John Antoniadis,
Danai Antonopoulou,
Sarah Buchner,
Shi Dai,
Paul Demorest,
Paulo Freire,
Jason Hessels,
Jerome Margueron,
Micaela Oertel,
Alessandro Patruno,
Andrea Possenti,
Scott Ransom,
Ingrid Stairs,
Ben Stappers
Abstract:
With an average density higher than the nuclear density, neutron stars (NS) provide a unique test-ground for nuclear physics, quantum chromodynamics (QCD), and nuclear superfluidity. Determination of the fundamental interactions that govern matter under such extreme conditions is one of the major unsolved problems of modern physics, and -- since it is impossible to replicate these conditions on Ea…
▽ More
With an average density higher than the nuclear density, neutron stars (NS) provide a unique test-ground for nuclear physics, quantum chromodynamics (QCD), and nuclear superfluidity. Determination of the fundamental interactions that govern matter under such extreme conditions is one of the major unsolved problems of modern physics, and -- since it is impossible to replicate these conditions on Earth -- a major scientific motivation for SKA. The most stringent observational constraints come from measurements of NS bulk properties: each model for the microscopic behaviour of matter predicts a specific density-pressure relation (its `Equation of state', EOS). This generates a unique mass-radius relation which predicts a characteristic radius for a large range of masses and a maximum mass above which NS collapse to black holes. It also uniquely predicts other bulk quantities, like maximum spin frequency and moment of inertia. The SKA, in Phase 1 and particularly in Phase 2 will, thanks to the exquisite timing precision enabled by its raw sensitivity, and surveys that dramatically increase the number of sources: 1) Provide many more precise NS mass measurements (high mass NS measurements are particularly important for ruling out EOS models); 2) Allow the measurement of the NS moment of inertia in highly relativistic binaries such as the Double Pulsar; 3) Greatly increase the number of fast-spinning NS, with the potential discovery of spin frequencies above those allowed by some EOS models; 4) Improve our knowledge of new classes of binary pulsars such as black widows and redbacks (which may be massive as a class) through sensitive broad-band radio observations; and 5) Improve our understanding of dense matter superfluidity and the state of matter in the interior through the study of rotational glitches, provided that an ad-hoc campaign is developed.
△ Less
Submitted 30 December, 2014;
originally announced January 2015.
-
The unusual glitch recoveries of the high magnetic field pulsar J1119$-$6127
Authors:
D. Antonopoulou,
P. Weltevrede,
C. M. Espinoza,
A. L. Watts,
S. Johnston,
R. M. Shannon,
M. Kerr
Abstract:
Providing a link between magnetars and radio pulsars, high magnetic field neutron stars are ideal targets to investigate how bursting/magnetospheric activity and braking torque variations are connected to rotational glitches. The last spin-up glitch of the highly magnetised pulsar J1119$-$6127 back in 2007 was the first glitch in a rotationally powered radio pulsar to be accompanied by radiative c…
▽ More
Providing a link between magnetars and radio pulsars, high magnetic field neutron stars are ideal targets to investigate how bursting/magnetospheric activity and braking torque variations are connected to rotational glitches. The last spin-up glitch of the highly magnetised pulsar J1119$-$6127 back in 2007 was the first glitch in a rotationally powered radio pulsar to be accompanied by radiative changes. Moreover, it was followed by an uncommon glitch relaxation that resulted in a smaller spin-down rate relative to the prediction of the pre-glitch timing model. Here, we present 4 years of new radio timing observations and analyse the total of 16 years of timing data for this source. The new data uncover an ongoing evolution of the spin-down rate, thereby allowing us to exclude permanent changes in the external or internal torque as a stand-alone cause of the peculiar features of the glitch recovery. Furthermore, no additional variations of the radio pulse profile are detected, strengthening the association of the previously observed transient emission features with the glitching activity. A self-consistent measurement of the braking index yields a value $n\simeq2.7$, indicating a trajectory in the $P-\dot{P}$ plane inclined towards the magnetars. Such a potential evolutionary link might be strengthened by a, possibly permanent, reduction of $\sim15\%$ in $n$ at the epoch of the 2007 glitch.
△ Less
Submitted 18 December, 2014;
originally announced December 2014.
-
Magnetically-driven crustquakes in neutron stars
Authors:
S. K. Lander,
N. Andersson,
D. Antonopoulou,
A. L. Watts
Abstract:
Crustquake events may be connected with both rapid spin-up `glitches' within the regular slowdown of neutron stars, and high-energy magnetar flares. We argue that magnetic field decay builds up stresses in a neutron star's crust, as the elastic shear force resists the Lorentz force's desire to rearrange the global magnetic-field equilibrium. We derive a criterion for crust-breaking induced by a ch…
▽ More
Crustquake events may be connected with both rapid spin-up `glitches' within the regular slowdown of neutron stars, and high-energy magnetar flares. We argue that magnetic field decay builds up stresses in a neutron star's crust, as the elastic shear force resists the Lorentz force's desire to rearrange the global magnetic-field equilibrium. We derive a criterion for crust-breaking induced by a changing magnetic-field configuration, and use this to investigate strain patterns in a neutron star's crust for a variety of different magnetic-field models. Universally, we find that the crust is most liable to break if the magnetic field has a strong toroidal component, in which case the epicentre of the crustquake is around the equator. We calculate the energy released in a crustquake as a function of the fracture depth, finding that it is independent of field strength. Crust-breaking is, however, associated with a characteristic local field strength of $2.4\times 10^{14}$ G for a breaking strain of $0.001$, or $2.4\times 10^{15}$ G at a breaking strain of $0.1$. We find that even the most luminous magnetar giant flare could have been powered by crustal energy release alone.
△ Less
Submitted 8 March, 2015; v1 submitted 18 December, 2014;
originally announced December 2014.
-
Neutron star glitches have a substantial minimum size
Authors:
C. M. Espinoza,
D. Antonopoulou,
B. W. Stappers,
A. Watts,
A. G. Lyne
Abstract:
Glitches are sudden spin-up events that punctuate the steady spin down of pulsars and are thought to be due to the presence of a superfluid component within neutron stars. The precise glitch mechanism and its trigger, however, remain unknown. The size of glitches is a key diagnostic for models of the underlying physics. While the largest glitches have long been taken into account by theoretical mo…
▽ More
Glitches are sudden spin-up events that punctuate the steady spin down of pulsars and are thought to be due to the presence of a superfluid component within neutron stars. The precise glitch mechanism and its trigger, however, remain unknown. The size of glitches is a key diagnostic for models of the underlying physics. While the largest glitches have long been taken into account by theoretical models, it has always been assumed that the minimum size lay below the detectability limit of the measurements. In this paper we define general glitch detectability limits and use them on 29 years of daily observations of the Crab pulsar, carried out at Jodrell Bank Observatory. We find that all glitches lie well above the detectability limits and by using an automated method to search for small events we are able to uncover the full glitch size distribution, with no biases. Contrary to the prediction of most models, the distribution presents a rapid decrease of the number of glitches below ~0.05 $μ$Hz. This substantial minimum size indicates that a glitch must involve the motion of at least several billion superfluid vortices and provides an extra observable which can greatly help the identification of the trigger mechanism. Our study also shows that glitches are clearly separated from all the other rotation irregularities. This supports the idea that the origin of glitches is different to that of timing noise, which comprises the unmodelled random fluctuations in the rotation rates of pulsars.
△ Less
Submitted 28 February, 2014;
originally announced February 2014.
-
Existence and regularity of solution for a Stochastic Cahn-Hilliard/Allen-Cahn equation with unbounded noise diffusion
Authors:
Dimitra C. Antonopoulou,
Geogia Karali,
Annie Millet
Abstract:
The Cahn-Hilliard/Allen-Cahn equation with noise is a simplified mean field model of stochastic microscopic dynamics associated with adsorption and desorption-spin flip mechanisms in the context of surface processes. For such an equation we consider a multiplicative space-time white noise with diffusion coefficient of sub-linear growth. Using technics from semigroup theory, we prove existence, and…
▽ More
The Cahn-Hilliard/Allen-Cahn equation with noise is a simplified mean field model of stochastic microscopic dynamics associated with adsorption and desorption-spin flip mechanisms in the context of surface processes. For such an equation we consider a multiplicative space-time white noise with diffusion coefficient of sub-linear growth. Using technics from semigroup theory, we prove existence, and path regularity of stochastic solution depending on that of the initial condition. Our results are also valid for the stochastic Cahn-Hilliard equation with unbounded noise diffusion, for which previous results were established only in the framework of a bounded diffusion coefficient. We prove that the path regularity of stochastic solution depends on that of the initial condition, and are identical to those proved for the stochastic Cahn-Hilliard equation and a bounded noise diffusion coefficient. If the initial condition vanishes, they are strictly less than 2-d/2 in space and 1/2-d/8 in time. As expected from the theory of parabolic operators in the sense of Petrovski, the bi-Laplacian operator seems to be dominant in the combined model.
△ Less
Submitted 4 October, 2013;
originally announced October 2013.
-
Glitch recoveries in radio-pulsars and magnetars
Authors:
Brynmor Haskell,
Danai Antonopoulou
Abstract:
Pulsar glitches are sudden increases in the spin frequency of an otherwise steadily spinning down neutron star. These events are thought to represent a direct probe of the dynamics of the superfluid interior of the star. However glitches can differ significantly from one another, not only in size and frequency, but also in the post-glitch response of the star. Some appear as simple steps in freque…
▽ More
Pulsar glitches are sudden increases in the spin frequency of an otherwise steadily spinning down neutron star. These events are thought to represent a direct probe of the dynamics of the superfluid interior of the star. However glitches can differ significantly from one another, not only in size and frequency, but also in the post-glitch response of the star. Some appear as simple steps in frequency, while others also display an increase in spin-down rate after the glitch. Others still show several exponentially relaxing components in the post-glitch recovery. We show that if glitches are indeed due to large scale unpinning of superfluid vortices, the different regions in which this occurs and respective timescales on which they recouple can lead to the various observed signatures. Furthermore we show that this framework naturally accounts for the peculiar relaxations of glitches in Anomalous X-ray Pulsars.
△ Less
Submitted 21 June, 2013;
originally announced June 2013.
-
A Hilbert expansions method for the rigorous sharp interface limit of the generalized Cahn-Hilliard Equation
Authors:
D. C. Antonopoulou,
G. D. Karali,
E. Orlandi
Abstract:
We consider Cahn-Hilliard equations with external forcing terms. Energy decreasing and mass conservation might not hold. We show that level surfaces of the solutions of such generalized Cahn-Hilliard equations tend to the solutions of a moving boundary problem under the assumption that classical solutions of the latter exist. Our strategy is to construct approximate solutions of the generalized Ca…
▽ More
We consider Cahn-Hilliard equations with external forcing terms. Energy decreasing and mass conservation might not hold. We show that level surfaces of the solutions of such generalized Cahn-Hilliard equations tend to the solutions of a moving boundary problem under the assumption that classical solutions of the latter exist. Our strategy is to construct approximate solutions of the generalized Cahn-Hilliard equation by the Hilbert expansion method used in kinetic theory and proposed for the standard Cahn-Hilliard equation, by Carlen, Carvalho and Orlandi, \cite {CCO}. The constructed approximate solutions allow to derive rigorously the sharp interface limit of the generalized Cahn-Hilliard equations. We then estimate the difference between the true solutions and the approximate solutions by spectral analysis, as in \cite {A-B-C}
△ Less
Submitted 5 January, 2013;
originally announced January 2013.
-
Crank-Nicolson Finite Element Discretizations for a 2D Linear Schrödinger-Type Equation Posed in a Noncylindrical Domain
Authors:
D. C. Antonopoulou,
G. D. Karali,
M. Plexousakis,
G. E. Zouraris
Abstract:
Motivated by the paraxial narrow-angle approximation of the Helmholtz equation in domains of variable topography that appears as an important application in Underwater Acoustics, we analyze a general Schrödinger-type equation posed on two-dimensional variable domains with mixed boundary conditions. The resulting initial- and boundary-value problem is transformed into an equivalent one posed on a r…
▽ More
Motivated by the paraxial narrow-angle approximation of the Helmholtz equation in domains of variable topography that appears as an important application in Underwater Acoustics, we analyze a general Schrödinger-type equation posed on two-dimensional variable domains with mixed boundary conditions. The resulting initial- and boundary-value problem is transformed into an equivalent one posed on a rectangular domain and is approximated by fully discrete, $L^2$-stable, finite element, Crank--Nicolson type schemes. We prove a global elliptic regularity theorem for complex elliptic boundary value problems with mixed conditions and derive $L^2$-error estimates of optimal order. Numerical experiments are presented which verify the optimal rate of convergence.
△ Less
Submitted 18 October, 2011;
originally announced October 2011.
-
A Finite Difference method for the Wide-Angle `Parabolic' equation in a waveguide with downsloping bottom
Authors:
D. C. Antonopoulou,
V. A. Dougalis,
G. E. Zouraris
Abstract:
We consider the third-order wide-angle `parabolic' equation of underwater acoustics in a cylindrically symmetric fluid medium over a bottom of range-dependent bathymetry. It is known that the initial-boundary-value problem for this equation may not be well posed in the case of (smooth) bottom profiles of arbitrary shape if it is just posed e.g. with a homogeneous Dirichlet bottom boundary conditio…
▽ More
We consider the third-order wide-angle `parabolic' equation of underwater acoustics in a cylindrically symmetric fluid medium over a bottom of range-dependent bathymetry. It is known that the initial-boundary-value problem for this equation may not be well posed in the case of (smooth) bottom profiles of arbitrary shape if it is just posed e.g. with a homogeneous Dirichlet bottom boundary condition. In this paper we concentrate on downsloping bottom profiles and propose an additional boundary condition that yields a well posed problem, in fact making it $L^2$-conservative in the case of appropriate real parameters. We solve the problem numerically by a Crank-Nicolson-type finite difference scheme, which is proved to be unconditionally stable and second-order accurate, and simulates accurately realistic underwater acoustic problems.
△ Less
Submitted 15 October, 2011;
originally announced October 2011.
-
Galerkin Methods for Parabolic and SCHR{Ö}DINGER Equations with Dynamical Boundary Conditions and Applications to Underwater Acoustics
Authors:
D. C. Antonopoulou,
V. A. Dougalis,
G. E. Zouraris
Abstract:
In this paper we consider Galerkin-finite element methods that approximate the solutions of initial-boundary-value problems in one space dimension for parabolic and Schrödinger evolution equations with dynamical boundary conditions. Error estimates of optimal rates of convergence in $L^2$ and $H^1$ are proved for the accociated semidiscrete and fully discrete Crank-Nicolson-Galerkin approximatio…
▽ More
In this paper we consider Galerkin-finite element methods that approximate the solutions of initial-boundary-value problems in one space dimension for parabolic and Schrödinger evolution equations with dynamical boundary conditions. Error estimates of optimal rates of convergence in $L^2$ and $H^1$ are proved for the accociated semidiscrete and fully discrete Crank-Nicolson-Galerkin approximations. The problem involving the Schrödinger equation is motivated by considering the standard `parabolic' (paraxial) approximation to the Helmholtz equation, used in underwater acoustics to model long-range sound propagation in the sea, in the specific case of a domain with a rigid bottom of variable topography. This model is contrasted with alternative ones that avoid the dynamical bottom boundary condition and are shown to yield qualitatively better approximations. In the (real) parabolic case, numerical approximations are considered for dynamical boundary conditions of reactive and dissipative type.
△ Less
Submitted 24 April, 2009;
originally announced April 2009.