Dynamically screened strongly quantized electron transport in binary neutron-star merger

We examine electron-transport coefficients in magnetized hot and dense electron-ion plasma relevant in binary neutron star merger simulation. We calculate electrical and thermal conductivities in low density, high temperature, highly magnetized plasma of binary neutron star mergers where quantum oscillatory behavior of electrons emerge. For pronounced thermodynamic effects, we consider zeroth Landau level population of electrons for the calculation of conductivity. We solve Boltzmann equation in presence of magnetic field to obtain the dissipative components of electrical and thermal conductivities. The dissipative coefficients are formulated considering frequency dependent dynamical screening in the quantized electron-ion scattering rate. Numerical estimations show that the effect of dynamical screening of photon propagator on electrical and thermal conductivities is pronounced. We observe that dynamical screening reduces the maxima of both the electrical and thermal conductivities by factors of thirty one and twenty respectively leading to a reduction in the corresponding time scales of these coefficients. The common scaling factor between electrical and thermal conductivity is also observed to follow cubic relationship with temperature violating Wiedemann–Franz law.


Introduction
Binary neutron star mergers and collider experiments are the sources of most extreme states of matter in a sreemoyee.sarkar@nmims.edub souvik.adhya@ifj.edu.pl the universe with densities several times nuclear saturation density and temperatures upto hundreds of MeV [1].The recent detection of gravitational wave signal GW170817 originating from binary neutron star (BNS) merger by the LIGO and Virgo detectors have opened up a new era in multi-messenger astronomy [2][3][4].Additionally, short gamma-ray bursts (SGRBs) were also detected by the Fermi satellite GRB170817A indicating the presence of huge magnetic field in the merging event [5][6][7].These mergers are unique astrophysical objects of significant sources of gravitational radiation, electromagnetic as well as neutrino emission [8].They offer a novel avenue to study highly non-linear gravitational effects blended with complex micro-physical processes; serving as Einstein's richest natural laboratory [9].
In the event post merging, a remnant neutron star is created and if the remnant possess a mass beyond Tolman-Oppenheimer-Volkoff (TOV) limiting mass, the merged object survives for 10-100 milliseconds and collapses thereafter.The description of neutron star mergers requires the knowledge of General Relativistic Magneto Hydro-dynamics (GRMHD) [10][11][12][13][14][15][16][17][18].Most of these general-relativistic simulations account for ideal Magneto Hydro-dynamics (MHD) limit.In a recent work Ref. [19], the authors have pointed out that Hall effect plays an important role in magnetic field decay of the merged object, hence, the Hall effect needs to be incorporated in the ideal MHD merger simulation.In this Ref.[19] the authors have considered the electrical conductivity (σ) of warm neutron star crust in the nonquantizing scenario calculated in Ref. [20].Apart from these studies, in Ref. [21] authors have analysed the relevance of thermal conductivity (κ), viscous coefficients in BNS mergers in the non-magnetic scenario.In view of these recent studies [19][20][21], we analyse the importance of dissipative mechanism in the merger simulation by evaluating quantized electrical and thermal conductivity coefficients with magnetically modified many-body effects in the present paper.We calculate quantized dissipative coefficients with frequency dependent screening in both hot and dense plasma relevant in binary neutron star merger.The results can be implemented in analysing the magnetothermal evolution [22] of the merged compact object as well.
We consider fully ionized plasma of electrons and ions.Heat and charge in this medium are transported by electrons.The dominant electron transport mechanism is scattering on ions in the liquid phase.In presence of extreme magnetic field (B) (B∼ 10 16 G) and density (ρ) (ρ ∼ 10 12 g cm −3 ), the classical description of electrons breaks down.Therefore, one should incorporate Landau quantization of energy levels in the formalism.This quantization occurs for a particular set of temperature, density and magnetic field in case of neutron star.Thus, the inclusion of Landau quantization eventually modifies the non-magnetic electrical conductivity to great extent [23][24][25][26].In the present paper, we focus mainly on the strongly quantizing case, since in this domain, the transport coefficients receive major modification due to the magnetic field.
The calculation of electrical and thermal conductivities by solving Boltzmannn equation in ultra-compressed plasma have been studied by several authors over the last few decades [27][28][29][30], see for a review [31].This requires the information of scattering rate of plasma constituents.The calculations of scattering rates considering screened Coulomb potential have already been observed in different Refs.[23][24][25][26]32].In all these calculations it has been assumed that ions are static scatterers.This formulation can not be easily transported to the relativistic domain where dynamical effects are important for reliable description of transport coefficients.Medium modified Hard-Thermal-Loop (HTL) propagators for hot and Hard-Dense-Loop (HDL) for dense plasma include dynamical effects of relativistic medium through frequency dependent screening [33][34][35][36][37][38].While Debye screening in plasma is related to the static/longitudinal photon exchange, the exchange of magnetic/ transverse photons contribute to dynamical frequency dependent screening of the plasma particles.It is observed in different studies [36,37,[39][40][41][42][43] that for ultra-degenerate case, both in Quantum Chromodynamics (QCD) and Quantum Electrodynamics (QED) plasmas, the transverse interactions not only become important but they dominate over their longitudinal interaction.In a recent calculation [20], the authors have included many-body effects through the HTL modified propagator in the calculation of non-quantized elec-trical conductivity in magnetized, warm neutron star crust.Motivated by all these calculations of dynamical screening in different coefficients, we include HDL modified photon propagator in quantized σ and κ in the context of BNS merger in the present paper.We incorporate plasma screening through magnetic Debye mass.The current calculation is important in two ways.First, in this paper we consider the quantized transport coefficients in estimating dissipation coefficients in BNS merger.Second, this formulation includes Landau damping in quantized electron-ion interaction rate in both hot and dense relativistic, magnetized plasma.Here, we perform the calculations of σ and κ in an extreme scenario of temperature ∼ 12 MeV, density ∼ 10 12 g cm −3 and magnetic field B ∼ 10 16 G.Finally, with the strongly quantized σ and κ, we estimate the relevant dissipative time scales and compare it with the survival time period of the post-merger object.Validity of Wiedemann-Franz law has also been discussed.
The paper is organised as follows.In section II, we derive the longitudinal electrical and thermal conductivities in a dynamically screened QED plasma.Next, in section III we present the constraints on temperature, magnetic field and density of hot and dense plasma to become relativistic and strongly quantized in the BNS merger scenario.We present numerical variations for σ and κ with temperature, magnetic field and density for the dense, relativistic plasma along with estimation of corresponding dissipative time scales in section III.Finally, we summarize and discuss the impact of dynamical screening on both the coefficients and decay time scales in section IV.Throughout the manuscript we will use following notation for four vectors p = ( p , p z , p ⊥ ) and k = ( k , k). 1

Electrical conduction from Transport Theory
In the current paper, we consider fully ionized plasma of two components: electrons (e) and positive ions of charge Ze (Z atomic number of the nucleus).In compact objects, the huge magnetic field quantizes the motion of electrons in the QED plasma.In this section we derive the coefficients for electrical conduction, electrical conductivity and coefficients for thermal conduction, thermal conductivity, in magnetized electron-ion plasma from transport theory.In presence of magnetic field both the coefficients are anisotropic and the conductivity tensor is given below, In the above matrix equation, σ /κ and σ ⊥ /κ ⊥ are the parallel and perpendicular components of σ and κ respectively in presence of external magnetic field along z direction.σ H / κ H , are the Hall components of the conductivities.In the current paper we present the calculation of quantized σ /κ in electron ion plasma.For the rest of the paper, we re-define σ , κ as σ and κ respectively.The σ is related to the electric current density (j) and satisfies the constitutive relation j = σE where E is the electric field.For thermal conductivity Q = −κ∇T , where, Q is the amount of charge transported through the plasma under the temperature gradient ∇T .j and Q are related to displacement of the electronic distributions from their equilibrium configuration due to the presence of electric field and temperature gradient in the plasma respectively and can be calculated from kinetic theory, In the above equation e is the charge of an electron, v p is the velocity of the electrons, p is the energy of electrons, µ is the chemical potential of electrons.Φ is the off-equilibrium distribution function which arises due to the presence of electromagnetic field in plasma.Φ is obtained by solving Boltzmann equation in presence of magnetic field.In presence of perturbation along z direction, the distribution function (f n,pz,s ) evolves ac-cording to the magnetically modified Boltzmann equation as given below [24], In the above equation, f n,pz,s describes the population of electrons defined by the quantum state n, s, p z , n is the number of the Landau level, s is the spin and p z is the z component of electron's momentum.v z is the z component of the velocity of the particle.The third term in the LHS of eq.( 4) arises from the Lorentz force term F = e( E+ v p × B).In absence of external magnetic field, the Lorentz force term vanishes.
The RHS of Eq.( 4) contains the information of scattering rate of electrons with the ions present in the medium, where, sum is over final state quantum numbers n , p z , s .
where, f 0 ( p ) is the equilibrium Fermi-distribution functions and Φ n,pz,s is the perturbation due to the electromagnetic field.
We now proceed to calculate the collision integral considering strongly quantizing magnetic field.To calculate the interaction rate, we consider an electron with momentum p = ( p , p z , p ⊥ ) and mass m exchanges a virtual photon of momentum q = (q 0 , q) with an inmedium ion of momentum k = ( k , k) and mass M .The electron emerges with momentum p = ( p , p z , p ⊥ ) and ion with momentum k = ( k , k ) (Fig. (1)).In order to obtain finite interaction rate, we use the HDL re-summed photon propagator with transverse and longitudinal components.
We start by re-visiting the expression of the interaction rate (I f i ) without electromagnetic modification [24], where, M f i is the electron-ion scattering matrix.M f i can be written as [20], where, are the components of currents.e = √ 4πe, v k is the velocity of ion with momentum k.The ∆ T and ∆ L are transverse and longitudinal HDL effective photon propagators respectively.The form of the electronic spinors (u(p), u(p )) are given in the Appendix.A.
To proceed further, here, we describe the screening mechanism of electron-ion plasma.In earlier calculations [27][28][29][30], the authors have implemented longitudinal component of photon propagator to screen the Coulomb potential, In the above equation m 2 D = e 2 dn e /dµ, where, n e is the number density of electrons.Following linear response theory, there exists additional weak dynamical screening, related to the energy transfer to the constituents of plasma known as Landau damping.This arises because of non-zero frequency of the plasma.We implement the effects of non-zero frequency in both the electric and magnetic components of the photon propagator computed within the HDL formalism, where, P µν and Q µν are the transverse and longitudinal projectors, respectively, and the effective longitudinal (∆ L ) and transverse (∆ T ) propagators are given by the following expressions, In the above equations Π T and Π L are the transverse and longitudinal HDL photon self-energies and are given by, In presence of strong magnetic field, electron density present in the plasma changes, leading to a modification in the screening.At low temperature and strong magnetic field, the presence of sharp Fermi surface modifies the nature of screening.In the relativistic domain m µ , the Debye mass is given by [44], One obtains eqn.( 16) using m 2 D = e 2 dn e /dµ, where, n e in the relativistic domain (m µ) is given by eBµ/2π 2 .In the screening mass, we have ignored the finite temperature correction.This assumption is valid as we consider only the degenerate domain (T < µ) of the plasma.In addition, the conductivities are weakly dependent on the screening mass as discussed in [24].Hence, inclusion of finite temperature Debye mass would mark negligible correction in the final results of the conductivities.
In order to proceed further, we compute the phase space factor in the interaction rate given in eq.( 7).We do not consider the change of momentum of ions in the phase space factor.Hence, the linearized phase space factor can be written as, The final expression for the interaction rate thus becomes (the details of the derivation are given in Appendix A), n i is the number density of ions and u, ζ and F are given in the Appendix A. In order to find the transport coefficients, it is useful to define a dimensionless scattering rate a and a dimensionless perturbation to the distribution function (Ψ ) defined as, where, σ 0 = πZ 2 e 4 /ω 2 B and ω B = eB/m.The electric charge e is related to the fine structure constant by α = e 2 /(4π) = 1/137.Using the above two equations, we obtain the dimensionless form of the linearized Boltzmann equation as given below [45], γ = ± denotes the scattering channel for forward (+) and backward reactions (−).In the current paper, we present the results for the strongly quantizing scenario (i.e.zeroth Landau level) which provides the maximum effect with finite magnetic field in contrast to the nonmagnetic scenario.For zeroth Landau level, n = n = 0, spin degeneracy is absent and backward scattering (γ = −1) is the only allowed channel for scattering.Hence, after solving the dimensionless Boltzmann equation the off equilibrium distribution function is obtained as, where, E = p /m, Ψ ≡ Ψ 0,−1 and We obtain the expressions for σ and κ with only zeroth Landau level population by inserting the value of Φ in the equations below, We then insert j and q in the constitutive relations j = σE and Q = −κ∇T to obtain both σ and κ.The final form of electrical conductivity is given below, here, δ 0 is a constant and is given by The thermal conductivity coefficient can be obtained by integrating following expression, From the above two equations it is evident that at low temperature satisfying ( p − µ) ∼ T , ∂f 0 /∂ p = δ( p − µ).After integration thus σ becomes temperature independent and κ varies linearly with temperature T .Hence, κ/σ varies linearly with temperature till ( p − µ) ∼ T is satisfied.This linear relationship of σ/κ with temperature is known as Wiedemann-Franz law.
In the next section, we study the variation of σ and κ with different parameters and the effects of the frequency dependent screening derived in the current section.

Results
In this section we describe the behaviour of σ and κ with density, temperature, magnetic field in the hot and dense QED plasma.First we present the physical conditions of the plasma for the calculation of the transport coefficients.

Physical conditions
Physical properties of the BNS merger, which forms an unstable configuration are different from isolated neutron stars.We consider simplest possible constituents of post-merger object of electron-ion plasma with fully ionized ions and free mobile electrons in the low density (up to 10 12 g cm −3 ), high magnetic field (up to 10 17 G) and high temperature (T ∼ 15 MeV) regime of BNS merger.Electron density n e is related to ion density n i via n e = Zn i .We consider the magnetic field present along the z direction.Scattering of electrons with ions only contribute in electrical conductivity.In the absence of magnetic field the electron density can be written as, where, f 0 ( p ) = 1/(exp( p −µ T ) + 1), T is the temperature.In absence of magnetic field energy of the electrons are given by p = p 2 f + m 2 , where, p f is the Fermi momentum.
The magnetically modified electronic energy states are obtained as solutions of Dirac equations in presence of finite magnetic field [46,47].The positive energy states are denoted by quantum numbers , p z , n, s. p z is the electron momentum along the field which we consider along z direction, s = ±1 is the helicity, and n = 0, 1, 2 enumerates the Landau levels.For nonzero B, µ ≡ µ(B) = (2π 2 n e )/mω B .The energy of the relativistic electrons in presence of magnetic field is p = p 2 z + m 2 + 2nω B m.The ground Landau level is non-degenerate with respect to spin while the higher levels are doubly degenerate.The number density of electrons in presence of magnetic field is written as, where, the sums are over n, s.The magnetic field strongly quantizes the motion of electrons and different transport coefficients receive significant contribution when the electrons are confined to the zeroth Landau level.We do not consider the situation when ions receive quantum modifications due to the magnetic field.
Parameters which determine zeroth level population are as follows [48], In the above equation, B 12 = B/10 12 , ω ce is the cyclotron frequency for electrons.B is strongly quantizing if ρ < ρ B and T T ce .It is convenient to introduce the relativistic parameters x r = p F /m ∼ 1.008( ρ6Z A ) 1/3 (where ρ 6 = ρ/10 6 ), T r ∼ 5.930×10 9 K.The electron-ion plasma is relativistic for x r 1 and T T r .Thus the electrons become relativistic when T > 5 × 10 9 K and density ρ ∼ 10 6 g cm −3 .
The momentum of an electron is related to the energy via the relation From this expression one can obtain the maximum Landau level that the electrons can populate and is given by the integer part of ν = (E 2 − 1)/2b.The energy of the electrons is constrained by the relation (E 2 − 2b) < 1 to meet the condition of lowest Landau level.This is an important condition for the plots of σ and κ as we describe later in this section.The parameters for density, temperature and magnetic field are appropriately chosen for relativistic quantized electrons to simultaneously meet the physical conditions applicable for the merging scenario also.We consider Fe and Mo for the numerical analysis of both the coefficients.The reason behind choosing these elements is given in the Appendix.B.

Variation with density
Fig. 2 shows the variation of σ and κ given in Eqs.(24) and (25) with ρ for different temperatures and different magnetic fields for different elements Fe and Mo.In order to consider electrons to be relativistic, the density and temperatures are chosen as ρ 10 6 g cm −3 and T 5 × 10 9 K respectively.For fixed B, µ increases with ρ and electrons start to populate higher Landau levels.Since, we are interested in population of the zeroth Landau level, the density and temperature should also satisfy ρ < ρ B and T T ce as given in Eq. (28).With these two conditions, both the coefficients have been obtained by numerically integrating the expressions in Eqns.24

Effect of HDL modified propagator
In the Fig. 5 we explicitly estimate the modification induced in σ and κ when electron-ion interact through HDL photon propagator and also compare the case when these plasma constituents interact through static, longitudinal photon propagator.The upper panel is for variation of σ and the lower panel is for κ.The left panel considers magnetic field of B=10 17  ρ [g cm -3 ] ρ [g cm -3 ] ficients when the interaction rate is screened through static, longitudinal photon propagator relevant in nonrelativistic plasma and black curves when the constituents interact through frequency dependent HDL photon propagator required in the relativistic plasma.We find that the inclusion of HDL propagator reduces the values of both the coefficients substantially.For σ the maximum value of the curve decreases by a factor of 31 in comparison to the case when the interaction rate is Debye screened, whereas, for κ the factor of reduction is 20.
It is known that modifications to different equilibrium and non-equilibrium properties of plasma due to inclusion of HDL propagators emerge from frequency dependent photon propagator.The reduction in both the transport coefficients shown above arises due to the increase in the interaction rate caused by non-zero frequency in both the longitudinal and transverse components of the photon propagator.The electromagnetic perturbation to the distribution function (Ψ ) is inversely proportional to I f i and σ, κ are directly proportional to Ψ .Hence, increment in Ψ decreases both the coefficients.The increment in I f i with the inclusion of HDL photon propagator can be seen from Eq.( 22).The second term in Eq.( 22) is dependent on v k < 1. v k reduces the numerator of second term in Eq.( 22).In the denominator of the same equation q 2 ⊥ < 2mω B and m 2 D < 2mω B lead to u < 1, ζ < 1 respectively.Hence, the reduction due to small value of v k gets compensated and interaction rate increases due to frequency dependent dynamical screening.

Estimation of dissipation time scales
In this paper we give quantitative estimations of both the time scales related to electrical and thermal transport coefficients.The magnetic field decay (or diffusion) timescale due to Ohmic dissipation is given by the wellknown expression [19,49], τ σ = 4πσλ 2 B /c 2 .The time scale is dependent on the two factors magnetic field scale height (λ B ) and σ.For a typical value of λ B ∼ 10 −4 c.m ( λ B ≥ λ mf p (electron mean free path) ∼ 10 −4 cm, λ B < λ lim (limiting magnetic field scale height) ∼ 10 −3 cm) we obtain τ σ ∼ 15 ms which is well within the range of survival time period of neutron star merger.Choice of λ B satisfies both the conditions specified by Ref. [19], λ mf p < λ B > λ lim .With the same choice of λ B , τ σ ∼ 600 ms in the non-relativistic scenario.We emphasize on the reduction of time scale due to inclusion of dynamical screening in the calculation in contrast to the scenario when electrons-ions interact through longitudinal plasmon.
Thermal conduction time scale is defined as τ κ = c v ∆z 2 /6κ, where, c v is the specific heat and ∆z is the region which is hotter than surrounding by a temperature difference ∆T .The contribution of degenerate relativistic electrons in specific heat is, c v = 5.4 × T 9 , n 0 is the nuclear saturation density and T 9 = T /10 9 .For n e = 0.57 × 10 35 cm −3 , temperature T = 2×10 10 K and ∆z ∼ 25 c.m, thermal equilibration time scale becomes ∼ 50 ms which is of the order of time scale of the merged object.For ∆z ∼ 1 km, the time scale is much greater than the survival time period of the merged compact star.However, the estimation of time scales presented here demands a more detailed hydrodynamic calculation to assess the inclusion of dissipation and resistivity in magneto-hydrodynamic simulation of BNS mergers.

Validity of Wiedemann-Franz law
In a typical degenerate plasma, electrons' conduction dominates system's dissipation coefficients through electrical and thermal coefficients.This suggests linear increment in σ/κ with T known as Wiedemann-Franz law.From the plot of σ in fig.

Summary and conclusions
In this work, we have calculated the quantized longitudinal electronic transport coefficients in dynamically screened hot and dense magnetized QED plasma involved in binary neutron star merger simulations.The calculation considers scattering of electrons with ions through screened electromagnetic force in electron-ion plasma.We have presented the plots for the variation of transport coefficients with density ( ∼ 10 12 gm cm −3 ), temperature (∼ 10 10 K) and magnetic field (∼ 10 17 G) for two elements Mo and Fe.The scales for the generation of the plots are chosen so that they obey the conditions of relativistic electrons at density ρ > 10 6 g cm −3 , temperature T > 5.93 × 10 9 K and zeroth Landau level population by obeying √ E 2 − 1/2b 1.These two constraints indicate, high magnetic field and low density regime of BNS mergers as the valid physical domain of our calculation.
For the calculation of both electrical and thermal conductivities, we have assumed particles are slightly out of equilibrium which allows us to solve the Boltzmann equation numerically.We calculate electron-ion scattering amplitude with screened electromagnetic interaction of magnetically modified spinors.The off-equilibrium distribution function has been obtained by solving the Boltzmann kinetic equation in relaxation time approximation.However, we have not considered the finite size of the nuclei and ion structure function for the calculation of the relaxation rate.ρ [g cm -3 ] ρ [g cm -3 ] The electromagnetic interaction between electrons and ions have been incorporated through HDL propagator in the calculation.The calculation should account for magnetically modified anisotropic photon propagator; however in the present paper, we have considered only isotropic HDL propagator.For the screening mass zero temperature magnetically modified Debye mass has been incorporated but in the hot and dense plasma, inclusion of finite temperature screening mass would be more relevant.
We have found that the inclusion of the HDL propagator in the relativistic plasma reduces the values of both σ and κ in contrast to static screening in the nonrelativistic plasma.The frequency dependent screening in the HDL propagator enhances the interaction rate.The off equilibrium distribution function is inversely proportional to the interaction rate and the transport coefficients are directly proportional to Φ. Hence, enhancement in the interaction rate decrease both the transport coefficients at high density in the degenerate regime.We also observe cubic temperature scaling between κ and σ violating the Wiedemann-Franz law.
A rough estimation of diffusion and thermal equilibration time scales from the expressions of quantized σ, κ including frequency dependent screening have also been attempted.The obtained Ohmic time scale has found out to be of the same order as the survival time of the merged compact object due to reduction in the value of σ (for λ B cm).The thermal equilibration time scale can also be found to match with the time scale of the merged object when ∆z ∼ cm.One can infer from the analysis that many-body effects play an important role in determining the dissipative time scales relevant in the neutron star merger in the low density, high temperature and high magnetic field regime if the length-scales are of the order mentioned above.This estimation of dissipation time scales of the transport coefficients requires rigorous hydrodynamic formulation for more realistic values of λ B and ∆z.The current calculation of σ and κ can be implemented in modeling the magneto-thermal evolution of the merged compact object as well.
In the current paper, we have considered only longitudinal component of the transport coefficients (and ignored all the other tensorial components) in the equation for magnetic-field evolution in plasma.The realistic estimate of τ σ and τ κ can only be obtained if all the components of conductivity tensor are known in the background of relevant equation of state for neutron star merger.However, our calculations presented in this paper provide a significant step towards conjoining the complex microphysical effects in plasma with the GRMHD simulations.The present formalism of trans-port theory can be extended in relativistic, magnetised QCD matter with certain modifications like inclusion of QCD coupling constant, diagrams of strong interaction and proper vertex corrections.
where, we have used |p | = |p−q| = p − p. q.The angular integrals of k can be expressed as, Eq.( 30) can further be written as, where, we have changed the variable dy B = dq x /mω B following momentum conservation p x − p x = q x .In the above equation, n i is the number density of ions which can be expressed in terms of electron number density as n i = n e /Z.Further, n e can be expressed in terms of the Debye mass as n e = µm 2 D /3e 2 .Next we introduce a dimensionless variable y = q z /| q|, In order to calculate |M| 2 , we use following electronic spinor in presence of magnetic field, Using above spinors and the expression for photon propagator (Eq.( 15)) in Eq. (8) in Section.2 one obtains, where, s and s are ±, and exp(iq y y) Hn (ξ ) Hn (ξ)dy.
H n (ξ) is the Hermite polynomial, ξ = √ mω B .Inserting I n ,n in Eq.( 36) and performing the y integration we obtain, where, η = ±, F n ,n (u x + q 2 y ).The functions L n−n n (u) are Laguerre polynomials and F n ,n (u) are normalized to ∞ 0 F 2 n ,n du = 1.Now, to perform the integration in y in Eq.( 34) we use the following sum rule, .
In the limiting case, Re Π T,L (y=∞) = m 2 D /3, Re Π T (y=0) = 0, Re Π L (y=0) = m 2 D .Using the above relations the interaction rate becomes, where, We change the variable q y to u and perform the integration as follows, Finally, the particle scattering rate becomes, where, ζ = m 2 D /2mω B .

Appendix B
In this appendix, we show the compariosn of the conductivities considering two equation of states (EOSs), BPS (Baym, Pethick and Surtherland) [50] and magnetic BPS model [51] and comment on considering Fe and Mo for numerical analysis in the paper.In ref. [50] the equation of state of zero-temperature matter in complete nuclear equilibrium is given for mass densities below 5 × 10 14 g cm −3 .In Ref. [51], the BPS [50] has been extended to include the physical parameters for a low density plasma in presence of high magnetic field relevant for neutron star crust.In the fig.7, we plot σ and κ with ρ for these two different EOSs.From the plot fig.7 it is evident that the BPS curve matches well with the Mo curve for both σ and κ.We also include Fe for reference.ρ [g cm -3 ] ρ [g cm -3 ] Fig. 7: The variation of σ (left) and κ (right) with ρ for different elements.The magnetic field is chosen as 10 17 G and temperature as 2 × 10 10 T. We show the behaviour of the transport coefficients for two different EOS s given in Refs.[50,51]

Fig. 1 :
Fig. 1: Feynman diagram describing the electron-ion scattering amplitude in a static medium (left panel).Feynman diagram contributing to the electron-ion scattering amplitude via exchange of resummed HDL photon propagator (right panel).
Fig.2shows the variation of σ and κ given in Eqs.(24) and(25) with ρ for different temperatures and different magnetic fields for different elements Fe and Mo.In order to consider electrons to be relativistic, the density and temperatures are chosen as ρ 10 6 g cm −3 and T 5 × 10 9 K respectively.For fixed B, µ increases with ρ and electrons start to populate higher Landau levels.Since, we are interested in population of the zeroth Landau level, the density and temperature should also satisfy ρ < ρ B and T T ce as given in Eq.(28).With these two conditions, both the coefficients have been obtained by numerically integrating the expressions in Eqns.24 and 25.In fig.2, the upper and lower panel plots are for the variation of σ and κ with ρ respectively.The left panel plots have been drawn considering magnetic field B=10 17 G and 7.5 × 10 16 G in the right panel.Each plot has three curves corresponding to three different temperatures 1.4×10 11 K, 2×10 11 K and 5 × 10 11 K for each element.This is observed from the

Fig. 2 :
Fig. 2: The comparison of σ (upper plots) and κ (lower plots) with ρ for Fe (solid lines) and Mo (dotted lines).The left panel plots are for magnetic field 10 17 G and the right ones are for B = 7.5 × 10 16 G.In each plot we have three curves for each element corresponding to three different temperatures 1.4 × 10 11 K, 2 × 10 11 K, 5 × 10 11 K.

s - 1 ]Fig. 5 : 4 )Fig. 6 :
Fig. 5: Plots of σ (upper panel) and κ (lower panel) with ρ comparing static, Coulomb screening in non-relativistic plasma (blue) and dynamical screening in relativistic plasma (black) at fixed temperature of 2 × 10 10 K.The left panel curves are for B= 10 17 G and the right curves are for B= 5 × 10 16 G.The choice of elements are Mo (dotted lines) and Fe (solid lines).