@scopatz I just noticed a bug here - could you throw in this line at the top of the NewDecay function:

```c++
  atom(); // force evaluation of atom-composition if not calculated already
```