A new solution of the pulsar equation

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.


The Pulsar Equation
Right after the discovery of pulsars in 1968 [1], Goldreich and Julian were the first to sketch the general structure of the pulsar magnetosphere [2].The detailed structure of the axisymmetric pulsar magnetosphere was obtained much later by Contopoulos, Kazanas and Fendt, who clearly identified the presence of a global magnetospheric poloidal electric current sheet [3].Since then, several authors have addressed the pulsar magnetosphere problem with force-free electrodynamic (FFE) simulations ( [4,5], etc.) magnetohydrodynamic (MHD) simulations ( [6], etc.) 'ab initio' particle-in-cell (PIC) simulations ( [7][8][9], etc.) and most recently with Machine Learning [10,11].Several inconsistencies in recent stateof-the-art numerical simulations (e.g., the extent of the co-rotating closed-line region, the microscopic treatment and the thickness of the current sheet, magnetospheric dissipation, etc.) led us to believe that one can only trust them qualitatively (not quantitatively) to make meaningful comparisons with observations.This is why, in [11], we proposed to return to the basics and obtain the reference ideal steady-state force-free magnetosphere in a novel independent way, namely with Machine Learning.
The magnetospheres of neutron stars are dominated by the electric and magnetic fields.The physical conditions allow us to neglect gravity, thermal pressure and particle inertia because they are several orders of magnitude smaller than the electromagnetic forces.Therefore, force balance in the bulk of the pulsar magnetosphere is reduced to where E and B are the electric and magnetic fields, respectively, ρ e = ∇ • E/(4π) is the electric charge density, and J is the electric current density.Under axisymmetric conditions, one may define the magnetic flux function Ψ such that B r = ∂Ψ/∂θ/(r 2 sin θ) and B θ = −∂Ψ/∂r/(r sin θ), which automatically satisfies the condition ∇ • B = 0.Under steady-state axisymmetric force-free ideal conditions, Equation (1) becomes the well-known pulsar equation, namely in spherical coordinates r, θ [12].Here, I = I(Ψ) is the distribution of poloidal electric current from which one obtains the azimuthal component of the magnetic field B ϕ = I/(r sin θ), and I ′ ≡ dI/dΨ.Equation ( 2) is an elliptical partial differential equation with a singularity along the light cylinder, the nominal distance where the co-rotational velocity becomes equal to the speed of light (R LC ≡ c/Ω, where Ω is the angular velocity of stellar rotation).This equation was first solved with a special numerical technique developed by [3] that was used again since then by several others (e.g., [13][14][15]).
In Section 2, we present a new type of solution of Equation ( 2) that is very different from the canonical solutions found in the literature.In Section 3, we obtain new solutions of the ideal force-free pulsar magnetosphere for various amounts of open magnetic flux using the novel Machine Learning methodology of [11] (hereafter Paper I; see description in Appendix A).We show that the new type of solution is the limit of a continuous sequence of standard magnetospheres with open and closed field lines in which the closed line region approaches closer and closer to the light cylinder and the amount of open field lines gradually drops to zero.In Section 4 we discuss the possibility that an abrupt transition from our new magnetospheric solution to a standard magnetosphere with open and closed field lines may generate a fast radio burst.We conclude in Section 5 with a discussion of our results.We consider only axisymmetric magnetospheres, and in future work, we will investigate whether our results apply also to an oblique rotator.

A Confined Magnetosphere
Let us consider here a theoretical limit that has never been considered before, namely one with no open field lines and no poloidal electric current (I = 0).In that limit, the whole magnetosphere of the central dipole is contained within the light cylinder, with boundary conditions Ψ(r, θ = 0) = Ψ(r, θ = π) = 0, a dipole field Ψ(r * , θ) = Ψ max sin 2 θ along the surface of the central star (Ψ max is the total dipolar magnetic flux that emanates from each stellar hemisphere), and Ψ(r = R LC / sin θ, θ) = 0 along the light cylinder.The solution of Equation ( 2) with these boundary conditions is obtained with the method described in [16] and is shown in Figure 1.This configuration co-rotates with the central star and is wrapped around by an electrically charged layer at the light cylinder.It is a unique feature of special relativity that allows such a spatially confined electromagnetic field solution, namely the presence of electric charges and electric fields that are generated by the rigid co-rotation of the magnetosphere.In the region next to the light cylinder in particular, the electric charge density is positive/negative and generates (via its co-rotation) an azimuthal electric current along/opposite to the direction of rotation for an aligned/counter-aligned rotator, respectively.The magnetic field points downwards/upwards and the electric field outwards/inwards, respectively, so it is obvious that the inward Lorentz force J × B/c is balanced by the outward electrostatic force ρ e E.
It is, however, interesting that such a confined solution with a vacuum outside is not possible without a light cylinder, as is the case in non-relativistic MHD.In particular, pressure balance across the electrically charged current sheet along the light cylinder may be written as or equivalently The requirement for the continuity of B 2 − E 2 across an infinitely thin relativistic current sheet, Equation (3), stems from the integration of the force-free equation ρ e E + J × B/c = 0 across it [17,18].Here, x ≡ R/R LC is the cylindrical radius in units of the light cylinder radius R LC , and x| IN → 1 − .It is assumed that the whole magnetosphere is confined inside the light cylinder, it remains untwisted (i.e., B ϕ = 0), it co-rotates with the central star, and there is nothing outside.Obviously, in the limit x| IN = 1 − , pressure balance (i.e., Equation ( 4)) is satisfied; thus, solutions spatially confined inside the light cylinder are possible.In this new solution, the value of B z | IN is non-zero and can be determined numerically.We emphasize once again that such spatially self-confined solutions with a vacuum outside are possible only along the light cylinder, not inside or outside.We propose that this singular configuration may be realized in nature under special conditions during the stellar collapse that led to the formation of a spinning neutron star in a supernova explosion.The initial neutron star magnetosphere will be squeezed to a thin layer above the surface of the neutron star and will rotate with the neutron star.The solutions of the pulsar equation shown in Figure 2 were obtained with the same method by setting Ψ = 0 along a sphere of radius r ejecta , and along the light cylinder wherever r ejecta sin θ > R LC .Thus, the whole magnetosphere is confined between the neutron star surface and a spherical shell of supernova ejecta.As the supernova explosion proceeds and the ejecta expand, the size r ejecta of the magnetosphere will increase.During that stage, the whole magnetosphere co-rotates with the central neutron star; thus, it will be confined inside the light cylinder wherever r ejecta sin θ > R LC .It is implied here that, during the supernova explosion, as r ejecta grows slowly at subluminal velocities, the magnetosphere evolves fast into a sequence of steady-state solutions.When r ejecta grows beyond r ejecta = R LC , any part of the magnetosphere that crosses the light cylinder cannot open up to infinity as in the magnetosphere of an isolated pulsar.Therefore, it will be wound very fast by the stellar rotation, the Lorentz force will squeeze it vertically, and it will eventually detach from the rest of the magnetosphere that remains confined inside the light cylinder as in the solutions shown in Figure 2. It is obvious that this is an idealized situation in which a neutron star and its magnetosphere formed and remained at the center of the collapse (in general, the crust of a young neutron star may not have formed/solidified yet, and the neutron star would also potentially experience a kick that will move it away from the center of the collapse).It is also conceivable that such a confined configuration is attained when a neutron star accretes from a surrounding disk or from external winds from a donor star (e.g., [19,20]).It is also obvious that the solutions shown in Figures 1 and  2 are force-free approximations since they were obtained without taking into consideration the baryon-rich and dense environment of a newborn neutron star.The idealized solution shown in Figure 1 is most probably unstable.We have not performed a formal stability analysis of the light cylinder boundary, but we understand this instability as follows: Whenever a part of the confined magnetosphere extends beyond the light cylinder, it cannot keep co-rotating with the central star without an azimuthal component of the magnetic field.In other words, when an initially poloidal magnetic field line crosses the light cylinder, it is swept back so that the particle drift velocity always remains subluminal.Without the azimuthal component B ϕ , v D > c outside the light cylinder, where x > 1.Therefore, any part of the confined magnetosphere that crosses the light cylinder will develop an azimuthal magnetic field component that will push it outwards toward infinity.Thus, the magnetosphere will transition to a standard solution with open and closed field lines, emitting an electromagnetic pulse of azimuthal magnetic field along the way.A similar effect is seen in numerical simulations during the evolution of initially poloidal magnetic field configurations when the stellar rotation is initiated (e.g., [4]).
We will return to a further discussion of this burst in Section 3 below.Notice that using Equation ( 5) with B ϕ = 0, the drift velocity is exactly equal to c at the light cylinder, and the neglect of inertial terms is more questionable for the confined solution than for the standard one near the light cylinder.Inertial effects are even more pronounced near the light cylinder because, due to the (1 − x 2 ) term in the denominator of the Goldreich-Julian charge density , the volume integral of ρ e and the corresponding mass content of the idealized confined magnetosphere diverge.This suggests that, if realized in nature, the confined solution will be limited to within some short distance from the light cylinder.

A Spectrum of Solutions
We will now show that the new solution obtained in the previous section is the limit of a continuous spectrum of solutions of the standard type, namely solutions that consist of a co-rotating closed-line region and an open-line region that extends to infinity.Figure 1 is the mathematical limit of a sequence of independent standard solutions, and does not imply a physical evolution of standard solutions.Nevertheless, each solution contains a different amount of electromagnetic energy, so magnetospheric transitions from a higher to a lower energy solution are possible (see below).
We obtain here magnetospheric solutions with the methodology of Paper I (see Appendix A) in which the separatrix between open and closed field lines is considered to be a mathematical contact discontinuity, namely a surface of zero thickness.The tip of the closed-line region is called the Y-point.We consider in particular Y-point positions from x Y ≈ r * /R LC ≪ 1 to x Y ≈ 1 very close to the light cylinder.For reasons of computational convenience, we chose a rather large star with radius r * = 0.25R LC (PINN methods cannot handle well large scale differences in their domain of application).We have also tried different (smaller) values of r * and obtained similar results.As is described in Appendix A, we first choose the angular opening of the polar cap θ pc from which a certain amount of magnetic flux Ψ open ≡ Ψ S = Ψ max sin 2 θ pc emanates and extends to infinity.This forms the so-called open-line region.In the solutions shown in Figure 3, one particular value of θ pc , namely θ pc = 1.176(r * /R LC ) 1/2 , was specially chosen with the following in mind: it is straightforward to calculate that, in a dipolar magnetic field configuration, the magnetic field line that crosses the light cylinder corresponds to * Ω/(2c) (here, B * is the polar value of the surface magnetic field), and In previous high-resolution solutions of the pulsar equation (e.g., [14]), the magnetic field line that crosses the light cylinder corresponds to Ψ LC Timokhin = 1.23Ψ dipole LC and θ pc Timokhin ≈ 1.176(r * /R LC ) 1/2 .Therefore, this particular value of θ pc was chosen so that the Y-point lies very close to the light cylinder.Instead, in Paper I, it was found that R Y = 0.88R LC , which is closer to R dipole (Ψ S ) = r * / sin 2 θ pc = 0.81R LC .This unexpected result is due to our clear treatment of the separatrix surface as a surface of zero thickness which allows the manifestation of the T-point at its tip (see Appendix B).This is not the case in all previous MHD, FFE, and PIC simulations where the separatrix surface has an unphysically large finite thickness which at its tip leads to an exterior Y-point, and a T-point some distance inside it (see Figure 3 of [21] for details).
We also obtained the solution for Ψ open = 0.866Ψ dipole LC that corresponds to θ pc = 0.97(r * /R LC ) 1/2 and R Y = 0.94R LC .The reader can check that this value of Ψ open yields a pulsar spindown rate equal to Ė = 0.75 Ėvacuum (90 • ).Here, Ėvacuum (λ) is the spin- down rate of a vacuum dipole rotator with inclination angle λ, and Ėvacuum (90 One may tentatively generalize our result for nonzero pulsar inclination angles according to [4] as Ė(λ) ≈ 0.75 Ėvacuum (90 • )(1 + sin 2 λ), and since Ėvacuum (λ) = Ėvacuum (90 • ) sin 2 λ, we obtain that It is interesting that in all previous solutions of the FFE pulsar magnetosphere, the above ratio was found to be greater than 3 (e.g., [22]).This value is significantly larger than the ratio of spindown rates ĖON / ĖOFF observed in the intermittent pulsars PSR B1931+24, PSR J1832+0029 and PSR J1841-0500 for their corresponding 'ON' and 'OFF' states (1.5, 1.7 and 2.5, respectively, e.g., [23,24]).The inability to account for observed values lower than 3 was the reason that led to the development of resistive magnetospheric solutions (e.g., [22,25]).With our new solutions, it seems that there is no need for magnetospheric resistivity anymore.This result certainly merits further investigation.We were able to apply our method to even smaller values of Ψ open and θ pc .These solutions are very hard to obtain with the methodology of Paper I because in these, the Y-point approaches very close to the light cylinder and the convergence of the solution in the open-line region between the Y-point and the light cylinder is very difficult.In the limit Ψ open = θ pc = 0, one reaches the solution shown in Figure 1, which was obtained independently as a solution of the pulsar equation (eq.2) with standard methods (see [16]).This spectrum of solutions with smaller and smaller spindown rates is new and has never been obtained before in previous MHD, FFE, nor PIC simulations.Our results are summarized in Table 1.In Figure 4, the distribution of B z is plotted along the equator inside the light cylinder and shows its clear increase right inside the Y-point at its tip, as is needed in order to satisfy Equation (A2).Although this effect was known, it has never been observed so clearly before, and Figure 4 improves Figure 11 of [14].It is now seen clearly that, with the improved methodology of [11], the minimum B z value of this sequence of solutions lies at around x Y ∼ 0.8, which is near the position of the Y-point in the high-resolution PIC simulations of [26][27][28][29].It may be that the Y-point naturally relaxes to that position (not at the light cylinder), as this corresponds to a minimum energy solution.This result may answer an important question in the recent literature on the pulsar magnetosphere.Note, however, that the numerical experiments in [30] lead to a different interpretation, namely that the Y-point can be pushed arbitrarily close to the stellar surface depending on the number of PIC particles that are either injected or reach the dissipative region beyond the Y-point.This is an interesting interpretation since it suggests that x Y may decrease as the magnetosphere evolves with pulsar spindown.This will lead to a value of the braking index smaller than its canonical value of 3. Otherwise, if x Y remains unchanged, all solutions obtained in Table 1 will have braking indices equal to 3.
In Figure 5, the variation in the open magnetic flux Ψ open is plotted in units of Ψ dipole LC with the position of the Y-point.Our present result that Ψ open → 0 as x Y → 1 is in tension with all previous solutions of the pulsar magnetosphere.The dotted blue line corresponds to the fit 1.1 Ψ open /Ψ dipole LC /x Y , while the grey band corresponds to the range of values obtained from the simulations of [14].For example, solution No 4 from Table 1 contains the same amount of open magnetic flux as the one obtained by [27], yet their Y-point lies at x Y = 0.85, while ours lies at x Y = 0.71.We believe that this pronounced difference is mainly due to our treatment of the separatrix current sheet as a contact discontinuity.In all previous solutions (e.g., [3,4,14,27,28]), the separatrix contains a finite amount of poloidal magnetic flux and has a pronounced nonzero thickness.This affects the pressure balance across it which significantly modifies the position of the magnetospheric Y-point.The difference may also be partially due to our choice of a rather large stellar radius r * = R LC /4 compared to previous solutions.Finally, in Figure 6 the variation in the pulsar spindown energy loss Ė is plotted in units of Ėvacuum (90 • ) with the position of the Y-point.Our present result that Ė → 0 as x Y → 1 is also in tension with all previous solutions of the pulsar magnetosphere.The grey band corresponds to the spindown energy losses obtained from the simulations of [3,14,31].As is acknowledged above, the difference between the dotted line and the grey band may be partially due to our choice of a rather large stellar radius.The value shown with the horizontal line corresponds to the particular value Ė = 0.75 Ėvacuum (90 • ) mentioned in Equation ( 6).  1. B z is in arbitrary units.The dotted line represents the standard 1/r 3 dipolar dependence.One clearly infers both the outward stretching of the initially dipolar field (at small r distances the solid lines lie below the dotted one), and the inwards squeezing of the field right behind the Y-point that is required to satisfy pressure balance.Never before has anyone seen such a large increase in B z inside the Y-point, even in the solutions in the literature where x Y → 1.Our results improve Figure 11 of [14].

A Possible Connection with Fast Radio Bursts
Time-dependent numerical simulations of the pulsar magnetosphere start with an initial configuration B(r, θ, ϕ; t = 0) = B dipole (r, θ, ϕ; t = 0) that corresponds to a magnetostatic dipole in vacuum, and E(r, θ, ϕ; t = 0) = 0 everywhere.At time t = 0, the programmer introduces electric fields E(r * , θ, ϕ; t) = −r sin θ Ω/c φ × B dipole (r * , θ, ϕ; t) along the surface of the star, thus effectively setting the central star in rotation with angular velocity Ω along the z-axis.Here, r * is the radius of the star, and B dipole (r * , θ, ϕ; t) corresponds to a dipole field that rotates with angular velocity Ω along the z-axis.This is the procedure followed in the pioneering paper of Spitkovsky [4] and in all time-dependent numerical simulations since then.The abrupt introduction of a non-zero electric field along the surface of the star generates a transient blast wave of poloidal electric and azimuthal magnetic fields that sweeps through the static magnetosphere.This is analogous to the transient blast wave of electric and magnetic fields generated when an electric charge is abruptly set in motion [32] or an electric current is introduced [16].The aim of such simulations is to obtain the steady-state solution; thus, the programmer is not interested in the blast wave and waits for it to leave the inner magnetosphere.Indeed, after that wave sweeps through the computational grid that corresponds to the inner magnetosphere, a breathing steady-state solution is established with closed and open magnetic field lines.That initial artificial polarized blast wave has the characteristics of a fast radio burst (hereafter FRB).
FRBs represent one of the most intriguing enigmas in the field of astrophysics today.First discovered in 2007 [33], FRBs are transient radio pulses that originate from distant galaxies, characterized by their extremely high energy and incredibly short duration, typically lasting just a few milliseconds.Despite their fleeting nature, these cosmic phenomena release more energy in a fraction of a second than the sun in an entire day.The exact mechanisms that produce FRBs remain a subject of intense debate and speculation among scientists.Hypotheses range from highly magnetized neutron stars, known as magnetars, to more exotic theories involving cosmic strings or even extraterrestrial intelligence.However, despite over a decade of observation and research, no single theory has been able to fully explain all observed properties of FRBs.
We propose here that FRBs may be related to the electromagnetic burst of the toroidal magnetic field that follows an abrupt transition from the new solution that was obtained in this paper to one of the standard magnetospheric solutions of Figure 1 and Table 1.The duration of the burst will be on the order of the light crossing time of the light cylinder, namely R LC /c = P/(2π), where P is the spin period of the pulsar.The total energy of the burst is estimated to be equal to the difference between the electromagnetic energy of the new solution and the electromagnetic energy inside the light cylinder of a standard magnetospheric solution.An estimate of the former is obtained as follows: according to Figures 1 and 4, the dipolar magnetic flux Ψ dipole LC that would have extended beyond the light cylinder, is now confined to close within an annular region of area ∼ R LC δ where δ ∼ 0.1R LC .That region has volume ∼ R 2 LC δ ∼ 0.1R 3 LC , and in it, An estimate of the latter is Hence, the energy of the burst will be equal to the difference between the two, namely Notice that in the energy estimates of this section, we use a physical neutron star radius of r * ≈ 10 km.
As we will now see, this is more than enough to power an FRB.In the case of an accreting recycled millisecond pulsar with a surface magnetic field on the order of B * ∼ 10 8−9 G, where L ⊙ is the solar luminosity.This value is below the lowest limit of the low-energy FRBs presented by [34]; hence, it does not seem consistent with actual FRB observations.For a newborn neutron star with a millisecond period formed via a core collapse supernova, Actual FRB energies lie between the above extreme limits, namely between 10 37 and 10 41 erg (e.g., [35,36]).Therefore, we propose that FRBs correspond to millisecond pulsars with surface magnetic fields between 10 11 and 10 12 G.

Discussion and Conclusions
We have improved the solutions of the pulsar equation when the tip of the closedline region, the so-called Y-point, approaches the light cylinder.The new solutions were obtained under the assumption that the separatrix surface between open and closed field lines is a zero-thickness contact discontinuity.The new solutions are similar to the old ones of [3,14,31] when the Y-point lies inside about 80% of the light cylinder radius.As the Y-point approaches closer to the light cylinder, however, the solutions diverge from the standard ones found in the literature, and in the extreme limit that the closed-line region touches the light cylinder, the open-line region disappears completely.As is seen in Equation (A1) in Appendix B and in Figure 4, the closer the Y-point to the light cylinder, the stronger the divergence of B p right inside the Y-point at the tip of the closed-line region needs to be.We see no other way to implement this divergence of B p (x Y → 1 − ) except with more and more field lines to enter the closed-line region and compress the magnetic field at its tip.This effect is clearly seen in the sequence of solutions shown in Figure 3.In the limit where B p → ∞ as x Y = 1 − , all field lines enter the closed-line region confined inside the light cylinder, and the open line region disappears completely (Figure 1).
One may legitimately ask why these solutions were never before obtained in the literature.There are several reasons.Firstly, in all previous solutions, the separatrix surface along which the magnetospheric electric current circuit closes on the star had non-zero thickness due to the finite resolution of the code.In particular, at its tip, the outer part of the separatrix forms a Y-point that touches the light cylinder, while the inner part forms a clear T-point some distance inside it.This is clearly explained in Figure 4 of [21].Unfortunately, the ideal force-free conditions of Equation ( 2) are not valid inside current sheets (as, e.g., in Figure 1b of [4]).Nevertheless, the fact that PIC simulations also form a separatrix that contains a nonzero amount of poloidal magnetic flux (see high-resolution detail in Figure 9 of [27]) led us to consider the possibility that this effect is indeed physical.It is interesting that the separatrix current sheet contains a nonzero guide field and is thus different from the equatorial current sheet which is an ordinary Harris-type current sheet with 180 • field inversion across it.It also seems that the width of the former type of current sheet grows with time as is seen in Figure 1c of [4].This effect merits further investigation with hybrid PIC-FFE simulations (Soudais, Cerutti & Contopoulos 2024, submitted).Now that the separatrix has zero thickness, it forms a clear T-point at its tip, which indeed lies some distance inside the light cylinder.Secondly, time-dependent simulations always start from a vacuum dipole configuration that extends to all space and therefore, some fraction of it already crosses the light cylinder.After the central star is set into rotation, the magnetosphere is divided into open and closed field lines, i.e., there is no chance to obtain a configuration totally enclosed inside the light cylinder as the solution shown in Figure 1.Thirdly, after the time-dependent evolution relaxes to a steady-state solution, it is reasonable to expect that it will further relax to the lowest energy solution.Therefore, time-dependent simulations cannot yield the full sequence of solutions up to the latter singular solution as our Machine Learning methodology can.It is interesting that, now that the full sequence of solutions is obtained, it is shown that the solution where x Y ∼ 80% R LC is a minimum energy solution, and this is why previous time-dependent high-resolution PIC simulations relax to that solution.We have thus improved upon the result of [21] which placed the minimum energy solution around x Y = 92% R LC .
We argue that the confined singular solution of Figure 1 and all solutions between that one and the minimum energy solution are unstable; thus, if the magnetospheric configuration of Figure 1 forms during a supernova explosion or during a high accretion event from a surrounding disk or wind from a donor star, it will transition to the minimum energy solution by emitting a burst of energy equal to the energy difference between these two solutions.We propose that this burst may be related to an FRB.We would like to notice here that, as long as the surrounding material is there, the confined solution cannot transition to a standard solution in which part of the magnetosphere extends to infinity, and thus the accumulated electromagnetic energy of the confined solution cannot be released.What happens when the surrounding material dissolves and how fast electromagnetic energy is released remain to be seen in a time-dependent numerical simulation of this transition.In case the pulsar magnetosphere undergoes a later accretion event from an external disk or wind from a donor star that will re-confine it to the solution of Figure 1, the FRB will repeat.We must acknowledge of course that a few FRBs have already been observed from a magnetar, and that the cosmological sky rate of FRBs exceeds that of core-collapse supernovae; therefore, our model may just be a rare kind of FRB, if feasible at all.Moreover, we have not answered why would the emission of this state transition be in the MHz to GHz radio regime, and how it would escape the very dense accretion environment that formed the confined solution of Figure 1.Our model certainly merits further investigation with detailed MHD numerical simulations.
What is also interesting to investigate is whether a confined solution exists for nonaligned (oblique) pulsars in which ∂B/∂t ̸ = 0. We expect that oblique rotators would radiate electromagnetic radiation to infinity; thus, magnetic field lines would extend to infinity.Therefore, it seems impossible for an oblique magnetosphere to be totally confined inside the light cylinder.It may thus be that the confined solution of the type shown in Figure 1 is only possible in aligned pulsars.

Appendix A
In the first paper of this series ( [11], Paper I), a new methodology was presented for solving the pulsar equation which essentially consists of the following steps: 1.
The angular size θ pc of the polar cap that contains the so-called open magnetic field line is chosen a priori.The rest of the stellar surface contains closed field lines.

2.
The two regions of open and closed magnetic field lines are considered independently.An 'infinitely thin' physical discontinuity exists between the two regions in the form of a separatrix current sheet.In reality, the separatrix current sheet has some finite thickness which is orders of magnitude smaller than the characteristic dimensions of the system such as the radius of the star and the radius of the light cylinder.

3.
It is very helpful to deal with the equatorial current sheet in the open line region by mathematically reversing the polarity of the magnetic field that emanates from the south pole of the star.By carrying this out, the mathematical discontinuity in the equatorial region disappears; thus, one does not need to worry about it computationally.At the end of the calculation, we will return to the original polarity of the magnetic field, and the equatorial current sheet will reappear.

4.
The force-free problem is solved in the open and closed line regions independently.
Our method of choice is with Neural Networks so-called Physics Inspired (PINNs).One may also use classical grid methods in the two regions.5.
The shape r S (θ) of the separatrix is readjusted so that continuity of B 2 − E 2 is achieved at all points across it.One way to carry that out is to move the separatrix radially at each point with an amount that is proportional to the difference (B Repeat the procedure until pressure balance is achieved over the whole separatrix and its shape stabilizes.
In [11], Machine Learning was implemented in the PyTorch Deep Learning library.Code development was implemented in Python on Anaconda Jupyter Notebooks running on local GPUs, while production runs were performed in the Cloud in Google's Colab.Furthermore, Adam optimizers, SiLU activation functions, and ReduceLROnPlateau schedulers from the PyTorch library were used.Three independent NNs were trained: 1.
One PINN with two entries (the r, θ spherical coordinates), three hidden layers with 64 nodes each, and two exits Ψ, I.This PINN solves the pulsar equation in the open line region outside the separatrix.2.
One PINN with two entries r, θ, three hidden layers with 64 nodes each, and one exit Ψ.This PINN solves the pulsar equation in the closed line region inside the separatrix where I = 0.

3.
A third NN with one entry θ, two hidden layers with 128 nodes each, and one exit r S is trained to yield the shape r S (θ) of the separatrix at all angles θ.The number of internal nodes and layers in the two PINNs that solve the pulsar equation was chosen by trial and error to be able to reproduce known features of the solution.In practice, the two main PINNs are initially trained for 50,000 steps after which, satisfactory convergence is achieved in the two regions.After that initial training stage, the separatrix is displaced based on the resulting pressure differences between the closed and open line regions, and then the third NN is trained.The two main PINNs are re-trained for another 50,000 steps, and the process is repeated 10 times.In total, we run 500,000 training steps for the two main PINNs.In the resulting configuration, pressure balance is achieved across the separatrix to less than 1%.
As it is acknowledged in [11], the numerical accuracy of our method is not on par yet with that of standard methods such as finite difference, finite volume or spectral methods.This is also the case in the solution of [37].On the other hand, the numerical treatment of current sheets by standard methods it too is problematic (numerical current sheets have unphysical thicknesses and the numerical dissipation there is questionable).In order to improve our results in future iterations of this work, we will try to combine standard methods outside current sheets with our methodology of moving the separatrix.This is still a work in progress.Modulo the above disclaimers, we applied this method to obtain the sequence of solutions shown in Figure 3. satisfy Equation (A1) is to have B p IN ̸ = 0.This can only be satisfied inside a T-point, not inside a Y-point.Nevertheless, we will continue calling the tip of the closed-line region the Y-point because of its overall bifurcating shape, and we will denote the rescaled radial distance of the tip of the closed-line region with x Y .
Furthermore, it is obvious that Equation (A1) becomes singular when x Y approaches very close to unity.In that limit, Equation (A1) yields As is shown in this paper, contrary to what was obtained in previous works, it is not possible to reach this limit and have a significant open-line region outside.The closer the Y-point approaches the light cylinder, the more field lines must be squeezed inside the closed-line region to support the growth of B p at its tip required by Equation (A2), the fewer field lines are left in the open-line region.As is shown with the new confined magnetospheric solution of Figure 1, in the limit x Y → 1 − the open-line region disappears completely.

Figure 2 .
Figure 2. Sequence of co-rotating magnetospheric solutions confined by an expanding shell of plasma material (grey region).As the ejecta expand to infinity, the solution approaches the solution shown in Figure 1.

Figure 4 .
Figure 4.The distribution of B z (r, θ = π/2) inside the Y-point for the various solutions listed in Table1.B z is in arbitrary units.The dotted line represents the standard 1/r 3 dipolar dependence.One clearly infers both the outward stretching of the initially dipolar field (at small r distances the solid lines lie below the dotted one), and the inwards squeezing of the field right behind the Y-point that is required to satisfy pressure balance.Never before has anyone seen such a large increase in B z inside the Y-point, even in the solutions in the literature where x Y → 1.Our results improve Figure11of[14].