Magnetospheric Model Research Articles

Coherence of Multidimensional Pair Production Discharges in Polar Caps of Pulsars

We report on the first self-consistent multidimensional particle-in-cell numerical simulations of nonhomogeneous pair discharges in polar caps of rotation-powered pulsars. By introducing strong inhomogeneities in the initial plasma distribution in our simulations, we analyze the degree of self-consistently emerging coherence of discharges across magnetic field lines. In 2D, we study discharge evolution for a wide range of physical parameters and boundary conditions corresponding to both the absent and free escape of charged particles from the surface of a neutron star. We also present the results of the first 3D simulations of discharges in a polar cap with a distribution of the global magnetospheric current appropriate for a pulsar with 60° inclination angle. For all parameters, we find the coherence scale of pair discharges across magnetic field lines to be of the order of the gap height. We also demonstrate that the popular “spark” model of pair discharges is incompatible with the universally adopted force-free magnetosphere model: intermittent discharges fill the entire zone of the polar cap that allows pair cascades, leaving no space for discharge-free regions. Our findings disprove the key assumption of the spark model about the existence of isolated distinct discharge columns.

The Origins of Narrow Spectra of Fast Radio Bursts

Observations find that some fast radio bursts (FRBs) have extremely narrowband spectra, i.e., Δν/ν 0 ≪ 1. We show that, when the angular size of the emission region is larger than the Doppler beaming angle, the observed spectral width (Δν/ν 0) exceeds 0.58 due to the high-latitude effects for a source outside the magnetosphere, even when the spectrum in the source’s comoving frame is monochromatic. The angular size of the source for magnetospheric models of FRBs can be smaller than the Doppler beaming angle, in which case this geometric effect does not influence the observed bandwidth. We discuss various propagation effects to determine if any could transform a broad-spectrum radio pulse into a narrow spectrum signal at the observer’s location. We find that plasma lensing and scintillation can result in a narrow bandwidth in the observed spectrum. However, the likelihood of these phenomena being responsible for the narrow observed spectra with Δν/ν 0 < 0.58 in the fairly large observed sample of FRBs is exceedingly small.

Morphologies of Bright Complex Fast Radio Bursts with CHIME/FRB Voltage Data

We present the discovery of 12 apparently nonrepeating fast radio burst (FRB) sources, detected by the Canadian Hydrogen Intensity Mapping Experiment (CHIME) telescope. These sources, only one of which has been presented previously in the first CHIME/FRB catalog, were selected from a database comprising O(103) CHIME/FRB full-array raw voltage data recordings, based on their large signal-to-noise ratios and complex morphologies. Our study examines the time-frequency characteristics of these bursts, including drifting, microstructure, and periodicities. The events in this sample display a variety of unique drifting phenomenologies that deviate from the linear negative drifting phenomenon seen in many repeating FRBs, and motivate a possible new framework for classifying drifting archetypes. Additionally, we detect microstructure features of duration ≲50 μs in seven events, with some as narrow as ≃7 μs. We find no evidence of significant periodicities between subburst components. Furthermore, we report the polarization characteristics of seven events, including their polarization fractions and Faraday rotation measures (RMs). The observed ∣RM∣ values span a wide range of 17.24(2)–328.06(2) rad m−2, with apparent linear polarization fractions between 0.340(1) and 0.946(3). The morphological properties of the bursts in our sample appear broadly consistent with predictions from both relativistic shock and magnetospheric models of FRB emission, as well as propagation through discrete ionized plasma structures. We address these models and discuss how they can be tested using our improved understanding of morphological archetypes.

Multi-wavelength pulse profiles from the force-free neutron star magnetosphere

Context. The last two decades have witnessed dramatic progress in our understanding of neutron star magnetospheres thanks to force-free and particle-in-cell simulations. However, the associated particle dynamics and its emission mechanisms and locations have not been fully constrained, notably in X-rays. Aims. In this paper, we compute a full atlas of radio, X-ray, and γ-ray pulse profiles, relying on the force-free magnetosphere model. Our goal is to use such a data bank of multi-wavelength profiles to fit a substantial number of radio-loud γ-ray pulsars that have also been detected in non-thermal X-rays to decipher the X-ray radiation mechanism and sites. Using results from the third γ-ray pulsar catalogue (3PC), we investigate the statistical properties of this population. Methods. We assume that radio emission emanates from field lines rooted to the polar caps, at varying height above the surface, close to the surface, at an altitude about 5–10% of the light cylinder radius, r L. The X-ray photons are produced in the separatrix region within the magnetosphere; that is, the current sheet formed by the jump from closed to open magnetic field lines. We allow for substantial variations in emission height. The γ-rays are produced within the current sheet of the striped wind, outside the light cylinder. Results. A comprehensive set of radio, X-ray, and γ-ray light curves was computed. Based on only geometric considerations about magnetic obliquity, line-of-sight inclination, and the radio beam cone opening angle, pulsars can be classified as radio-loud or quiet and as γ-ray-loud or quiet. We found that the 3PC sample is compatible with an isotropic distribution of obliquity and line of sight. Conclusions. The atlases constructed in this work are the fundamental tools with which to explore individual pulsars and fit their multi-wavelength pulse profiles in order to constrain their magnetic topology, the emission sites, and the observer’s line of sight.

A phase-resolved Fermi-LAT analysis of the mode-changing pulsar PSR J2021+4026 shows hints of a multipolar magnetosphere

Context. The radio-quiet γ-ray pulsar PSR J2021+4026 is a peculiar Fermi-LAT pulsar showing repeated and quasi-periodic mode changes. Its γ-ray flux shows repeated variations between two states at intervals of ∼3.5 years. These events occur over timescales &lt; 100 days and are correlated with sudden changes in the spin-down rate. Multiwavelength observations also revealed an X-ray phase shift relative to the γ-ray profile for one of the events. PSR J2021+4026 is currently the only known isolated γ-ray pulsar showing significant variability, and thus it has been the object of thorough investigations. Aims. The goal of our work is to study the mode changes of PSR J2021+4026 with improved detail. By accurately characterizing variations in the γ-ray spectrum and pulse profile, we aim to relate the Fermi-LAT observations to theoretical models. We also aim to interpret the mode changes in terms of variations in the structure of a multipolar dissipative magnetosphere. Methods. We continually monitored the rotational evolution and the γ-ray flux of PSR J2021+4026 using more than 13 years of Fermi-LAT data with a binned likelihood approach. We investigated the features of the phase-resolved spectrum and pulse profile, and from these we inferred the macroscopic conductivity, the electric field parallel to the magnetic field, and the curvature radiation cutoff energy. These physical quantities are related to the spin-down rate and the γ-ray flux and therefore are relevant to the theoretical interpretation of the mode changes. We introduced a simple magnetosphere model that combines a dipole field with a strong quadrupole component. We simulated magnetic field configurations to determine the positions of the polar caps for different sets of parameters. Results. We clearly detect the previous mode changes and confirm a more recent mode change that occurred around June 2020. We provide a full set of best-fit parameters for the phase-resolved γ-ray spectrum and the pulse profile obtained in five distinct time intervals. We computed the relative variations in the best-fit parameters, finding typical flux changes between 13% and 20%. Correlations appear between the γ-ray flux and the spectral parameters, as the peak of the spectrum shifts by ∼10% toward lower energies when the flux decreases. The analysis of the pulse profile reveals that the pulsed fraction of the light curve is larger when the flux is low. Finally, the magnetosphere simulations show that some configurations could explain the observed multiwavelength variability. However, self-consistent models are required to reproduce the observed magnitudes of the mode changes.

Radio Polarization of Millisecond Pulsars with Multipolar Magnetic Fields

NICER has observed a few millisecond pulsars where the geometry of the X-ray-emitting hotspots on the neutron star have been analyzed in order to constrain the mass and radius from X-ray light-curve modeling. One example, PSR J0030 + 0451, has been shown to possibly have significant multipolar magnetic fields at the stellar surface. Using force-free simulations of the magnetosphere structure, it has been shown that the radio, X-ray, and γ-ray light curves can be modeled simultaneously with an appropriate field configuration. An even more stringent test is to compare predictions of the force-free magnetosphere model with observations of radio polarization. This paper attempts to reproduce the radio polarization of PSR J0030 + 0451 using a force-free magnetospheric solution. As a result of our modeling, we can reproduce certain features of the polarization well.

Penetrating electric field during the Nov 3–4, 2021 geomagnetic storm

We simulated the Nov 3–4, 2021 geomagnetic storm event penetrating electric field using the Multiscale Atmosphere-Geospace Environment (MAGE) model and compared with the NASA ICON observation. The ICON observation showed sudden enhancement of the vertical ion drift when the penetrating electric field arrived at the equatorial region. The MAGE model simulated vertical ion drifts have the similarly fast enhancement that shown in the ICON data at the same UT time and satellite location. Hence, ICON ion drift data was able to verify MAGE simulation, which couples the magnetospheric model was able to simulate the penetrating electric field very well.

Tests of Disk-locking in T Tauri Stars of the Orion Nebula Cluster

We look for specific correlations predicted by magnetospheric accretion models for young stars that assume disk-locking using stellar and accretion parameters derived from low-resolution HST-STIS spectra of 33 T Tauri stars in the Orion Nebula Cluster. Our results provide strong support for the correlation predicted by a model that does not constrain the stellar magnetic field to a specific geometry, while little support is found for the correlation that assumes a dipolar stellar field geometry. These results support the conclusions from similar studies of older T Tauri stars in Taurus and NGC 2264 and underscore the role that trapped flux plays in regulating a young star’s angular momentum as detailed by Ostriker & Shu. While our targets were all selected to be accreting young stars based on photometric indices, approximately half of the observed stars show no significant signs of accretion in our spectra, illustrating the difficulty in using photometric indices to find accreting stars and the possible role that variability has in the appearance of spectra in young stars. Although values of accretion luminosity derived from our models agree well with those derived from Hβ luminosity for strongly accreting stars, we find that accretion luminosity derived from Hβ luminosity is not a reliable parameter for discriminating between weakly accreting and nonaccreting T Tauri stars due to chromospheric emission that is likely present in all T Tauri stars.

Revised Magnetospheric Model Reveals Signatures of Field‐Aligned Current Systems at Mercury

AbstractMercury is the smallest and innermost planet of our solar system and has a dipole‐dominated internal magnetic field that is relatively weak, very axisymmetric and significantly offset toward north. Through the interaction with the solar wind, a magnetosphere is created. Compared to the magnetosphere of Earth, Mercury's magnetosphere is smaller and more dynamic. To understand the magnetospheric structures and processes we use in situ MESSENGER data to develop further a semi‐empiric model of the magnetospheric magnetic field, which can explain the observations and help to improve the mission planning for the BepiColombo mission en‐route to Mercury. We present this semi‐empiric KTH22‐model, a modular model to calculate the magnetic field inside the Hermean magnetosphere. Korth et al. (2015, https://doi.org/10.1002/2015JA021022, 2017, https://doi.org/10.1002/2017gl074699) published a model, which is the basis for the KTH22‐model. In this new version, the representation of the neutral sheet current magnetic field is more realistic, because it is now based on observations rather than ad‐hoc assumptions. Furthermore, a new module is added to depict the eastward ring shaped current magnetic field. These enhancements offer the possibility to improve the main field determination. In addition, analyzing the magnetic field residuals allows us to investigate the field‐aligned currents and their possible dependencies on external drivers. We see increasing currents under more disturbed conditions inside the magnetosphere, but no clear dependence on the z‐component of the interplanetary magnetic field nor on the magnetosheath plasma β.

On the Two Approaches to Incorporate Wave‐Particle Resonant Effects Into Global Test Particle Simulations

AbstractEnergetic electron dynamics in the Earth's radiation belts and near‐Earth plasma sheet are controlled by multiple processes operating on very different time scales: from storm‐time magnetic field reconfiguration on a timescale of hours to individual resonant wave‐particle interactions on a timescale of milliseconds. The most advanced models for such dynamics either include test particle simulations in electromagnetic fields from global magnetospheric models, or those that solve the Fokker‐Plank equation for long‐term effects of wave‐particle resonant interactions. The most prospective method, however, would be to combine these two classes of models, to allow the inclusion of resonant electron scattering into simulations of electron motion in global magnetospheric fields. However, there are still significant outstanding challenges that remain regarding how to incorporate the long term effects of wave‐particle interactions in test‐particle simulations. In this paper, we describe in details two approaches that incorporate electron scattering in test particle simulations: stochastic differential equation (SDE) approach and the mapping technique. Both approaches assume that wave‐particle interactions can be described as a probabilistic process that changes electron energy, pitch‐angle, and thus modifies the test particle dynamics. To compare these approaches, we model electron resonant interactions with field‐aligned whistler‐mode waves in dipole magnetic fields. This comparison shows advantages of the mapping technique in simulating the nonlinear resonant effects, but also underlines that more significant computational resources are needed for this technique in comparison with the SDE approach. We further discuss applications of both approaches in improving existing models of energetic electron dynamics.

Blandford-Znajek jets in MOdified Gravity

General relativity (GR) will be imminently challenged by upcoming experiments in the strong gravity regime, including those testing the energy extraction mechanisms for black holes. Motivated by this, we explore magnetospheric models and black hole jet emissions in Modified Gravity (MOG) scenarios. Specifically, we construct new power emitting magnetospheres in a Kerr-MOG background which are found to depend non-trivially on the MOG deformation parameter. This may allow for high-precision tests of GR. In addition, a complete set of analytic solutions for vacuum magnetic field configurations around static MOG black holes are explicitly derived, and found to comprise exclusively Heun's polynomials.

Local bow shock environment during magnetosheath jet formation: results from a hybrid-Vlasov simulation

Abstract. Magnetosheath jets are plasma structures that are characterised by enhanced dynamic pressure and/or plasma velocity. In this study, we investigate the formation of magnetosheath jets in four two-dimensional simulation runs of the global magnetospheric hybrid-Vlasov model Vlasiator. We focus on jets whose origins were not clearly determined in a previous study using the same simulations (Suni et al., 2021) to have been associated with foreshock structures of enhanced dynamic pressure and magnetic field. We find that these jets can be divided into two categories based on their direction of propagation, either predominantly antisunward or predominantly toward the flanks of the magnetosphere. As antisunward-propagating jets can potentially impact the magnetopause and have effects on the magnetosphere, understanding which foreshock and bow shock phenomena are associated with them is important. The antisunward-propagating jets have properties indistinguishable from those of the jets found in the previous study. This indicates that the antisunward jets investigated in this paper belong to the same continuum as the jets previously found to be caused by foreshock structures; however, due to the criteria used in the previous study, they did not appear in this category before. These jets together make up 86 % of all jets in this study. The flankward-propagating jets make up 14 % of all jets and are different, showing no clear association with foreshock structures and exhibiting temperature anisotropy unlike the other jets. We suggest that they could consist of quasi-perpendicular magnetosheath plasma, indicating that these jets could be associated with local turning of the shock geometry from quasi-parallel to quasi-perpendicular. This turning could be due to bow shock reformation at the oblique shock caused by foreshock ultralow-frequency (ULF) wave activity.

A global 3D simulation of magnetospheric accretion – I. Magnetically disrupted discs and surface accretion

ABSTRACT We present a 3D ideal MHD simulation of magnetospheric accretion on to a non-rotating star. The accretion process unfolds with intricate 3D structures driven by various mechanisms. First, the disc develops filaments at the magnetospheric truncation radius (RT) due to magnetic interchange instability. These filaments penetrate deep into the magnetosphere, form multiple accretion columns, and eventually impact the star at ∼30o from the poles at nearly the free-fall speed. Over 50 per cent (90 per cent) of accretion occurs on just 5 per cent (20 per cent) of the stellar surface. Secondly, the disc region outside RT develops large-scale magnetically dominated bubbles, again due to magnetic interchange instability. These bubbles orbit at a sub-Keplerian speed, persisting for a few orbits while leading to asymmetric mass ejection. The disc outflow is overall weak because of mostly closed field lines. Thirdly, magnetically supported surface accretion regions appear above the disc, resembling a magnetized disc threaded by net vertical fields, a departure from traditional magnetospheric accretion models. Stellar fields are efficiently transported into the disc region due to above instabilities, contrasting with the ‘X-wind’ model. The accretion rate on to the star remains relatively steady with a 23 per cent standard deviation. The periodogram reveals variability occurring at around 0.2 times the Keplerian frequency at RT, linked to the large-scale magnetic bubbles. The ratio of the spin-up torque to $\dot{M}(GM_*R_T)^{1/2}$ is around 0.8. Finally, after scaling the simulation, we investigate planet migration in the inner protoplanetary disc. The disc driven migration is slow in the MHD turbulent disc beyond RT, while aerodynamic drag plays a significant role in migration within RT.

Unusual Magnetospheric Dynamics During Intense Substorm Initiated by Strong Magnetospheric Compression

AbstractWe investigate an unusual sequence and peculiar features of magnetotail changes during a storm‐range substorm initiated by the interplanetary shock. Auroral observations and measurements at several favorably distributed magnetospheric spacecraft allowed the construction of an adaptive time‐dependent magnetospheric model to quantitatively characterize the configurational changes and mapping variations. Several passages of low‐altitude spacecraft in polar orbits near midnight help reveal the magnetic configuration of the nightside tail‐dipole transition region. In this event, an intense auroral and convection activity (accompanied by an up to 1,500 nT increase in the SuperMag AL‐index index) emerged in the highly compressed magnetosphere after the passage of interplanetary shock followed by strongly southward interplanetary magnetic field (IMF). This directly driven phase of the activity continued for an hour and resulted in the formation of a hybrid magnetic configuration with dipolarized midtail and stretched field lines in the transition region. Observations of energetic particle isotropy boundary latitudes near midnight are consistent with the modeled magnetic configuration. In concert with a downward turn of the solar wind (SW) flow, and weakening of the IMF driver and convection, an unusual stretching signature of the inner magnetosphere magnetic field was observed as close as at r ∼ 5 to 7 Re; which resulted mostly from an increasing downward tilt of the thin azimuthal current. Classic substorm breakup signatures commenced at fairly low, ∼60° magnetic latitude, deep in the closed field line region, in association with the current sheet upward motion. It was followed by strong stepwise poleward auroral expansion. We discuss how these signatures deviate from standard substorm scenarios and may be potentially imparted by the aforementioned changes in SW flow direction and pressure.

Plasma Properties in the Earth's Magnetosheath Near the Subsolar Magnetopause: Implications for Geocoronal Density Estimates

AbstractCombined in situ ion measurements and remote sensing of energetic neutral atoms are used to determine the geocoronal Hydrogen density at large (∼10 RE) distances from the Earth. This method for determining the geocoronal density requires global magnetospheric modeling. Observations in the Earth's subsolar magnetosheath from the Magnetospheric Multiscale mission are used to determine the accuracy of using global models to predict the geocoronal density. On average, gas dynamic and magnetohydrodynamic (MHD) models and observations are in reasonable agreement, with differences &lt;25%. In addition, the MHD model subsolar magnetopause is about 0.5 RE sunward of the observed location. However, variations around averages are large (up to a factor of 2), indicating that global models introduce relatively large uncertainties in geocoronal density estimates. Finally, the critical ion flux in the Interstellar Boundary Explorer IBEX‐Hi energy range is often minimally affected by fluctuations of a factor of 2 in the density.

Constraining the FRB mechanism from scintillation in the host galaxy

ABSTRACT Most fast radio burst (FRB) models can be divided into two groups based on the distance of the radio emission region from the central engine. The first group of models, the so-called ‘nearby’ or magnetospheric models, invoke FRB emission at distances of 109 cm or less from the central engine, while the second ‘far-away’ models involve emission from distances of 1011 cm or greater. The lateral size for the emission region for the former class of models (≲107 cm) is much smaller than the second class of models (≳109 cm). We propose that an interstellar scattering screen in the host galaxy is well-suited to differentiate between the two classes of models, particularly based on the level of modulations in the observed intensity with frequency, in the regime of strong diffractive scintillation. This is because the diffractive length scale for the host galaxy’s interstellar medium scattering screen is expected to lie between the transverse emission-region sizes for the ‘nearby’ and the ‘far-away’ class of models. Determining the strength of flux modulation caused by scintillation (scintillation modulation index) across the scintillation bandwidth (∼1/2πδts) would provide a strong constraint on the FRB radiation mechanism when the scatter broadening (δts) is shown to be from the FRB host galaxy. The scaling of the scintillation bandwidth as ∼ν4.4 may make it easier to determine the modulation index at ≳ 1 GHz.

Solving the pulsar equation using physics-informed neural networks

ABSTRACT In this study, Physics-Informed Neural Networks (PINNs) are skilfully applied to explore a diverse range of pulsar magnetospheric models, specifically focusing on axisymmetric cases. The study successfully reproduced various axisymmetric models found in the literature, including those with non-dipolar configurations, while effectively characterizing current sheet features. Energy losses in all studied models were found to exhibit reasonable similarity, differing by no more than a factor of three from the classical dipole case. This research lays the groundwork for a reliable elliptic Partial Differential Equation solver tailored for astrophysical problems. Based on these findings, we foresee that the utilization of PINNs will become the most efficient approach in modelling three-dimensional magnetospheres. This methodology shows significant potential and facilitates an effortless generalization, contributing to the advancement of our understanding of pulsar magnetospheres.

Monitoring Hα Emission from the Wide-orbit Brown-dwarf Companion FU Tau B

Monitoring mass accretion onto substellar objects provides insights into the geometry of the accretion flows. We use the Lulin One-meter Telescope to monitor Hα emission from FU Tau B, a ∼19 M Jup brown-dwarf companion at 5.″7 (719 au) from the host star, for six consecutive nights. This is the longest continuous Hα monitoring for a substellar companion near the deuterium-burning limit. We aim to investigate if accretion near the planetary regime could be rotationally modulated as suggested by magnetospheric accretion models. We find tentative evidence that Hα mildly varies on hourly and daily timescales, though our sensitivity is not sufficient to definitively establish any rotational modulation. No burst-like events are detected, implying that accretion onto FU Tau B is overall stable during the time baseline and sampling windows over which it was observed. The primary star FU Tau A also exhibits Hα variations over timescales from minutes to days. This program highlights the potential of monitoring accretion onto substellar objects with small telescopes.

The Gamma-Ray Pulsar Phenomenology in View of 3D Kinetic Global Magnetosphere Models

We develop kinetic plasma models of pulsar magnetospheres with magnetic-field-line-dependent plasma injection that reveal the importance of various magnetosphere regions in regulating the γ-ray emission. We set different particle injection rates for the so-called open, closed, and separatrix zones. Moderate particle injection rates in open and closed zones ensure a global field structure close to the force-free structure, while the dissipation occurs mainly in and around the equatorial current sheet. The particles that are injected into the separatrix zone affect the particle populations that enter the equatorial current sheet region, and therefore the corresponding accelerating electric fields, particle energies, the spectral cutoff energy, and γ-ray efficiency. The separatrix zone models reproduce the recently discovered fundamental plane of γ-ray pulsars consistent with curvature radiation emission, the γ-ray light-curve shapes, and the radio-lag versus peak-separation correlation reported in the Fermi second pulsar catalog. The model beaming factors indicate that the pulsar total γ-ray luminosities listed in the Fermi catalogs are overestimations of the actual ones. We find that the radiation reaction limited regime starts ceasing to govern the high-energy emission for . Our results also indicate that toward high magnetic inclination angles, the Y point around the rotational equator migrates well inside the light cylinder, sparking additional peaks in the γ-ray pulse profiles. We find that an equivalent enhanced particle injection beyond the Y point strengthens these features, making the model γ-ray light curves inconsistent with those observed.

Energetic Electron Flux Dropouts Measured by ELFIN in the Ionospheric Projection of the Plasma Sheet

AbstractLow‐altitude observations of magnetospheric particles provide a unique opportunity for remote probing of the magnetospheric and plasma states during active times. We present the first statistical analysis of a specific pattern in such observations, energetic electron flux dropouts in the low‐altitude projection of the plasma sheet. Using 3.5 years of data from the ELFIN CubeSats we report the occurrence distribution of 145 energetic electron flux dropout events and identify characteristics, including their prevalence in the dusk and premidnight sectors, their association with substorms and enhanced auroral activities, and their correlation with the region‐1 (R1) field‐aligned current region. We also investigate three representative dropout events which benefit from satellite conjunctions between ELFIN, GOES, and THEMIS, to better understand the magnetospheric drivers and magnetic field conditions that lead to such dropouts as viewed by ELFIN. One class of dropouts may be associated with magnetic field mapping distortions due to local enhancements and thinning of cross‐tail current sheets and amplification of R1 field‐aligned currents. The other class may be associated with the increase in perpendicular anisotropy of magnetospheric electrons due to magnetic field dipolarizations near premidnight. These plasma sheet flux dropouts at ELFIN provide a valuable tool for refining magnetospheric models, thereby improving the accuracy of field‐line mapping during substorms.