Revisiting spin Hamiltonian parameters in a Kitaev material via Bayesian optimization of magnetization curves
Phys. Rev. B 114, 274420 – Published 12 November, 2026
DOI: https://doi.org/10.1103/s9sc-4m5y
Abstract
Determining the spin Hamiltonian of a magnetic compound is crucial for understanding its magnetic properties. A standard approach is to derive model parameters from ab initio calculations based on the crystal structure. However, the resulting Hamiltonian can depend sensitively on methodological details of the ab initio procedure. This issue is particularly evident in , a candidate Kitaev material. Here we present an alternative, data-driven approach to determine the spin Hamiltonian parameters of by Bayesian optimization of experimental magnetization curves along the - and -axis directions. We optimize five parameters, namely the Kitaev interaction , off-diagonal interactions and , the Heisenberg interaction , and the -axis -factor . The parameter set that minimizes the cost function is , where the exchange couplings are in millielectronvolts. We find that the cost function is insensitive to the absolute value of the Kitaev coupling . Thus, the magnetization data alone do not determine its energy scale. The cost function also depends only weakly on and , while the optimization favors a large positive . By computing the static spin structure factor, magnetic susceptibility, and specific heat, we show that these quantities favor the large- scenario over the small- scenario and that the parameter set that minimizes the cost function yields good agreement with the experimental structure factor and susceptibility, and qualitative agreement with the specific heat. The combination of Bayesian optimization and accurate low-energy solvers provides an effective approach for determining parameters of spin Hamiltonians. This methodology opens a systematic route to determining spin Hamiltonians in quantum magnets from experimental data.