In this work, we propose a data-driven method to discover the latent space and learn the
corresponding latent dynamics for a collisional-radiative (CR) model in radiative plasma
simulations. The CR model, consisting of high-dimensional stiff ordinary differential
equations (ODEs), must be solved at each grid point in the configuration space, leading
to significant computational costs in plasma simulations. Our method employs a physicsassisted
autoencoder to extract a low-dimensional latent representation of the original CR
system. A flow map neural network is then used to learn the latent dynamics. Once trained,
the reduced surrogate model predicts the entire latent dynamics given only the initial condition
by iteratively applying the flow map. The radiative power loss is then reconstructed
using a decoder. Numerical experiments demonstrate that the proposed architecture can
accurately predict both the full-order CR dynamics and the radiative power loss rate.