I can't understand how authors of the article solve the system of rate equations with Runge-Kutta method.(adsbygoogle = window.adsbygoogle || []).push({});

Pulsed laser light is transmited through nonlinear absorbing sample. So, here is the system of rate equations which governs the populations of the levels in the three-level system.

[itex] \frac{dN_{0}}{dt} = -\frac{\sigma_{0} I N_{0}}{h \omega} - \frac {\beta I^{2} } {2 h \omega} - \frac{\gamma I^{3} }{3 h \omega} + \frac{N_{1}}{\tau_{s1}} [/itex]

[itex] \frac{dN_{1}}{dt} = \frac{\sigma_{0} I N_{0}}{h \omega} - \frac{\sigma_{1} I N_{1}}{h \omega} - \frac{N_{1}}{\tau_{s1}} +\frac{N_{2}}{\tau_{sn}} [/itex]

[itex] \frac{dN_{2}}{dt} = \frac{\sigma_{1} I N_{1}}{h \omega} + \frac {\beta I^{2} } {2 h \omega} + \frac{\gamma I^{3} }{3 h \omega} - \frac{N_{2}}{\tau_{sn}} [/itex]

Intensity transmitted through the sample is given by

[itex] \frac{dI}{dz} = -\sigma_{0} I N_{0} - \sigma_{1} I N_{1} - \beta I^{2} - \gamma I^{3}[/itex]

where [itex] I(r,z,t) = I_{00} \left( \frac{w_{0}^{2}}{w^{2}(z)} \right) \ \exp \left( - \frac{t^{2}}{\tau_{p}^{2}} \right) \ \exp \left( - \frac{2r^{2}}{w^{2}(z)} \right)[/itex]

Here [itex] \sigma_{i} [/itex] is the cross-section of the corresponding state, [itex] \beta[/itex] and [itex] \gamma [/itex] - two- and three-photon absoption coefficients, [itex] w(z) [/itex] - radius of the laser (assumed to be Gaussian) beam, [itex]N_{i}(t)[/itex] - population of the i-th state, [itex] \tau_{i}[/itex] - life-time of the corresponding state.

Here we come to the point where my incomprehension is rising. It is stated in the article that the differential equations were first decoupled and then integrated over time, length, and along the radial direction before solving them numerically with Runge-Kutta fourth order method. Limits of integration of [itex]r, t, z[/itex] are [itex](0,\infty), (-\infty, \infty)[/itex] and [itex](0,L)[/itex],respectively, where L is the sample length.

So, i wonder how i should do this "de-coupling"?

I thought that the RK method is aplicable for the systems where the derivatives are taken with respect to the same variable. But in this exaple the state population is the function of time, and in the equation for the intensity there is the coordinate derivative. May be i'm wrong?

**Physics Forums | Science Articles, Homework Help, Discussion**

Join Physics Forums Today!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# Solving rate equations with RK method. Don't understand how.

Loading...

**Physics Forums | Science Articles, Homework Help, Discussion**