Search NASA⌕ Search

SEARCH · Search NASA

Results for “MAGNETOHYDRODYNAMICS”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12

Numerical assessment of the impact of the guiding-centre approximation on fast ion simulations in NSTX

Guiding-centre (GC) and full-orbit (FO) simulations of the beam-injected fast ion distribution and the corresponding neutron emissivity have been carried out for magnetohydrodynamics-quiescent National Spherical Torus eXperiment (NSTX) plasmas, using a combination of ASCOT5 and DRESS, to assess the suitability of the GC approximation for fast ions in NSTX. It was found that GC and FO simulations predicted substantially different steady-state distributions in both position and velocity space and different neutron emissivity profiles, leading to a 15% reduction in the predicted global neutron rate for FO relative to GC. These changes accompany a higher magnetic moment in FO, and correspond to a change in particle orbits from co-passing to trapped and stagnation orbits. ASCOT5 was also benchmarked against TRANSP/NUBEAM with input loaded entirely from TRANSP/NUBEAM output files, with agreement found between the GC simulations when finite Larmor radius (FLR) corrections were omitted. ASCOT5 FO and TRANSP/NUBEAM with FLR produced fast ion distributions which differed in localised regions, but predicted global neutron rates which agree within 3%.

ASCOT↗

Effects of plasma squareness on RMP response in MAST-U

The effect of the squareness of the plasma boundary shape on the plasma response to the applied n = 2 (n is the toroidal mode number) resonant magnetic perturbation, and consequently on controlling the edge-localized mode (ELM) in MAST-U in the double-null configuration, is numerically systematically investigated. The modeling is performed utilizing three different plasma response models: the single-fluid model, the two-fluid model and the magnetohydrodynamic-kinetic hybrid model. The effect of the plasma squareness on ELM control is found to depend on the equilbrium kinetic profiles. Fixing the equilibrium radial profiles to that from a recent MAST-U discharge, the experimental plasma shape is found to be close to optimal for ELM control in terms of squareness. The optimal coil phasing does not significantly vary with varying plasma squareness. These findings robustly hold with (i) two figures of merit (FoMs) defined for predicting the ELM control feasibility, and (ii) three plasma response models employed for evaluating these FoMs.

MAST-U↗

Variability of MHD instabilities in benign termination of high-current runaway electron beams in the JET and DIII-D tokamaks

Benign termination, in which magnetohydrodynamic (MHD) instabilities deconfine runaway electrons (REs) following hydrogenic injections, is a promising strategy for mitigating dangerous RE loads after disruptions. Recent experiments on the Joint European Torus (JET) have explored this scenario at higher pre-disruptive plasma currents than are achievable on other devices, revealing challenges in obtaining benign terminations at I p ≥2.5 MA. This work analyzes the evolution of these high-current RE beams and their terminating MHD events using fast magnetic sensor measurements and EFIT equilibrium reconstructions for approximately 40 JET and 20 DIII-D tokamak discharges. On JET, unsuccessful non-benign terminations occur at low edge safety factor (q edge ≈ 2), and are preceded by intermittent, non-terminating MHD events at higher rational qedge. Trends in the internal inductance I i indicate more peaked RE current profiles in the high-I p non-benign population, which may hinder successful recombination through re-ionization of the companion plasma. In contrast, benign terminations on JET typically occur at higher q edge ≥3 and exhibit less peaked RE current profiles. DIII-D displays a broader range of terminating edge safety factors, again correlated with the measured values. Across both tokamaks, the RE current peaking is therefore found to determine which MHD instability boundary is encountered, a result confirmed by linear resistive MHD modeling with the CASTOR3D code. Measured growth rates are similar for benign and non-benign cases, indicating that ideal MHD timescales at low density after hydrogenic injection do not alone explain efficient RE deconfinement. Instead, non-benign cases are most readily characterized by their comparably lower overall MHD perturbation amplitudes δB. These observations suggest that the interplay between ideal and resistive dynamics governs the termination process, with implications for extrapolating benign RE termination to high-I p reactor scenarios.

MHD instabilities↗

Physics basis for the reference flat-top plasma scenario in the ST–E1 fusion power plant

As part of the U.S. Department of Energy’s Milestone-Based Fusion Energy Development Program, Tokamak Energy has completed the pre-concept design of the ST–E1 fusion power plant. ST–E1 is envisaged to operate in two phases: a pilot plant phase, targeting sustained net power production of 300 - 500 MWe for a duration >1 hr, followed by a commercial power plant phase targeting steady-state operations and a normalised overnight capital cost of ⩽12 000 $\$$/kWe. The design process adopted was highly iterative, integrating all major plant systems and progressing in a phased fidelity approach. At the pre-conceptual stage, the emphasis has been on exploring the design space, identifying the main system-level trade-offs, and making the key decisions that define the overall plant concept, rather than optimising a single operating point. This paper, part of a focused collection detailing the ST–E1 pre-concept design, addresses the development of a series of reference flat-top plasma operating points for the pilot plant phase. A modelling workflow was established to develop and assess candidate plasma design points and explore key dependencies. The workflow includes integrated core plasma modelling, magnetohydrodynamic (MHD) stability assessment, equilibrium generation, scrape-off-layer and exhaust modelling, heating & current drive design and optimisation, and turbulent transport modelling. Using this framework, the impact of several key parameters on the flat-top operating space was investigated, including the density limit, core radiation fraction and divertor power loading, level of external heating and curent drive power and assumed pedestal characteristics. The MHD stability, controllability and micro-stability characteristics of these plasmas were also analysed. These investigations informed the definition of a set of fully non-inductive, flat-top reference operating points that satisfy the high-level ST–E1 mission, including a low and high density case, a case that is stable to resistive wall modes and a case with reduced divertor power loading.

ST–E1↗

Recent advances in plasma control and physics research in the Large Helical Device

The Large Helical Device (LHD), the largest superconducting helical system in the world, is equipped with advanced heating and diagnostic tools, facilitating plasma control and physics research. Data assimilation was employed for electron temperature control using a real-time Thomson scattering system and real time prediction code. A virtual LHD environment enabled visualization of escaping high-energy tritium ions and demonstrated that these ions impact the rear side of the divertor plate. Pioneering results crucial to plasma control have also been achieved. Real-time wall conditioning using Lithium granule dropping improved bulk ion energy and particle transport while simultaneously enhancing the heavy impurity transport. Progress has also been made in the investigation of turbulence-driven transport. At the confinement bifurcation, ion-scale turbulence decreased, while electron-scale turbulence increased. A change in the anisotropy of turbulent eddies was also observed at the confinement bifurcation. Coexistence of local and non-local turbulence was identified in electron-scale turbulence. Non-local turbulence exhibited the rapid spatial propagation of perturbations throughout the plasma, while local turbulence followed the temperature gradient. A transition between drift-wave turbulence and magnetohydrodynamics (MHD) turbulence was observed with the turbulence minimized at the transition condition. Machine learning analysis was employed to evaluate the temperate and density conditions of this turbulence transition. Then, real-time control of fueling and heating was applied to maintain the turbulence transition condition, improving the energy confinement enhancement factor by 20%. In addition, evidence was obtained for collisionless ion heating by energetic-ion-driven geodesic acoustic modes and MHD bursts. These achievements represent unique contributions to the development of fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Toroidal phase of post-disruption runaway electron loss to wall in presence of applied 3-D fields in DIII-D

The toroidal distribution of runaway electrons (REs) striking the centre post in DIII-D post-disruption RE final loss events is simulated utilizing the MARS-F code combined with the REORBIT module, by tracing the guiding-centre drift orbits of test REs in the presence of external (applied) plus internal (intrinsic instability) 3-D fields. To better recover the experimental results, three different equilibria with different safety factor profiles are adopted. Dominant resistive magnetohydrodynamic instabilities are found to be m/n=2/1, 3/2 and 1/1 modes (where m and n are the poloidal and toroidal mode numbers). Externally applied n=1 resonant magnetic perturbations (RMPs), with both even and odd parity configuration, are considered in each equilibrium. Locking of MHD instability to the applied RMPs is simulated by minimizing the total perturbed magnetic energy in the plasma. Here, the simulated toroidal locations, where most REs strike in the presence of an edge m/n= 2/1 instability and the n=1 RMPs in even parity, are found to agree well with experiments. Modeling also captures the measured n = 2 RE loss pattern on the wall for experiments where an edge m/n=3/2 instability is dominant instead of the m/n=2/1 mode. The key overall findings include (a) the measured global RE toroidal impact pattern in experiments can be explained by a locked resistive kink instability, (b) the toroidal phase of the locked mode is well predicted by minimizing the total perturbed magnetic energy, and (c) this edge locked mode determines the toroidal phase of the peak RE wall impact.

Bai, Xue [University of California San Diego, La J↗

Transfer learning nonlinear plasma dynamic transitions in low dimensional embeddings via deep neural networks

Deep learning algorithms provide a new paradigm to study high-dimensional dynamical behaviors, such as those in fusion plasma systems. Development of novel, data-driven model reduction methods, coupled with detection of abnormal modes with plasma physics, opens a unique opportunity to identify plasma instabilities through automated construction of parsimonious models that can be tuned to balance accuracy and cost. Our fusion transfer learning (FTL) model demonstrates success in rapidly reconstructing nonlinear kink mode structures by learning from a limited amount of nonlinear simulation data. The knowledge transfer process leverages a pre-trained neural encoder–decoder network, initially trained on linear simulations, to effectively capture nonlinear dynamics. The low-dimensional embeddings extract the coherent structures of interest, while preserving the inherent dynamics of the complex system. Experimental results highlight FTL’s capacity to capture transitional behaviors and dynamical features in plasma dynamics—a task often challenging for conventional methods. The model developed in this study is generalizable and can be extended broadly through transfer learning to address various magnetohydrodynamics modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Stellar wind impact on early atmospheres around unmagnetized Earth-like planets

ABSTRACT Stellar rotation at early ages plays a crucial role in the survival of primordial atmospheres around Earth-mass exoplanets. Earth-like planets orbiting fast-rotating stars may undergo complete photoevaporation within the first few hundred Myr driven by the enhanced stellar XUV [X-rays and extreme ultraviolet (EUV)] radiation, while planets orbiting slow-rotating stars are expected to experience difficulty in losing their primordial envelopes. Besides the action of stellar radiation, stellar winds induce additional erosion on these primordial atmospheres, altering their morphology, extent, and causing supplementary atmospheric losses. In this paper, we study the impact of activity-dependent stellar winds on primordial atmospheres to evaluate the extent to which the action of these winds can be significant in the whole planetary evolution at early evolutionary stages. We performed 3D magnetohydrodynamical (MHD) simulations of the interaction of photoevaporating atmospheres around unmagnetized Earth-mass planets in the time span between 50 and 500 Myr, analysing the joint evolution of stellar winds and atmospheres for both fast- and slow-rotating stars. Our results reveal substantial changes in the evolution of primordial atmospheres when influenced by fast-rotating stars, with a significant reduction in extent at early ages. In contrast, atmospheres embedded in the stellar winds from slow-rotating stars remain largely unaltered. The interaction of the magnetized stellar winds with the ionized upper atmospheres of these planets allows us to evaluate the formation and evolution of different MHD structures, such as double bow shocks and induced magnetospheres. This work will shed light on the first evolutionary stages of Earth-like exoplanets, which are of crucial relevance in terms of planet habitability.

Astronomy & Astrophysics↗

AuriDESI: mock catalogues for the DESI Milky Way Survey

The Dark Energy Spectroscopic Instrument Milky Way Survey (DESI MWS) will explore the assembly history of the Milky Way by characterizing remnants of ancient dwarf galaxy accretion events and improving constraints on the distribution of dark matter in the outer halo. We present mock catalogues that reproduce the selection criteria of MWS and the format of the final MWS data set. These catalogues can be used to test methods for quantifying the properties of stellar halo substructure and reconstructing the Milky Way’s accretion history with the MWS data, including the effects of halo-to-halo variance. The mock catalogues are based on a phase-space kernel expansion technique applied to star particles in the Auriga suite of six high-resolution lambda-cold dark matter magnetohydrodynamic zoom-in simulations. They include photometric properties (and associated errors) used in DESI target selection and the outputs of the MWS spectral analysis pipeline (radial velocity, metallicity, surface gravity, and temperature). They also include information from the underlying simulation, such as the total gravitational potential and information on the progenitors of accreted halo stars. We discuss how the subset of halo stars observable by MWS in these simulations corresponds to their true content and properties. These mock Milky Ways have rich accretion histories, resulting in a large number of substructures that span the whole stellar halo out to large distances and have substantial overlap in the space of orbital energy and angular momentum.

dynamics↗

Magnetorotational dynamo can generate large-scale vertical magnetic fields in 3D GRMHD simulations of accreting black holes

ABSTRACT Jetted astrophysical phenomena with black hole engines, including binary mergers, jetted tidal disruption events, and X-ray binaries, require a large-scale vertical magnetic field for efficient jet formation. However, a dynamo mechanism that could generate these crucial large-scale magnetic fields has not been identified and characterized. We have employed three-dimensional global general relativistic magnetohydrodynamical simulations of accretion discs to quantify, for the first time, a dynamo mechanism that generates large-scale magnetic fields. This dynamo mechanism primarily arises from the non-linear evolution of the magnetorotational instability (MRI). In this mechanism, large non-axisymmetric MRI-amplified shearing wave modes, mediated by the axisymmetric azimuthal magnetic field, generate and sustain the large-scale vertical magnetic field through their non-linear interactions. We identify the advection of magnetic loops as a crucial feature, transporting the large-scale vertical magnetic field from the outer regions to the inner regions of the accretion disc. This leads to a larger characteristic size of the, now advected, magnetic field when compared to the local disc height. We characterize the complete dynamo mechanism with two time-scales: one for the local magnetic field generation, $t_{\rm gen}$, and one for the large-scale scale advection, $t_{\rm adv}$. Whereas the dynamo we describe is non-linear, we explore the potential of linear mean field models to replicate its core features. Our findings indicate that traditional $\alpha$-dynamo models, often computed in stratified shearing box simulations, are inadequate and that the effective large-scale dynamics is better described by the shear current effects or stochastic $\alpha$-dynamos.

Astronomy & Astrophysics↗

The three-dimensional structure of black hole accretion flows within the plunging region

ABSTRACT We analyse, using new analytical models and numerical general relativistic magnetohydrodynamic simulations, the three-dimensional properties of accretion flows inside the plunging region of black hole spacetimes (i.e. at radii smaller than the innermost stable circular orbit). These simulations are of thick discs, with aspect ratios of order unity $h/r \sim 1$, and with a magnetic field geometry given by the standard low-magnetization ‘SANE’ configuration. This work represents the first step in a wider analysis of this highly relativistic region. We show that analytical expressions derived in the ‘thin disc’ limit describe the numerical results remarkably well, despite the large aspect ratio of the flow. We further demonstrate that accretion within this region is typically mediated via spiral arms, and that the geometric properties of these spiral structures can be understood with a simple analytical model. These results highlight how accretion within the plunging region is fundamentally two-dimensional in character, which may have a number of observational implications. We derive a modified theoretical description of the pressure within the plunging region which accounts for turbulent heating and may be of use to black hole image modelling.

Astronomy & Astrophysics↗

The influence of the cloud virial parameter on the initial mass function

ABSTRACT Crucial for star formation is the interplay between gravity and turbulence. The observed cloud virial parameter, $\alpha _{\mathrm{vir}}$, which is the ratio of twice the turbulent kinetic energy to the gravitational energy, is found to vary significantly in different environments, where the scatter among individual star-forming clouds can exceed an order of magnitude. Therefore, a strong dependence of the initial mass function (IMF) on $\alpha _{\mathrm{vir}}$ may challenge the notion of a universal IMF. To determine the role of $\alpha _{\mathrm{vir}}$ on the IMF, we compare the star-particle mass functions obtained in high-resolution magnetohydrodynamical simulations including jet and heating feedback, with $\alpha _{\mathrm{vir}}=0.0625$, 0.125, and 0.5. We find that varying $\alpha _{\mathrm{vir}}$ from $\alpha _{\mathrm{vir}}\sim 0.5$ to $\alpha _{\mathrm{vir}}< 0.1$ shifts the peak of the IMF to lower masses by a factor of $\sim 2$ and increases the star formation rate by a similar factor. The dependence of the IMF and star formation rate on $\alpha _{\mathrm{vir}}$ is non-linear, with the dependence subsiding at $\alpha _{\mathrm{vir}}< 0.1$. Our study shows a systematic dependence of the IMF on $\alpha _{\mathrm{vir}}$. Yet, it may not be measurable easily in observations, considering the uncertainties, and the relatively weak dependence found in this study.

Mathew, Sajay Sunny (ORCID:0000000283818195)↗

Auriga Streams – I: disrupting satellites surrounding Milky Way-mass haloes at multiple resolutions

In a hierarchically formed Universe, galaxies accrete smaller systems that tidally disrupt as they evolve in the host’s potential. We present a complete catalogue of disrupting galaxies accreted onto Milky Way-mass haloes from the Auriga suite of cosmological magnetohydrodynamic zoom-in simulations. We classify accretion events as intact satellites, stellar streams, or phase-mixed systems based on automated criteria calibrated to a visually classified sample, and match accretions to their counterparts in haloes re-simulated at higher resolution. Most satellites at the present day have lost substantial amounts of stellar mass – 67 per cent have $f_\text{bound} < 0.97$ (our threshold of lost stellar mass to no longer be considered intact), while 53 per cent satisfy a more stringent $f_\text{bound} < 0.8$. Streams typically outnumber intact systems, contribute a smaller fraction of overall accreted stars, and are substantial contributors at intermediate distances from the host centre ($\sim$0.1 to $\sim 0.7R_\text{200m}$, or $\sim$35 to $\sim$250 kpc for the Milky Way). We also identify accretion events that disrupt to form streams around massive intact satellites instead of the main host. Streams are more likely than intact or phase-mixed systems to have experienced pre-processing, suggesting this mechanism is important for setting disruption rates around Milky Way-mass haloes. All of these results are preserved across different simulation resolutions, though we do find some hints that satellites disrupt more readily at lower resolution. The Auriga haloes suggest that disrupting satellites surrounding Milky Way-mass galaxies are the norm and that a wealth of tidal features waits to be uncovered in upcoming surveys.

galaxies: haloes↗

2D end-to-end modelling of kilonovae from binary neutron star merger remnants

We investigate the kilonova emission resulting from outflows produced in a 3D general-relativistic magnetohydrodynamic (GRMHD) simulation of a hypermassive neutron star (HMNS) remnant. We map the outflows into the flash hydrodynamics code to model their expansion in axisymmetry, and study the effects of employing different r-process heating rates. Except for the highest heating rate prescription, we find no significant differences with respect to overall ejecta dynamics and morphology compared to the simulation without heating. Once homologous expansion is attained, typically after ~2s for these ejecta, we map the outflows to the sedona radiative transfer code and compute the spectral evolution of the kilonova and broad-band light curves in various Legacy Survey of Space and Time (LSST) bands. The kilonova properties depend on the remnant lifetime, with peak luminosities and peak time-scales increasing for longer lived remnants that produce more massive ejecta. For all models, there is a strong dependence of both the bolometric and broad-band light curves on the viewing angle. While the short-lived (12 ms) remnant produces higher luminosities when viewed from angles closer to the pole, longer lived remnants (240 ms and 2.5 s) are more luminous when viewed from angles closer to the equator. Our results highlight the importance of self-consistent, long-term modelling of merger ejecta, and taking viewing-angle dependence into account when interpreting observed kilonova light curves. We find that magnetized outflows from an HMNS – if it survives long enough – could explain blue kilonovae, such as the blue emission seen in AT2017gfo.

79 ASTRONOMY AND ASTROPHYSICS↗

Collapse of magnetized white dwarfs as site of heavy-element formation and kilonova signal

We present the first end-to-end calculation connecting the accretion-induced collapse (AIC) of a magnetized, rapidly rotating white dwarf to observable kilonova signatures, combining two-dimensional (2D) general-relativistic neutrino-magnetohydrodynamic simulations, followed by radiation hydrodynamics with in-situ nuclear network and 2D Monte Carlo radiative transfer with spatially resolved heating rates. Unlike all previous unmagnetized AIC models – which predicted proton-rich, $^{56}$Ni-dominated ejecta – strong magnetic fields eject ${\approx }\, 0.2\, \mathrm{ M}_\odot$ of neutron-rich material ($\langle Y_e \rangle \sim 0.24$) on dynamical time-scales, before neutrino irradiation can raise the electron fraction, enabling strong r-process nucleosynthesis up to and beyond the third peak. The resulting kilonova is lanthanide-rich ($X_{\rm lan} \approx 8~{{\ \rm per\ cent}}$) and dominated by near-infrared emission. We compute synthetic light curves in the Large Synoptic Survey Telescope and J ames Webb Space Telescope bands and find striking agreement, without parameter tuning, between the observations of AT 2023vfi/GRB 230307A and our broadband light curves for polar viewing angles. These results establish magnetized AIC as a viable channel for heavy r-process element production and a compelling progenitor candidate for long-duration gamma-ray bursts with kilonova signatures.

MHD↗

Probing the axion-electron coupling with NuSTAR observations of galaxies

We search for the existence of ultralight axions coupling to electrons and photons using data from the NuSTAR telescope directed toward the galaxies M82, M87, and M31. We focus on electron bremsstrahlung and Compton scattering for axion production in stars, summing over the stellar populations found in the target galaxies when computing the axion luminosity. We then compute the hard x-ray signal that arises from the conversion of these axions to photons in each galaxy’s magnetic fields, inferred from analog galaxies in cosmological magnetohydrodynamic simulations. Analyzing NuSTAR data toward these galaxies between roughly 20 to 70 keV, we find no evidence for axions and set leading constraints on the combined axion-electron and axion-photon coupling at the level of |𝑔 𝑎⁢𝑒⁢𝑒 ×𝑔 𝑎⁢𝛾⁢𝛾 | ≲ 8.3 ×10 −27 GeV −1 for 𝑚 𝑎 ≲ 10 −10 eV at 95% confidence, with M82 providing the most stringent constraints.

Axion-like particles↗

Gamma-ray signatures of 𝑟-process radioactivity from the collapse of magnetized white dwarfs

We predict the gamma-ray line emission from 𝑟-process nuclei synthesized in the ejecta of the accretion-induced collapse (AIC) of a magnetized, rapidly rotating white dwarf. Using ejecta from a two-dimensional general-relativistic neutrino-magnetohydrodynamic simulation, further evolved with a radiation-hydrodynamics code coupled to an in situ nuclear reaction network, we construct angle-dependent gamma-ray spectra in the 0.01–10 MeV band via composition-dependent ray tracing through the ejecta. The emission between ∼1 and 10 d is dominated by 132 I (𝑡 1/2 = 2.3 h), continuously replenished by the decay of its parent 132 Te (𝑡 1/2 = 3.2 d), with additional contributions from 131 I, 133 Xe, and 132 Te . At 𝑡 ≳ 20 d, 56 Co (from 56 Ni decay) becomes the primary emitter. The simultaneous presence of 𝑟 process and iron-peak gamma-ray lines is distinctive of AIC ejecta and absent in binary neutron star mergers, where iron-peak nuclei are generally not synthesized. Comparing with the 3⁢𝜎 continuum sensitivities of planned MeV gamma-ray telescopes (COSI, AMEGO-X, e-ASTROGAM, GRAMS, GammaTPC), we find the brightest 𝑟-process lines detectable to ∼10 Mpc by GammaTPC and GRAMS, with the signal approaching their sensitivity threshold at 30 Mpc. As a result, the 𝑟-process spectral features survive time integration over ∼30 d exposures, demonstrating robustness against the long observation times required by gamma-ray detectors.

Nuclear reactions↗

Electrical conductivity of a warm neutron star crust in magnetic fields: Neutron-drip regime

We compute the anisotropic electrical conductivity tensor of the inner crust of a compact star at nonzero temperature by extending a previous work on the conductivity of the outer crust. The physical scenarios, where such crust is formed, involve protoneutron stars born in supernova explosions, binary neutron star mergers, and accreting neutron stars. The temperature-density range studied covers the transition from a semidegenerate to a highly degenerate electron gas and assumes that the nuclei form a liquid, i.e., the temperature is above the melting temperature of the lattice of nuclei. The electronic transition probabilities include (i) the screening of electron-ion interaction in the hard-thermal-loop approximation for the QED plasma, (ii) the correlations of the ionic component in a one-component plasma, and (iii) finite nuclear size effects. The conductivity tensor is obtained from the Boltzmann kinetic equation in relaxation time approximation accounting for the anisotropy introduced by a magnetic field. The sensitivity of the results towards the matter composition of the inner crust is explored by using several compositions of the inner crust, which were obtained using different nuclear interactions and methods of solving the many-body problem. The standard deviations of relaxation time and components of the conductivity tensor from the average are below ≤25% except close to crust-core transition, where nonspherical nuclear structures are expected. Finally, our results can be used in dissipative magnetohydrodynamics simulations of warm compact stars.

Physics↗