This project contains the code for a two-level quantum Bloch solver for arbitrary magnetic fields. Using the unitary integration method, this code can be easily extended to arbitrary potentials in the semiclassical limit Currently, the code base is designed for simple fields and simulations of the radio-frequency neutron spin flipper developed by Indiana University.
Thing to add in the future
- Methods for integrating orbital angular momentum
- Add ray-tracing option for 2D and 3D potentials
- Expand code to include infinite-dimensional algebras
- Enhance visualization code
Other algorithms/methods may also prove useful:
- Wei-Normann method (factorization)
- Higher-order Magnus method (enhancing convergence)
- Fer method (for seperable potentials)
- Wilcox method (simplier version of the Fer method)