Search

We present a new inference framework for neutron star astrophysics based on conditional variational autoencoders. Once trained, the generator block of the model reconstructs the neutron star equation of state from a given set of mass-radius observations. While the pressure of dense matter is the focus of the present study, the proposed model is flexible enough to accommodate the reconstructing of any other quantity related to dense matter equation of state. Our results show robust reconstructing performance of the model, allowing to make instantaneous inference from any given observation set. The present framework contrasts with computationally expensive evaluation of the equation-of-state posterior probability distribution based on Markov chain Monte Carlo methods., Comment: 9 pages, 7 figures

We explore the anomaly detection framework based on Normalizing Flows (NF) models introduced in \cite{PhysRevC.106.065802} to detect the presence of a large (destabilising) dense matter phase transition in neutron star (NS) observations of masses and radii, and relate the feasibility of detection with parameters of the underlying mass-radius sequence, which is a functional of the dense matter equation of state. Once trained on simulated data featuring continuous $M(R)$ solutions (i.e., no phase transitions), NF is used to determine the likelihood of a first-order phase transition in a given set of $M(R)$ observations featuring a discontinuity, i.e., perform the anomaly detection. Different mock test sets, featuring two branch solutions in the $M(R)$ diagram, were parameterized by the NS mass at which the phase transition occurs, $M_c$, and the radius difference between the heaviest hadronic star and lightest hybrid star, $\Delta R$. We analyze the impact of these parameters on the NF performance in detecting the presence of a first-order phase transition. Among the results, we report that given a set of 15 stars with radius uncertainty of $0.2$ km, a detection of a two-branch solution is possible with 95\% accuracy if $\Delta R > 0.4$ km., Comment: 9 pages, 8 figures

In this paper, we present a new method to search for a short, a few tens of milliseconds long, postmerger gravitational-wave signal following the merger of two neutron stars. Such a signal could follow the event GW170817 observed by LIGO and Virgo detectors. Our method is based on a matched filtering statistic and an approximate template of the postmerger signal in the form of a damped sinusoid. We test and validate our method using postmerger numerical simulations from the CoRe database. We find no evidence of the short postmerger signal in the LIGO data following the GW170817 event and we obtain upper limits. For short postmerger signals investigated, our best upper limit on the root sum square of the gravitational-wave strain emitted from 1.15 kHz to 4 kHz is $h_{\text{rss}}^{50\%}=1.8\times 10^{-22}/\sqrt{\text{Hz}}$ at 50% detection efficiency. The distance corresponding to this best upper limit is 4.64 Mpc., Comment: LaTeX, 30 pages, 15 figures, identical to the published version

Observations of neutron stars may be used to study aspects of extremely dense matter, specifically a possibility of phase transitions to exotic states, such as de-confined quarks. We present a novel data analysis method for detecting signatures of dense-matter phase transitions in sets of mass-radius measurements, and study its sensitivity with respect to the size of observational errors and the number of observations. The method is based on machine learning anomaly detection coupled with normalizing flows technique: the algorithm trained on samples of astrophysical observations featuring no phase transition signatures interprets a phase transition sample as an ''anomaly''. For the sake of this study, we focus on dense-matter equations of state leading to detached branches of mass-radius sequences (strong phase transitions), use an astrophysically-informed neutron-star mass function, and various magnitudes of observational errors and sample sizes. The method is shown to reliably detect cases of mass-radius relations with phase transition signatures, while increasing its sensitivity with decreasing measurement errors and increasing number of observations. We discuss marginal cases, when the phase transition mass is located near the edges of the mass function range. Evaluated on the current state-of-art selection of real measurements of electromagnetic and gravitational-wave observations, the method gives inconclusive results, which we interpret as due to small available sample size, large observational errors and complex systematics., Comment: 9 pages, 6 figures, accepted for publication in Physical Review C

It is currently unknown if neutron stars (NSs) are composed of nucleons only or are hybrid stars, i.e., in addition to nucleonic crusts and outer cores, they also possess quark cores. Quantum chromodynamics allows for such a possibility, but accurate calculations relevant for compact stars are still elusive. Here we investigate some crust-breaking aspects of hybrid stars. We show that the crust-breaking frequency and maximum fiducial ellipticity are sensitive to the quark-hadron density jump and equation of state stiffness. Remarkably, the crust-breaking frequency related to static tides scales linearly with the mass of the star (for a given companion's mass), and its slope encompasses information about the microphysics of the star. However, for precise crust-breaking frequency predictions, relativistic corrections to Kepler's third law and the Newtonian tidal field should not be ignored. When a liquid quark core touches an elastic hadronic phase (the result of a significant energy-density jump), the maximum ellipticity can increase around an order of magnitude when compared to a liquid quark core touching a liquid hadronic phase. That is relevant because it would increase the odds of detecting continuous gravitational waves from NSs. Our order-of-magnitude analysis also suggests that a given upper limit to the ellipticity(crust-breaking frequency) could have representatives in stars with either small or intermediate(large) energy-density jumps. Therefore, when upper limits to the ellipticity for isolated stars are better constrained or electromagnetic radiation (e.g., gamma-ray precursors) is detected along with gravitational waves in inspiraling binary systems, they may help constrain some aspects of phase transitions in NSs., Comment: 12 pages, 6 figures. Improved discussions on boundary conditions, maximum fiducial ellipticities, crust-breaking frequencies, error estimation due to relativistic corrections to some Newtonian quantities used (orbital frequency and tidal field). Accepted for publication in ApJ

3XMM J185246.6+003317 is a transient magnetar located in the vicinity of the supernova remnant Kes\,79. So far, observations have only set upper limits to its surface magnetic field and spindown, and there is no estimate for its mass and radius. Using ray-tracing modelling and Bayesian inference for the analysis of several light curves spanning a period of around three weeks, we have found that it may be one of the most massive neutron stars to date. In addition, our analysis suggests a multipolar magnetic field structure with a subcritical field strength and a carbon atmosphere composition. Due to the time-resolution limitation of the available light curves, we estimate the surface magnetic field and the mass to be $\log_{10} (B/{\rm G}) = 11.89^{+0.19}_{-0.93}$ and $M=2.09^{+0.16}_{-0.09}$~$M_{\odot}$ at $1\sigma$ confidence level, while the radius is estimated to be $R=12.02^{+1.44}_{-1.42}$ km at $2\sigma$ confidence level. They were verified by simulations, i.e., data injections with known model parameters, and their subsequent recovery. The best-fitting model has three small hot spots, two of them in the southern hemisphere. These are, however, just first estimates and conclusions, based on a simple ray-tracing model with anisotropic emission; we also estimate the impact of modelling on the parameter uncertainties and the relevant phenomena on which to focus in more precise analyses. We interpret the above best-fitting results as due to accretion of supernova layers/interstellar medium onto 3XMM J185246.6+003317 leading to burying and a subsequent re-emergence of the magnetic field, and a carbon atmosphere being formed possibly due to hydrogen/helium diffusive nuclear burning. Finally, we briefly discuss some consequences of our findings for superdense matter constraints., Comment: 14 pages, 5 figures, 4 tables. Accepted for publication in Journal of High Energy Astrophysics (JHEAP)

If very low mass primordial black holes (PBH) within the asteroid/moon-mass range indeed reside in galactic dark matter halos, they must necessarily collide with galactic neutron stars (NSs). These collisions must, again necessarily, form light black holes (LBHs) with masses of typical NSs, $M_{\rm LBH} \approx \,1-2\,M_{\odot}$. LBHs may be behind events already detected by ground-based gravitational-wave detectors (GW170817, GW190425, and others such as a mixed stellar black hole-neutron star mass event GW191219_163120), and most recently by microlensing (OGLE-BLG-2011-0462). Although the status of these observations as containing LBHs is not confirmed, there is no question that gravitational-wave detectors and microlensing are in principle and in practice capable of detecting LBHs. We have calculated the creation rate of LBHs resulting from these light primordial black hole collisions with neutron stars. On this basis, we claim that if improved gravitational-wave detectors and microlensing statistics of the LBH events would indicate that the number of LBHs is significantly lower that what follows from the calculated creation rate, then this would be an unambiguous proof that there is no significant light PBH contribution to the galactic dark matter halos. Otherwise, if observed and calculated numbers of LBHs roughly agree, then the hypothesis of primordial black hole existence gets strong observational support, and in addition their collisions with neutron stars may be considered a natural creation channel for the LBHs, solving the problem of their origin, as it is known that they cannot be a product of standard stellar evolution.

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., Comment: 18 pages, 17 figures; accepted for publication in ApJ; very minor edits and no changes to numbers

The internal composition of neutron stars is still an open issue in astrophysics. Their innermost regions are impervious to light propagation and gravitational waves mostly carry global aspects of stars, meaning that only indirect inferences of their interiors could be obtained. Here we assume a hypothetical future scenario in which an equation of state softening due to a phase transition is identified and estimate the observational accuracy to differentiate a sharp phase transition from a smoother one (associated with a mixed phase/state due to the unknown value of the surface tension of dense matter) in a region of a hybrid star by means of some electromagnetic and gravitational wave observables. We show that different transition constructions lead to similar sequences of stellar configurations due to their shared thermodynamic properties. In the most optimistic case - a strong quark-hadron density jump phase transition - radius observations require fractional uncertainties smaller than $1\%-2\%$ to differentiate mixed states from sharp phase transitions. For tidal deformabilities, relative uncertainties should be smaller than $5\%-10\%$. However, for masses around the onset of stable quark cores, relative tidal deformability differences associated with strong sharp phase transitions and mixed states could be much larger (up to around $20\%-30\%$). All the above suggests that 2.5- and 3rd generation gravitational wave detectors and near-term electromagnetic missions may be able to start assessing some particular aspects of phase transitions in neutron stars. In addition, it points to some limitations on the equation of state recovery using typical neutron star observables and the impact of systematic uncertainties on modellings of the equation of state of hybrid stars. Finally, we briefly discuss other observables that may also be relevant for the probe of mixed states in stars., Comment: 17 pages, 15 figures (some updated); Improved abstract, presentation and discussions (assumptions, phase transitions, EOS recovery, predictions, etc); New references added; Accepted for publication in PRD

As of this moment, fifty gravitational waves (GW) detections have been announced, thanks to the observational efforts of the LIGO-Virgo Collaboration, working with the Advanced LIGO and the Advanced Virgo interferometers. The detection of signals is complicated by the noise-dominated nature of the data. Conventional approaches in GW detection procedures require either precise knowledge of the GW waveform in the context of matched filtering searches or coincident analysis of data from multiple detectors. Furthermore, the analysis is prone to contamination by instrumental or environmental artifacts called glitches which either mimic astrophysical signals or reduce the overall quality of data. In this paper, we propose an alternative generic method of studying GW data based on detecting anomalies. The anomalies we study are transient signals, different from the slow non-stationary noise of the detector. Presented in the manuscript anomalies are mostly based on the GW emitted by the mergers of binary black hole systems. However, the presented study of anomalies is not limited only to GW alone, but also includes glitches occurring in the real LIGO/Virgo dataset available at the Gravitational Waves Open Science Center., Comment: Accepted to the Machine Learning: Science and Technology; 29 pages, 18 figures

Aims. We estimate the number of pulsars, detectable as continuous gravitational wave sources with the current and future gravitational-wave detectors, assuming a simple phenomenological model of evolving non-axisymmetry of the rotating neutron star. Methods. We employ a numerical model of the Galactic neutron star population, with the properties established by comparison with radio observations of isolated Galactic pulsars. We generate an arbitrarily large synthetic population of neutron stars and evolve their period, magnetic field, and position in space. We use a gravitational wave emission model based on exponentially decaying ellipticity - a non-axisymmetry of the star, with no assumption of the origin of a given ellipticity. We calculate the expected signal in a given detector for a 1 year observations and assume a detection criterion of the signal-to-noise ratio of 11.4 - comparable to a targeted continous wave search. We analyze the population detectable separately in each detector: Advanced LIGO, Advanced Virgo, and the planned Einstein Telescope. In the calculation of the expected signal we neglect signals frequency change due to the source spindown and the Earth motion with respect to the Solar barycentre. Results. With conservative values for the neutron stars evolution: supernova rate once per 100 years, initial ellipticity $\epsilon_{0}$ = 1e-5 with no decay of the ellipticity $\eta$ = $t_\rm{hub}$ = 1e4 Myr, the expected number of detected neutron stars is below one: 0.15 (based on a simulation of 10 M stars) for the Advanced LIGO detector. A broader study of the parameter space ($\epsilon_{0}$ , $\eta$) is presented. With the planned sensitivity for the Einstein Telescope, and assuming the same ellipiticity model, the expected detection number is: 26.4 pulsars during a 1-year long observing run., Comment: 10 pages, submitted to A&A

The internal composition of neutron stars is currently largely unknown. Due to the possibility of phase transitions in quantum chromodynamics, stars could be hybrid and have quark cores. We investigate some imprints of elastic quark phases (only when perturbed) on the dynamical stability of hybrid stars. We show that they increase the dynamical stability window of hybrid stars in the sense that the onset of instabilities happen at larger central densities than the ones for maximum masses. In particular, when the shear modulus of a crystalline quark phase is taken at face value, the relative radius differences between elastic and perfect-fluid hybrid stars with null radial frequencies (onset of instability) would be up to $1-2\%$. Roughly, this would imply a maximum relative radius dispersion (on top of the perfect-fluid predictions) of $2-4\%$ for stars in a given mass range exclusively due to the elasticity of the quark phase. In the more agnostic approach where the estimates for the quark shear modulus only suggest its possible order of magnitude (due to the many approximations taken in its calculation), the relative radius dispersion uniquely due to a quark phase elasticity might be as large as $5-10\%$. Finally, we discuss possible implications of the above dispersion of radii for the constraint of the elasticity of a quark phase with electromagnetic missions such as NICER, eXTP and ATHENA., Comment: 15 pages, 2 figures. New discussions (EOS uncertainties) inserted. Accepted for publication in ApJ

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., Comment: 11 pages, 10 figures; accepted for publication in MNRAS

PSR J1813-1749 is one of the most energetic rotation-powered pulsars known, producing a pulsar wind nebula (PWN) and gamma-ray and TeV emission, but whose spin period is only measurable in X-ray. We present analysis of two Chandra datasets that are separated by more than ten years and recent NICER data. The long baseline of the Chandra data allows us to derive a pulsar proper motion mu_R.A.=-(0.067"+/-0.010") yr^-1 and mu_decl.=-(0.014"+/-0.007") yr^-1 and velocity v_perp~900-1600 km/s (assuming a distance d=3-5 kpc), although we cannot exclude a contribution to the change in measured pulsar position due to a change in brightness structure of the PWN very near the pulsar. We model the PWN and pulsar spectra using an absorbed power law and obtain best-fit absorption NH=(13.1+/-0.9)x10^22 cm^-2, photon index Gamma=1.5+/-0.1, and 0.3-10 keV luminosity Lx~5.4x10^34 erg/s (d/5 kpc)^2 for the PWN and Gamma=1.2+/-0.1 and Lx~9.3x10^33 erg/s (d/5 kpc)^2 for PSR J1813-1749. These values do not change between the 2006 and 2016 observations. We use NICER observations from 2019 to obtain a timing model of PSR J1813-1749, with spin frequency nu=22.35 Hz and spin frequency time derivative nudot=(-6.428+/-0.003)x10^-11 Hz/s. We also fit nu measurements from 2009-2012 and our 2019 value and find a long-term spin-down rate nudot=(-6.3445+/-0.0004)x10^-11 Hz/s. We speculate that the difference in spin-down rates is due to glitch activity or emission mode switching., Comment: 9 pages, 5 figures; accepted for publication in MNRAS

Aims: The aim of this work is to study the application of the artificial neural networks guided by the autoencoder architecture as a method for precise reconstruction of the neutron star equation of state, using their observable parameters: masses, radii and tidal deformabilities. In addition we study how well the neutron star radius can be reconstructed using the gravitational-wave only observations of tidal deformability, i.e. quantities which are not related in a straightforward way. Methods: Application of artificial neural network in the equation of state reconstruction exploits the non-linear potential of this machine learning model. Since each neuron in the network is basically a non-linear function, it is possible to create a complex mapping between the input sets of observations and the output equation of state table. Within the supervised training paradigm, we construct a few hidden layer deep neural network on a generated data set, consisting of a realistic equation of state for the neutron star crust connected with a piecewise relativistic polytropes dense core, with parameters representative to the state-of-the art realistic equations of state. Results: We demonstrate the performance of our machine learning implementation with respect to the simulated cases with varying number of observations and measurement uncertainties. Furthermore we study the impact of the neutron star mass distributions on the results. Finally, we test the reconstruction of the equation of state trained on parametric polytropic training set using the simulated mass--radius and mass--tidal-deformability sequences based on realistic equations of state. Neural networks trained with a limited data set are able to generalize the mapping between global parameters and equation of state input tables for realistic models., Comment: 8, pages, 7 figures, accepted in Astronomy and Astrophysics

Machine learning has emerged as a popular and powerful approach for solving problems in astrophysics. We review applications of machine learning techniques for the analysis of ground-based gravitational-wave detector data. Examples include techniques for improving the sensitivity of Advanced LIGO and Advanced Virgo gravitational-wave searches, methods for fast measurements of the astrophysical parameters of gravitational-wave sources, and algorithms for reduction and characterization of non-astrophysical detector noise. These applications demonstrate how machine learning techniques may be harnessed to enhance the science that is possible with current and future gravitational-wave detectors.

Gravitational wave astronomy is expected to provide independent constraints on neutron star properties, such as their dense matter equation of state. This is possible with the measurements of binary components' tidal deformability, which alter the point-particle gravitational waveforms of the late inspiral phase of neutron-star binaries. Although current gravitational wave detectors are not sensitive enough for a precise determination of the individual tidal deformations of the components, a large number of combined observations with future detectors will decrease uncertainties in this quantity. Here we provide a first study of the tidal deformability effects due to the elasticity/solidity of the crust (hadronic phase) in a hybrid neutron star, as well as the influence of a quark-hadronic phase density jump on tidal deformations. We employ the framework of nonradial perturbations with zero frequency and study hadronic phases presenting elastic aspects when perturbed (with the shear modulus approximately $1\%$ of the pressure). We find that the relative tidal deformation change in a hybrid star with a perfect-fluid quark phase and a hadronic phase presenting an elastic part is never larger than about $2-4\%$ (with respect to a perfect-fluid counterpart). These maximum changes occur when the elastic region of a hybrid star is larger than approximately $60\%$ of the star's radius, which may happen when its quark phase is small and the density jump is large enough, or even when a hybrid star has an elastic mixed phase. For other cases, the relative tidal deformation changes due to an elastic crust are negligible ($10^{-5}-10^{-1}\%$), therefore unlikely to be measured even with third generation detectors. Thus, only when the size of the elastic hadronic region of a hybrid star is over half of its radius, the effects of elasticity could have a noticeable impact on tidal deformations., Comment: 11 pages, 3 figures. Minor changes to the abstract. Accepted for publication in ApJ

Gravitational waves astronomy allows us to study objects and events invisible in electromagnetic waves. It is crucial to validate the theories and models of the most mysterious and extreme matter in the Universe: the neutron stars. In addition to inspirals and mergers of neutrons stars, there are currently a few proposed mechanisms that can trigger radiation of long-lasting gravitational radiation from neutron stars, such as e.g., elastically and/or magnetically driven deformations: mountains on the stellar surface supported by the elastic strain or magnetic field, free precession, or unstable oscillation modes (e.g., the r-modes). The astrophysical motivation for continuous gravitational waves searches, current LIGO and Virgo strategies of data analysis and prospects are reviewed in this work., Comment: 47 pages; 7 figures; published Universe 'Neutron Star Astrophysics' special issue

Among astrophysical sources in the Advanced LIGO and Advanced Virgo detectors' frequency band are rotating non-axisymmetric neutron stars emitting long-lasting, almost-monochromatic gravitational waves. Searches for these continuous gravitational-wave signals are usually performed in long stretches of data in a matched-filter framework e.g., the F-statistic method. In an all-sky search for a priori unknown sources, large number of templates is matched against the data using a pre-defined grid of variables (the gravitational-wave frequency and its derivatives, sky coordinates), subsequently producing a collection of candidate signals, corresponding to the grid points at which the signal reaches a pre-defined signal-to-noise threshold. An astrophysical signature of the signal is encoded in the multi-dimensional vector distribution of the candidate signals. In the first work of this kind, we apply a deep learning approach to classify the distributions. We consider three basic classes: Gaussian noise, astrophysical gravitational-wave signal, and a constant-frequency detector artifact ("stationary line"), the two latter injected into the Gaussian noise. 1D and 2D versions of a convolutional neural network classifier are implemented, trained and tested on a broad range of signal frequencies. We demonstrate that these implementations correctly classify the instances of data at various signal-to-noise ratios and signal frequencies, while also showing concept generalization i.e., satisfactory performance at previously unseen frequencies. In addition we discuss the deficiencies, computational requirements and possible applications of these implementations., Comment: 22 pages, 11 figures, accepted by MLST

Among promising sources of gravitational waves are long-lived nearly periodic signals produced by rotating, asymmetric neutron stars. Depending on the astrophysical scenario, the sources of asymmetry may have thermal, viscous, elastic and/or magnetic origin. In this work we introduce a follow-up procedure for an all-sky search for gravitational wave signals from rotating neutron stars. The procedure denoted as Followup implements matched-filtering $\mathcal{F}$-statistic method. We describe data analysis methods and algorithms used in the procedure. We present tests of the Followup for artificial signals added to white, Gaussian noise. The tests show a good agreement with the theoretical predictions. The Followup will become part of the Time-Domain $\mathcal{F}$-statistic pipeline that is routinely used for all-sky searches of LIGO and Virgo detector data., Comment: Published in CQG; 23 pages, 11 figures

Our existence in the Universe resulted from a rare combination of circumstances. The same must hold for any highly developed extraterrestrial civilisation, and if they have ever existed in the Milky Way, they would likely be scattered over large distances in space and time. However, all technologically advanced species must be aware of the unique property of the galactic centre: it hosts Sagittarius A* (Sgr A*), the closest supermassive black hole to anyone in the Galaxy. A civilisation with sufficient technical know-how may have placed material in orbit around Sgr A* for research, energy extraction, and communication purposes. In either case, its orbital motion will necessarily be a source of gravitational waves. We show that a Jupiter-mass probe on the retrograde innermost stable circular orbit around Sgr A* emits, depending on the black hole spin, at a frequency of $f_{GW} = 0.63 - 1.07$ mHz and with a power of $P_{GW}=2.7 \times\, 10^{36} - 2.0 \times\, 10^{37}$ erg/s. We discuss that the energy output of a single star is sufficient to stabilise the location of an orbiting probe for a billion years against gravitational wave induced orbital decay. Placing and sustaining a device near Sgr A* is therefore astrophysically possible. Such a probe will emit an unambiguously artificial continuous gravitational wave signal that is observable with LISA-type detectors., Comment: 7 pages, 2 figures, published in Scientific Reports

The global network of ground-based gravitational-wave detectors (the Advanced LIGO and the Advanced Virgo) is sensitive at the frequency range corresponding to relativistic stellar-mass compact objects. Among the promising types of gravitational-wave sources are binary systems and rotating, deformed neutron stars. I will describe these sources and present predictions of how their observations will contribute to modern astrophysics in the near future., Comment: 6 pages; contribution to XXXVIII Polish Astronomical Society meeting, Zielona G\'ora, 11-14 September 2017

Potentially interesting gravitational-wave candidates (outliers) from the blind all-sky searches have to be confirmed or rejected by studying their origin and precisely estimating their parameters. We present the design and first results for the followup procedure of the {\tt Polgraw} all-sky search pipeline: a coherent search for almost-monochromatic gravitational-wave signals in several-day long time segments using the $F$-statistic method followed by the coincidences between the candidate signals. Approximate parameters resulting in these two initial steps are improved in the final followup step, in which the signals from detectors are studied separately, together with the network combination of them, and the true parameters and signal-to-noise values are established., Comment: 4 pages, 4 figures, published in Proceedings of the Polish Astronomical Society, vol. 7, 37-40 (2018)

Using parametric equations of state (relativistic polytropes and a simple quark bag model) to model dense-matter phase transitions, we study global, measurable astrophysical parameters of compact stars such as their allowed radii and tidal deformabilities. We also investigate the influence of stiffness of matter before the onset of the phase transitions on the parameters of the possible exotic dense phase. The aim of our study is to compare the parameter space of the dense matter equation of state permitting phase transitions to a sub-space compatible with current observational constraints such as the maximum observable mass, tidal deformabilities of neutron star mergers, radii of configurations before the onset of the phase transition, and to give predictions for future observations. We studied solutions of the Tolman-Oppenheimer-Volkoff equations for a flexible set of parametric equations of state. We compare tidal deformabilities of stars with weak and strong phase transitions with the results of the GW170817 neutron star merger. Specifically, we study characteristic phase transition features in the $\Lambda_1-\Lambda_2$ relation, and estimate the deviations of our results from the approximate formul{\ae} for $\tilde{\Lambda}-R(M_1)$ and $\Lambda$-compactness proposed in the literature. We find constraints on the hybrid equations of state to produce stable neutron stars on the twin branch. For the exemplary equations of state most of the high-mass twins occur for the minimum values of the density jump $\lambda = 1.33-1.54$; corresponding values of the square of the speed of sound are $\alpha = 0.7-0.37$. We compare results with observations of gravitational waves and with the theoretical causal limit and find that the minimum radius of a twin branch is 9.5 - 10.5 km. For these solutions, the phase transition occurs below 0.56 [fm$^{-3}$]., Comment: 11 pages, 21 figures, published in A&A

The formal division of the distance indicators into primary and secondary leads to difficulties in description of methods which can actually be used in two ways: with, and without the support of the other methods for scaling. Thus instead of concentrating on the scaling requirement we concentrate on all methods of distance determination to extragalactic sources which are designated, at least formally, to use for individual sources. Among those, the Supernovae Ia is clearly the leader due to its enormous success in determination of the expansion rate of the Universe. However, new methods are rapidly developing, and there is also a progress in more traditional methods. We give a general overview of the methods but we mostly concentrate on the most recent developments in each field, and future expectations., Comment: Space Science Rewiev (in press)

Periodic (almost monochromatic) gravitational waves emitted by rotating, asymmetric neutron stars are intriguing potential signals in the sensitivity band of Advanced LIGO and Advanced Virgo detectors. These signals are related to elastic and magnetic stresses in the neutron-star interior, as well as to various possible instabilities, and thus are interesting from the point of view of the largely-unknown neutron star structure. I will describe the main challenges related to these searches, the current state of the data-analysis methods and plans for the future., Comment: 8 pages, 3 figures (summary talk at the Rencontres de Moriond 2017 Gravitation session)

If primordial black holes with masses of $10^{25}\,\mbox{g}\gtrsim m \gtrsim 10^{17}\,\mbox{g}$ constitute a non-negligible fraction of the galactic dark-matter haloes, their existence should have observable consequences: they necessarily collide with galactic neutron stars, nest in their centers and accrete the dense matter, eventually converting them to neutron-star mass black holes while releasing the neutron-star magnetic field energy. Such processes may explain the fast radio bursts phenomenology, in particular their millisecond durations, large luminosities ${\sim}10^{43}$ erg/s, high rate of occurrence $\gtrsim 1000/\mbox{day}$, as well as high brightness temperatures, polarized emission and Faraday rotation. Longer than the dynamical timescale of the Bondi-like accretion for light primordial black holes allows for the repeating fast radio bursts. This explanation follows naturally from (assumed) existence of the dark matter primordial black holes and requires no additional unusual phenomena, in particular no unacceptably large magnetic fields of neutron stars. In our model, the observed rate of fast radio bursts throughout the Universe follows from the presently known number of neutron stars in the Galaxy., Comment: 9 pages, 6 figures. Accepted by ApJ

We study the structure and evolution of the hyperaccreting disks and outflows in the gamma ray bursts central engines. The torus around a stellar mass black hole is composed of free nucleons, Helium, electron-positron pairs, and is cooled by neutrino emission. Accretion of matter powers the relativistic jets, responsible for the gamma ray prompt emission. The significant number density of neutrons in the disk and outflowing material will cause subsequent formation of heavier nuclei. We study the process of nucleosynthesis and its possible observational consequences. We also apply our scenario to the recent observation of the gravitational wave signal, detected on September 14th, 2015 by the two Advanced LIGO detectors, and related to an inspiral and merger of a binary black hole system. A gamma ray burst that could possibly be related with the GW150914 event was observed by the Fermi satellite. It had a duration of about 1 second and appeared about 0.4 seconds after the gravitational-wave signal. We propose that a collapsing massive star and a black hole in a close binary could lead to the event. The gamma ray burst was powered by a weak neutrino flux produced in the star remnant's matter. Low spin and kick velocity of the merged black hole are reproduced in our simulations. Coincident gravitational-wave emission originates from the merger of the collapsed core and the companion black hole., Comment: 11 pages, 4 figures. Accepted to Special Issue of Galaxies. Gamma-Ray Bursts. Recent Theoretical Models and Observations. Based on the contributions presented during Symposium 11 at the European Week of Astronomy and Space Science (4-5 July, 2016, Athens, Greece). Eds. I. Contopoulos, S. Oates

Rapidly rotating neutron stars are promising sources of continuous gravitational wave radiation for the LIGO and Virgo interferometers. The majority of neutron stars in our galaxy have not been identified with electromagnetic observations. All-sky searches for isolated neutron stars offer the potential to detect gravitational waves from these unidentified sources. The parameter space of these blind all-sky searches, which also cover a large range of frequencies and frequency derivatives, presents a significant computational challenge. Different methods have been designed to perform these searches within acceptable computational limits. Here we describe the first benchmark in a project to compare the search methods currently available for the detection of unknown isolated neutron stars. We employ a mock data challenge to compare the ability of each search method to recover signals simulated assuming a standard signal model. We find similar performance among the short duration search methods, while the long duration search method achieves up to a factor of two higher sensitivity. We find the absence of second derivative frequency in the search parameter space does not degrade search sensivity for signals with physically plausible second derivative frequencies. We also report on the parameter estimation accuracy of each search method, and the stability of the sensitivity in frequency, frequency derivative and in the presence of detector noise., Comment: 16 pages, 11 figures

Data from the Fermi Gamma-ray Burst Monitor satellite observatory suggested that the recently discovered gravitational wave source, a pair of two coalescing black holes, was related to a gamma-ray burst. The observed high-energy electromagnetic radiation (above 50 keV) originated from a weak transient source and lasted for about 1 second. Its localization is consistent with the direction to GW150914. We speculate about the possible scenario for the formation of a gamma-ray burst accompanied by the gravitational-wave signal. Our model invokes a tight binary system consisting of a massive star and a black hole which leads to the triggering of a collapse of the star's nucleus, the formation of a second black hole, and finally to the binary black hole merger. For the most-likely configuration of the binary spin vectors with respect to the orbital angular momentum in the GW150914 event, the recoil speed (kick velocity) acquired by the final black hole through gravitational wave emission is of the order of a few hundred km/s and this might be sufficient to get it closer to the envelope of surrounding material and capture a small fraction of matter from the remnant of the host star. The gamma-ray burst is produced by the accretion of this remnant matter onto the final black hole. The moderate spin of the final black hole suggests that the gamma-ray burst jet is powered by weak neutrino emission rather than the Blandford-Znajek mechanism, and hence explains the low power available for the observed GRB signal., Comment: 20 pages, 6 figures; to be published in New Astronomy

The SageManifolds project aims at extending the mathematics software system Sage towards differential geometry and tensor calculus. Like Sage, SageManifolds is free, open-source and is based on the Python programming language. We discuss here some details of the implementation, which relies on Sage's parent/element framework, and present a concrete example of use., Comment: 17 pages, 6 figures, many small changes w.r.t. v1; to appear in the Proceedings of "ERE2014: Almost 100 years after Einstein's revolution" (Journal of Physics: Conference Series)

The parallelization, design and scalability of the \sky code to search for periodic gravitational waves from rotating neutron stars is discussed. The code is based on an efficient implementation of the F-statistic using the Fast Fourier Transform algorithm. To perform an analysis of data from the advanced LIGO and Virgo gravitational wave detectors' network, which will start operating in 2015, hundreds of millions of CPU hours will be required - the code utilizing the potential of massively parallel supercomputers is therefore mandatory. We have parallelized the code using the Message Passing Interface standard, implemented a mechanism for combining the searches at different sky-positions and frequency bands into one extremely scalable program. The parallel I/O interface is used to escape bottlenecks, when writing the generated data into file system. This allowed to develop a highly scalable computation code, which would enable the data analysis at large scales on acceptable time scales. Benchmarking of the code on a Cray XE6 system was performed to show efficiency of our parallelization concept and to demonstrate scaling up to 50 thousand cores in parallel., Comment: 11 pages, 9 figures. Submitted to Computer Physics Communications

We consider a scenario for the longest duration gamma ray bursts, resulting from the collapse of a massive rotating star in a close binary system with a companion black hole. The primary black hole born during the core collapse is first being spun up and increases its mass during the fallback of the stellar envelope just after its birth. As the companion black hole enters the outer envelope, it provides an additional angular momentum to the gas. After the infall and spiral-in towards the primary, the two black holes merge inside the circumbinary disk. The second episode of mass accretion and high final spin of the post-merger black hole prolongs the gamma ray burst central engine activity. The observed events should have two distinct peaks in the electromagnetic signal, separated by the gravitational wave emission. The gravitational recoil of the burst engine is also possible., Comment: 9 pages, 4 figures, accepted for publication in A&A

Collisions of particles in black holes' ergospheres may result in an arbitrarily large center of mass energy. This led recently to the suggestion (Banados et al., 2009) that black holes can act as ultimate particle accelerators. If the energy of an outgoing particle is larger than the total energy of the infalling particles the energy excess must come from the rotational energy of the black hole and hence this must involve a Penrose process. However, while the center of mass energy diverges the position of the collision makes it impossible for energetic particles to escape to infinity. Following an earlier work on collisional Penrose processes (Piran & Shaham 1977) we show that even under the most favorable idealized conditions the maximal energy of an escaping particle is only a modest factor above the total initial energy of the colliding particles. This implies that one shouldn't expect collisions around a black hole to act as spectacular cosmic accelerators., Comment: 4 pages, 5 figures (minor corrections, matches PRL accepted version)

Using numerical simulations based on solving the general relativistic hydrodynamic equations, we study the dynamics of a phase transition in the dense core of isolated rotating neutron stars, triggered by the back bending instability reached via angular momentum loss. In particular, we investigate the dynamics of a migration from an unstable configuration into a stable one, which leads to a mini-collapse of the neutron star and excites sizeable pulsations in its bulk until it acquires a new stable equilibrium state. We consider equations of state with softening at high densities, a simple analytic one with a mixed hadron-quark phase in an intermediate pressure interval and pure quark matter at very high densities, and a microphysical one that has a first-order phase transition, originating from kaon condensation. Although the marginally stable initial models are rigidly rotating, we observe that during the collapse (albeit little) differential rotation is created. We analyze the emission of gravitational radiation, which in some models is amplified by mode resonance effects, and assess its prospective detectability by interferometric detectors. We expect that the most favorable conditions for dynamic migration exist in very young magnetars. We find that the damping of the post-migration pulsations strongly depends on the character of the equation of state softening. The damping of pulsations in the models with the microphysical equation of state is caused by dissipation associated with matter flowing through the density jump at the edge of the dense core. If at work, this mechanism dominates over all other types of dissipation, like bulk viscosity in the exotic-phase core, gravitational radiation damping, or numerical viscosity., Comment: 23 pages, 18 figures, minor modifications

Coalescing compact star binaries are expected to be among the strongest sources of gravitational radiation to be seen by laser interferometers. We present calculations of the final phase of inspiral of equal mass irrotational neutron star binaries and strange quark star binaries. Six types of equations of state at zero temperature are used - three realistic nuclear equations of state of various softness and three different MIT bag models of strange quark matter. We study the precoalescing stage within the Isenberg-Wilson-Mathews approximation of general relativity using a multidomain spectral method. The gravitational-radiation driven evolution of the binary system is approximated by a sequence of quasi-equilibrium configurations at fixed baryon number and decreasing separation. We find that the innermost stable circular orbit (ISCO) is given by an orbital instability for binary strange quark stars and by the mass-shedding limit for neutron star binaries. The gravitational wave frequency at the ISCO, which marks the end of the inspiral phase, is found to be around 1100-1460 Hz for two 1.35 solar masses irrotational strange stars described by the MIT bag model and between 800 Hz and 1230 Hz for neutron stars., Comment: 6 pages, 2 figures, Advances in Space Research, in press,doi:10.1016/j.asr.2006.09.021, available online www.sciencedirect.com, paper updated according to the reviewer's suggestions (conclusions unchanged)

We present results of star variability analysis in OGLE II first bulge field. Photometric database was derived by means of image subtraction method (Wozniak 2000) and contains 4597 objects pre-classified as variables. We analyzed all the light curves in order to find periodic variables, non periodic but all-time variables and stars showing episodic changes in their brightness, e.g. gravitational lenses. Variability of 3969 stars was confirmed, among them we detected 12 lensing event candidates whose light curves are shown. We found 762 periodic variables. Algorithmic methods let us identify 71 RR Lyrae and 110 W UMa stars. Almost all luminous red giants are found to be variable. Most of red clump giants are not variable. Classification of all 4597 objects is presented online., Comment: 24 pages, 19 figures, results can be obtained at ftp://ftp.astrouw.edu.pl/pub/mizerski/