We currently have an issue in how we handle ices. We do not supply a time dependent dust temperature to it. Doing so causes instabilities in the solver. Long term, the best way forward would be to get dust distribution-independent freeze out and evaporation rates from KIDA (Reboussin et al. 2014: https://kida.astrochem-tools.org/networks.html)
Overall, it does not make a big difference for problems where getting the temperature correct is key (e.g., typical star formation and ISM simulations).
We currently have an issue in how we handle ices. We do not supply a time dependent dust temperature to it. Doing so causes instabilities in the solver. Long term, the best way forward would be to get dust distribution-independent freeze out and evaporation rates from KIDA (Reboussin et al. 2014: https://kida.astrochem-tools.org/networks.html)
Overall, it does not make a big difference for problems where getting the temperature correct is key (e.g., typical star formation and ISM simulations).