spinoct

Optimal control paths and energy-efficient switching pulses for classical spin dynamics (Landau-Lifshitz-Gilbert).

spinoct computes the control (an applied magnetic field, an electric current, or both) that drives a magnetic moment from one state to another in a given time for the least dissipated energy. It is the reusable engine behind Espira, and it is deliberately independent of any material database, so it works on any spin Hamiltonian.

Why this exists

The optimal control of magnetization switching has a small, rigorous literature (Kwiatkowski, Badarneh, Berkov and Bessarab, Phys. Rev. Lett. 126, 177206 (2021), and the papers that follow it), but no open implementation. The analytic solutions live only as equations in journal articles, and the group's own code is not public. spinoct is that implementation, with the closed-form results built in as positive controls that every numerical solver is validated against before it is trusted.

Install

pip install spinoct              # core: numpy + scipy
pip install "spinoct[torch]"     # add the batched GPU solvers

The dimensional contract

Read spinoct.units first. The same physical quantity is written four different ways across this literature, and the central quantity of the package, the switching cost Phi = int |b|^2 dt, is in tesla-squared-seconds, not joules. It becomes an energy only through an explicit circuit model. The package refuses to hide that assumption: spinoct.units.CircuitModel is a required, described object, never a buried constant.

Status

Pre-1.0, under active development. The analytic uniaxial optimal control path, its closed-form pulse and cost, and the negative-parameter Jacobi elliptic machinery it needs are complete and validated. The numerical image-based solver, the GRAPE and CRAB constrained solvers, and the batched GPU lane are in progress.

License

MIT. See LICENSE.