ABSTRACT
The magnetic field is believed to play an important role in at least some core-collapse supernovae (CCSN) if its magnitude reaches $10^{15}\, \rm {G}$, which is a typical value for a ...magnetar. In the presence of fast rotation, such a strong magnetic field can drive powerful jet-like explosions if it has the large-scale coherence of a dipole. The topology of the magnetic field is, however, probably much more complex with strong multipolar and small-scale components and the consequences for the explosion are so far unclear. We investigate the effects of the magnetic field topology on the dynamics of CCSN and the properties of the forming proto-neutron star (PNS) by comparing pre-collapse fields of different multipolar orders and radial profiles. Using axisymmetric special relativistic MHD simulations and a two-moment neutrino transport, we find that higher multipolar magnetic configurations lead to generally less energetic explosions, slower expanding shocks, and less collimated outflows. Models with a low order multipolar configuration tend to produce more oblate PNS, surrounded in some cases by a rotationally supported toroidal structure of neutron-rich material. Moreover, magnetic fields which are distributed on smaller angular scales produce more massive and faster rotating central PNS, suggesting that higher order multipolar configurations tend to decrease the efficiency of the magnetorotational launching mechanism. Even if our dipolar models systematically display a far more efficient extraction of the rotational energy of the PNS, fields distributed on smaller angular scales are still capable of powering magnetorotational explosions and shape the evolution of the central compact object.
Abstract
We study the long-term quasi-steady evolution of the force-free magnetosphere of a magnetar coupled to its internal magnetic field. We find that magnetospheric currents can be maintained on ...long time-scales of the order of thousands of years. Meanwhile, the energy, helicity and twist stored in the magnetosphere all gradually increase over the course of this evolution, until a critical point is reached, beyond which a force-free magnetosphere cannot be constructed. At this point, some large-scale magnetospheric rearrangement, possibly resulting in an outburst or a flare, must occur, releasing a large fraction of the stored energy, helicity and twist. After that, the quasi-steady evolution should continue in a similar manner from the new initial conditions. The time-scale for reaching this critical point depends on the overall magnetic field strength and on the relative fraction of the toroidal field. The energy stored in the force-free magnetosphere is found to be up to ∼30 per cent larger than the corresponding vacuum energy. This implies that for a 1014 G field at the pole, the energy budget available for fast magnetospheric events is of the order of a few 1044 erg. The spin-down rate is estimated to increase by up to ∼60 per cent, since the dipole content in the magnetosphere is enhanced by the currents present there. A rough estimate of the braking index n reveals that it is systematically n < 3 for the most part of the evolution, consistent with actual measurements for pulsars and early estimates for several magnetars.
We perform numerical evolutions of the fully nonlinear Einstein (complex, massive) Klein-Gordon and Einstein (complex) Proca systems, to assess the formation and stability of spinning bosonic stars. ...In the scalar (vector) case these are known as boson (Proca) stars. Firstly, we consider the formation scenario. Starting with constraint-obeying initial data, describing a dilute, axisymmetric cloud of spinning scalar or Proca field, gravitational collapse toward a spinning star occurs, via gravitational cooling. In the scalar case the formation is transient, even for a nonperturbed initial cloud; a nonaxisymmetric instability always develops ejecting all the angular momentum from the scalar star. In the Proca case, by contrast, no instability is observed and the evolutions are compatible with the formation of a spinning Proca star. Secondly, we address the stability of an existing star, a stationary solution of the field equations. In the scalar case, a nonaxisymmetric perturbation develops, collapsing the star to a spinning black hole. No such instability is found in the Proca case, where the star survives large amplitude perturbations; moreover, some excited Proca stars decay to, and remain as, fundamental states. Our analysis suggests bosonic stars have different stability properties in the scalar (vector) case, which we tentatively relate to its toroidal (spheroidal) morphology. A parallelism with instabilities of spinning fluid stars is briefly discussed.
The detection of gravitational waves from core-collapse supernova (CCSN) explosions is a challenging task, yet to be achieved, in which it is key the connection between multiple messengers, including ...neutrinos and electromagnetic signals. In this work, we present a method for detecting these kind of signals based on machine learning techniques. We tested its robustness by injecting signals in the real noise data taken by the Advanced LIGO-Virgo network during the second observing run, O2. We trained a newly developed Mini-Inception Resnet neural network using time-frequency images corresponding to injections of simulated phenomenological signals, which mimic the waveforms obtained in 3D numerical simulations of CCSNe. With this algorithm we were able to identify signals from both our phenomenological template bank and from actual numerical 3D simulations of CCSNe. We computed the detection efficiency versus the source distance, obtaining that, for signal to noise ratio higher than 15, the detection efficiency is 70% at a false alarm rate lower than 5%. We notice also that, in the case of the O2 run, it would have been possible to detect signals emitted at 1 kpc of distance, while lowering down the efficiency to 60%, the event distance reaches values up to 14 kpc.
While gravitational waves have been detected from mergers of binary black holes and binary neutron stars, signals from core collapse supernovae, the most energetic explosions in the modern Universe, ...have not been detected yet. Here we present a new method to analyse the data of the LIGO, Virgo, and KAGRA network to enhance the detection efficiency of this category of signals. The method takes advantage of a peculiarity of the gravitational wave signal emitted in the core collapse supernova and it is based on a classification procedure of the time-frequency images of the network data performed by a convolutional neural network trained to perform the task to recognize the signal. We validate the method using phenomenological waveforms injected in Gaussian noise whose spectral properties are those of the LIGO and Virgo advanced detectors and we conclude that this method can identify the signal better than the present algorithm devoted to select gravitational wave transient signal.
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.
► We overlap the global distribution of protected areas (PAs) and legal mining activities for four metals. ► A smaller number of mines occurred within PAs, but more mines within the proximity of PAs ...than expected by chance. ► 14% of the global protected coverage holds mines within its bounds and in close proximity. ► Our results indicate metal mines as a potential threat to global PA system.
The global protected area (PA) system has a key role to play in biological conservation, and it is thus vital to understand the factors that are likely to limit this potential. Attention to date has focused foremost on the consequences of biases in the spatial distribution of PAs for their effectiveness and efficiency in representing biodiversity. What is less clear is the extent to which these biases may also have affected the likelihood with which PAs coincide with or are influenced by particular kinds of threatening processes, further undermining their role. An obvious candidate for such concerns is metal mining activities. Here we demonstrate that approximately 7% of mines for four key metals directly overlap with PAs and a further 27% lie within 10km of a PA boundary. Moreover, those PAs with mining activity within their boundaries constitute around 6% of the total areal coverage of the global terrestrial PA system, and those with mining activity within or up to 10km from their boundary constitute nearly 14% of the total area. Given the distances over which mining activities can have influences, the persistence of their effects (often long after actual operations have closed down), and the rapidly growing demand for metals, there is an urgent need to limit or mitigate such conflicts.
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 properly 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 simulations of a region close to the surface of a differentially rotating protoneutron star in non-ideal magnetohydrodynamics 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 protoneutron stars, a more realistic value for the magnetic 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 protoneutron 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.
Abstract
The study of the electrodynamics of static, axisymmetric, and force-free Kerr magnetospheres relies vastly on solutions of the so-called relativistic Grad–Shafranov equation (GSE). Different ...numerical approaches to the solution of the GSE have been introduced in the literature, but none of them has been fully assessed from the numerical point of view in terms of efficiency and quality of the solutions found. We present a generalization of these algorithms and give a detailed background on the algorithmic implementation. We assess the numerical stability of the implemented algorithms and quantify the convergence of the presented methodology for the most established set-ups (split-monopole, paraboloidal, BH disc, uniform).