pith. sign in

arxiv: 2606.10750 · v1 · pith:CXV3T2RLnew · submitted 2026-06-09 · 🌌 astro-ph.IM · cond-mat.mtrl-sci

UnReal-B : Real-Space DFT Solver for Matter in Extreme Magnetic Fields

Bhalchandra S. Pujari , Andrey Tokarev , Dipanjan Mitra This is my paper

Pith reviewed 2026-06-27 11:48 UTC · model grok-4.3

classification 🌌 astro-ph.IM cond-mat.mtrl-sci
keywords density functional theoryneutron starsstrong magnetic fieldsadiabatic approximationelectronic structurereal-space methodscondensed matter
0
0 comments X

The pith

UnReal-B is a real-space DFT solver that calculates the electronic structure of one-dimensional atomic chains in magnetic fields from 10^12 to 10^15 gauss using the adiabatic approximation.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper introduces UnReal-B as an open-source real-space density functional theory solver designed for one-dimensional chains of matter under the extreme magnetic fields found on neutron stars. It relies on the adiabatic approximation to simplify the numerical treatment of electronic structure. The implementation is benchmarked on astrophysically relevant elements and reproduces published results with good accuracy. The code is released openly to support reproducible work on neutron-star surface modeling.

Core claim

UnReal-B provides a streamlined numerical framework for calculating the electronic structure of strongly magnetized condensed matter by employing the adiabatic approximation in a real-space Density Functional Theory solver for one-dimensional chains. It demonstrates excellent agreement with published results for several astrophysically relevant elements while maintaining a comparatively simple and transparent implementation.

What carries the argument

The adiabatic approximation applied inside a real-space density functional theory solver for one-dimensional atomic chains.

If this is right

  • The solver enables reproducible calculations of neutron-star surface matter for elements already benchmarked.
  • Open-source release supports community extensions driven by new observational data on neutron-star surfaces.
  • The framework supplies a transparent base for adding further numerical capabilities while preserving the adiabatic simplification.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same code structure could be applied to other one-dimensional systems outside astrophysics if the magnetic-field strengths remain comparable.
  • Direct tests against full three-dimensional calculations at selected field strengths would map the approximation's breakdown points.
  • Integration with radiative-transfer codes could link the computed electronic structures to predicted spectra from neutron-star atmospheres.

Load-bearing premise

The adiabatic approximation remains valid and accurate for the electronic structure calculations across the full range of magnetic fields considered.

What would settle it

A side-by-side comparison of energy levels or charge densities obtained from UnReal-B and from a calculation that drops the adiabatic approximation at B = 10^15 G would show large differences if the approximation fails.

Figures

Figures reproduced from arXiv: 2606.10750 by Andrey Tokarev, Bhalchandra S. Pujari, Dipanjan Mitra.

Figure 1
Figure 1. Figure 1: a) Convergence against k-points. b) Computational scaling c) Convergence against z-grid. (B_field) is specified in Gauss, which determines the mag￾netic length and Landau level spacing. The nuclear charge (atomic_Z) and num_electrons define the atomic species. The lattice_constant set the periodicity of the 1D chain. The model’s complexity is further controlled by M_max (the num￾ber of magnetic orbitals) a… view at source ↗
Figure 2
Figure 2. Figure 2: Example of EOS data for Fe, B12= 500. Points and solid line represent the data while the dotted line is for the quadratic fit. 0.0 0.2 0.4 0.6 0.8 1.0 k / ( /a) 1000 800 600 400 200 0 200 400 Energy (eV) Bands & DOS C chain | B = 1.0×10¹² G 0.00 0.02 0.04 DOS [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Band Structure of carbon chain. All bands corre￾sponds to 𝜈 = 0 state. The colors indicate the 𝑚 values, with red being larger value of 𝑚. The fermi level is set at 0 eV. strictly linear runtime dependence on 𝑁𝑘 due to the indepen￾dent nature of the 𝑘-point loops, the runtime scales super￾linearly with 𝑁𝑧 . This behavior reflects the increasing com￾putational complexity of the real-space operations—such as… view at source ↗
Figure 4
Figure 4. Figure 4: Band Structure of iron chain under 𝐵12 = 500. The solid lines indicate the 𝜈 = 0 states and the dashed lines are for 𝜈 = 1 states. The colors indicate the 𝑚 values, with yellow being largest value of 𝑚. The fermi level is set at 0 eV and the number of bands with 𝜈 = 0 are restricted below fermi level. Beyond the electronic quantum numbers, the macro￾scopic periodicity of the 1D chain—represented by the lat… view at source ↗
Figure 6
Figure 6. Figure 6: C B12=1. Fe B12 = 500. density The Density of States for the Fe chain underscores the extreme degeneracy of the system, with a massive density peak pinned at the Fermi level. The ability of UnReal-B to resolve this dense manifold of states while maintaining numerical stability based linear FFT convolution approach. This degeneracy can further be seen via band occupations. Fig (5) compares the occupation fr… view at source ↗
Figure 7
Figure 7. Figure 7: Scaling of the UnReal-B total energy per cell with magnetic field strength. (a) |𝐸∞| versus 𝐵12 on logarithmic axes for H, He, C and Fe; symbols are computed values (FFT solver), grey lines are guides of slope 2∕5, and fitted exponents 𝛼 (from |𝐸∞| ∝ 𝐵 𝛼 12) are listed in the legend. (b) Compensated coefficient |𝐸∞|∕(𝑍9∕5𝑏 2∕5); its near-constancy confirms the Landau–Thomas–Fermi law of Eq. (14). The dotte… view at source ↗
read the original abstract

As new observational technologies reveal increasingly detailed properties of neutron star surfaces, the demand for accessible and extensible theoretical modeling tools continues to grow. We present UnReal-B, a real-space Density Functional Theory solver for one-dimensional chains of matter in extreme magnetic fields $B \approx 10^{12} - 10^{15},\mathrm{G}$. By employing the adiabatic approximation, UnReal-B, provides a streamlined numerical framework for calculating the electronic structure of strongly magnetized condensed matter. The solver is benchmarked against published results for several astrophysically relevant elements, demonstrating excellent agreement while maintaining a comparatively simple and transparent implementation. Released as open-source software, UnReal-B facilitates reproducible and community-driven investigations of neutron-star surface matter and provides a foundation for future developments motivated by emerging observational constraints.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The manuscript introduces UnReal-B, an open-source real-space Density Functional Theory solver for one-dimensional chains of matter in extreme magnetic fields (B ≈ 10^12–10^15 G). It employs the adiabatic approximation to compute electronic structure, benchmarks the code against published results for several astrophysically relevant elements, and claims excellent agreement while emphasizing a simple and transparent implementation.

Significance. If the benchmarks and implementation details hold, the open-source release of a streamlined DFT tool for neutron-star surface modeling would be a useful contribution to the field, supporting reproducible investigations motivated by new observational data.

major comments (2)
  1. [Methods and Results sections] The central claim of accuracy across B ≈ 10^12–10^15 G rests on the adiabatic approximation, yet no section quantifies its error (e.g., via residuals against full 3D Landau-level calculations or higher-order transverse corrections) or identifies the B threshold where the approximation breaks down. This omission is load-bearing because the abstract asserts 'excellent agreement' without supporting sensitivity tests at the upper end of the range.
  2. [Results section] Benchmark comparisons are described only qualitatively ('excellent agreement'); the manuscript provides neither tabulated residuals, convergence metrics, nor explicit sensitivity tests at B = 10^15 G, which prevents independent assessment of whether the reported agreement actually validates the approximation over the full claimed interval.
minor comments (1)
  1. [Abstract] The abstract would be strengthened by naming the specific elements used in the benchmarks and the quantitative metrics (e.g., energy differences or density residuals) that define 'excellent agreement'.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for the constructive comments, which have helped us improve the clarity and rigor of the manuscript. We address each major point below and have made revisions to strengthen the presentation of the adiabatic approximation and the benchmark results.

read point-by-point responses

  1. Referee: [Methods and Results sections] The central claim of accuracy across B ≈ 10^12–10^15 G rests on the adiabatic approximation, yet no section quantifies its error (e.g., via residuals against full 3D Landau-level calculations or higher-order transverse corrections) or identifies the B threshold where the approximation breaks down. This omission is load-bearing because the abstract asserts 'excellent agreement' without supporting sensitivity tests at the upper end of the range.

    Authors: We agree that an explicit discussion of the adiabatic approximation's validity range would strengthen the manuscript. The approximation is standard in the literature for the B range considered (e.g., as used in prior works on magnetized neutron-star crusts), but we acknowledge the value of identifying its limitations. In the revised manuscript we have added a dedicated paragraph in the Methods section that cites existing comparisons between adiabatic and full 3D treatments, states the expected breakdown threshold near or above 10^15 G for the densities of interest, and includes new sensitivity tests at B = 10^15 G. Direct residuals against independent full 3D Landau-level calculations are not feasible within the scope of this work, which focuses on releasing an accessible real-space 1D solver; such comparisons would require a separate computational framework. revision: partial

  2. Referee: [Results section] Benchmark comparisons are described only qualitatively ('excellent agreement'); the manuscript provides neither tabulated residuals, convergence metrics, nor explicit sensitivity tests at B = 10^15 G, which prevents independent assessment of whether the reported agreement actually validates the approximation over the full claimed interval.

    Authors: We concur that quantitative presentation of the benchmarks is necessary. The revised Results section now includes tabulated residuals (absolute and relative differences) for each element and field strength against the published reference values, convergence metrics with respect to spatial grid spacing and simulation-box size, and explicit numerical results at B = 10^15 G demonstrating that agreement remains at the level reported for lower fields. revision: yes

standing simulated objections not resolved
  • Direct quantitative residuals against independent full 3D Landau-level calculations, as this would require implementing and validating an entirely separate 3D solver outside the scope and design goals of the UnReal-B 1D real-space framework.

Circularity Check

0 steps flagged

No significant circularity

full rationale

The paper describes a numerical real-space DFT solver (UnReal-B) that implements the adiabatic approximation for electronic structure in extreme magnetic fields and validates it by direct benchmarking against independently published results for astrophysically relevant elements. No derivation chain reduces a claimed prediction or result to its own fitted inputs, self-citations, or definitional tautologies; the central output is a computational framework whose accuracy is assessed externally rather than by internal consistency alone. The adiabatic approximation is adopted as a modeling choice whose validity range is asserted but not derived from the solver itself.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the validity of the adiabatic approximation for the stated field range; no free parameters or new physical entities are mentioned.

axioms (1)
  • domain assumption The adiabatic approximation is valid for magnetic fields B ≈ 10^12 - 10^15 G
    Explicitly employed by the solver as stated in the abstract.

pith-pipeline@v0.9.1-grok · 5671 in / 1048 out tokens · 24289 ms · 2026-06-27T11:48:39.304594+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

42 extracted references · 2 canonical work pages

  1. [1]

    Density-functional-theory calculations of matter in strong magnetic fields. II. Infinite chains and condensed matter , author=. Physical Review A—Atomic, Molecular, and Optical Physics , volume=. 2006 , publisher=

    2006

  2. [2]

    Density-functional-theory calculations of matter in strong magnetic fields. I. Atoms and molecules , author=. Physical Review A—Atomic, Molecular, and Optical Physics , volume=. 2006 , publisher=

    2006

  3. [3]

    Matter and radiation in the strong magnetic fields of neutron stars , author=. Ph. D. Thesis , year=

  4. [4]

    The Astrophysical Journal , volume=

    Structure of solid iron in superstrong magnetic fields , author=. The Astrophysical Journal , volume=. 1974 , publisher=

    1974

  5. [5]

    The Astrophysical Journal , volume=

    Theory of pulsars: polar gaps, sparks, and coherent microwave radiation , author=. The Astrophysical Journal , volume=. 1975 , publisher=

    1975

  6. [6]

    Part II: Including Magnetic Fields , author=

    The Structure of the Density-Potential Mapping. Part II: Including Magnetic Fields , author=. ACS Physical Chemistry Au , volume=. 2023 , publisher=

    2023

  7. [7]

    Universe , volume=

    Dense Matter and Compact Stars in Strong Magnetic Fields , author=. Universe , volume=. 2026 , publisher=. doi:10.3390/universe12050122 , note=

    work page doi:10.3390/universe12050122 2026

  8. [8]

    Zane, Silvia and Taverna, Roberto and González–Caniulef, Denis and Muleri, Fabio and Turolla, Roberto and Heyl, Jeremy and Uchiyama, Keisuke and Ng, Mason and Tamagawa, Toru and Caiazzo, Ilaria and Di Lalla, Niccolò and Marshall, Herman L. and Bachetti, Matteo and La Monaca, Fabio and Gau, Ephraim and Di Marco, Alessandro and Baldini, Luca and Negro, Mich...

    work page doi:10.3847/2041-8213/acb703 2023

  9. [9]

    Physical Review Letters , volume=

    Matter in Superstrong Magnetic Fields: The Surface of a Neutron Star , author=. Physical Review Letters , volume=. 1971 , publisher=

    1971

  10. [10]

    Reviews of Modern Physics , volume=

    Matter in strong magnetic fields , author=. Reviews of Modern Physics , volume=. 2001 , publisher=

    2001

  11. [11]

    The Astrophysical Journal , volume=

    Variational calculation of ground-state energy of iron atoms and condensed matter in strong magnetic fields , author=. The Astrophysical Journal , volume=. 1977 , publisher=

    1977

  12. [12]

    Monthly Notices of the Royal Astronomical Society , volume=

    Density functional calculations of the ground-state energies of atoms and infinite linear molecules in very strong magnetic fields , author=. Monthly Notices of the Royal Astronomical Society , volume=. 1985 , publisher=

    1985

  13. [13]

    Physical Review A , volume=

    Hartree-Fock calculations of atoms and molecular chains in strong magnetic fields , author=. Physical Review A , volume=. 1986 , publisher=

    1986

  14. [14]

    Monthly Notices of the Royal Astronomical Society , volume=

    Condensed surfaces of magnetic neutron stars, thermal surface emission, and particle acceleration above pulsar polar caps , author=. Monthly Notices of the Royal Astronomical Society , volume=. 2007 , publisher=

    2007

  15. [15]

    Computer Physics Communications , volume=

    Multi-electron systems in strong magnetic fields I: The 2D Landau--Hartree--Fock--Roothaan method , author=. Computer Physics Communications , volume=. 2014 , publisher=

    2014

  16. [16]

    Physical Review Letters , volume=

    Density-functional theory in strong magnetic fields , author=. Physical Review Letters , volume=. 1987 , publisher=

    1987

  17. [17]

    Physics-Uspekhi , volume=

    Atmospheres and radiating surfaces of neutron stars , author=. Physics-Uspekhi , volume=. 2014 , publisher=

    2014

  18. [18]

    A review , author=

    Magnetars: the physics behind observations. A review , author=. Reports on Progress in Physics , volume=. 2015 , publisher=

    2015

  19. [19]

    and Watts, Anna L

    Riley, Thomas E. and Watts, Anna L. and Bogdanov, Slavko and others , journal=. A. 2019 , publisher=

    2019

  20. [20]

    Soviet Physics JETP , volume=

    Heavy atom in an ultrastrong magnetic field , author=. Soviet Physics JETP , volume=. 1970 , note=

    1970

  21. [21]

    Physical Review Letters , volume=

    Statistical Model of Atoms in Intense Magnetic Fields , author=. Physical Review Letters , volume=. 1971 , publisher=

    1971

  22. [22]

    Semiclassical regions , author=

    Asymptotics of heavy atoms in high magnetic fields: II. Semiclassical regions , author=. Communications in Mathematical Physics , volume=. 1994 , publisher=

    1994

  23. [23]

    The Astrophysical Journal , volume=

    Surface structure of neutron stars with high magnetic fields , author=. The Astrophysical Journal , volume=. 1989 , publisher=

    1989

  24. [24]

    Annual Review of Astronomy and Astrophysics , volume=

    Magnetars , author=. Annual Review of Astronomy and Astrophysics , volume=. 2017 , doi=

    2017

  25. [25]

    Space Science Reviews , volume=

    Physics of Neutron Star Crusts , author=. Space Science Reviews , volume=. 2015 , doi=

    2015

  26. [26]

    The Astrophysical Journal Letters , volume=

    PSR J0740+6620 Mass and Radius from NICER Data and Implications for the Properties of Neutron Star Matter , author=. The Astrophysical Journal Letters , volume=. 2021 , doi=

    2021

  27. [27]

    The Astrophysical Journal Letters , volume=

    A NICER View of the Massive Pulsar PSR J0740+6620 Informed by Radio Timing and XMM-Newton Spectroscopy , author=. The Astrophysical Journal Letters , volume=. 2021 , doi=

    2021

  28. [28]

    The Astrophysical Journal Letters , volume=

    The Radius of PSR J0030+0451 from NICER and XMM-Newton Observations , author=. The Astrophysical Journal Letters , volume=

  29. [29]

    The Astrophysical Journal Letters , volume=

    A Joint Analysis of NICER and XMM-Newton Observations of PSR J0030+0451 , author=. The Astrophysical Journal Letters , volume=

  30. [30]

    Science , volume=

    Polarized X-rays from a Magnetar , author=. Science , volume=. 2022 , doi=

    2022

  31. [31]

    Physical Review Letters , volume=

    GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral , author=. Physical Review Letters , volume=. 2017 , doi=

    2017

  32. [32]

    The Astrophysical Journal Letters , volume=

    Multi-messenger Observations of a Binary Neutron Star Merger , author=. The Astrophysical Journal Letters , volume=. 2017 , doi=

    2017

  33. [33]

    Nature , volume=

    A Variable Absorption Feature in the X-ray Spectrum of a Magnetar , author=. Nature , volume=. 2013 , doi=

    2013

  34. [34]

    Monthly Notices of the Royal Astronomical Society , volume=

    Unifying the Observational Diversity of Isolated Neutron Stars via Magneto-thermal Evolution Models , author=. Monthly Notices of the Royal Astronomical Society , volume=. 2013 , doi=

    2013

  35. [35]

    Monthly Notices of the Royal Astronomical Society , volume=

    Strong Crustal Magnetic Fields in Neutron Stars , author=. Monthly Notices of the Royal Astronomical Society , volume=. 2013 , doi=

    2013

  36. [36]

    Physical Review A , volume=

    Ground-state Properties of Matter in Strong Magnetic Fields , author=. Physical Review A , volume=. 1987 , doi=

    1987

  37. [37]

    The Astrophysical Journal , volume=

    Hydrogen Phases on the Surface of a Strongly Magnetized Neutron Star , author=. The Astrophysical Journal , volume=. 1997 , doi=

    1997

  38. [38]

    Density-functional-theory Calculations of Matter in Strong Magnetic Fields. I. Atoms and Molecules , author=. Physical Review A , volume=. 2006 , doi=

    2006

  39. [39]

    Density-functional-theory Calculations of Matter in Strong Magnetic Fields. II. Infinite Chains and Condensed Matter , author=. Physical Review A , volume=. 2006 , doi=

    2006

  40. [40]

    Universe , volume=

    Decoding the Nature of Coherent Radio Emission in Pulsars I: Observational Constraints , author=. Universe , volume=. 2024 , doi=

    2024

  41. [41]

    Advances in Space Research , volume=

    Cohesive Property of Magnetized Neutron Star Surfaces: Computations and Implications , author=. Advances in Space Research , volume=. 2007 , doi=

    2007

  42. [42]

    Astronomy and Astrophysics , volume=

    Partially Screened Gap Model of Inner Acceleration Region in Pulsars , author=. Astronomy and Astrophysics , volume=. 2003 , doi=

    2003