Can you just use powsimp instead to put powers back together?

```
>>> collect( sol, exp(x/3))
sqrt(exp(x)) + 1
>>> powsimp(_)
exp(x/2) + 1
```