Recent experimental studies of the spin-1 honeycomb antiferromagnet Na3Ni2BiO6 have revealed a pronounced one-third magnetization plateau under applied magnetic fields, highlighting the presence of strong magnetic frustration and anisotropy in this material. Such behavior has been attributed to substantial bond-dependent Kitaev interactions in combination with single-ion anisotropy, placing Na3Ni2BiO6 among honeycomb compounds of interest for unconventional magnetic phases. Motivated by these observations, we present a first-principles–based analysis of the magnetic interactions in Na3Ni2BiO6. By combining density-functional calculations with microscopic modeling, we extract the relevant exchange parameters and construct an effective spin model that quantitatively reproduces both the elastic neutron-scattering spectra and the magnetization curve. The model captures the experimentally observed zero-field zigzag magnetic order, and proposes a double-zigzag state at intermediate magnetic fields, realizing the one-third magnetization plateau in a simpler way than suggested in previous works. Crucially, we show that the one-third magnetization plateau does not require Kitaev interactions; instead, it arises from the interplay of strong out-of-plane single-ion anisotropy and competing ferromagnetic nearest-neighbor J1 and antiferromagnetic third-neighbor J3 Heisenberg couplings. These results establish a consistent microscopic description of Na3Ni2BiO6 and clarify the origin of its field-induced plateau phase.

Datasets obtained from classical Monte Carlo (Metropolis-Hastings) for magnetization plots and static structure factor data.