Showing 1–50 of 62 results for author: Gourgouliatos, K N

A new solution of the pulsar equation

Authors: Ioannis Contopoulos, Ioannis Dimitropoulos, Dimitris Ntotsikas, Konstantinos N. Gourgouliatos

Abstract: We present the first new type of solution of the pulsar equation since 1999. In it, the whole magnetosphere is confined inside the light cylinder and an electrically charged layer wraps around it and holds it together. The reason this new solution has never been obtained before is that all current time-dependent simulations are initialized with a vacuum dipole configuration that extends to infinit… ▽ More We present the first new type of solution of the pulsar equation since 1999. In it, the whole magnetosphere is confined inside the light cylinder and an electrically charged layer wraps around it and holds it together. The reason this new solution has never been obtained before is that all current time-dependent simulations are initialized with a vacuum dipole configuration that extends to infinity, thus their final steady-state solution also extends to infinity. Under special conditions, such a confined configuration may be attained when the neutron star first forms in the interior of a collapsing star during a supernova explosion, or when it accretes from an external wind or disk from a donor star. It is shown that this new maximally closed non-decelerating solution is the limit of a continuous sequence of standard magnetospheres with open and closed field lines when the amount of open field lines gradually drops to zero. The minimum energy solution in this sequence is a standard magnetosphere in which the closed field line region extends up to about 80% of the light cylinder. We estimate that the released energy when the new solution transitions to the minimum energy one is enough to power a fast radio burst. △ Less

Submitted 11 April, 2024; v1 submitted 8 April, 2024; originally announced April 2024.

Comments: 15 pages, 6 figures, in press in Journal Universe Special Volume (MDPI) "A New Horizon of Pulsar and Neutron Star: The 55-Year Anniversary" (Eds. C. Zhang, D. Blaschke, J. Petri)

Revisiting thermoelectric effects in the crust of neutron stars

Authors: Dionysios Gakis, Konstantinos N. Gourgouliatos

Abstract: Context. Large thermal variations have been observed in neutron stars that typically are not aligned with density gradients. Such terms may activate the Biermann battery effect, leading to thermo-electric interactions and the generation of electromotive force. Aims. We aim to identify the impact a temperature anisotropy on a neutron star's crust can have in the evolution of its magnetic field, thr… ▽ More Context. Large thermal variations have been observed in neutron stars that typically are not aligned with density gradients. Such terms may activate the Biermann battery effect, leading to thermo-electric interactions and the generation of electromotive force. Aims. We aim to identify the impact a temperature anisotropy on a neutron star's crust can have in the evolution of its magnetic field, through the thermo-electric terms. Methods. We consider a neutron star crust with large temperature gradients, associated with long-lived hot spots, described by a Gaussian-type function localized. We simulate the interplay between the battery term and the Hall and Ohmic evolution numerically, for axisymmetric systems. Results. The results indicate that for crust temperatures of $\sim$$10^9$ K the toroidal field can be amplified up to $\sim$$10^{14}$ - $10^{15}$ G near the points of maximum temperature gradients, and it locally changes the architecture of the poloidal field lines. For internal crustal temperatures around $\sim$$10^8$ K, the temperature gradient generates fields of about two orders of magnitude lower. In such cases, saturation is achieved after some hundred thousand years, after which the battery and Ohmic dissipation balance each other. Conclusions. We conclude that the thermoelectric effect can impact the overall magnetic field evolution, provided that the thermal gradient is maintained for a sufficiently long time. Neutron stars endowned with moderate strength magnetic fields may be affected by the thermoelectric effect if the hotspots survive for timescales of a few kiloyears. △ Less

Submitted 22 February, 2024; originally announced February 2024.

Comments: 13 pages, 11 figures, submitted to A&A

Combined magnetic field evolution in neutron star cores and crusts: Ambipolar diffusion, Hall effect and Ohmic dissipation

Authors: Dimitrios Skiathas, Konstantinos N. Gourgouliatos

Abstract: Neutron star magnetic field evolution is mediated through the Hall effect and Ohmic dissipation in the crust while ambipolar diffusion is taking place in the core. These effects have been studied in detail in either part of the star, however, their combined, simultaneous evolution and interplay has not been explored in detail yet. Here, we present simulation results of the simultaneous evolution o… ▽ More Neutron star magnetic field evolution is mediated through the Hall effect and Ohmic dissipation in the crust while ambipolar diffusion is taking place in the core. These effects have been studied in detail in either part of the star, however, their combined, simultaneous evolution and interplay has not been explored in detail yet. Here, we present simulation results of the simultaneous evolution of the magnetic field in the core due to ambipolar diffusion and the crust due to Hall effect and Ohmic decay, under the assumption of axial symmetry. We find that a purely poloidal field generates a toroidal field in the crust, due to the Hall effect, that sinks into the core. A purely toroidal field remains toroidal and spreads into the core and the crust. Finally, for a mixed poloidal-toroidal field, the north-south symmetry is broken due to the Hall effect in the crust, however, ambipolar diffusion, tends to restore it. We examine the role of ambipolar diffusion to the magnetic field decay and we compare the rate of the conversion of magnetic field energy into heat, finding that it enhances the magnetic field decay in neutron stars. △ Less

Submitted 26 February, 2024; v1 submitted 17 January, 2024; originally announced January 2024.

Comments: 11 pages, 8 figures, accepted by MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 528, Issue 3, March 2024, Pages 5178-5188

Twisted magnetar magnetospheres

Authors: D. Ntotsikas, K. N. Gourgouliatos, I. Contopoulos, S. K. Lander

Abstract: Magnetar magnetospheres are strongly twisted, and are able to power sudden energetic events through the rapid release of stored electromagnetic energy. In this paper, we investigate twisted relativistic force-free axisymmetric magnetospheres of rotating neutron stars. We obtain numerical solutions of such configurations using the method of simultaneous relaxation for the magnetic field inside and… ▽ More Magnetar magnetospheres are strongly twisted, and are able to power sudden energetic events through the rapid release of stored electromagnetic energy. In this paper, we investigate twisted relativistic force-free axisymmetric magnetospheres of rotating neutron stars. We obtain numerical solutions of such configurations using the method of simultaneous relaxation for the magnetic field inside and outside the light-cylinder. We introduce a toroidal magnetic field in the region of closed field-lines that is associated with a poloidal electric current distribution in that region, and explore various mathematical expressions for that distribution. We find that, by increasing the twist, a larger fraction of magnetic field-lines crosses the light-cylinder and opens up to infinity, thus increasing the size of the polar caps and enhancing the spin-down rate. We also find that, for moderately to strongly twisted magnetospheres, the region of closed field-lines ends at some distance inside the light-cylinder. We discuss the implications of these solutions on the variation of magnetar spin-down rates, moding and nulling of pulsars, the relation between the angular shear and the twist and the overall shape of the magnetosphere. △ Less

Submitted 14 November, 2023; originally announced November 2023.

Comments: Accepted by MNRAS

On the pulsar Y-point

Authors: I. Contopoulos, D. Ntotsikas, K. N. Gourgouliatos

Abstract: The pulsar magnetosphere is divided into a corotating region of closed field lines surrounded by open field lines that emanate from the two poles of the star, extend to infinity and are separated by an equatorial current sheet. The three regions meet at a magnetospheric Y-point. In steady-state solutions of the ideal force-free magnetosphere, the Y-point may lie at any distance inside the light cy… ▽ More The pulsar magnetosphere is divided into a corotating region of closed field lines surrounded by open field lines that emanate from the two poles of the star, extend to infinity and are separated by an equatorial current sheet. The three regions meet at a magnetospheric Y-point. In steady-state solutions of the ideal force-free magnetosphere, the Y-point may lie at any distance inside the light cylinder. Time-dependent force-free simulations, however, develop closed-line regions that extend all the way to the light cylinder. On the other hand, particle (PIC) solutions consistently develop smaller closed-line regions. In order to understand this effect, we solve the pulsar equation with an improved numerical method. We show that the total electromagnetic energy stored in the ideal force-free magnetosphere manifests a subtle minimum when the closed-line region extends to only 90% of the light cylinder, and thus argue that the system will spontaneously choose this particular configuration. Furthermore, we argue that the intersection of the corotating region with the equatorial current sheet is at right angles, literally leading to a T-point. △ Less

Submitted 19 September, 2023; originally announced September 2023.

Constraining supermassive primordial black holes with magnetically induced gravitational waves

Authors: Theodoros Papanikolaou, Konstantinos N. Gourgouliatos

Abstract: Primordial black holes (PBHs) can answer a plethora of cosmic conundra, among which the origin of the cosmic magnetic fields. In particular, supermassive PBHs with masses $M_\mathrm{PBH}>10^{10} M_\odot$ and furnished with a plasma-disk moving around them can generate through the Biermann battery mechanism a seed primordial magnetic field which can later be amplified so as to provide the magnetic… ▽ More Primordial black holes (PBHs) can answer a plethora of cosmic conundra, among which the origin of the cosmic magnetic fields. In particular, supermassive PBHs with masses $M_\mathrm{PBH}>10^{10} M_\odot$ and furnished with a plasma-disk moving around them can generate through the Biermann battery mechanism a seed primordial magnetic field which can later be amplified so as to provide the magnetic field threading the intergalactic medium. In this work, we derive the gravitational wave (GW) signal induced by the magnetic anisotropic stress of such a population of magnetised PBHs. Interestingly enough, by using GW constraints from Big Bang Nucleosynthesis (BBN) and an effective model for the galactic/turbulent dynamo amplification of the magnetic field, we set a conservative upper bound constraint on the abundances of supermassive PBHs at formation time, $Ω_\mathrm{PBH,f}$ as a function of the their masses, namely that $Ω_\mathrm{PBH,f}\leq 2.5\times 10^{-10}\left(\frac{M}{10^{10}M_\odot}\right)^{45/22}$. Remarkably, these constraints are comparable, and, in some mass ranges, even tighter compared to the constraints on $Ω_\mathrm{PBH,f}$ from large-scale structure (LSS) probes; hence promoting the portal of magnetically induced GWs as a new probe to explore the enigmatic nature of supermassive PBHs. △ Less

Submitted 12 September, 2023; v1 submitted 8 June, 2023; originally announced June 2023.

Comments: Accepted in Physical Review D

The rocket effect mechanism in neutron stars in supernova remnants

Authors: V. Agalianou, K. N. Gourgouliatos

Abstract: While the dipole magnetic field axis of neutron stars is usually postulated to cross the star's centre, it may be displaced from this location, as it has been recently indicated in the millisecond pulsar J0030+0451. Under these conditions, the electromagnetic rocket effect may be activated, where the magnetic field exerts a net force, accelerating the star. This post-natal kick mechanism relies on… ▽ More While the dipole magnetic field axis of neutron stars is usually postulated to cross the star's centre, it may be displaced from this location, as it has been recently indicated in the millisecond pulsar J0030+0451. Under these conditions, the electromagnetic rocket effect may be activated, where the magnetic field exerts a net force, accelerating the star. This post-natal kick mechanism relies on asymmetric electromagnetic radiation from an off-centre dipole may be relevant to the high spatial velocities of pulsars $\sim 10^{3}$ $\rm{km/s}$. Here, we explore its impact in young pulsars associated with supernova remnants and we compare the observational data on characteristic quantities, such as the braking index and proper motion, with results obtained from the rocket effect. Using a Markov Chain Monte Carlo analysis, we explore the required conditions, for the initial spin periods and the distance between the magnetic axis and the star's center, so that the velocity kick due to the rocket effect approaches the present velocity. We find that the electromagnetic rocket effect can account for typical pulsar transverse velocities assuming an initial spin period of 3.8 $\rm{ms}$ and a dipole field whose distance from the centre of the star is approximately 7 $\rm{km}$. We also explore the influence of the rocket effect on the braking index of a neutron star, and we find that for the sample studied this impact is minimal. Finally, we apply the rocket effect model on the pulsars J0030+0451 and J0538+2817, which are likely candidates for this mechanism. △ Less

Submitted 28 April, 2023; originally announced May 2023.

Primordial magnetic field generation via primordial black hole disks

Authors: Theodoros Papanikolaou, Konstantinos N. Gourgouliatos

Abstract: Large scale primordial magnetic fields (PMFs) threading the intergalactic medium are observed ubiquitously in the Universe playing a key role in the cosmic evolution. Their origin is still debated constituting a very active field of research. In the present article, we propose a novel natural ab initio mechanism for the origin of such PMFs through the portal of supermassive primordial black holes… ▽ More Large scale primordial magnetic fields (PMFs) threading the intergalactic medium are observed ubiquitously in the Universe playing a key role in the cosmic evolution. Their origin is still debated constituting a very active field of research. In the present article, we propose a novel natural ab initio mechanism for the origin of such PMFs through the portal of supermassive primordial black holes (PBHs) forming between the Big Bang Nucleosynthesis and the recombination era. In particular, by considering PBHs furnished with a locally isothermal disk we study the generation of a Biermann battery induced seed magnetic field (MF) due to the vortexlike motion of the primordial plasma around the black hole. Finally, by considering monochromatic PBH mass distributions and deriving the relevant MF power spectrum we make a conservative estimate for the seed PMF in intergalactic scales and at redshift $z=30$, when typical galaxies are considered to form, which reads as $B\simeq 10^{-30}\mathrm{G}\left(\frac{\ell_\mathrm{R}}{10^6}\right)^2\left(\frac{M_\mathrm{PBH}}{10^{14}M_\odot}\right)^{5/2}$, where $M_\mathrm{PBH}$ is the PBH mass and $\ell_\mathrm{R}\equiv R_\mathrm{d}/R_\mathrm{ISCO}$, is the ratio of the radius of the disk, $R_\mathrm{d}$ over the radius of the innermost stable circular orbit, $R_\mathrm{ISCO}$. Interestingly enough, by requiring to seed a PMF of the order of $10^{-30}\mathrm{G}$ necessary to give rise to a present day $10^{-18}\mathrm{G}$ in intergalactic scales, we find a lower bound on the PBH mass within the range $[10^{10}- 10^{16}]M_\odot$ depending on the radius of the PBH disk. △ Less

Submitted 16 May, 2023; v1 submitted 24 January, 2023; originally announced January 2023.

Comments: 1)To appear in Physical Review D. 2)Major revisions regarding the derivation of the magnetic field power spectrum

Orbital analysis of the Pluto-Charon moon system's mutual interactions and forced frequencies

Authors: Dionysios Gakis, Konstantinos N. Gourgouliatos

Abstract: Context. The orbits of the four small moons in the Pluto-Charon system, Styx, Nix, Kerberos and Hydra, are circumbinary, as the former form a binary dwarf planet. Consequently, the orbit of each one of them is characterized by a number of frequencies, arising by the central binary and the mutual gravitational interactions. Aims. In this work, we identify the most prominent of these forced frequenc… ▽ More Context. The orbits of the four small moons in the Pluto-Charon system, Styx, Nix, Kerberos and Hydra, are circumbinary, as the former form a binary dwarf planet. Consequently, the orbit of each one of them is characterized by a number of frequencies, arising by the central binary and the mutual gravitational interactions. Aims. In this work, we identify the most prominent of these forced frequencies using Fast Fourier Transformations. Methods. Two methods are implemented, a semi-analytic and a numerical one, and comparisons are being made. Results. The results indicate that as a first approximation, moon orbits may well be modelled as the superposition of a series of inevitable oscillations, induced by Pluto and Charon, deviating from circular ones, even if the eccentricity is set to zero. Moreover, the mutual gravitational effects are significant in their long term evolution, especially for the lighter moons Styx and Kerberos, activating modes that dominate the low-frequency region of the power spectrum. This becomes evident through the comparison of simulations where only one moon is included along with the binary dwarf planet and simulations of the entire six-body system. These modes become noticeable over long integration times and may affect the orbits of the lighter moons of the system. △ Less

Submitted 22 February, 2023; v1 submitted 5 January, 2023; originally announced January 2023.

Comments: 11 pages, 7 figures, accepted for publication in A&A

Journal ref: A&A 670, A152 (2023)

Orbits of the TOI-1338 and TIC-172900988 systems

Authors: Dionysios Gakis, Konstantinos N. Gourgouliatos

Abstract: Recent observations by TESS revealed the existence of circumbinary planets in the systems of TOI-1338 and TIC-172900988. The purpose of this work is to model the planetary orbits in these two systems and study them under the perspective of previous theoretical models. Each planet's distance from the barycenter through time is simulated using n-body integrations and is compared with outcomes from a… ▽ More Recent observations by TESS revealed the existence of circumbinary planets in the systems of TOI-1338 and TIC-172900988. The purpose of this work is to model the planetary orbits in these two systems and study them under the perspective of previous theoretical models. Each planet's distance from the barycenter through time is simulated using n-body integrations and is compared with outcomes from a semi-analytic, a geometric and a Keplerian-based approach. Furthermore, we infer the most prominent frequencies of both planets' orbits induced by the central binaries. We confirm that both systems appear to be stable. Lastly, we examine the implications of an additional candidate planet in TOI-1338 system finding that an extra, 48 M$_\oplus$ planet that has been hinted from observations could be located at 0.8 AU without generating any radical changes to the orbits of the other members of the system. △ Less

Submitted 20 December, 2022; v1 submitted 9 December, 2022; originally announced December 2022.

Comments: 11 pages, 9 figures, accepted for publication in MNRAS

Dynamical Friction due to fuzzy dark matter on satellites described by axisymmetric logarithmic potentials

Authors: Andreas Vitsos, Konstantinos N. Gourgouliatos

Abstract: A plausible dark matter candidate is an ultralight bosonic particle referred to as fuzzy dark matter. The equivalent mass-energy of the fuzzy dark matter boson is $\sim 10^{-22}$eV and has a corresponding de Broglie wavelength of kiloparsec scale, thus exhibiting wave behaviour in scales comparable to a galactic core, which could not appear in conventional cold dark matter models. The presence of… ▽ More A plausible dark matter candidate is an ultralight bosonic particle referred to as fuzzy dark matter. The equivalent mass-energy of the fuzzy dark matter boson is $\sim 10^{-22}$eV and has a corresponding de Broglie wavelength of kiloparsec scale, thus exhibiting wave behaviour in scales comparable to a galactic core, which could not appear in conventional cold dark matter models. The presence of fuzzy dark matter in galactic clusters will impact the motion of their members through dynamical friction. In this work, we present simulations of the dynamical friction on satellites traversing an initially uniform fuzzy dark matter halo. We focus on the satellites whose shapes are beyond spherical symmetry described by ellipsoidal and logarithmic potentials. We find that the wakes created on the fuzzy dark matter halo due to the passage of such satellites are qualitatively different from those generated by spherically symmetric ones. Furthermore, we quantify the dynamical friction coefficient for such systems, finding that the same satellite may experience a drag differing by a factor of $5$ depending on its ellipticity and the direction of motion. Finally, we find that the dynamical friction time-scale is close to Hubble time, assuming a satellite of $10^{11}$M$_{\odot}$ traversing at $10^{3}$km/s a FDM halo whose mean density is $\sim 10^6$M$_{\odot}$kpc$^{-3}$. △ Less

Submitted 10 April, 2023; v1 submitted 12 November, 2022; originally announced November 2022.

Comments: 17 pages, 7 figures, Accepted by JCAP

Journal ref: JCAP04(2023)071

The Laser Interferometer Space Antenna mission in Greece White Paper

Authors: Nikolaos Karnesis, Nikolaos Stergioulas, Georgios Pappas, Charis Anastopoulos, John Antoniadis, Theocharis Apostolatos, Spyros Basilakos, Kyriakos Destounis, Areti Eleni, Georgios Lukes-Gerakopoulos, Konstantinos N. Gourgouliatos, Kostas D. Kokkotas, George Kottaras, V K Oikonomou, Theodoros Papanikolaou, Leandros Perivolaropoulos, Manolis Plionis, Emmanuel N. Saridakis, Theodoros Sarris, Elias C. Vagenas, Wolf von Klitzing

Abstract: The Laser Interferometer Space Antenna (LISA) mission, scheduled for launch in the mid-2030s, is a gravitational wave observatory in space designed to detect sources emitting in the millihertz band. LISA is an ESA flagship mission, currently entering the Phase B development phase. It is expected to help us improve our understanding about our Universe by measuring gravitational wave sources of diff… ▽ More The Laser Interferometer Space Antenna (LISA) mission, scheduled for launch in the mid-2030s, is a gravitational wave observatory in space designed to detect sources emitting in the millihertz band. LISA is an ESA flagship mission, currently entering the Phase B development phase. It is expected to help us improve our understanding about our Universe by measuring gravitational wave sources of different types, with some of the sources being at very high redshifts $z\sim 20$. On the 23rd of February 2022 we organized the 1$^\mathrm{st}$ {\it LISA in Greece Workshop}. This workshop aimed to inform the Greek scientific and tech industry community about the possibilities of participating in LISA science and LISA mission, with the support of the Hellenic Space Center (HSC). In this white paper, we summarize the outcome of the workshop, the most important aspect of it being the inclusion of $15$ Greek researchers to the LISA Consortium, raising our total number to $22$. At the same time, we present a road-map with the future steps and actions of the Greek Gravitational Wave community with respect to the future LISA mission. △ Less

Submitted 21 September, 2023; v1 submitted 9 September, 2022; originally announced September 2022.

Comments: 11 pages, 2 figures

Horizons: Nuclear Astrophysics in the 2020s and Beyond

Authors: H. Schatz, A. D. Becerril Reyes, A. Best, E. F. Brown, K. Chatziioannou, K. A. Chipps, C. M. Deibel, R. Ezzeddine, D. K. Galloway, C. J. Hansen, F. Herwig, A. P. Ji, M. Lugaro, Z. Meisel, D. Norman, J. S. Read, L. F. Roberts, A. Spyrou, I. Tews, F. X. Timmes, C. Travaglio, N. Vassh, C. Abia, P. Adsley, S. Agarwal , et al. (140 additional authors not shown)

Abstract: Nuclear Astrophysics is a field at the intersection of nuclear physics and astrophysics, which seeks to understand the nuclear engines of astronomical objects and the origin of the chemical elements. This white paper summarizes progress and status of the field, the new open questions that have emerged, and the tremendous scientific opportunities that have opened up with major advances in capabilit… ▽ More Nuclear Astrophysics is a field at the intersection of nuclear physics and astrophysics, which seeks to understand the nuclear engines of astronomical objects and the origin of the chemical elements. This white paper summarizes progress and status of the field, the new open questions that have emerged, and the tremendous scientific opportunities that have opened up with major advances in capabilities across an ever growing number of disciplines and subfields that need to be integrated. We take a holistic view of the field discussing the unique challenges and opportunities in nuclear astrophysics in regards to science, diversity, education, and the interdisciplinarity and breadth of the field. Clearly nuclear astrophysics is a dynamic field with a bright future that is entering a new era of discovery opportunities. △ Less

Submitted 16 May, 2022; originally announced May 2022.

Comments: 96 pages. Submitted to Journal of Physics G

Report number: LA-UR-22-23997

Initial periods and magnetic fields of neutron stars

Authors: Andrei P. Igoshev, Anastasia Frantsuzova, Konstantinos N. Gourgouliatos, Savina Tsichli, Lydia Konstantinou, Sergei B. Popov

Abstract: Initial distributions of pulsar periods and magnetic fields are essential components of multiple modern astrophysical models. Not enough work has been done to properly constrain these distributions using direct measurements. Here we aim to fill this gap by rigorously analysing properties of young neutron stars associated to supernova remnants. In order to perform this task, we compile a catalogue… ▽ More Initial distributions of pulsar periods and magnetic fields are essential components of multiple modern astrophysical models. Not enough work has been done to properly constrain these distributions using direct measurements. Here we aim to fill this gap by rigorously analysing properties of young neutron stars associated to supernova remnants. In order to perform this task, we compile a catalogue of 56 neutron stars uniquely paired to supernova remnants with known age estimate. Further, we analyse this catalogue using multiple statistical techniques. We found that distribution of magnetic fields and periods for radio pulsars are both well described using the log-normal distribution. The mean magnetic field is $\log_{10} [B/\mathrm{G}] = 12.44$ and standard deviation is $σ_B = 0.44$. Magnetars and central compact objects do not follow the same distribution. The mean initial period is $\log_{10} P_0 [P / \mathrm{s}] = -1.04_{-0.2}^{+0.15}$ and standard deviation is $σ_p = 0.53_{-0.08}^{+0.12}$. We show that the normal distribution does not describe the initial periods of neutron stars sufficiently well. Parameters of the initial period distribution are not sensitive to the exact value of the braking index. △ Less

Submitted 13 May, 2022; originally announced May 2022.

Comments: Submitted to MNRAS, 15 pages, 21 figures

Relativistic magnetic explosions

Authors: Maxim V. Barkov, Praveen Sharma, Konstantinos N. Gourgouliatos, Maxim Lyutikov

Abstract: Many explosive astrophysical events, like magnetars' bursts and flares, are magnetically driven. We consider dynamics of such magnetic explosions - relativistic expansion of highly magnetized and highly magnetically over-pressurized clouds. The corresponding dynamics is qualitatively different from fluid explosions due to the topological constraint of the conservation of the magnetic flux. Using a… ▽ More Many explosive astrophysical events, like magnetars' bursts and flares, are magnetically driven. We consider dynamics of such magnetic explosions - relativistic expansion of highly magnetized and highly magnetically over-pressurized clouds. The corresponding dynamics is qualitatively different from fluid explosions due to the topological constraint of the conservation of the magnetic flux. Using analytical, relativistic MHD as well as force-free calculations, we find that the creation of a relativistically expanding, causally disconnected flow obeys a threshold condition: it requires sufficiently high initial over-pressure and sufficiently quick decrease of the pressure in the external medium (the pre-explosion wind). In the subcritical case the magnetic cloud just "puffs-up" and quietly expands with the pre-flare wind. We also find a compact analytical solution to the Prendergast's problem - expansion of force-free plasma into vacuum. △ Less

Submitted 2 May, 2022; originally announced May 2022.

Comments: 31 pages, 23 figures, submitted to ApJ

Force-free and twisted, relativistic neutron star magnetospheres

Authors: Dimitris Ntotsikas, Konstantinos N. Gourgouliatos

Abstract: In this poster we present the structure of an axisymmetric, force-free magnetosphere of a twisted, aligned rotating dipole within a corotating plasma-filled magnetosphere. We explore various profiles for the twist. We find that as the current increases more field lines cross the light cylinder leading to more efficient spin-down. Moreover, we notice that the twist cannot be increased indefinitely… ▽ More In this poster we present the structure of an axisymmetric, force-free magnetosphere of a twisted, aligned rotating dipole within a corotating plasma-filled magnetosphere. We explore various profiles for the twist. We find that as the current increases more field lines cross the light cylinder leading to more efficient spin-down. Moreover, we notice that the twist cannot be increased indefinitely and after a finite twist of about $π/2$ the field becomes approximately radial. This could have implications for torque variations of magnetars related to outbursts. △ Less

Submitted 18 March, 2022; v1 submitted 17 March, 2022; originally announced March 2022.

Comments: 2 pages, 2 figures, IAU SYMPOSIUM International Astronomical Union Proceedings Series

Orbit determination of the moons of the Pluto-Charon system

Authors: Dionysios Gakis, Konstantinos N. Gourgouliatos

Abstract: Four small moons (Styx, Nix, Kerberos and Hydra) are at present known to orbit around the barycenter of Pluto and Charon, which are themselves considered a binary dwarf-planet due to their relatively high mass ratio. The central, non-axisymmetric potential induces moon orbits inconvenient to be described by Keplerian osculating elements. Here, we report that observed orbital variations, may not be… ▽ More Four small moons (Styx, Nix, Kerberos and Hydra) are at present known to orbit around the barycenter of Pluto and Charon, which are themselves considered a binary dwarf-planet due to their relatively high mass ratio. The central, non-axisymmetric potential induces moon orbits inconvenient to be described by Keplerian osculating elements. Here, we report that observed orbital variations, may not be the result of orbital eccentricities or observational uncertainties, but may be due to forced oscillations caused by the central binary. We show, using numerical integration and analytical considerations, that the differences reported on their orbital elements, may well arise from this intrinsic behavior rather than limitations on our instruments. △ Less

Submitted 27 February, 2022; originally announced February 2022.

Comments: 20 pages, 8 figures, accepted for publication in Celestial Mechanics and Dynamical Astronomy

Journal ref: Celestial Mechanics and Dynamical Astronomy 134, 14 (2022), 1-18

Magnetic Field Evolution in the Crust of Neutron Stars: Crust Failure and Plastic Flow

Authors: Konstantinos N. Gourgouliatos

Abstract: The evolution of the magnetic field in neutron star crusts because of the Hall effect has received significant attention over the last two decades, which is strongly justified because of the dominance of this effect in highly magnetised neutron stars. However, the applicability of the Hall effect is based on the assumption that the crust does not fail and sustains its rigidity. This assumption can… ▽ More The evolution of the magnetic field in neutron star crusts because of the Hall effect has received significant attention over the last two decades, which is strongly justified because of the dominance of this effect in highly magnetised neutron stars. However, the applicability of the Hall effect is based on the assumption that the crust does not fail and sustains its rigidity. This assumption can be violated for substantially strong magnetic fields. If this is the case, the evolution of the magnetic field is described by a different set of equations, which include the effects of a non-rigid crust. In this talk, after a brief review of the main characteristic of the Hall evolution, I will discuss the impact a plastic flow of the crust has on the magnetic field, studying axisymmetric models. Moreover, the way the crust fails impacts the overall evolution, with major differences appearing if the failure is local, intermediate or global. Quite remarkably, crustal failure and plasticity do not annul the Hall effect, and under certain circumstances they may even lead to a more dramatic evolution. I will discuss the impact of these effects in the context of neutron star timing behaviour, with special focus on timing noise, outbursts and glitches. △ Less

Submitted 14 February, 2022; originally announced February 2022.

Comments: 11 pages, 5 figures, IAU Symposium 363, Neutron Star Astrophysics at the Crossroads: Magnetars and the Multimessenger Revolution

Application of an Upwind Integration Method to Plane Parallel Hall-MHD

Authors: Georgios Chouliaras, K. N Gourgouliatos

Abstract: We study the impact of an Upwind scheme on the numerical convergence of simulations of the Hall and Ohmic effect in neutron stars crusts. While simulations of these effects have explored a variety of geometries and wide ranges of physical parameters, they are limited to relatively low values of the Hall parameter, playing the role of the magnetic Reynolds number, which should be not exceed a few h… ▽ More We study the impact of an Upwind scheme on the numerical convergence of simulations of the Hall and Ohmic effect in neutron stars crusts. While simulations of these effects have explored a variety of geometries and wide ranges of physical parameters, they are limited to relatively low values of the Hall parameter, playing the role of the magnetic Reynolds number, which should be not exceed a few hundred for numerical convergence. We study the evolution of the magnetic field in a plane-parallel Cartesian geometry. We discretise the induction equation using a finite difference scheme and then integrate it via the Euler forward method. Two different approaches are used for the integration of the advective terms appearing in the equation: a Forward Time and Central in Space (FTCS) and an Upwind scheme. We compare them in terms of accuracy and performance. We explore the impact of the Upwind method on convergence according to the ratio of planar to vertical field and the Hall parameter. In the limit of a low strength planar field the use of an Upwind scheme provides a vast improvement leading to the convergence of simulations where the Hall parameter is 2 orders of magnitude higher than that of the FTCS. Upwind is still better if the planar field is stronger, yet, the difference of the maximum value of the Hall parameter reached is within a factor of 10 or a few. Moreover, we notice if the schemes diverge their behaviour is very different, with FTCS producing infinite energy, while the Upwind scheme only temporarily increasing the overall magnetic field energy. Overall, the Upwind scheme enhances the efficiency of the simulations allowing the exploration of environments with higher value of electric conductivity getting us closer than before to realistic environmental conditions of magnetars. △ Less

Submitted 3 February, 2022; originally announced February 2022.

Magnetic Field Evolution in Neutron Star Crusts: Beyond the Hall Effect

Authors: Konstantinos N. Gourgouliatos, Davide De Grandis, Andrei Igoshev

Abstract: Neutron stars host the strongest magnetic fields that we know of in the Universe. Their magnetic fields are the main means of generating their radiation, either magnetospheric or through the crust. Moreover, the evolution of the magnetic field has been intimately related to explosive events of magnetars, which host strong magnetic fields, and their persistent thermal emission. The evolution of the… ▽ More Neutron stars host the strongest magnetic fields that we know of in the Universe. Their magnetic fields are the main means of generating their radiation, either magnetospheric or through the crust. Moreover, the evolution of the magnetic field has been intimately related to explosive events of magnetars, which host strong magnetic fields, and their persistent thermal emission. The evolution of the magnetic field in the crusts of neutron stars has been described within the framework of the Hall effect and Ohmic dissipation. Yet, this description is limited by the fact that the Maxwell stresses exerted on the crusts of strongly magnetised neutron stars may lead to failure and temperature variations. In the former case, a failed crust does not completely fulfil the necessary conditions for the Hall effect. In the latter, the variations of temperature are strongly related to the magnetic field evolution. Finally, sharp gradients of the star's temperature may activate battery terms and alter the magnetic field structure, especially in weakly magnetised neutron stars. In this review, we discuss the recent progress made on these effects. We argue that these phenomena are likely to provide novel insight into our understanding of neutron stars and their observable properties. △ Less

Submitted 21 January, 2022; v1 submitted 20 January, 2022; originally announced January 2022.

Comments: 18 Pages, 6 Figures, Symmetry MDPI

Polytropic Wind Solutions via the Complex Plane Strategy

Authors: Vasileios Karageorgopoulos, Konstantinos N. Gourgouliatos, Vassilis Geroyannis

Abstract: Solar-type stars generate spherical winds, which are pressure driven flows, that start subsonic, reach the sound speed at the sonic point and transition to supersonic flows. The sonic point, mathematically corresponds to a singularity of the system of differential equations describing the flow. In the problem of an isothermal wind, the Parker solution provides an exact analytical expression tuned… ▽ More Solar-type stars generate spherical winds, which are pressure driven flows, that start subsonic, reach the sound speed at the sonic point and transition to supersonic flows. The sonic point, mathematically corresponds to a singularity of the system of differential equations describing the flow. In the problem of an isothermal wind, the Parker solution provides an exact analytical expression tuned appropriately so that the singularity does not affect the solution. However, if the wind is polytropic it is not possible to find an analytical solution and a numerical approach needs to be followed. We study solutions of spherical winds that are driven by pressure within a gravitational field. The solutions pass smoothly from the critical point and allow us to study the impact of the changes of the polytropic index to these winds. We explore the properties of these solutions as a function of the polytropic index and the boundary conditions used. We apply the Complex Plane Strategy (CPS) and we obtain numerically solutions of polytropic winds. This allows us to avoid the singularity appearing in the equations through the introduction complex variables and integration on the complex plane. Applying this method, we obtain solutions with physical behaviour at the stellar surface, the sonic point and at large distances from the star. We further explore the role of the polytropic index in the flow and the effect of mass--loss rate and temperature on the solution. We find that the increase of polytropic index as well as the decrease of flow parameter both yield to a smoother velocity profile and lower velocities and shifts the transition point from subsonic to supersonic behaviour further from the star. Finally, we verify that the increasing of coronal temperature yields higher wind velocities and a weaker dependence on polytropic index. △ Less

Submitted 28 November, 2021; originally announced November 2021.

Comments: Article history to Astronomy and Computing: Submitted 26 December 2020, Accepted 5 August 2021, Available online 11 August 2021

Journal ref: Astronomy and Computing, Year: 2021, Volume: 36, Article id 100491

Axisymmetric magneto-plastic evolution of neutron-star crusts

Authors: Konstantinos N. Gourgouliatos, Samuel K. Lander

Abstract: Magnetic field evolution in neutron-star crusts is driven by the Hall effect and Ohmic dissipation, for as long as the crust is sufficiently strong to absorb Maxwell stresses exerted by the field and thus make the momentum equation redundant. For the strongest neutron-star fields, however, stresses build to the point of crustal failure, at which point the standard evolution equations are no longer… ▽ More Magnetic field evolution in neutron-star crusts is driven by the Hall effect and Ohmic dissipation, for as long as the crust is sufficiently strong to absorb Maxwell stresses exerted by the field and thus make the momentum equation redundant. For the strongest neutron-star fields, however, stresses build to the point of crustal failure, at which point the standard evolution equations are no longer valid. Here, we study the evolution of the magnetic field of the crust up to and beyond crustal failure, whence the crust begins to flow plastically. We perform global axisymmetric evolutions, exploring different types of failure affecting a limited region of the crust. We find that a plastic flow does not simply suppress the Hall effect even in the regime of a low plastic viscosity, but it rather leads to non-trivial evolution -- in some cases even overreacting and enhancing the impact of the Hall effect. Its impact is more pronouced in the toroidal field, with the differences on the poloidal field being less substantial. We argue that both the nature of magnetar bursts and their spindown evolution will be affected by plastic flow, so that observations of these phenomena may help to constrain the way the crust fails. △ Less

Submitted 1 July, 2021; v1 submitted 7 June, 2021; originally announced June 2021.

Comments: 12 page, 8 figures, Accepted by MNRAS

3D Magneto-thermal Simulations of Tangled Crustal Magnetic Field in Central Compact Objects

Authors: Andrei P. Igoshev, Konstantinos N. Gourgouliatos, Rainer Hollerbach, Toby S. Wood

Abstract: Central compact objects are young neutron stars emitting thermal X-rays with bolometric luminosities $L_X$ in the range $10^{32}$-$10^{34}$ erg/s. Gourgouliatos, Hollerbach and Igoshev recently suggested that peculiar emission properties of central compact objects can be explained by tangled magnetic field configurations formed in a stochastic dynamo during the proto-neutron star stage. In this ca… ▽ More Central compact objects are young neutron stars emitting thermal X-rays with bolometric luminosities $L_X$ in the range $10^{32}$-$10^{34}$ erg/s. Gourgouliatos, Hollerbach and Igoshev recently suggested that peculiar emission properties of central compact objects can be explained by tangled magnetic field configurations formed in a stochastic dynamo during the proto-neutron star stage. In this case the magnetic field consists of multiple small-scale components with negligible contribution of global dipolar field. We study numerically three-dimensional magneto-thermal evolution of tangled crustal magnetic fields in neutron stars. We find that all configurations produce complicated surface thermal patterns which consist of multiple small hot regions located at significant separations from each other. The configurations with initial magnetic energy of $2.5-10\times 10^{47}$ erg have temperatures of hot regions that reach $\approx 0.2$ keV, to be compared with the bulk temperature of $\approx 0.1$ keV in our simulations with no cooling. A factor of two in temperature is also seen in observations of central compact objects. The hot spots produce periodic modulations in light curve with typical amplitudes of $\leq 9-11$ %. Therefore, the tangled magnetic field configuration can explain thermal emission properties of some central compact objects. △ Less

Submitted 20 January, 2021; originally announced January 2021.

Comments: 12 pages, 25 figures, accepted for publication ApJ on January 19, 2021

Magnetic inhibition of the recollimation instability in relativistic jets

Authors: Jin Matsumoto, Serguei S. Komissarov, Konstantinos N. Gourgouliatos

Abstract: In this paper, we describe the results of three-dimensional relativistic magnetohydrodynamic simulations aimed at probing the role of regular magnetic field on the development of the instability that accompanies recollimation of relativistic jets. In particular, we studied the recollimation driven by the reconfinement of jets from active galactic nuclei (AGN) by the thermal pressure of galactic co… ▽ More In this paper, we describe the results of three-dimensional relativistic magnetohydrodynamic simulations aimed at probing the role of regular magnetic field on the development of the instability that accompanies recollimation of relativistic jets. In particular, we studied the recollimation driven by the reconfinement of jets from active galactic nuclei (AGN) by the thermal pressure of galactic coronas. We find that a relatively weak azimuthal magnetic field can completely suppress the recollimation instability in such jets, with the critical magnetisation parameter $σ_{cr} < 0.01$. We argue that the recollimation instability is a variant of the centrifugal instability (CFI) and show that our results are consistent with the predictions based on the study of magnetic CFI in rotating fluids. The results are discussed in the context of AGN jets in general and the nature of the Fanaroff-Riley morphological division of extragalactic radio sources in particular. △ Less

Submitted 18 March, 2021; v1 submitted 21 October, 2020; originally announced October 2020.

Comments: 12 pages, 8 figures, accepted for publication in MNRAS

Sleeping beasts: strong toroidal magnetic field in quiescent magnetars explains their large pulsed fraction

Authors: Andrei P. Igoshev, Rainer Hollerbach, Toby Wood, Konstantinos N. Gourgouliatos

Abstract: Magnetars are neutron stars (NSs) with extreme magnetic fields of strength $5 \times 10^{13}$ - $10^{15}$ G. They exhibit transient, highly energetic events, such as short X-ray flashes, bursts and giant flares, all of which are powered by their enormous magnetic energy. Quiescent magnetars have X-ray luminosities between $10^{29}$ and $10^{35}$ erg/s, and are further classified as either persiste… ▽ More Magnetars are neutron stars (NSs) with extreme magnetic fields of strength $5 \times 10^{13}$ - $10^{15}$ G. They exhibit transient, highly energetic events, such as short X-ray flashes, bursts and giant flares, all of which are powered by their enormous magnetic energy. Quiescent magnetars have X-ray luminosities between $10^{29}$ and $10^{35}$ erg/s, and are further classified as either persistent or transient magnetars. Their X-ray emission is modulated with the rotational period of the NS, with a typical relative amplitude (so-called pulsed fraction) between 10-58 per cent, implying that the surface temperature is significantly non-uniform despite the high thermal conductivity of the star's crust. Here, we present the first 3D magneto-thermal MHD simulations of magnetars with strong toroidal magnetic fields. We show that these models, combined with ray propagation in curved space-time, accurately describe the light-curves of most transient magnetars in quiescence and allow us to further constrain their rotational orientation. We find that the presence of a strong toroidal magnetic field explains the observed asymmetry in the surface temperature, and is the main cause of the strong modulation of thermal X-ray emission in quiescence. △ Less

Submitted 16 October, 2020; originally announced October 2020.

Comments: Original version of manuscript submitted to the Nature Astronomy. 11 pages, 6 figures

3D modelling of magneto-thermal evolution of neutron stars: method and test cases

Authors: Davide De Grandis, Roberto Turolla, Toby S. Wood, Silvia Zane, Roberto Taverna, Konstantinos N. Gourgouliatos

Abstract: Neutron stars harbour extremely strong magnetic fields within their solid outer crust. The topology of this field strongly influences the surface temperature distribution, and hence the star's observational properties. In this work, we present the first realistic simulations of the coupled crustal magneto-thermal evolution of isolated neutron stars in three dimensions with account for neutrino emi… ▽ More Neutron stars harbour extremely strong magnetic fields within their solid outer crust. The topology of this field strongly influences the surface temperature distribution, and hence the star's observational properties. In this work, we present the first realistic simulations of the coupled crustal magneto-thermal evolution of isolated neutron stars in three dimensions with account for neutrino emission, obtained with the pseudo-spectral code Parody. We investigate both the secular evolution, especially in connection with the onset of instabilities during the Hall phase, and the short-term evolution following episodes of localised energy injection. Simulations show that a resistive tearing instability develops in about a Hall time if the initial toroidal field exceeds ~$10^{15}$ G. This leads to crustal failures because of the huge magnetic stresses coupled with the local temperature enhancement produced by dissipation. Localised heat deposition in the crust results in the appearance of hot spots on the star surface which can exhibit a variety of patterns. Since the transport properties are strongly influenced by the magnetic field, the hot regions tend to drift away and get deformed following the magnetic field lines while cooling. The shapes obtained with our simulations are reminiscent of those recently derived from NICER X-ray observations of the millisecond pulsar PSR J0030+0451. △ Less

Submitted 9 September, 2020; originally announced September 2020.

Comments: 15 pages, 14 figures, accepted for publication in ApJ

Powering Central Compact Objects with a Tangled Crustal Magnetic Field

Authors: Konstantinos N. Gourgouliatos, Rainer Hollerbach, Andrei P. Igoshev

Abstract: Central Compact Objects (CCOs) are X-ray sources with luminosity ranging between $10^{32}$-$10^{34}$ erg~s$^{-1}$, located at the centres of supernova remnants. Some of them have been confirmed to be neutron stars. Timing observations have allowed the estimation of their dipole magnetic field, placing them in the range $\sim10^{10}$-$10^{11}$ G. The decay of their weak dipole fields, mediated by t… ▽ More Central Compact Objects (CCOs) are X-ray sources with luminosity ranging between $10^{32}$-$10^{34}$ erg~s$^{-1}$, located at the centres of supernova remnants. Some of them have been confirmed to be neutron stars. Timing observations have allowed the estimation of their dipole magnetic field, placing them in the range $\sim10^{10}$-$10^{11}$ G. The decay of their weak dipole fields, mediated by the Hall effect and Ohmic dissipation, cannot provide sufficient thermal energy to power their X-ray luminosity, as opposed to magnetars whose X-ray luminosities are comparable. Motivated by the question of producing high X-ray power through magnetic field decay while maintaining a weak dipole field, we explore the evolution of a crustal magnetic field that does not consist of an ordered axisymmetric structure, but rather comprises a tangled configuration. This can be the outcome of a non-self-excited dynamo, buried inside the crust by fallback material following the supernova explosion. We find that such initial conditions lead to the emergence of the magnetic field from the surface of the star and the formation of a dipolar magnetic field component. An internal tangled magnetic field of the order of $10^{14}$ G can provide sufficient Ohmic heating to the crust and power CCOs, while the dipole field it forms is approximately $10^{10}$ G, as observed in CCOs. △ Less

Submitted 5 May, 2020; originally announced May 2020.

Comments: 9 pages, Accepted by MNRAS

Nonaxisymmetric Hall instability: A key to understanding magnetars

Authors: Konstantinos N. Gourgouliatos, José A. Pons

Abstract: It is generally accepted that the non-linear, dynamical evolution of magnetic fields in the interior of neutron stars plays a key role in the explanation of the observed phenomenology. Understanding the transfer of energy between toroidal and poloidal components, or between different scales, is of particular relevance. In this letter, we present the first 3D simulations of the Hall instability in… ▽ More It is generally accepted that the non-linear, dynamical evolution of magnetic fields in the interior of neutron stars plays a key role in the explanation of the observed phenomenology. Understanding the transfer of energy between toroidal and poloidal components, or between different scales, is of particular relevance. In this letter, we present the first 3D simulations of the Hall instability in a neutron star crust, confirming its existence for typical magnetar conditions. We confront our results to estimates obtained by a linear perturbation analysis, which discards any interpretation as numerical instabilities and confirms its physical origin. Interestingly, the Hall instability creates locally strong magnetic structures that occasionally can make the crust yield to the magnetic stresses and generates coronal loops, similarly as solar coronal loops find their way out through the photosphere. This supports the viability of the mechanism, which has been proposed to explain magnetar outbursts. △ Less

Submitted 10 January, 2020; originally announced January 2020.

Comments: 5 pages, 3 Figures, Published in Physical Review Research, Rapid Communication

Journal ref: Phys. Rev. Research 1, 032049(R) (2019)

Magnetic Inhibition of Centrifugal Instability in Astrophysical Jets

Authors: Serguei S. Komissarov, Konstantinos N. Gourgouliatos, Jin Matsumoto

Abstract: Recently it was shown that the centrifugal instability may be important in the dynamics of astrophysical jets undergoing reconfinement by external pressure. However, these studies were limited to the case of unmagnetised flows. Here we explore the role of the magnetic field within both the Newtonian and relativistic frameworks. Since the jet problem is rather complicated, we focus instead on the s… ▽ More Recently it was shown that the centrifugal instability may be important in the dynamics of astrophysical jets undergoing reconfinement by external pressure. However, these studies were limited to the case of unmagnetised flows. Here we explore the role of the magnetic field within both the Newtonian and relativistic frameworks. Since the jet problem is rather complicated, we focus instead on the simpler problem of cylindrical rotation and axial magnetic field, which shares significant similarity with the jet problem, and consider only axisymmetric perturbations. The studied equilibrium configurations involve a cylindrical interface and they are stable to non-magnetic centrifugal and magneto-rotational instabilities everywhere except this interface. We use a heuristic approach to derive the local stability criterion for the interface in the magnetic case and numerical simulations to verify the role of the magnetic field. The theory and simulations agree quite well for Newtonian models but indicate a potential discrepancy for the relativistic models in the limit of high Lorentz factor of the rotational motion at the interface. In general, the magnetic field sets a critical wavelength below which the centrifugal modes are stabilised. We discuss the implication of our findings for the astrophysical jets, which suggest that the centrifugal instability develops only in jets with relatively low magnetisation. Namely, the magnetic pressure has to be below the thermal one and for the relativistic case the jets have to be kinetic-energy dominated. △ Less

Submitted 20 June, 2019; v1 submitted 28 May, 2019; originally announced May 2019.

Comments: submitted to MNRAS

Monte Carlo studies for the optimisation of the Cherenkov Telescope Array layout

Authors: A. Acharyya, I. Agudo, E. O. Angüner, R. Alfaro, J. Alfaro, C. Alispach, R. Aloisio, R. Alves Batista, J. -P. Amans, L. Amati, E. Amato, G. Ambrosi, L. A. Antonelli, C. Aramo, T. Armstrong, F. Arqueros, L. Arrabito, K. Asano, H. Ashkar, C. Balazs, M. Balbo, B. Balmaverde, P. Barai, A. Barbano, M. Barkov , et al. (445 additional authors not shown)

Abstract: The Cherenkov Telescope Array (CTA) is the major next-generation observatory for ground-based very-high-energy gamma-ray astronomy. It will improve the sensitivity of current ground-based instruments by a factor of five to twenty, depending on the energy, greatly improving both their angular and energy resolutions over four decades in energy (from 20 GeV to 300 TeV). This achievement will be possi… ▽ More The Cherenkov Telescope Array (CTA) is the major next-generation observatory for ground-based very-high-energy gamma-ray astronomy. It will improve the sensitivity of current ground-based instruments by a factor of five to twenty, depending on the energy, greatly improving both their angular and energy resolutions over four decades in energy (from 20 GeV to 300 TeV). This achievement will be possible by using tens of imaging Cherenkov telescopes of three successive sizes. They will be arranged into two arrays, one per hemisphere, located on the La Palma island (Spain) and in Paranal (Chile). We present here the optimised and final telescope arrays for both CTA sites, as well as their foreseen performance, resulting from the analysis of three different large-scale Monte Carlo productions. △ Less

Submitted 2 April, 2019; originally announced April 2019.

Comments: 48 pages, 16 figures, accepted for publication in Astroparticle Physics

Current closure through the neutron star crust

Authors: V. Karageorgopoulos, K. N. Gourgouliatos, I. Contopoulos

Abstract: Force-free pulsar magnetospheres develop a large scale poloidal electric current circuit that flows along open magnetic field lines from the neutron star to the termination shock. The electric current closes through the interior of the neutron star where it provides the torque that spins-down the star. In the present work, we study the internal electric current in an axisymmetric rotator. We evalu… ▽ More Force-free pulsar magnetospheres develop a large scale poloidal electric current circuit that flows along open magnetic field lines from the neutron star to the termination shock. The electric current closes through the interior of the neutron star where it provides the torque that spins-down the star. In the present work, we study the internal electric current in an axisymmetric rotator. We evaluate the path of the electric current by requiring the minimization of internal Ohmic losses. We find that, in millisecond pulsars, the current reaches the base of the crust, while in pulsars with periods of a few seconds, the bulk of the electric current does not penetrate deeper than about $100$ m. The region of maximum spin-down torque in millisecond pulsars is the base of the crust, while in slowly spinning ones it is the outer crust. We evaluate the corresponding Maxwell stresses and find that, in typical rotation-powered radio pulsars, they are well below the critical stress that can be sustained by the crust. For magnetar-level fields, the Maxwell stresses near the surface are comparable to the critical stress and may lead to the decoupling of the crust from the rest of the stellar rotation. △ Less

Submitted 2 June, 2019; v1 submitted 12 March, 2019; originally announced March 2019.

Comments: 10 pages, 6 Figures, Accepted by MNRAS

Magnetic-field evolution in a plastically-failing neutron-star crust

Authors: S. K. Lander, K. N. Gourgouliatos

Abstract: Under normal conditions in a neutron-star crust, ions are locked in place in the crustal lattice and only electrons are mobile, and magnetic-field evolution is thus directly related to the electron velocity. The evolution, however, builds magnetic stresses that can become sufficiently large for the crust to exceed its elastic limit, and to flow plastically. We consider the nature of this plastic f… ▽ More Under normal conditions in a neutron-star crust, ions are locked in place in the crustal lattice and only electrons are mobile, and magnetic-field evolution is thus directly related to the electron velocity. The evolution, however, builds magnetic stresses that can become sufficiently large for the crust to exceed its elastic limit, and to flow plastically. We consider the nature of this plastic flow and the back-reaction on the crustal magnetic field evolution. We formulate a plane-parallel model for the local failure, showing that surface motions are inevitable once the crust yields, in the absence of extra dissipative mechanisms. We perform numerical evolutions of the crustal magnetic field under the joint effect of Hall drift and Ohmic decay, tracking the build-up of magnetic stresses, and diagnosing crustal failure with the von Mises criterion. Beyond this point we solve for the coupled evolution of the plastic velocity and magnetic field. Our results suggest that to have a coexistence of a magnetar corona with small-scale magnetic features, the viscosity of the plastic flow must be roughly $10^{36}-10^{37}\ \textrm{g cm}^{-1}\textrm{s}^{-1}$. We find significant motion at the surface at a rate of $10-100$ centimetres per year, and that the localised magnetic field is weaker than in evolutions without plastic flow. We discuss astrophysical implications, and how our local simulations could be used to build a global model of field evolution in the neutron-star crust. △ Less

Submitted 14 June, 2019; v1 submitted 6 February, 2019; originally announced February 2019.

Comments: 15 pages, 7 figures. Published version; modest edits from v1 for clarity

Journal ref: MNRAS 486 (2019) 4130

Observational signatures of magnetic field structure in relativistic AGN jets

Authors: Christopher Prior, Konstantinos N. Gourgouliatos

Abstract: Context. Active Galactic Nuclei (AGN) launch highly energetic jets sometimes outshining their host galaxy. These jets are collimated outflows that have been accelerated near a supermassive black hole located at the centre of the galaxy. Their, virtually indispensable, energy reservoir is either due to gravitational energy released from accretion or due to the extraction of kinetic energy from the… ▽ More Context. Active Galactic Nuclei (AGN) launch highly energetic jets sometimes outshining their host galaxy. These jets are collimated outflows that have been accelerated near a supermassive black hole located at the centre of the galaxy. Their, virtually indispensable, energy reservoir is either due to gravitational energy released from accretion or due to the extraction of kinetic energy from the rotating supermassive black hole itself. In order to channel part of this energy to the jet, though, the presence of magnetic fields is necessary. The extent to which these magnetic fields survive in the jet further from the launching region is under debate. Nevertheless, observations of polarised emission and Faraday rotation measure confirm the existence of large scale magnetic fields in jets. Results. Radio emission distributions from the jets with dominant large-scale helical fields show asymmetry across their width. The Faraday rotation asymmetry is the same for fields with opposing chirality (handedness). For jets which are tilted towards the observer the synchrotron emission and fractional polarization can distinguish the fields chirality. When viewed either side-on or at a Blazar type angle only the fractional polarization can make this distinction. Further this distinction can only be made if the direction of the jet propagation velocity is known, along with the location of the jet's origin. The complex structure of the braided field is found not to be observable due to a combination of line of sight integration and limited resolution of observation. This raises the possibility that, even if asymmetric radio emission signatures are present, the true structure of the field may still be obscure. △ Less

Submitted 15 January, 2019; originally announced January 2019.

Comments: Accepted by Astronomy and Astrophysics (abstract shortened to fit within arXiv character limit)

Coupled axisymmetric pulsar magnetospheres

Authors: Konstantinos N. Gourgouliatos, Donald Lynden-Bell

Abstract: We present solutions of force-free pulsar magnetospheres coupled with a uniform external magnetic field aligned with the dipole magnetic moment of the pulsar. The inclusion of the uniform magnetic field has the following consequences: The equatorial current sheet is truncated to a finite length, the fraction of field lines that are open increases, and the open field lines are confined within a cyl… ▽ More We present solutions of force-free pulsar magnetospheres coupled with a uniform external magnetic field aligned with the dipole magnetic moment of the pulsar. The inclusion of the uniform magnetic field has the following consequences: The equatorial current sheet is truncated to a finite length, the fraction of field lines that are open increases, and the open field lines are confined within a cylindrical surface instead of becoming radial. A strong external magnetic field antiparallel to the dipole allows for solutions where the pulsar magnetic field is fully enclosed within an ellipsoidal surface. Configurations of fully enclosed or confined magnetospheres may appear in a double neutron star (DNS) where one of the components is a slowly rotating, strongly magnetised pulsar and the other a weakly magnetised millisecond pulsar. Depending on the geometry, twisted field lines from the millisecond pulsar could form an asymmetric wind. More dramatic consequences could appear closer to merger: there, the neutron star with the weaker magnetic field may undergo a stage where it alternates between an open and a fully enclosed magnetosphere releasing up $10^{38}$erg. Under favourable combinations of magnetic fields DNSs can spend up to 10$^3$ yr in the coupled phase, implying a Galactic population of $0.02$ systems. Potential observational signatures of this interaction are discussed including the possibility of powering recurring Fast Radio Bursts (FRBs). We also note that the magnetic interaction cannot have any significant effect on the spin evolution of the two pulsars. △ Less

Submitted 9 October, 2018; originally announced October 2018.

Comments: Accepted by MNRAS, 13 pages, 7 Figures

Reconfinement and Loss of Stability in Jets from Active Galactic Nuclei

Authors: Konstantinos N. Gourgouliatos, Serguei S. Komissarov

Abstract: Jets powered by active galactic nuclei appear impressively stable compared with their terrestrial and laboratory counterparts-they can be traced from their origin to distances exceeding their injection radius by up to a billion times. However, some less energetic jets get disrupted and lose their coherence on the scale of their host galaxy. Quite remarkably, on the same scale, these jets are expec… ▽ More Jets powered by active galactic nuclei appear impressively stable compared with their terrestrial and laboratory counterparts-they can be traced from their origin to distances exceeding their injection radius by up to a billion times. However, some less energetic jets get disrupted and lose their coherence on the scale of their host galaxy. Quite remarkably, on the same scale, these jets are expected to become confined by the thermal pressure of the intra-galactic gas. Motivated by these observations, we have started a systematic study of active galactic nuclei jets undergoing reconfinement via computer simulations. Here, we show that in the case of unmagnetized relativistic jets, the reconfinement is accompanied by the development of an instability and transition to a turbulent state. During their initial growth, the perturbations have a highly organized streamwise-oriented structure, indicating that it is not the Kelvin-Helmholtz instability, the instability which has been the main focus of the jet stability studies so far. Instead, it is closely related to the centrifugal instability. This instability is likely to be behind the division of active galactic nuclei jets into two morphological types in the Fanaroff-Riley classification. △ Less

Submitted 14 June, 2018; originally announced June 2018.

Comments: Published in Nature Astronomy

Journal ref: Nature Astronomy volume 2, pages 167-171 (2018)

Strongly magnetized pulsars: explosive events and evolution

Authors: Konstantinos N. Gourgouliatos, Paolo Esposito

Abstract: Well before the radio discovery of pulsars offered the first observational confirmation for their existence (Hewish et al., 1968), it had been suggested that neutron stars might be endowed with very strong magnetic fields of $10^{10}$-$10^{14}$G (Hoyle et al., 1964; Pacini, 1967). It is because of their magnetic fields that these otherwise small ed inert, cooling dead stars emit radio pulses and s… ▽ More Well before the radio discovery of pulsars offered the first observational confirmation for their existence (Hewish et al., 1968), it had been suggested that neutron stars might be endowed with very strong magnetic fields of $10^{10}$-$10^{14}$G (Hoyle et al., 1964; Pacini, 1967). It is because of their magnetic fields that these otherwise small ed inert, cooling dead stars emit radio pulses and shine in various part of the electromagnetic spectrum. But the presence of a strong magnetic field has more subtle and sometimes dramatic consequences: In the last decades of observations indeed, evidence mounted that it is likely the magnetic field that makes of an isolated neutron star what it is among the different observational manifestations in which they come. The contribution of the magnetic field to the energy budget of the neutron star can be comparable or even exceed the available kinetic energy. The most magnetised neutron stars in particular, the magnetars, exhibit an amazing assortment of explosive events, underlining the importance of their magnetic field in their lives. In this chapter we review the recent observational and theoretical achievements, which not only confirmed the importance of the magnetic field in the evolution of neutron stars, but also provide a promising unification scheme for the different observational manifestations in which they appear. We focus on the role of their magnetic field as an energy source behind their persistent emission, but also its critical role in explosive events. △ Less

Submitted 4 May, 2018; originally announced May 2018.

Comments: Review commissioned for publication in the White Book of "NewCompStar" European COST Action MP1304, 43 pages, 8 figures

Journal ref: The Physics and Astrophysics of Neutron Stars, eds L. Rezzolla, P. M. Pizzochero, D. I. Jones, N. Rea, and I. Vidaña, Astrophysics and Space Science Library, Springer Nature Switzerland, Basel, Vol. 457, pp. 57-93 (2018)

Relativistic Centrifugal Instability

Authors: Konstantinos N. Gourgouliatos, Serguei S. Komissarov

Abstract: Near the central engine, many astrophysical jets are expected to rotate about their axis. Further out they are expected to go through the processes of reconfinement and recollimation. In both these cases, the flow streams along a concave surface and hence, it is subject to the centrifugal force. It is well known that such flows may experience the Centrifugal Instability (CFI), to which there are m… ▽ More Near the central engine, many astrophysical jets are expected to rotate about their axis. Further out they are expected to go through the processes of reconfinement and recollimation. In both these cases, the flow streams along a concave surface and hence, it is subject to the centrifugal force. It is well known that such flows may experience the Centrifugal Instability (CFI), to which there are many laboratory examples. The recent computer simulations of relativistic jets from Active Galactic Nuclei undergoing the process of reconfinement show that in such jets CFI may dominate over the Kelvin-Helmholtz instability associated with velocity shear (Gourgouliatos and Komissarov 2018). In this letter, we generalise the Rayleigh criterion for CFI in rotating fluids to relativistic flows using a heuristic analysis. We also present the results of computer simulations which support our analytic criterion for the case of an interface separating two uniformly-rotating cylindrical flows. We discuss the difference between CFI and the Rayleigh-Taylor instability in flows with curved streamlines. △ Less

Submitted 31 January, 2018; v1 submitted 3 October, 2017; originally announced October 2017.

Comments: Accepted by MNRAS Letters (submitted initially to PRL and JFM)

Magnetic Axis Drift and Magnetic Spot Formation in Neutron Stars with Toroidal Fields

Authors: Konstantinos N. Gourgouliatos, Rainer Hollerbach

Abstract: We explore magnetic field configurations that lead to the formation of magnetic spots on the surface of neutron stars, and to the displacement of the magnetic dipole axis. We find that a toroidally dominated magnetic field is essential for the generation of a single spot with a strong magnetic field. Once a spot forms, it survives for several million years, even after the total magnetic field has… ▽ More We explore magnetic field configurations that lead to the formation of magnetic spots on the surface of neutron stars, and to the displacement of the magnetic dipole axis. We find that a toroidally dominated magnetic field is essential for the generation of a single spot with a strong magnetic field. Once a spot forms, it survives for several million years, even after the total magnetic field has decayed significantly. We find that the dipole axis is not stationary with respect to the neutron star's surface and does not in general coincide with the location of the magnetic spot. This is due to non-axisymmetric instabilities of the toroidal field that displace the poloidal dipole axis at rates that may reach $0.4^{\circ}$ per century. A misaligned poloidal dipole axis with the toroidal field leads to more significant displacement of the dipole axis than the fully aligned case. Finally we discuss the evolution of neutron stars with such magnetic fields on the $P-\dot{P}$ diagram and the observational implications. We find that neutron stars spend a very short time before they cross the Death-Line of the $P-\dot{P}$ diagram, compared to their characteristic ages. Moreover, the maximum intensity of their surface magnetic field is substantially higher than the dipole component of the field. We argue that SGR 0418+5729 could be an example of this type of behaviour, having a weak dipole field, yet hosting a magnetic spot responsible for its magnetar behaviour. The evolution on the pulse profile and braking index of the Crab pulsar, which are attributed to an increase of its obliquity, are compatible with the anticipated drift of the magnetic axis. △ Less

Submitted 13 December, 2017; v1 submitted 3 October, 2017; originally announced October 2017.

Comments: 12 pages, 8 figures, Accepted by the Astrophysical Journal

Journal ref: 2018 ApJ 852 21

The glitches and rotational history of the highly energetic young pulsar PSR J0537$-$6910

Authors: Robert D. Ferdman, Robert F. Archibald, Konstantinos N. Gourgouliatos, Victoria M. Kaspi

Abstract: We present a timing and glitch analysis of the young X-ray pulsar PSR J0537$-$6910, located within the Large Magellanic Cloud, using 13 years of data from the now decommissioned Rossi X-ray Timing Explorer. Rotating with a spin period of 16 ms, PSR J0537$-$6910 is the fastest spinning and most energetic young pulsar known. It also displays the highest glitch activity of any known pulsar. We have f… ▽ More We present a timing and glitch analysis of the young X-ray pulsar PSR J0537$-$6910, located within the Large Magellanic Cloud, using 13 years of data from the now decommissioned Rossi X-ray Timing Explorer. Rotating with a spin period of 16 ms, PSR J0537$-$6910 is the fastest spinning and most energetic young pulsar known. It also displays the highest glitch activity of any known pulsar. We have found 42 glitches over the data span, corresponding to a glitch rate of 3.2 yr$^{-1}$, with an overall glitch activity rate of $8.8\times 10^{-7}\,$yr$^{-1}$. The high glitch frequency has allowed us to study the glitch behavior in ways that are inaccessible in other pulsars. We observe a strong linear correlation between spin frequency glitch magnitude and wait time to the following glitch. We also find that the post-glitch spin-down recovery is well described by a single two-component model fit to all glitches for which we have adequate input data. This consists of an exponential amplitude $A = (7.6 \pm 1.0)\times 10^{-14}\,$s$^{-2}$ and decay timescale $τ= 27_{-6}^{+7}\,$d, and linear slope $m = (4.1\pm 0.4)\times 10^{-16}\,$s$^{-2}\,$d$^{-1}$. The latter slope corresponds to a second frequency derivative $\ddotν = (4.7\pm 0.5) \times 10^{-22}\,$s$^{-3}$, from which we find an implied braking index $n=7.4 \pm 0.8$. We also present a maximum-likelihood technique for searching for periods in event-time data, which we used to both confirm previously published values and determine rotation frequencies in later observations. We discuss the implied constraints on glitch models from the observed behavior of this system, which we argue cannot be fully explained in the context of existing theories. △ Less

Submitted 11 December, 2017; v1 submitted 29 August, 2017; originally announced August 2017.

Comments: 17 pages, 2 tables, 8 figures. Accepted for publication in the Astrophysical Journal

Resistive Tearing Instability in Electron-MHD: Application to Neutron Star Crusts

Authors: Konstantinos N. Gourgouliatos, Rainer Hollerbach

Abstract: We study a resistive tearing instability developing in a system evolving through the combined effect of Hall drift in the Electron-MHD limit and Ohmic dissipation. We explore first the exponential growth of the instability in the linear case and we find the fastest growing mode, the corresponding eigenvalues and dispersion relation. The instability growth rate scales as… ▽ More We study a resistive tearing instability developing in a system evolving through the combined effect of Hall drift in the Electron-MHD limit and Ohmic dissipation. We explore first the exponential growth of the instability in the linear case and we find the fastest growing mode, the corresponding eigenvalues and dispersion relation. The instability growth rate scales as $γ\propto B^{2/3} σ^{-1/3}$ where $B$ is the magnetic field and $σ$ the electrical conductivity. We confirm the development of the tearing resistive instability in the fully non-linear case, in a plane parallel configuration where the magnetic field polarity reverses, through simulations of systems initiating in Hall equilibrium with some superimposed perturbation. Following a transient phase, during which there is some minor rearrangement of the magnetic field, the perturbation grows exponentially. Once the instability is fully developed the magnetic field forms the characteristic islands and X-type reconnection points, where Ohmic decay is enhanced. We discuss the implications of this instability for the local magnetic field evolution in neutron stars' crusts, proposing that it can contribute to heating near the surface of the star, as suggested by models of magnetar post-burst cooling. In particular, we find that a current sheet a few meters thick, covering as little as $1\%$ of the total surface can provide $10^{42}~$erg in thermal energy within a few days. We briefly discuss applications of this instability in other systems where the Hall effect operates such as protoplanetary discs and space plasmas. △ Less

Submitted 8 September, 2016; v1 submitted 26 July, 2016; originally announced July 2016.

Comments: Accepted by MNRAS, 10 pages, 10 Figures

Magnetic field evolution in magnetar crusts through three dimensional simulations

Authors: Konstantinos N. Gourgouliatos, Toby Wood, Rainer Hollerbach

Abstract: Current models of magnetars require extremely strong magnetic fields to explain their observed quiescent and bursting emission, implying that the field strength within the star's outer crust is orders of magnitude larger than the dipole component inferred from spin-down measurements. This presents a serious challenge to theories of magnetic field generation in a proto-neutron star. Here, we presen… ▽ More Current models of magnetars require extremely strong magnetic fields to explain their observed quiescent and bursting emission, implying that the field strength within the star's outer crust is orders of magnitude larger than the dipole component inferred from spin-down measurements. This presents a serious challenge to theories of magnetic field generation in a proto-neutron star. Here, we present detailed modelling of the evolution of the magnetic field in the crust of a neutron star through 3-D simulations. We find that, in the plausible scenario of equipartition of energy between global-scale poloidal and toroidal magnetic components, magnetic instabilities transfer energy to non-axisymmetric, kilometre-sized magnetic features, in which the local field strength can greatly exceed that of the global-scale field. These intense small-scale magnetic features can induce high energy bursts through local crust yielding, and the localised enhancement of Ohmic heating can power the star's persistent emission. Thus, the observed diversity in magnetar behaviour can be explained with mixed poloidal-toroidal fields of comparable energies. △ Less

Submitted 5 April, 2016; originally announced April 2016.

Comments: Published by PNAS

Magnetar Activity via the Density-Shear Instability in Hall-MHD

Authors: Konstantinos N. Gourgouliatos, Todor Kondic, Maxim Lyutikov, Rainer Hollerbach

Abstract: We investigate the density-shear instability in Hall-MHD via numerical simulation of the full non-linear problem, in the context of magnetar activity. We confirm the development of the instability of a plane-parallel magnetic field with an appropriate intensity and electron density profile, in accordance with analytic theory. We find that the instability also appears for a monotonically decreasing… ▽ More We investigate the density-shear instability in Hall-MHD via numerical simulation of the full non-linear problem, in the context of magnetar activity. We confirm the development of the instability of a plane-parallel magnetic field with an appropriate intensity and electron density profile, in accordance with analytic theory. We find that the instability also appears for a monotonically decreasing electron number density and magnetic field, a plane-parallel analogue of an azimuthal or meridional magnetic field in the crust of a magnetar. The growth rate of the instability depends on the Hall properties of the field (magnetic field intensity, electron number density and the corresponding scale-heights), while being insensitive to weak resistivity. Since the Hall effect is the driving process for the evolution of the crustal magnetic field of magnetars, we argue that this instability is critical for systems containing strong meridional or azimuthal fields. We find that this process mediates the formation of localised structures with much stronger magnetic field than the average, which can lead to magnetar activity and accelerate the dissipation of the field and consequently the production of Ohmic heating. Assuming a $5\times10^{14}$G magnetic field at the base of crust, we anticipate that magnetic field as strong as $10^{15}$G will easily develop in regions of typical size of a few $10^{2}$ meters, containing magnetic energy of $10^{43}$erg, sufficient to power magnetar bursts. These active regions are more likely to appear in the magnetic equator where the tangential magnetic field is stronger. △ Less

Submitted 27 July, 2015; originally announced July 2015.

Comments: Accepted by MNRAS Letters, 3 Figures, 1 Table

Hall drift and the braking indices of young pulsars

Authors: K. N. Gourgouliatos, A. Cumming

Abstract: Braking index measurements of young radio pulsars are all smaller than the value expected for spin down by magnetic dipole braking. We investigate magnetic field evolution in the neutron star crust due to Hall drift as an explanation for observed braking indices. Using numerical simulations and a semi-analytic model, we show that a $\approx 10^{14}\ {\rm G}$ quadrupolar toroidal field in the neutr… ▽ More Braking index measurements of young radio pulsars are all smaller than the value expected for spin down by magnetic dipole braking. We investigate magnetic field evolution in the neutron star crust due to Hall drift as an explanation for observed braking indices. Using numerical simulations and a semi-analytic model, we show that a $\approx 10^{14}\ {\rm G}$ quadrupolar toroidal field in the neutron star crust at birth leads to growth of the dipole moment at a rate large enough to agree with measured braking indices. A key factor is the density at which the crust yields to magnetic stresses that build up during the evolution, which sets a characteristic minimum Hall timescale. The observed braking indices of pulsars with inferred dipole fields of $\lesssim 10^{13}\ {\rm G}$ can be explained in this picture, although with a significant octupole component needed in some cases. For the stronger field pulsars, those with $B_d\gtrsim 10^{13}\ {\rm G}$, we find that the magnetic stresses in the crust exceed the maximum shear stress before the pulsar reaches its current age, likely quenching the Hall effect. This may have implications for the magnetar activity seen in the high magnetic field radio pulsar PSR~J1846-0258. Observations of braking indices may therefore be a new piece of evidence that neutron stars contain subsurface toroidal fields that are significantly stronger than the dipole field, and may indicate that the Hall effect is important in a wider range of neutron stars than previously thought. △ Less

Submitted 19 November, 2014; v1 submitted 13 June, 2014; originally announced June 2014.

Comments: 8 pages, 7 figures, Accepted by MNRAS

Journal ref: MNRAS (2015) Vol. 446 1121-1128

Hall Attractor in Axially Symmetric Magnetic Fields in Neutron Star Crusts

Authors: Konstantinos N. Gourgouliatos, Andrew Cumming

Abstract: We have found an attractor for an axially symmetric magnetic field evolving under the Hall effect and subdominant ohmic dissipation, resolving the question of the long term fate of the magnetic field in neutron star crusts. The electron fluid is in isorotation, analogous to Ferraro's law, with its angular velocity being approximately proportional to the poloidal magnetic flux, $Ω\propto Ψ$. This e… ▽ More We have found an attractor for an axially symmetric magnetic field evolving under the Hall effect and subdominant ohmic dissipation, resolving the question of the long term fate of the magnetic field in neutron star crusts. The electron fluid is in isorotation, analogous to Ferraro's law, with its angular velocity being approximately proportional to the poloidal magnetic flux, $Ω\propto Ψ$. This equilibrium is the long term configuration of a magnetic field evolving because of the Hall effect and ohmic dissipation. For an initial dipole dominated field the attractor consists mainly of a dipole and an octupole component accompanied by an energetically negligible quadrupole toroidal field. The field dissipates in a self-similar way: although higher multipoles should have been decaying faster, the toroidal field mediates transfer of energy into them from the lower ones, leading to an advection diffusion equilibrium and keeping the ratio of the poloidal multipoles almost constant. This has implications for the structure of the intermediate age neutron stars, suggesting that their poloidal field should consist of a dipole and a octupole component accompanied by a very weak toroidal quadrupole. For initial conditions that have a higher multipole $\ell$ structure the attractor consists mainly of $\ell$ and $\ell+2$ poloidal components. △ Less

Submitted 29 April, 2014; v1 submitted 28 November, 2013; originally announced November 2013.

Comments: Accepted by Physical Review Letters, 5 pages, 4 figures

Hall Effect in Neutron Star Crusts: Evolution, Endpoint and Dependence on Initial Conditions

Authors: K. N. Gourgouliatos, A. Cumming

Abstract: We present new simulations of the evolution of axially symmetric magnetic fields in neutron star crusts under the influence of the Hall effect and subdominant Ohmic dissipation. In the Hall effect, differential rotation of the electron fluid generates toroidal field by winding of the poloidal field. For this reason, we focus on the influence of the initial choice of the electron angular velocity p… ▽ More We present new simulations of the evolution of axially symmetric magnetic fields in neutron star crusts under the influence of the Hall effect and subdominant Ohmic dissipation. In the Hall effect, differential rotation of the electron fluid generates toroidal field by winding of the poloidal field. For this reason, we focus on the influence of the initial choice of the electron angular velocity profile on the subsequent and long term magnetic evolution. Whereas previous simulations have generally chosen angular velocities increasing outwards, corresponding to the lowest order Ohmic mode in the crust, a more realistic choice is an angular velocity decreasing outwards, corresponding to the MHD equilibrium field that is likely present at the time of crust formation. We find that the evolution passes through three basic phases. The early evolution is a response to the initial conditions. During the second phase the field consists of poloidal and toroidal components which eventually relax to an isorotation state in which the angular velocity of the electrons becomes constant along poloidal magnetic field lines, causing Hall evolution to saturate. In the third phase the field dissipates slowly while maintaining isorotation. We discuss the implications for the long term field structure and observable properties of isolated neutron stars. △ Less

Submitted 27 November, 2013; originally announced November 2013.

Comments: 13 pages, 8 figures, Accepted by MNRAS

An Anti-Glitch in a Magnetar

Authors: R. F. Archibald, V. M. Kaspi, C. -Y. Ng, K. N. Gourgouliatos, D. Tsang, P. Scholz, A. P. Beardmore, N. Gehrels, J. A. Kennea

Abstract: Magnetars are neutron stars showing dramatic X-ray and soft $γ$-ray outbursting behaviour that is thought to be powered by intense internal magnetic fields. Like conventional young neutron stars in the form of radio pulsars, magnetars exhibit "glitches" during which angular momentum is believed to be transferred between the solid outer crust and the superfluid component of the inner crust. Hithert… ▽ More Magnetars are neutron stars showing dramatic X-ray and soft $γ$-ray outbursting behaviour that is thought to be powered by intense internal magnetic fields. Like conventional young neutron stars in the form of radio pulsars, magnetars exhibit "glitches" during which angular momentum is believed to be transferred between the solid outer crust and the superfluid component of the inner crust. Hitherto, the several hundred observed glitches in radio pulsars and magnetars have involved a sudden spin-up of the star, due presumably to the interior superfluid rotating faster than the crust. Here we report on X-ray timing observations of the magnetar 1E 2259+586 which we show exhibited a clear "anti-glitch" -- a sudden spin down. We show that this event, like some previous magnetar spin-up glitches, was accompanied by multiple X-ray radiative changes and a significant spin-down rate change. This event, if of origin internal to the star, is unpredicted in models of neutron star spin-down and is suggestive of differential rotation in the neutron star, further supporting the need for a rethinking of glitch theory for all neutron stars. △ Less

Submitted 29 May, 2013; originally announced May 2013.

Journal ref: Nature, 497, 2013, p591-593

Hall equilibria with toroidal and poloidal fields: application to neutron stars

Authors: K. N. Gourgouliatos, A. Cumming, A. Reisenegger, C. Armaza, M. Lyutikov, J. A. Valdivia

Abstract: We present solutions for Hall equilibria applicable to neutron star crusts. Such magnetic configurations satisfy a Grad-Shafranov-type equation, which is solved analytically and numerically. The solutions presented cover a variety of configurations, from purely poloidal fields connected to an external dipole to poloidal-toroidal fields connected to an external vacuum field, or fully confined withi… ▽ More We present solutions for Hall equilibria applicable to neutron star crusts. Such magnetic configurations satisfy a Grad-Shafranov-type equation, which is solved analytically and numerically. The solutions presented cover a variety of configurations, from purely poloidal fields connected to an external dipole to poloidal-toroidal fields connected to an external vacuum field, or fully confined within the star. We find that a dipole external field should be supported by a uniformly rotating electron fluid. The energy of the toroidal magnetic field is generally found to be a few percent of the total magnetic field energy for the fields with an external component. We discuss the evolution due to Ohmic dissipation which leads to slowing down of the electron fluid. We also find that the transition from an MHD equilibrium to a state governed by Hall effect, generates spontaneously an additional toroidal field in regions where the electron fraction changes. △ Less

Submitted 26 June, 2013; v1 submitted 27 May, 2013; originally announced May 2013.

Comments: 12 pages, 9 figures, accepted by MNRAS

Hall Equilibria: Solutions with toroidal and poloidal magnetic fields in Neutron Star Crusts

Authors: K. N. Gourgouliatos, A. Cumming, M. Lyutikov, A. Reisenegger

Abstract: We present Hall equilibrium solutions for neutron stars crusts containing toroidal and poloidal magnetic field. Some simple cases are solved analytically while more complicated configurations are found numerically through a Gauss-Seidel elliptic partial differential equation solver. We present Hall equilibrium solutions for neutron stars crusts containing toroidal and poloidal magnetic field. Some simple cases are solved analytically while more complicated configurations are found numerically through a Gauss-Seidel elliptic partial differential equation solver. △ Less

Submitted 1 May, 2013; originally announced May 2013.

Comments: Conference Proceedings, Magnetic Fields in the Universe IV (2013)

Timing Noise in Pulsars and Magnetars and the Magnetospheric Moment of Inertia

Authors: David Tsang, Konstantinos N. Gourgouliatos

Abstract: We examine timing noise in both magnetars and regular pulsars, and find that there exists a component of the timing noise with strong magnetic field dependence above B ~ 10^{12.5} G. The dependence of the timing noise floor on the magnetic field is also reflected in the smallest observable glitch size. We find that magnetospheric torque variation cannot explain this component of timing noise. We c… ▽ More We examine timing noise in both magnetars and regular pulsars, and find that there exists a component of the timing noise with strong magnetic field dependence above B ~ 10^{12.5} G. The dependence of the timing noise floor on the magnetic field is also reflected in the smallest observable glitch size. We find that magnetospheric torque variation cannot explain this component of timing noise. We calculate the moment of inertia of the magnetic field outside of a neutron star and show that this timing noise component may be due to variation of this moment of inertia, and could be evidence of rapid global magnetospheric variability. △ Less

Submitted 17 July, 2013; v1 submitted 18 February, 2013; originally announced February 2013.

Comments: 5 pages, 3 figures, accepted by ApJL

Magnetic structure of Coronal Mass Ejections

Authors: Maxim Lyutikov, Konstantinos Nektarios Gourgouliatos

Abstract: We present several models of the magnetic structure of solar coronal mass ejections (CMEs). First, we model CMEs as expanding force-free magnetic structures. While keeping the internal magnetic field structure of the stationary solutions, expansion leads to complicated internal velocities and rotation, while the field structures remain force-free. Second, expansion of a CME can drive resistive d… ▽ More We present several models of the magnetic structure of solar coronal mass ejections (CMEs). First, we model CMEs as expanding force-free magnetic structures. While keeping the internal magnetic field structure of the stationary solutions, expansion leads to complicated internal velocities and rotation, while the field structures remain force-free. Second, expansion of a CME can drive resistive dissipation within the CME changing the ionization states of different ions. We fit in situ measurements of ion charge states to the resistive spheromak solutions. Finally, we consider magnetic field structures of fully confined stable magnetic clouds containing both toroidal and poloidal magnetic fields and having no surface current sheets. Expansion of such clouds may lead to sudden onset of reconnection events. △ Less

Submitted 2 October, 2012; originally announced October 2012.

Comments: proceedings, XIIth Hvar Astrophys. Colloquium