Linear Two-Dimensional MHD of Accretion Disks: Crystalline structure and Nernst coefficient

We analyse the two-dimensional MHD configurations characterising the steady state of the accretion disk on a highly magnetised neutron star. The model we describe has a local character and represents the extension of the crystalline structure outlined in Coppi (2005), dealing with a local model too, when a specific accretion rate is taken into account. We limit our attention to the linearised MHD formulation of the electromagnetic back-reaction characterising the equilibrium, by fixing the structure of the radial, vertical and azimuthal profiles. Since we deal with toroidal currents only, the consistency of the model is ensured by the presence of a small collisional effect, phenomenologically described by a non-zero constant Nernst coefficient (thermal power of the plasma). Such an effect provides a proper balance of the electron force equation via non zero temperature gradients, related directly to the radial and vertical velocity components. We show that the obtained profile has the typical oscillating feature of the crystalline structure, reconciled with the presence of viscosity, associated to the differential rotation of the disk, and with a net accretion rate. In fact, we provide a direct relation between the electromagnetic reaction of the disk and the (no longer zero) increasing of its mass per unit time. The radial accretion component of the velocity results to be few orders of magnitude below the equatorial sound velocity. Its oscillating-like character does not allow a real matter in-fall to the central object (an effect to be searched into non-linear MHD corrections), but it accounts for the out-coming of steady fluxes, favourable to the ring-like morphology of the disk.

💡 Research Summary

The paper presents a linear, two‑dimensional magnetohydrodynamic (MHD) analysis of a thin accretion disk surrounding a highly magnetised neutron star. Building on the “crystalline structure” concept introduced by Coppi (2005), the authors incorporate a finite accretion rate and examine how toroidal electric currents interact with the disk’s magnetic field. Because only toroidal currents are considered, the model would be inconsistent without a small collisional effect; this is introduced phenomenologically as a constant Nernst coefficient, which represents the thermal power of the plasma. The Nernst term generates non‑zero temperature gradients that balance the electron force equation, linking the radial (v_r) and vertical (v_z) velocity components to the thermal structure of the disk.

Solving the linearised MHD equations under these assumptions yields an oscillatory solution for the magnetic perturbation, density, and pressure—exactly the “crystalline” pattern previously described. The oscillations are superimposed on a viscous background associated with the differential rotation (Ω ∝ r^‑3/2) of the disk. Viscosity, modeled in the spirit of an α‑disk, transports angular momentum outward and produces a net inward mass flux. However, the resulting radial inflow speed is only a few thousandths of the local sound speed, i.e., v_r ≈ 10⁻³–10⁻⁴ c_s. This very slow, quasi‑periodic radial motion does not correspond to a steady, monotonic infall of matter; instead, the flow oscillates back and forth, a behaviour that would be smoothed out only by non‑linear MHD effects not captured in the present linear treatment.

The temperature gradients required by the Nernst term are directly proportional to v_r and v_z, confirming that the collisional term is essential for maintaining equilibrium. Without it, the electron pressure gradient could not balance the Lorentz force, and the linear solution would break down. The authors therefore argue that the Nernst coefficient, though small, is a crucial ingredient for a self‑consistent description of the disk’s electromagnetic back‑reaction.

The paper also discusses the physical implications of the crystalline structure. The periodic density enhancements naturally lead to a ring‑like morphology, a feature that has been observed in several high‑resolution images of accretion disks. The oscillatory character of the solution suggests that, at the linear level, the disk can sustain steady fluxes of mass and angular momentum without a large-scale, continuous inflow. The authors acknowledge that to obtain a genuine accretion flow capable of delivering matter to the neutron star, one must go beyond the linear approximation and include non‑linear instabilities such as the magnetorotational instability (MRI) or shear‑driven turbulence.

In summary, the study provides (1) a clear demonstration that toroidal currents combined with a constant Nernst coefficient can produce a self‑consistent, oscillatory MHD equilibrium; (2) quantitative estimates showing that the induced radial inflow is orders of magnitude below the sound speed, consistent with a quasi‑steady, ring‑forming disk; and (3) a roadmap for future work, emphasizing the need for non‑linear simulations to capture realistic accretion onto highly magnetised neutron stars. The results enrich our theoretical understanding of how magnetic fields, thermal effects, and viscosity cooperate to shape the structure and evolution of compact‑object accretion disks.