An Efficient Implementation of Flux Formulae in Multidimensional Relativistic Hydrodynamical Codes

We derive and analyze a simplified formulation of the numerical viscosity terms appearing in the expression of the numerical fluxes associated to several High-Resolution Shock-Capturing schemes. After some algebraic pre-processing, we give explicit expressions for the numerical viscosity terms of two of the most widely used flux formulae, which implementation saves computational time in multidimensional simulations of relativistic flows. Additionally, such treatment explicitely cancells and factorizes a number of terms helping to amortiguate the growing of round-off errors. We have checked the performance of our formulation running a 3D relativistic hydrodynamical code to solve a standard test-bed problem and found that the improvement in efficiency is of high practical interest in numerical simulations of relativistic flows in Astrophysics.Comment: 10 pages, accepted for publication in Computer Physics Communication

Mirror World, Supersymmetric Axion and Gamma Ray Bursts

A modification of the relation between axion mass and the PQ constant permits a relaxation of the astrophysical constraints, considerably enlarging the allowed axion parameter space. We develop this idea in this paper, discussing a model for an {\it ultramassive} axion, which essentially represents a supersymmetric Weinberg-Wilczek axion of the mirror world. The experimental and astrophysical limits allow a PQ scale f_a ~ 10^4-10^6 GeV and a mass m_a ~ 1MeV, which can be accessible for future experiments. On a phenomenological ground, such an {\it ultramassive} axion turns out to be quite interesting. It can be produced during the gravitational collapse or during the merging of two compact objects, and its subsequent decay into e+e- provides an efficient mechanism for the transfer of the gravitational energy of the collapsing system to the electron-positron plasma. This could resolve the energy budget problem in the Gamma Ray Bursts and also help in understanding the SN type II explosion phenomena.Comment: 20 pages, 5 eps figures, added footnote and reference

Does the plasma composition affect the long term evolution of relativistic jets?

We study the influence of the matter content of extragalactic jets on their morphology, dynamics and emission properties. For this purpose we consider jets of extremely different compositions including pure leptonic and baryonic plasmas. Our work is based on two-dimensional relativistic hydrodynamic simulations of the long-term evolution of powerful extragalactic jets propagating into a homogeneous environment. The equation of state used in the simulations accounts for an arbitrary mixture of electrons, protons and electron-positron pairs. Using the hydrodynamic models we have also computed synthetic radio maps and the thermal Bremsstrahlung X-ray emission from their cavities. Although there is a difference of about three orders of magnitude in the temperatures of the cavities inflated by the simulated jets, we find that both the morphology and the dynamic behaviour are almost independent on the assumed composition of the jets. Their evolution proceeds in two distinct epochs. During the first one multidimensional effects are unimportant and the jets propagate ballistically. The second epoch starts when the first larger vortices are produced near the jet head causing the beam cross section to increase and the jet to decelerate. The evolution of the cocoon and cavity is in agreement with a simple theoretical model. The beam velocities are relativistic ($\Gamma \simeq 4$) at kiloparsec scales supporting the idea that the X-ray emission of several extragalactic jets may be due to relativistically boosted CMB photons. The radio emission of all models is dominated by the contribution of the hot spots. All models exhibit a depression in the X-rays surface brightness of the cavity interior in agreement with recent observations.Comment: 41 pages, 12 png figures. See also high quality versions of the figures and several movies under http://www.mpa-garching.mpg.de/~maa/JETS/LongTermEv/index.htm

Population properties of compact objects from the second LIGO-Virgo gravitational-wave transient catalog

We report on the population of 47 compact binary mergers detected with a false-alarm rate of (BBH) population not discernible until now. First, the primary mass spectrum contains structure beyond a power law with a sharp high-mass cutoff; it is more consistent with a broken power law with a break at 39.7-+9.120.3 M? or a power law with a Gaussian feature peaking at 33.1-+5.64.0 M? (90% credible interval). While the primary mass distribution must extend to ~65 M? or beyond, only 2.9-+1.73.5% of systems have primary masses greater than 45 M?. Second, we find that a fraction of BBH systems have component spins misaligned with the orbital angular momentum, giving rise to precession of the orbital plane. Moreover,12%-44% of BBH systems have spins tilted by more than 90°, giving rise to a negative effective inspiral spin parameter, ceff. Under the assumption that such systems can only be formed by dynamical interactions, we infer that between 25% and 93% of BBHs with nonvanishing ceff| \u3e 0.01 are dynamically assembled. Third, we estimate merger rates, finding RBBH = 23.9-+8.614.3 Gpc-3 yr-1 for BBHs and RBNS = 320-+240490 Gpc-3 yr-1 for binary neutron stars. We find that the BBH rate likely increases with redshift (85% credibility) but not faster than the star formation rate (86% credibility). Additionally, we examine recent exceptional events in the context of our population models, finding that the asymmetric masses of GW190412 and the high component masses of GW190521 are consistent with our models, but the low secondary mass of GW190814 makes it an outlier

On the equivalence between the Scheduled Relaxation Jacobi method and Richardson's non-stationary method

The Scheduled Relaxation Jacobi (SRJ) method is an extension of the classical Jacobi iterative method to solve linear systems of equations (Au=b) associated with elliptic problems. It inherits its robustness and accelerates its convergence rate computing a set of P relaxation factors that result from a minimization problem. In a typical SRJ scheme, the former set of factors is employed in cycles of M consecutive iterations until a prescribed tolerance is reached. We present the analytic form for the optimal set of relaxation factors for the case in which all of them are strictly different, and find that the resulting algorithm is equivalent to a non-stationary generalized Richardson's method where the matrix of the system of equations is preconditioned multiplying it by D=diag(A). Our method to estimate the weights has the advantage that the explicit computation of the maximum and minimum eigenvalues of the matrix A (or the corresponding iteration matrix of the underlying weighted Jacobi scheme) is replaced by the (much easier) calculation of the maximum and minimum frequencies derived from a von Neumann analysis of the continuous elliptic operator. This set of weights is also the optimal one for the general problem, resulting in the fastest convergence of all possible SRJ schemes for a given grid structure. The amplification factor of the method can be found analytically and allows for the exact estimation of the number of iterations needed to achieve a desired tolerance. We also show that with the set of weights computed for the optimal SRJ scheme for a fixed cycle size it is possible to estimate numerically the optimal value of the parameter ω in the Successive Overrelaxation (SOR) method in some cases. Finally, we demonstrate with practical examples that our method also works very well for Poisson-like problems in which a high-order discretization of the Laplacian operator is employed (e.g., a 9- or 17-points discretization). This is of interest since the former discretizations do not yield consistently ordered A matrices and, hence, the theory of Young cannot be used to predict the optimal value of the SOR parameter. Furthermore, the optimal SRJ schemes deduced here are advantageous over existing SOR implementations for high-order discretizations of the Laplacian operator in as much as they do not need to resort to multi-coloring schemes for their parallel implementation

A guide to LIGO-Virgo detector noise and extraction of transient gravitational-wave signals

The LIGO Scientific Collaboration and the Virgo Collaboration have cataloged eleven confidently detected gravitational-wave events during the first two observing runs of the advanced detector era. All eleven events were consistent with being from well-modeled mergers between compact stellar mass objects: black holes or neutron stars. The data around the time of each of these events have been made publicly available through the gravitational wave open science center. The entirety of the gravitational-wave strain data from the first and second observing runs have also now been made publicly available. There is considerable interest among the broad scientific community in understanding the data and methods used in the analyses. In this paper, we provide an overview of the detector noise properties and the data analysis techniques used to detect gravitational-wave signals and infer the source properties. We describe some of the checks that are performed to validate the analyses and results from the observations of gravitational-wave events. We also address concerns that have been raised about various properties of LIGO–Virgo detector noise and the correctness of our analyses as applied to the resulting data.B P Abbott, R Abbott ... Miftar Ganija … Craig Ingram … David J. Ottaway … Peter J. Veitch ... et al

Population properties of compact objects from the second ligo-virgo gravitational-wave transient catalog

We report on the population of 47 compact binary mergers detected with a false-alarm rate of 0.01 eff are dynamically assembled. Third, we estimate merger rates, finding = - + - - BBH 23.9 Gpc yr 8.6 14.3 3 1 for BBHs and = - + - - BNS 320 Gpc yr 240 490 3 1 for binary neutron stars. We find that the BBH rate likely increases with redshift (85% credibility) but not faster than the star formation rate (86% credibility). Additionally, we examine recent exceptional events in the context of our population models, finding that the asymmetric masses of GW190412 and the high component masses of GW190521 are consistent with our models, but the low secondary mass of GW190814 makes it an outlier.R. Abbott … David Ottaway … Alexei Ciobanu … Peter Veitch … Craig Ingram … Daniel Brown … et al

GW190425: Observation of a Compact Binary Coalescence with Total Mass ∼ 3.4 M o

On 2019 April 25, the LIGO Livingston detector observed a compact binary coalescence with signal-to-noise ratio 12.9. The Virgo detector was also taking data that did not contribute to detection due to a low signal-to-noise ratio, but were used for subsequent parameter estimation. The 90% credible intervals for the component masses range from to if we restrict the dimensionless component spin magnitudes to be smaller than 0.05). These mass parameters are consistent with the individual binary components being neutron stars. However, both the source-frame chirp mass and the total mass of this system are significantly larger than those of any other known binary neutron star (BNS) system. The possibility that one or both binary components of the system are black holes cannot be ruled out from gravitational-wave data. We discuss possible origins of the system based on its inconsistency with the known Galactic BNS population. Under the assumption that the signal was produced by a BNS coalescence, the local rate of neutron star mergers is updated to 250-2810

Search for gravitational-wave signals associated with gamma-ray bursts during the second observing run of advanced LIGO and advanced Virgo

We present the results of targeted searches for gravitational-wave transients associated with gamma-ray bursts during the second observing run of Advanced LIGO and Advanced Virgo, which took place from 2016 November to 2017 August. We have analyzed 98 gamma-ray bursts using an unmodeled search method that searches for generic transient gravitational waves and 42 with a modeled search method that targets compact-binary mergers as progenitors of short gamma-ray bursts. Both methods clearly detect the previously reported binary merger signal GW170817, with p-values of <9.38 × 10−6 (modeled) and 3.1 × 10−4 (unmodeled). We do not find any significant evidence for gravitational-wave signals associated with the other gamma-ray bursts analyzed, and therefore we report lower bounds on the distance to each of these, assuming various source types and signal morphologies. Using our final modeled search results, short gamma-ray burst observations, and assuming binary neutron star progenitors, we place bounds on the rate of short gamma-ray bursts as a function of redshift for z ≤ 1. We estimate 0.07–1.80 joint detections with Fermi-GBM per year for the 2019–20 LIGO-Virgo observing run and 0.15–3.90 per year when current gravitational-wave detectors are operating at their design sensitivities

Non-motor symptom burden in patients with Parkinson's disease with impulse control disorders and compulsive behaviours : results from the COPPADIS cohort

The study was aimed at analysing the frequency of impulse control disorders (ICDs) and compulsive behaviours (CBs) in patients with Parkinson's disease (PD) and in control subjects (CS) as well as the relationship between ICDs/CBs and motor, nonmotor features and dopaminergic treatment in PD patients. Data came from COPPADIS-2015, an observational, descriptive, nationwide (Spain) study. We used the validated Questionnaire for Impulsive-Compulsive Disorders in Parkinson's Disease-Rating Scale (QUIP-RS) for ICD/CB screening. The association between demographic data and ICDs/CBs was analyzed in both groups. In PD, this relationship was evaluated using clinical features and treatment-related data. As result, 613 PD patients (mean age 62.47 ± 9.09 years, 59.87% men) and 179 CS (mean age 60.84 ± 8.33 years, 47.48% men) were included. ICDs and CBs were more frequent in PD (ICDs 12.7% vs. 1.6%, p < 0.001; CBs 7.18% vs. 1.67%, p = 0.01). PD patients had more frequent previous ICDs history, premorbid impulsive personality and antidepressant treatment (p < 0.05) compared with CS. In PD, patients with ICDs/CBs presented younger age at disease onset, more frequent history of previous ICDs and premorbid personality (p < 0.05), as well as higher comorbidity with nonmotor symptoms, including depression and poor quality of life. Treatment with dopamine agonists increased the risk of ICDs/CBs, being dose dependent (p < 0.05). As conclusions, ICDs and CBs were more frequent in patients with PD than in CS. More nonmotor symptoms were present in patients with PD who had ICDs/CBs compared with those without. Dopamine agonists have a prominent effect on ICDs/CBs, which could be influenced by dose