src.Atomic.SEsolver.setMatrixR

src.Atomic.SEsolver.setMatrixR(_Rmat, _Rji_spon, _Rji_stim, _Rij, _idxI, _idxJ)

Compute the radiative rate matrix.

Parameters

• _Rmat (numpy.2darray of np.double) – radiative rate matrix, a 2D Array to store computed results, [\(s^{-1}\)]

• _Rji_spon (numpy.1darray of np.double) – spontaneous radiative transition rate, [\(s^{-1}\)]

• _Rji_stim (numpy.1darray of np.double) – spontaneous radiative transition rate, [\(s^{-1}\)]

• _Rij (numpy.1darray of np.double) – upward radiative transition rate, [\(s^{-1}\)]

• _idxI (numpy.1darray of np.uint16) – level index of lower level i, [-]

• _idxJ (numpy.1darray of np.uint16) – level index of upper level j, [-]

Notes

Refer to 1 Equation(9.80).

\[\sum_{j \neq i} n_j (R_{ji}+C_{ji}) - n_i \sum_{j \neq i}(R_{ij}+C_{ij}) = 0\]

References

1

Ivan Hubeny, Dimitri Mihalas, “Theory of Stellar Atmosphere: An Introduction to Astrophysical Non-equilibrium Quantitative Spectroscopic Analysis”, Princeton University Press, pp. 282, 2015.