ElliotB256

There is a recent repo for implementing this algorithmn, wouldn't hurt to try it and see how it performs. I expect it may need a bit of profiling https://docs.rs/nbody_barnes_hut/0.1.0/nbody_barnes_hut/

Good overview in [Leyuan Wu, Geophys. J. Int.(2018)215,1865–1886](https://watermark.silverchair.com/ggy379.pdf?token=AQECAHi208BE49Ooan9kkhW_Ercy7Dm3ZL_9Cf3qfKAc485ysgAAAsEwggK9BgkqhkiG9w0BBwagggKuMIICqgIBADCCAqMGCSqGSIb3DQEHATAeBglghkgBZQMEAS4wEQQMoZSD6UJDpP6O2_DJAgEQgIICdDxn4yXn5lOPKPVJhfBYRZmbJgIDzJGQFPlmthn_PvvNA1JSj3AjrcEdodzevs0dbPJil_kyga-P_NXwoCfI2YjmGe1UKVWRP6O8eQfPOxB_vghW9tHveHyManHzvhcBT6t1tle6zVH6XVAded5OHEOLHdIPAyX12qrFFSZ3MmDoKxc6Yhk2hq6akC8gPwIuPFntP8PxbwF_Abczr49RodsVCtjGiurfHPrWBFuJeJxkQX_5UdixW-S09MXz-BEV7vLY-H85ZMX0Y6jtQ-jSnhHu7BX5hUBofOZBNF1dVql43OYr-PVIwU5KKt37PjOrIqxjAaeGNyJAn98P7x7yjFwBtR0ocSvsjAv2Cp3taXXYZ34HmScEkcz7XIZQ8SqhMWnvwGVw-kNNrofNyIStfeqmKXDGdc_4TWNEt--rBdY3NP-rZiCA6fSGmcZGlCtjhdJt8sKcGZ7hy8awtDI5sVcKVNf3jPBiMbLO32p4O8sKhNWgnIe4oG9cDxeNCWNu02DEDTClr-c-dbYU4rNx3WtfzloWStOoptfTP6NQmenKyfo4TyDGP6IafeoRifiJd5LATGbcL4gpB4BF3inp7XNh0xuYR-B_8qHv1yTGQYUgR5wgb0rfJyPni4wbtZUxGdSzQc6_v4kPcXvmV5jXMVlbx_oNFtSqkyJps82cpE_wZgiB4gVq3W2WwZuzEbGU0i5lpP3UYLCdF43vetfYDU1aRKcGvdc8IUG0CoY5CAdmJ6JpgyGAniEwuwB7Mr7-80E8RhzLg7kivONYiTQ9vpaN43RSrpa-xhMHVYqk7YFjDVBfhoxwgH-1h5N4aahcM6eFUZo)

That's a good idea! I don't think there is at the moment, but it would be quite simple to ingest a file and populate the simulation from that. A straightforward...

How did you get on with the iterator approach?

That's great to hear! If you want to contribute the code, consider making a pull request and then we can merge it into atomECS

# Work at Oxford * **Sr cold atom source:** simulations to guide design of a source of ultracold strontium atoms for use in the AION project, using the 461nm cooling...

I think the doc is incorrect in this case. The 'polarisation' is either +1 or -1 depending on whether the beam is sigma_+ or sigma_- polarised. The `scatter1,2,3` components have...

> but it doesn't seem to be the most satisfactory solution because the user must ensure that the polarizatino is orthogonal to the wavector you could define private fields with...

I think we probably need to distinguish between two regimes: * far-detuned/conversative dipole trap potentials. In this case the beams should interfere and the combined intensity should dictate the trapping...

Side note - Also worth bearing in mind the current approach, which is intensity based and does not have interference, has a 'one-to-N' concept in the sampler, because each atom...