# Zero Order Extinction Reaction

**URL:** <https://discourse.julialang.org/t/zero-order-extinction-reaction/108749>\
**Category:** Chemistry\
**Tags:** question, catalyst\
**Created:** [January 12, 2024, 10:53pm UTC](https://discourse.julialang.org/t/zero-order-extinction-reaction/108749 "2024-01-12T22:53:58Z")\
**Posts on this page:** 1\
**Showing post:** 2

<div class="post-metadata">

**Author:** ![ChrisRackauckas](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/chrisrackauckas/32/77_2.png) [@ChrisRackauckas](https://discourse.julialang.org/u/ChrisRackauckas)\
**Post date:** [January 13, 2024, 11:58pm UTC](https://discourse.julialang.org/t/zero-order-extinction-reaction/108749/2 "2024-01-13T23:58:34Z")

</div>

Yes this is something that will be very numerically sensitive and using some kind of callback is generally a good idea. The `isoutofdomain` solver option and `PositiveDomain` callback in the DiffEqCallbacks.jl library were designed for this kind of application. The one that you have there introduces error so I would recommend the two standard approaches over that.

---

_[View the full topic](https://discourse.julialang.org/t/zero-order-extinction-reaction/108749)._