BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Pablo Cerdá-Durán

Pablo Cerdá-Durán appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
gr-qcastro-ph.SRastro-ph.HEmath-phmath.MP

Catalog footprint

What is connected

16works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

16 published item(s)

preprint2022arXiv

Impact of ultralight bosonic dark matter on the dynamical bar-mode instability of rotating neutron stars

We investigate the effects ultralight bosonic field dark matter may have on the dynamics of unstable differentially-rotating neutron stars prone to the bar-mode instability. To this aim we perform numerical simulations in general relativity of rotating neutron stars accreting an initial spherically symmetric bosonic field cloud, solving the Einstein-(complex, massive) Klein-Gordon-Euler and the Einstein-(complex) Proca-Euler systems. We find that the presence of the bosonic field can critically modify the development of the bar-mode instability of neutron stars, depending on the total mass of the bosonic field and on the boson particle mass. In some cases, the accreting bosonic field can even quench the dominant $\ell=m=2$ mode of the bar-deformation by dynamically forming a mixed (fermion-boson) star that retains part of the angular momentum of the original neutron star. However, the mixed star undergoes the development of a mixed bar that leads to significant gravitational-wave emission, substantially different to that of the isolated neutron star. Our results indicate that dark-matter accretion in neutron stars could change the frequency of the expected emission of the bar-mode instability, which would have an important impact on ongoing searches for continuous gravitational waves.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Can fermion-boson stars reconcile multi-messenger observations of compact stars?

Mixed fermion-boson stars are stable, horizonless, everywhere regular solutions of the coupled Einstein-(complex, massive) Klein-Gordon-Euler system. While isolated neutron stars and boson stars are uniquely determined by their central energy density, mixed configurations conform an extended parameter space that depends on the combination of the number of fermions and (ultra-light) bosons. The wider possibilities offered by fermion-boson stars could help explain the tension in the measurements of neutron star masses and radii reported in recent multi-messenger observations and nuclear-physics experiments. In this work we construct equilibrium configurations of mixed fermion-boson stars with realistic equations of state for the fermionic component and different percentages of bosonic matter. We show that our solutions are in excellent agreement with multi-messenger data, including gravitational-wave events GW170817 and GW190814 and X-ray pulsars PSR J0030+0451 and PSR J0740+6620, as well as with nuclear physics constraints from the PREX-2 experiment.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Coherent magneto-elastic oscillations in superfluid magnetars

We study the effect of superfluidity on torsional oscillations of highly magnetised neutron stars (magnetars) with a microphysical equation of state by means of two-dimensional, magnetohydrodynamical- elastic simulations. The superfluid properties of the neutrons in the neutron star core are treated in a parametric way in which we effectively decouple part of the core matter from the oscillations. Our simulations confirm the existence of two groups of oscillations, namely continuum oscillations that are confined to the neutron star core and are of Alfvénic character, and global oscillations with constant phase and that are of mixed magneto-elastic type. The latter might explain the quasi-periodic oscillations observed in magnetar giant flares, since they do not suffer from the additional damping mechanism due to phase mixing, contrary to what happens for continuum oscillations. However, we cannot prove rigorously that the coherent oscillations with constant phase are normal modes. Moreover, we find no crustal shear modes for the magnetic field strengths typical for magnetars.We provide fits to our numerical simulations that give the oscillation frequencies as functions of magnetic field strength and proton fraction in the core.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the maximum magnetic field amplification by the magnetorotational instability in core-collapse supernovae

Whether the magnetorotational instability (MRI) can amplify initially weak magnetic fields to dynamically relevant strengths in core collapse supernovae is still a matter of active scientific debate. Recent numerical studies have shown that the first phase of MRI growth dominated by channel flows is terminated by parasitic instabilities of the Kelvin-Helmholtz type that disrupt MRI channel flows and quench further magnetic field growth. However, it remains to be prop- erly assessed by what factor the initial magnetic field can be amplified and how it depends on the initial field strength and the amplitude of the perturbations. Different termination criteria leading to different estimates of the amplification factor were proposed within the parasitic model. To determine the amplification factor and test which criterion is a better predictor of the MRI termination, we perform three-dimensional shearing-disc and shearing-box simula- tions of a region close to the surface of a differentially rotating proto-neutron star in non-ideal MHD with two different numerical codes. We find that independently of the initial magnetic field strength, the MRI channel modes can amplify the magnetic field by, at most, a factor of 100. Under the conditions found in proto-neutron stars a more realistic value for the mag- netic field amplification is of the order of 10. This severely limits the role of the MRI channel modes as an agent amplifying the magnetic field in proto-neutron stars starting from small seed fields. A further amplification should therefore rely on other physical processes, such as for example an MRI-driven turbulent dynamo.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Partially implicit Runge-Kutta methods for wave-like equations

In this work we present a new class of Runge-Kutta (RK) methods for solving systems of hyperbolic equations with a particular structure, generalization of a wave-equation. The new methods are {\it partially implicit} in the sense that a proper subset of the equations of the system contains some terms which are treated implicitly. These methods can be viewed as a particular case of the implicit-explicit (IMEX) RK methods for systems of equations with wave-like structure. For these systems, the optimal methods with the new structure are easier to derive than the IMEX ones, specially when aiming at higher-order (up to fourth-order in this work). The methods are constructed considering the classical strong-stability-preserving optimal explicit RK methods for the purely explicit part. The resulting partially implicit RK methods do not require any inversion of operators and hence their computational cost per iteration is similar to those of explicit RK methods. We analyse the stability and convergence properties and show their practical applicability in several numerical examples. Our results show that, compared with explicit RK methods, the new methods have better stability properties (larger steps are allowed) and in general show smaller discretization error.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Termination of the magnetorotational instability via parasitic instabilities in core-collapse supernovae

The magnetorotational instability (MRI) can be a powerful mechanism amplifying the magnetic field in core-collapse supernovae. Whether initially weak magnetic fields can be amplified by this instability to dynamically relevant strengths is still a matter of debate. One of the main uncertainties concerns the process that terminates the growth of the instability. Parasitic instabilities of both Kelvin-Helmholtz and tearing-mode type have been suggested to play a crucial role in this process, disrupting MRI channel flows and quenching magnetic field amplification. We perform two-dimensional and three-dimensional sheering-disc simulations of a differentially rotating protoneutron star layer in non-ideal magnetohydrodynamics with unprecedented high numerical accuracy, finding that Kelvin-Helmholtz parasitic modes dominate tearing modes in the regime of large hydrodynamic and magnetic Reynolds numbers, as encountered close to the surface of protoneutron stars. They also determine the maximum magnetic field stress achievable during the exponential growth of the MRI. Our results are consistent with the theory of parasitic instabilities based on a local stability analysis. To simulate the Kelvin-Helmholtz instabilities properly, a very high numerical resolution is necessary. Using ninth-order spatial reconstruction schemes, we find that at least eight grid zones per MRI channel are necessary to simulate the growth phase of the MRI and reach an accuracy of ~10 per cent in the growth rate, while more than ~60 zones per channel are required to achieve convergent results for the value of the magnetic stress at MRI termination.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Are pulsars born with a hidden magnetic field?

The observation of several neutron stars in the center of supernova remnants and with significantly lower values of the dipolar magnetic field than the average radio-pulsar population has motivated a lively debate about their formation and origin, with controversial interpretations. A possible explanation requires the slow rotation of the proto-neutron star at birth, which is unable to amplify its magnetic field to typical pulsar levels. An alternative possibility, the hidden magnetic field scenario, considers the accretion of the fallback of the supernova debris onto the neutron star as responsible for the submergence (or screening) of the field and its apparently low value. In this paper we study under which conditions the magnetic field of a neutron star can be buried into the crust due to an accreting, conducting fluid. For this purpose, we consider a spherically symmetric calculation in general relativity to estimate the balance between the incoming accretion flow and the magnetosphere. Our study analyses several models with different specific entropy, composition, and neutron star masses. The main conclusion of our work is that typical magnetic fields of a few times 1e12 G can be buried by accreting only 1e-3 - 1e-2 solar mass, a relatively modest amount of mass. In view of this result, the Central Compact Object scenario should not be considered unusual, and we predict that anomalously weak magnetic fields should be common in very young (< few kyr) neutron stars.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Modulating the magnetosphere of magnetars by internal magneto-elastic oscillations

We couple internal torsional, magneto-elastic oscillations of highly magnetized neutron stars (magnetars) to their magnetospheres. The corresponding axisymmetric perturbations of the external magnetic field configuration evolve as a sequence of linear, force-free equilibria that are completely determined by the background magnetic field configuration and by the perturbations of the magnetic field at the surface. The perturbations are obtained from simulations of magneto-elastic oscillations in the interior of the magnetar. While such oscillations can excite travelling Alfvén waves in the exterior of the star only in a very limited region close to the poles, they still modulate the near magnetosphere by inducing a time-dependent twist between the foot-points of closed magnetic field lines that exit the star at a polar angle $\gtrsim 0.19\,$rad. Moreover, we find that for a dipole-like background magnetic field configuration the magnetic field modulations in the magnetosphere, driven by internal oscillations, can only be symmetric with respect to the equator. This is in agreement with our previous findings, where we interpreted the observed quasi-periodic oscillations in the X-ray tail of magnetar bursts as driven by the family of internal magneto-elastic oscillations with symmetric magnetic field perturbations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Gravitational wave signatures in black-hole-forming core collapse

We present numerical simulations in general relativity of collapsing stellar cores. Our initial model consists of a low metallicity rapidly-rotating progenitor which is evolved in axisymmetry with the latest version of our general relativistic code CoCoNuT, which allows for black hole formation and includes the effects of a microphysical equation of state (LS220) and a neutrino leakage scheme to account for radiative losses. The motivation of our study is to analyze in detail the emission of gravitational waves in the collapsar scenario of long gamma-ray bursts. Our simulations show that the phase during which the proto-neutron star (PNS) survives before ultimately collapsing to a black hole is particularly optimal for gravitational wave emission. The high-amplitude waves last for several seconds and show a remarkable quasi-periodicity associated with the violent PNS dynamics, namely during the episodes of convection and the subsequent nonlinear development of the standing-accretion shock instability (SASI). By analyzing the spectrogram of our simulations we are able to identify the frequencies associated with the presence of g-modes and with the SASI motions at the PNS surface. We note that the gravitational waves emitted reach large enough amplitudes to be detected with third-generation detectors as the Einstein Telescope within a Virgo cluster volume at rates ~< 0.1 /y.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Imprints of superfluidity on magneto-elastic QPOs of SGRs

Our numerical simulations show that axisymmetric, torsional, magneto-elastic oscillations of magnetars with a superfluid core can explain the whole range of observed quasi-periodic oscillations (QPOs) in the giant flares of soft gamma-ray repeaters. There exist constant phase, magneto-elastic QPOs at both low (f<150 Hz) and high frequencies (f>500 Hz), in full agreement with observations. The range of magnetic field strengths required to match the observed QPO frequencies agrees with that from spin-down estimates. These results strongly suggest that neutrons in magnetar cores are superfluid.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Magneto-elastic oscillations of neutron stars: exploring different magnetic field configurations

We study magneto-elastic oscillations of highly magnetized neutron stars (magnetars) which have been proposed as an explanation for the quasi-periodic oscillations (QPOs) appearing in the decaying tail of the giant flares of soft gamma-ray repeaters (SGRs). We extend previous studies by investigating various magnetic field configurations, computing the Alfvén spectrum in each case and performing magneto-elastic simulations for a selected number of models. By identifying the observed frequencies of 28 Hz (SGR 1900+14) and 30 Hz (SGR 1806-20) with the fundamental Alfvén QPOs, we estimate the required surface magnetic field strength. For the magnetic field configurations investigated (dipole-like poloidal, mixed toroidal-poloidal with a dipole-like poloidal component and a toroidal field confined to the region of field lines closing inside the star, and for poloidal fields with an additional quadrupole-like component) the estimated dipole spin-down magnetic fields are between 8x10^14 G and 4x10^15 G, in broad agreement with spin-down estimates for the SGR sources producing giant flares. A number of these models exhibit a rich Alfvén continuum revealing new turning points which can produce QPOs. This allows one to explain most of the observed QPO frequencies as associated with magneto-elastic QPOs. In particular, we construct a possible configuration with two turning points in the spectrum which can explain all observed QPOs of SGR 1900+14. Finally, we find that magnetic field configurations which are entirely confined in the crust (if the core is assumed to be a type I superconductor) are not favoured, due to difficulties in explaining the lowest observed QPO frequencies (f<30 Hz).

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Gravitational waves in dynamical spacetimes with matter content in the Fully Constrained Formulation

The Fully Constrained Formulation (FCF) of General Relativity is a novel framework introduced as an alternative to the hyperbolic formulations traditionally used in numerical relativity. The FCF equations form a hybrid elliptic-hyperbolic system of equations including explicitly the constraints. We present an implicit-explicit numerical algorithm to solve the hyperbolic part, whereas the elliptic sector shares the form and properties with the well known Conformally Flat Condition (CFC) approximation. We show the stability andconvergence properties of the numerical scheme with numerical simulations of vacuum solutions. We have performed the first numerical evolutions of the coupled system of hydrodynamics and Einstein equations within FCF. As a proof of principle of the viability of the formalism, we present 2D axisymmetric simulations of an oscillating neutron star. In order to simplify the analysis we have neglected the back-reaction of the gravitational waves into the dynamics, which is small (<2 %) for the system considered in this work. We use spherical coordinates grids which are well adapted for simulations of stars and allow for extended grids that marginally reach the wave zone. We have extracted the gravitational wave signature and compared to the Newtonian quadrupole and hexadecapole formulae. Both extraction methods show agreement within the numerical errors and the approximations used (~30 %).

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Dynamical spacetimes and gravitational radiation in a Fully Constrained Formulation

This contribution summarizes the recent work carried out to analyze the behavior of the hyperbolic sector of the Fully Constrained Formulation (FCF) derived in Bonazzola et al. 2004. The numerical experiments presented here allows one to be confident in the performances of the upgraded version of CoCoNuT's code by replacing the Conformally Flat Condition (CFC) approximation of the Einstein equations by the FCF.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Magneto-elastic torsional oscillations of magnetars

We extend a general-relativistic ideal magneto-hydrodynamical code to include the effects of elasticity. Using this numerical tool we analyse the magneto-elastic oscillations of highly magnetised neutron stars (magnetars). In simulations without magnetic field we are able to recover the purely crustal shear oscillations within an accuracy of about a few per cent. For dipole magnetic fields between 5 x 10^13 and 10^15 G the Alfvén oscillations become modified substantially by the presence of the crust. Those quasi-periodic oscillations (QPOs) split into three families: Lower QPOs near the equator, Edge QPOs related to the last open field line and Upper QPOs at larger distance from the equator. Edge QPOs are called so because they are related to an edge in the corresponding Alfvén continuum. The Upper QPOs are of the same kind, while the Lower QPOs are turning-point QPOs, related to a turning point in the continuous spectrum.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Numerical viscosity in hydrodynamics simulations in general relativity

We present an alternative method to estimate the numerical viscosity in simulations of astrophysical objects, which is based in the damping of fluid oscillations. We apply the method to general relativistic hydrodynamic simulations using a spherical coordinates. We perform 1D-spherical and 2D- axisymmetric simulations of radial oscillations in spherical systems. We calibrate first the method with simulations with added bulk viscosity and study the differences between different numerical schemes. We apply the method to radial oscillations of neutron stars and we conclude that the main source of numerical viscosity in this case is the surface of the star. We expect that this method could be useful to compute the resolution requirements and limitations of the numerical simulations in different astrophysical scenarios in the future.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Alfven QPOs in magnetars in the anelastic approximation

We perform two-dimensional simulations of Alfven oscillations in magnetars, modeled as relativistic stars with a dipolar magnetic field. We use the anelastic approximation to general relativistic magnetohydrodynamics, which allows for an effective suppression of fluid modes and an accurate description of Alfven waves. In addition, we compute Alfven oscillation frequencies along individual magnetic field lines with a semi-analytic approach, employing a short-wavelength approximation. Our main findings are as follows: a) we confirm the existence of two families of quasi-periodic oscillations (QPOs), with harmonics at integer multiples of the fundamental frequency, as was found in the linear study of Sotani, Kokkotas & Stergioulas; b) the QPOs appearing near the magnetic axis are split into two groups, depending on their symmetry across the equatorial plane. The antisymmetric QPOs have only odd integer-multiple harmonics; c) the continuum obtained with our semi-analytic approach agrees remarkably well with QPOs obtained via the two-dimensional simulations, allowing for a clear interpretation of the QPOs as corresponding to turning points of the continuum. This agreement will allow for a comprehensive study of Alfven QPOs for a larger number of different models, without the need for time-consuming simulations. Finally, we construct empirical relations for the QPO frequencies and compare them to observations of known Soft Gamma Repeaters. We find that, under the assumptions of our model and if the magnetic field of magnetars is characterized by a strong dipolar component, and QPOs are produced near the magnetic pole, then one can place an upper limit to the mean surface strength of the magnetic field of about 3-8 10^15 G.

0 citations0 reviewsNo code linkedOpen access