General features of the linear crystalline morphology of accretion disks

In this paper, we analyze the so-called Master Equation of the linear backreaction of a plasma disk in the central object magnetic field, when small scale ripples are considered. This study allows to single out two relevant physical properties of the linear disk backreaction: (i) the appearance of a vertical growth of the magnetic flux perturbations; (ii) the emergence of sequence of magnetic field O-points, crucial for the triggering of local plasma instabilities. We first analyze a general Fourier approach to the solution of the addressed linear partial differential problem. This technique allows to show how the vertical gradient of the backreaction is, in general, inverted with respect to the background one. Instead, the fundamental harmonic solution constitutes a specific exception for which the background and the perturbed profiles are both decaying. Then, we study the linear partial differential system from the point of view of a general variable separation method. The obtained profile describes the crystalline behavior of the disk. Using a simple rescaling, the governing equation is reduced to the second order differential Whittaker equation. The zeros of the radial magnetic field are found by using the solution written in terms Kummer functions. The possible implications of the obtained morphology of the disk magnetic profile are then discussed in view of the jet formation.


Introduction
The Standard Model of accretion in astrophysics [1,2] is based on the idea that the angular momentum transport across the plasma disk is ensured by effective viscosity, due to the turbulence phenomena triggered by the magnetorotational instability (MRI) [3], see also Refs. [4,5,6]. However, the necessity to account for a magnetic field, due to the central object, together with the long-wavelength response of the plasma disk, implies that the azimuthal component of the generalized Ohm law is consistent only postulating an anomalous plasma resistivity, see Refs. [7,8]. This effective resistivity makes the disk magnetic field limited in amplitude, because of its diffusive nature [7]. As a result, the possibility to deal with the high magnetic field strength required for the jet generation leads to investigate specific scenarios [9].
In Refs. [10,11] (see also Ref. [12]), it has been investigated the possibility for a small scale backreaction of the plasma disk to the magnetic field of the central object, outlining the emergence of a crystalline morphology of the magnetic field microscales, i.e., a local radial oscillation of the magnetic flux function. In particular, in the limit of a linear backreaction, in Ref. [10] it has been written down a Master Equation for the magnetic flux function, which is the starting point of the present investigation. The relevance of this new paradigm relies on the possibility to use the O-points of the resulting configuration as the sites where the jet formation could start, as discussed in Refs. [13,14].
Here, we develop a detailed general analysis of the linear Master Equation for the crystalline profile of the plasma disk, outlining two main significant features: i) performing a general Fourier analysis, we clarify how the decreasing behavior of the plasma backreaction with the vertical quote (discussed in Refs. [10,11], see also Ref. [15]) is actually a peculiar property of the fundamental harmonics in the radial oscillation; ii) by analyzing the Master Equation in the framework of a separable solution, we clarify the existence of regions of the parameter space, where a relevant sequence of O-points of the magnetic configuration can take place.
Both the results mentioned above have a precise physical implication. The point i) shows how the higher order harmonics of the radial oscillation of the flux function are actually associated to an inverse behavior of the vertical gradient between the background and the linear perturbation. This issue calls attention for its possible implications on the morphology of the nonlinear backreaction which is associated with the fragmentation of the disk in a ring-like morphology, see Ref. [11] and also Refs. [16,17,18]. The point ii) demonstrates that the topology of the induced radial component of the magnetic field can be characterized by a sequence of O-points, where such a radial component vanishes (noting that the dipole-like magnetic field of the central object is essentially vertical on the equatorial plane). This feature, absent in the analysis in Refs. [10,11], suggests that the magnetic micro-structure, emerging when the perturbation scale is sufficiently small (i.e., for sufficiently high values of the plasma β parameter), could be an interesting mechanism in triggering sites for the jet formation. It is worth noting that we also outline the existence of a parameter region for the model in which both the radial and the vertical dependence of the backreaction magnetic surfaces oscillate.

Standard model for accretion
The basic idea for accretion onto a compact object has been formulated in Ref. [1] and Ref. [2], where essentially thin disk configurations were considered. Actually, thin accretion disks are commonly present around stellar and black hole systems and they have the advantage to elucidate the accretion features in a one-dimensional model, i.e., essentially the radial dependence on the r coordinate is described, being the vertical one (z) averaged out of the problem and the azimuthal dependence on φ is forbidden by the axial symmetry of the steady disk configuration [7].
For such a steady axially symmetric one dimensional representation of a thin disk configuration, accreting material onto a compact object, having a mass M and an almost vertical magnetic field B = Bê z , is determined via the implementation of the mass conservation equation and of the momentum conservation ones.
If we define the disk accretion rate aṡ where Σ and u r denote the superficial density and the radial infalling velocity of the fluid, respectively, then the mass conservation equation reduces to the following condition dṀ /dr = 0 .
The radial component of the momentum conservation equation splits into two components, one stating that the disk angular velocity ω is essentially equal to the Keplerian one ω K and the advection and the pressure gradient balance each other separately, i.e., ρ 0 (r) and p 0 (r) denoting the values of the mass density and of the pressure, taken on the equatorial plane of the disk. The separation of the gravo-static equilibrium into Eqs. (3) and (4) is a natural assumption when the thinness of the disk is taken into account, see Ref. [7], and it relies on the idea that the advection and pressure gradient terms balance each other, on a different scale with respect to the Keplerian motion of the disk. In Ref. [19], it has been systematically investigated the equilibrium configuration of a disk, demonstrating that, only when the disk is thick enough, the differential angular velocity significantly deviates from the Keplerian value (a configuration very far from the applicability presented in Eqs. (3) and (4) is recovered in the so-called ADAF [20] and see also Refs. [21,22] for the hydrostatic equilibrium in strong gravitational field).
When, in the next Section, we will analyze the magnetic micro-structures generated at small spatial scales, we will see that a deviation from the Keplerian rotation of the disk can be also induced by the plasma backreaction and in that case, it is compensated by the Lorentz force. However, there, the advection term is not present, but, as it has been discussed in Ref. [8], such contributions, at high β values of the plasma, live on a macroscopical spatial scale, decoupling from the microscale balance of the Lorentz force with the correction to the centripetal Keplerian force (see also the discussion in Sect. 3.1).
The vertical component of the momentum equation determines the decay law of the mass density far away the equatorial plane. The details of such decay depend on the particular equation of state we adopt for the fluid, but a reliable expression, valid both for the general polytropic and isothermal case, reads as where H(r) denotes the half-width of the disk and if it is thin we must have H/r 1. The spatial scale of the disk vertical extension is fixed in the vertical gravostatic equilibrium via the ratio between the sound velocity to the Keplerian angular velocity, i.e., H ∼ v s /ω K . The sound velocity can be reliably estimated via the relation v s ∼ K B T /m, being T the disk temperature, K B the Boltzmann constant and m the proton mass.
Finally, the azimuthal component of the momentum conservation equation, which regulates the angular momentum transport across the disk, provides the following expression of the accretion rate: here η v denotes the shear viscosity coefficient, associated to the friction of the disk layers in differential rotation with ω K . This viscosity effect can not be due to the kinetic conditions of the plasma, which is actually quasiideal (see Ref. [23]), and in Ref. [1] η v has been justified via a turbulent behavior emerging in the disk spatial microscales. Such a coefficient can be estimated by comparing the exact expression of the (r, φ) component of the viscous stress tensor τ ij with its interpretation in terms of the correlation function of the turbulent radial and azimuthal velocity component v r and v φ , respectively, i.e., we have from which we get the basic Shakura expression The parameter α ≤ 1 has been introduced to phenomenologically account that all the supersonic velocity fluctuations are unavoidably damped in the turbulent regime. By combining this expression for η v with the relation (6), we arrive to the final expression for the accretion rate of the disk in terms of basic quantities (particle density, temperature, central body mass), namelẏ The possibility to deal with a turbulent behavior of the plasma requires that a sufficiently significant instability can affect the equilibrium configuration. Since in the case of a Keplerian disk the steady profile turns out to be stable [3], after some years of investigation, the source of turbulence has been identified in the MRI [24,25]. Firstly derived in Ref. [26] and Ref. [27], MRI is an Alfvénic instability taking place when the magnetic field intensity within the plasma is below a given threshold. In particular, for a Keplerian disk, the simplest condition to deal with MRI reads as follows v A being the Alfvén speed and k the perturbation wavenumber. In correspondence to an assigned value of the magnetic field, a minimal spatial scale exists λ = 2π/k for which MRI holds. Such a scale, say λ min is easily identified as β denoting the standard plasma parameter, which in astrophysics can take also rather large values. In Ref. [10], it has been shown that just on the small scale λ min , there stated as ∼ v A /ω K , the plasma backreaction to the central body can also take a very different nature with respect turbulence. In fact, a steady perturbation to the equilibrium is allowed for which the correction of the centripetal force associated to the differential rotation, is directly linked to the perturbed Lorentz force raised in the plasma by the emergence of current densities on very small scales. We observe that the background magnetic field is current-free since it is due to the central body magnetosphere and, within the disk, it is a vacuum field. In Ref. [11], this idea of a microscopic crystalline morphology of the perturbed magnetic field, i.e., its radial oscillating behavior on the scale λ min , is implemented also in the limit of a nonlinear beackreaction. As a result, it is possible to show that the disk morphology is decomposed in a series of microscopic ring-like structures, i.e., the mass density of the disk acquires periodic nodes, see also the global model developed in Ref. [12].
In the simplest case of a linear small scale backreaction within the disk, the magnetic flux function is currugated and this perturbed profile is associated to a linear two-dimensional partial differential equation, dubbed in what follows as Master Equation for the disk crystalline structure. We will discuss in detail the morphology of the solutions associated to this equation in order to extract information about the small scale properties of the disk steady profile on which MRI can develop. We conclude by observing that in Ref. [8], it was argued how the presence of a small scale crystalline disk backreaction and the associated relevant micro-current density can have significant implications on the problem of the "anomalous" disk resistivity, required to account for the observed accretion rates in systems like X-ray binary stars.

Physical model
We consider a steady and axialsymmetric thin disk configuration, characterized by a magnetic field B having the following poloidal form where we use the aforementioned standard cylindrical coordinates (r , φ , z) (ê r ,ê φ andê z being their versors) and ψ(r, z) is the magnetic flux function. The disk also possesses a purely azimuthal velocity field v, i.e., where the angular velocity ω is a function of ψ because we are in the range of validity of the so-called co-rotation theorem [28], i.e., Eqs. (12) and (13) holds together with the stationary and axial symmetry hypotheses. Since, we are assuming the plasma disk is quasi-ideal (actually it is true in many range of observed mass density and temperature), we neglect the poloidal velocities, especially the radial component, which are due to effective dissipation, according to the Shakura idea of accretion discussed in the previous Section. Furthermore, we observe that the co-rotation condition (13) prevents the emergence of an azimuthal magnetic field component via the dynamo effect.
We now split the magnetic flux function around a fiduciary radius r = r 0 , as follows where ψ 0 is the vacuum contribution of the central object around which the disk develops (essentially a vertical magnetic field comes out from the dipole-like nature of the field and from the thinness of the disk), while ψ 1 is a small (still steady) correction, here considered of very small scale with respect to the background quantities. By other words, we are studying a small backreaction of the plasma which is embedded in the central object magnetic field, whose spatial (radial and vertical) scales are sufficiently small to produce non-negligible local currents in the disk. According to Eqs. (14) and (12), also the magnetic field is expressed as B = B 0 + B 1 . According to the validity of the co-rotation theorem, at any order of perturbation of the steady configuration, we expand the angular velocity as follows In Ref. [10] (see also Refs. [11,12]), it was shown that, in the linear regime, i.e., |B 1 | |B 0 |, the equilibrium configuration near r 0 reduces to the radial equilibrium only, which, at the zeroth and first order in ψ, gives the following two equations respectively. We recall that H denotes the half-depth of the disk and k 0 ∼ 1/λ min is the typical wave-number of the small scale backreaction, taking the explicit form and, for a thin isothermal disk, it results k 0 H = √ 3β ≡ 1/ . As postulated above, in order to deal with small scale perturbations, we have to require that the value of β is sufficiently large.
In order to study the solutions of Eq. (17), it is convenient to introduce dimensionless quantities, as follows where δ 2 = H/k 0 . Hence, we get which, in what follows we dub the Master Equation for the crystalline structure of the plasma disk. We recall that, according to Eq.(5), the u variable takes a finite range across the thin disk configuration. Moreover, it is also easy to check the validity of the relations where the existence of linear perturbation regime requires |Y | 1. It is relevant to investigate the solutions of Eq. (20) in view of determining the physical conditions under which the crystalline structure, discussed in Refs. [10,29], can actually take place.

Reconciling accretion with magnetic micro-structures
Let us now clarify how the magnetic micro-structures can be relevant for the accretion picture of plasma disk. Actually, in Refs. [10,11,12] the profile of the backreaction on small scale has been analyzed in the presence of differential angular velocity only, i.e. without any poloidal matter flux and therefore accretion features are not addressed.
In the original idea, proposed in Refs. [30,31], the formation of magnetic micro-structures was intended to be an intriguing paradigm in which accretion could be sustained also in absence of dissipative effects, only on the basis of the ideal plasma morphology. The crucial point of this reformulation consists of the possibility that, in the case of nonlinear plasma backreaction, a large number of X-points can form. In fact, when the vertical magnetic field is dominated by the backreaction, it structure is intrinsically characterized by a radial oscillation and therefore X-points appear with the same periodicity of the magnetic surface oscillation. Clearly, near such points, the plasma disk manifests its porosity, simply because the prescription of the ideal electron force balance v r B z = 0 allow now for non-zero radial velocity also in the absence of turbulent viscosity. Of course, in order such a plasma porosity becomes an efficient tool for the disk accretion, a mechanism able to pump the plasma into the X-point is required. In Ref. [30], it was postulated that the pumping of plasma could be supported by intermittent ballooning instabilities, i.e. modulating a naturally mechanism observed in laboratory plasma physics to the astrophysical context scenario.
Furthermore, in Ref. [8], it was emphasized how, limiting attention to a one-dimensional model, a natural inconsistency appears between the Shakura idea of accretion and the emergence of magnetic micro-structures. The nature of this incoherent formulation of the small scale plasma backreaction into a visco-resistive accretion scenario can be easily recognized by observing that, in the Ohm law, the largest contribution to the radial infalling velocity comes from the induced oscillating current densities and therefore, it has, in turn, an oscillating character too (which contradicts the smooth Shakura radial infalling). In a more recent two-dimensional formulation [32] (see also Ref. [33]), it has been reconciled the disk magnetic micro-structure with the Shakura scenario of accretion, by making use of the smallness of the poloidal velocity field with respect to the toroidal differential rotation of the disk. Addressing the generalized Grad-Shafranov equation for an accretion disk (see Ref. [19]), we have overcome the inconsistency observing that the contribution of the microscale phenomenon must be averaged out on the disk macroscales. As a result, the effect of the current induced by the backreaction allows a proper balance of the Ohm law, even in the absence of resistivity high values, but the induced poloidal velocity is averaged to zero on a macroscopical scale.
This scenario suggests that the emergence of magnetic micro-structures across the accreting disk could provide a viable alternative paradigm to the necessity of the socalled disk anomalous resistivity, without rejecting the solid scenario of angular momentum transport as driven by the effective viscosity associated to the turbulence that MRI is able to trigger. However, the presence of the small scale backreaction requires very high values of the plasma β parameter and therefore there is a disk temperature interval in which MRI can still survive, but the disk is too cold to manifest magnetic micro-structures. In such a region of the disk parameters, the necessity of an effective large resistivity to account for accretion onto compact objects, like X-rays binary stars, can probably no way avoided.

A viable solution to the Master Equation
We investigate the solution of Eq. (20) for the crystalline morphology where we consider the regime 1, corresponding to large values of the β plasma parameter. Since the crystalline structure is associated to a periodicity in the x-dependence, we naturally search a solution for the dimensionless first-order magnetic flux function Y in the form where F n and ϕ n are amplitude and phase, respectively. Above, we included a u-dependent phase φ n (u) in order to properly weight both sin and cos functions in the Fourier expansion. Once set to zero the coefficients of the different trigonometric functions, this expression provides the two following equations Eq.(24) admits the solution with Θ = const., which reduces Eq.(25) to the following closed form in F n We observe that, if the periodicity is not associated to the fundamental wave-number k 0 but to a close one k = χk 0 (where χ 1), the equation above is simply mapped by replacing the integer n by χn. For 1, it is clear that we must have n = 1 and then we assume 1−χ 2 n 2 = . By this choice, Eq.(27) becomes independent of , reading as (noting F ≡ F 1 ) We underline how we set 1 − χ 2 n 2 = , thus removing a free parameter from the model, in order to properly recover, for Θ = 0, the same behavior discussed in Ref. [10], and it corresponds, from a physical point of view, to deal with a radial wave-number near k 0 . In the next Section, we will study the Master Equation by using the separation variable method and limiting therefore attention to linear behaviors only. The present case is recovered when the separation constant γ obeys the relation γ = 1 − , again ensuring that the radial wave-number, fixed by the γ value, remains close to k 0 . There we will provide a systematic analysis of the mathematical problem for different positive values of γ in order to span all the parameter space, without focusing attention to physical constraints. However, it remains clear that, for sufficiently low values of γ, the short wave-length approximation could break down and, on the contrary, for too large γ values, the backreaction scale would become smaller than the Debye length, invalidating the MHD scheme. Such peculiar situations take place in rather extreme conditions and they are fixed by the details of the disk configuration, so that we will not account for them in the next Section. Let us consider the problem of solving Eq.(28) with boundary condition F (0) =F , beingF a generic constant (we consider 0 <F 1 since such an equation is invariant for F → −F and because the linear regime must be preserved) and F (0) = 0 (here and in the following the prime denotes u derivatives). This last condition has been selected in order to guarantee the symmetry of the solution with respect to the equatorial plane. The plots of the numerical solution, in correspondence of different values of Θ at fixedF , are in Fig. 1. Such a solution is a Gaussian for Θ = 0 as derived in Ref. [10], while for Θ = 0 it starts to increase for increasing |u| values. It is worth recalling that, when F is greater than unity, the linear approximation fails and the Master Equation is no longer predictive.
We note that, only for Θ = 0, an exponentially decaying behavior of the linear perturbation above the equatorial plane is obtained, just like the background profile [7]. Such a natural similarity in the vertical behavior of the background and of the perturbations was also at the ground of the analyses in Refs. [10,11,29], but the solution derived above demonstrates how this feature is a peculiar property of the simplest case only, i.e., when the phases ϕ n are constant quantities.
We observe that in Eq. (26), in the linear regime for F n 1, we must require that the constant Θ be correspondingly small in order to avoid an exceedingly large value of the u-derivative of the phase φ n . Despite the ef-fective k z value is rather large in the considered scenario (since it is estimated as k z ∼ δ −1 ∼ 1/H √ ), however, we can explore only region of u-values u ∼ k z limited by the request that the linearity of the Master Equation is not broken. In this respect, the increasing behavior of F with large absolute values of u is physically acceptable. Such a behavior can be considered as predictive in the present configurational scheme, only up to the validity of the linearity Y 1 ⇒ F 1, so that such an increasing behavior is never associated to a diverging non-physical feature: our model is a local one, we are exploring a small region around a fiducial z 0 value and only in that neighborhood, our linear equilibrium makes sense. Similar considerations hold also for the radial dependence, when the oscillating behavior were mapped into an increasing one with the radial coordinate x, namely the trigonometric function could become hyperbolic ones (see the next Section).
We also that, from a physical point of view, the increasing behavior is not surprising. In fact, it can be explained because, while the background profile is fixed by a gravostatic equilibrium, the perturbation morphology has nothing to do with gravity and it is mainly produced by the balance of the centripetal (say centrifugal) and the Lorentz forces.
However, it is important to stress that for sufficiently small values of Θ a peak of the function F around the equatorial plane is still present. This feature takes place because for small values of Θ the effect of the nonlinear term Θ/F 3 is minimized.
From a physical point of view, a process able to induce a crystalline structure in the plasma disk (for instance a sound or a gravitational wave), i.e., a boundary condition able to select the solution associated to a radial corrugation of the magnetic flux function, should be vertically coherent in phase to produce the same vertical gradient of the background profile in the perturbation too. On the contrary, a vertical shear also in the phase of the underlying process is responsible for an inversion of the vertical gradient profile between the background and the perturbation.
The most important physical conclusion of the present analysis is to demonstrate that the crystalline structure is a very general feature of the linearized perturbed equilibrium, but, being associated to small values of (i.e., large values of the plasma β parameter), the radial oscillation can be associated to a basic wave-number only, very close to the natural one k 0 .

Separable solution of the Master Equation
Since the approach in the previous section gave us the solution Y = F (u) sin(x) as the only viable case for 1, it is natural to investigate solutions relying on the separation of variable method. In particular, we will be able to explore the situation in which the wave-number of the radial oscillation (x-dependence) is a generic constant times the fundamental wave-number k 0 .
In order to investigate solutions relying on the separation of variable method, the Master Equation (Eq.(20))   can be rewritten using Y (x, u) = Y 1 (x)Y 2 (u) in the following decoupled form: where we have introduced the arbitrary constant γ. The first equation has two possible kind of solutions, namely Eq. (30) can be rewritten in the form of the the Kummer equation of parabolic cylinder functions using the scaling u → u/ √ 2. In particular, we get where A = (γ − 1)/(2 ).
Let us now analyze the properties of the solution by numerically integrating Eq.(33) using initial conditions Y 2 (0) = 0.01 and Y 2 (0) = 0. We also fix the parameter = 0.01 and run the γ parameter in order investigate the morphology of the solution for different values of the constant A. From the numerical analysis, a range of the parameter A emerges (dependent on the choice of ) for which the behavior of the function Y 2 becomes to oscillate, as well sketched in Fig. 2. Since these values of A still corresponds to positive values of the constant γ (for our setup, γ 0.9), we obtain a range of the solution of the Master Equation where the radial and vertical dependence of the perturbed magnetic flux function both have an oscillating profile as represented in Fig. 3 and in the detailed cross sections in Figs. 4 and 5.
Clearly, when γ takes negative values, the x dependence of the function Y is no longer oscillating (while the u-dependence still remains oscillating) and the crystalline morphology of the disk is suppressed. As we shall discuss below, the oscillating nature of the vertical dependence of the perturbed magnetic flux function has the relevant physical implication that, according to the expression for the magnetic field in Eqs. (21) and (22), a series of O-points appears in the disk configuration, in which the radial component B r vanishes.

Whittaker equation and O-points
The many different behaviors of the solution of Eq. (33) can be better understood by investigating peculiar properties of the corresponding analytic solutions. In particular, using the following change of variables: y = u 2 /2 and Y 2 = (2y) −1/4 W (y), Eq.(33) rewrites now as which is a particular case of Whittaker's equation, i.e., with K = −A/2 and µ = ±1/4. The linearly independent solutions of this equation are given by K,µ = e −y/2 y µ+1/2 U (1/2 + µ − K, 1 + 2µ, y) , (37) where M and U are the Kummer and Tricomi functions [34], also known as confluent hypergeometric functions of the the first and second kind, respectively. We will use only the function M because the function U is multivalued, and we study the solutions only for finite y since the model holds for a limited range of u values. Moreover, for γ ≤ 1, the function W K,µ is bounded since M (1/2 + µ − K, 1 + 2µ, y) behaves like where Γ denotes the standard gamma function. Specifically, the Kummer function M (a, b, y) is formally defined as where the rising factorial denoted with the index (n) is defined as The sum is converging everywhere, while the derivative of M with respect to y is given by the following recurrent equation: dM (a, b, y)/dy = (a/b)M (a + 1, b + 1).
Using this formalism, we get the following solutions of Eq.(33): where c ± 1 are integration constants and we denoted with + /− the choice µ = ±1/4.
The zeros of the derivative of the solutions above correspond to that one of the radial component of the magnetic field from Eq. (22). Using recursive expressions, we get In this respect, we note that algorithms are present for systematically finding the roots of Kummer functions [34].

Asymptotic solutions
Taking into account the solutions for Y 2 (u) in Eqs.(38) and (39), let us now define the following functions where In the asymptotic limit −A u 2 (for A < 0) and finite u (we recall that the thin-disk hypothesis limit the range of u values), we obtain where p = √ −A, while θ r (u) and θ i (u) are assigned polynomials. Thus, the function Z(A, u), the solution Y 2 (u) and its derivative with respect to u oscillate giving rise to a sequence of O-points across the vertical axis.
For positive values of A, we get instead the following asymptotic behavior for A u 2 : where θ 1 is again an assigned polynomial. This case for A > 0 is not relevant for the present physical discussion.

Laguerre polynomials
We have shown how, in the case −A u 2 , there is a countable set of O-points, but it is also possible to find other zeros series in the correspondence of smaller |A|. It is easy to verify that the following relation takes place between the Kummer function M (a, b, y) and the Laguerre polynomials L  [1,10]. Standard color scheme from red (n = 1) to blue (n = 10). The number of oscillations increases with increasing n.
The equation above is a particular form of Eq.(41) and, since it is well-known that for each n there is a finite number of zeros of the Laguerre polynomials which increases with n, we have shown how a countable set of zeros on the z axis emerges. In particular, the zeros of Laguerre polynomial L (α) n belong to the interval (0, n+α+(n−1) √ n + α), and in Fig. 9 we plot, as an example, dY − 2 /du from Eq.(49) for the first 10 integer n values.

Discussion
We have shown how the z-dependence of the magnetic field is associated to O-points of the corresponding configuration. In such points, the radial component of the magnetic field, due to the plasma backreaction only, vanishes. Here, we outline a new feature induced by the crystalline morphology of the plasma response to the dipole field of the central object.
Let us now discuss the relevance of the O-points derived above, in view the physical properties of the plasma disk thought as an accreting structure. Our analysis has been performed for an ideal plasma disk in the absence of poloidal components of the velocity field, i.e., v r = v z ≡ 0. However, in the Standard Model of accretion [1], the plasma must have a finite electric conductivity and the azimuthal component of the generalized Ohm law reads where J φ is the corresponding azimuthal component of the current density and σ is the constant coefficient of electric conductivity. In the standard mechanism for accretion onto a compact object, both the quantities B r and v z are essentially negligible and, where σ takes very large values, the accretion is suppressed, i.e., since B z = 0, we must get v r 0.
In the spirit of the present analysis, in those regions where the plasma is quasi-ideal, we find the relation Clearly, nearby the numerous O-points of the magnetic profile, where B r is nearly vanishing, we get very high values of the vertical velocity, which suggests the existence of privileged plasma sites for the jet formation. This perspective has been investigated in some detail in Refs. [13,14]. It is worth observing that, when an anomalous resistivity is considered, the term v z B r can be actually safely neglected since the balance of the Ohm law is guaranteed by the relation v r B z = −J φ /σ. Of course, if we can state that v z B r ∼ v r B z , we naturally have large values of v z (due to the smallness of B r in a thin disk). However, this situation can not be regarded as the most general one, unless some other information on the velocity field is assigned by extra physics in the disk. On the contrary, when the values of sigma are very large, the balance between the two terms is mandatory and v z is very large near O-points. Clearly, also the solution v r ∼ v z = 0 is viable, but this situation is just that one discussed in Sect. 3, where differential disk angular velocity is included only. However, in order to get the equality v r B z = v z B r and the O-points, simultaneously, J φ /σ must be a negligible contribution, i.e., k 0 can not exceed a critical magnitude fixed by the value of the electric conductivity σ. If not, we can still have a large value of v z , but again in specific cases only.
Actually, the study of the crystalline micro-structures in the presence of poloidal velocity has not yet been fully investigated, see Ref. [32] for recent developments (the most important difficulty relies on the nonlinear advection terms in the poloidal velocity components). However, the picture associated to the Master Equation is reliable in view of the smallness characterizing the poloidal velocity components with respect to the toroidal rotation velocity, as expected in a thin accretion disk.
The main merit of the present analysis consists of having outlined the general character that the O-points takes in the linear profile of the perturbed magnetic field, as soon as the short wavelength backreaction of the plasma is excited in disk with high β values.

Conclusions
We analyzed general features of the small scale morphology of the backreaction that is generated in a thin plasma, embedded in the magnetic field of the gravitationally confining central object. In particular, we studied the Master Equation of the so-called crystalline structure [10], searching for a satisfactory characterization of the admissible profile in the linear backreaction limit.
We have followed two different, but complementary, approaches. The first one is based on a general Fourier expansion, while the second one relies on a separated variable procedure. This analysis offered a complete spectrum of the available solutions, outlining new behaviors in the plasma equilibrium, absent in the basic treatment in Refs. [10,11].
The main merit of this systematic study of the linear backreaction of the plasma disk consists of the relevant physical implications of the obtained solutions, as discussed above, i.e., the inverse behavior of the vertical gradient with respect to the background one (as the highest harmonics are considered) and the possible emergence of a consistent sequence of O-points, where the radial magnetic field component vanishes.
The first of these issues is of interest in view of the nonlinear scenario of the plasma backreaction discussed in Refs. [11,12], where, see also Ref. [29], the possibility for the radial fragmentation of the disk into a microscopic array of rings is demonstrated. In this respect, the linear inversion of the vertical gradients suggests that, in the nonlinear regime, a separation of the disk into two symmetric components, up and down the equatorial plane, could take place, with significant implications on the transport features across the disk.
The second result, concerning the appearance of Opoint series of the magnetic configuration, is relevant in view of the realization of conditions for jet emission. In fact, nearby the O-points, when the plasma resistivity is suppressed, the vertical component of the plasma velocity can take very large values, see the discussion in Refs. [13,14] about the seeds of a vertical matter flux from the disk. Using the solution given in terms of hypergeometric confluent functions, we have also analytically shown the existence of the O-point series and such a series critically depends on the parameter of the separation of variable technique.
The systematic analysis here pursued of the Master Equation for the crystalline structure of a stellar accretion disk, constitutes a well-grounded starting point for the investigation of nonlinear features of the magnetic confinement of the plasma disk. In fact, as discussed in Ref. [11], the extreme nonlinear regime is characterized by a dominant influence of the MHD-force in confining the plasma with respect to the radial gravitational field. Furthermore, in Ref. [35] it is outlined how the request that the Master Equation holds in the nonlinear case too turns out to be a proper choice to close the equilibrium system, preserving the crystalline morphology.

Authors contributions
All the authors were involved in the preparation of the manuscript. All the authors have read and approved the final manuscript.