- [ ] add a link to reference publication - [ ] check constant `k2` used in PhD thesis of Marta Korcyl https://arxiv.org/pdf/1410.5250.pdf  Current implementation uses a single value of `k2` for whole energy range.