Quick evaluation of electrical and thermal conductivities of degenerate fully ionized plasmas without magnetic field
|Department of Theoretical Astrophysics
conductivity codes have been updated in September 2019.
We provide here three programs: condegin19, condegsc19, condBSk.
In the first of them,
we assume that the charge distribution of an ion is steplike
(that is, uniform within the radius rnuc).
The characteristic dimensionless parameter that determines
importance of the finite size of the ion is the ratio
of rnuc to the ion-sphere radius
In the inner crust (where dripped neutrons are present), we use an interpolation formula by Itoh & Kohyama (1983, ApJ 275, 858):
The key subroutine of condegin19 is CONDEGIN, whose input parameters are temperature T, number density of ions ni, magnetic field B, ion charge number Z, ion mass number A, the number of neutrons per nucleus (in a Wigner-Seitz cell) A', and the effective impurity charge Zimp (which is the square root of the usual impurity parameter Q). In general, A' may be greater than A because of the presence of dripped neutrons. The output parameters are the longitudinal, transverse, and off-diagonal components of the electrical and thermal conductivity tensors. All dimensional input and output parameters of CONDEGIN are in relativistic units. The conversion from/to conventional units is exemplified in the auxiliary MAIN program.
The other two codes, condegsc19 and condBSk, take into account that the proton charge profile in a nucleus is not steplike, but goes from a finite value to zero continuously near rnuc. In "condegsc19", we use the charge profile parametrization worked out by Oyamatsu (1993, Nucl. Phys. A 561, 431). The corresponding nuclear shape parameters, as well as Z, A, and A', are interpolated as functions of baryon density throughout the neutron-star crust. This interpolation (named the smooth composition model by Kaminker et al. (1999)) is performed in the subroutine OYAFORM. The parameter xnuc and an auxiliary parameter x'nuc, which enter the fit for the effective Coulomb logarithm, are expressed through the nuclear shape parameters as explained in Appendix A of Ref.. Accordingly, the key subroutine CONDEGsc has smaller number of free input arguments than CONDEGIN described above: the input parameters of CONDEGsc are log10(T [K]), log10(rho [g/cc]), B12=B/1012 G, and Zimp. The output parameters are the longitudinal, transverse, and off-diagonal components of the electrical and thermal conductivity tensors in CGS units.
CONDBSK is analogous to CONDEGsc, but instead of Oyamatsu's profiles we use the nuclear form factors worked out recently in a series of works by N.Chamel and coauthors. A detailed paper on the used parametrizations, realized in subroutines FRACRUST and NUCPAR of the latter code, is in preparation.
Please communicate to us your opinion about this resource or send any suggestions or remarks.
Fortran codes to download (at a choice):