QED theory of multiphoton cascade transitions in atoms and ions is developed. This theory allows for the accurate description of the process important for astrophysical studies of the cosmological hydrogen recombination. In particular the $ 3s\rightarrow1s+2\gamma $, $ 4s\rightarrow1s+2\gamma $ and $ 3p\rightarrow1s+3\gamma $ processes are considered and some controversies existing in the literature are resolved.

A simple model based on QED is presented for the estimation of contribution of the excited level few-photon decays to the radiation escape from the matter in the epoch of the cosmological hydrogen recombination. It is shown that apart from the widely studied two-photon decays, some specific 3-photon decays can contribute on the level of 0.1% accuracy, required by the recent astrophysical observations.

The relativistic and nonrelativistic approaches for the calculations of the two-photon decay rates of highly excited states in hydrogen are compared. The dependence on the principal quantum number (n) of the ns, nd, and np initial states is investigated up to n = 100 for the nonresonant emissions. For the ns states together with the main E1E1 channel the contributions of higher multipoles (M1M1, E2E2, E1M2) are considered. For the np states the E1M1 and E1E2 channels are evaluated. Moreover, the simple analytical formula for the E1M1 decay is derived in the nonrelativistic limit.

A rigorous QED theory of the multiphoton decay of excited states in hydrogen atom is presented. The "two-photon" approximation is formulated which is limited by the one-photon and two-photon transitions including cascades transitions with two-photon links. This may be helpful for the strict description of the recombination process in hydrogen atom and, in principle, for the history of the hydrogen recombination in the early Universe.

The problem of the evaluation of the two-photon decay width of excited states in hydrogen is considered. Two different approaches to the evaluation of the width including cascades channels are employed: the summation of the transition probabilities for various decay channels and the evaluation of the imaginary part of the Lamb shift. As application the decay channels for the 3s level of the hydrogen atom are evaluated, including the cascade transition probability 3s-2p-1s as well as "pure" two-photon decay probability and the interference between both channels. An important role should be assigned to the "pure" two-photon probability in astrophysics context, since processes of this kind provide a possibility for the decoupling of radiation and matter in the early Universe. We demonstrate the ambiguity of separation of the "pure" two-photon contribution and criticize the existing attempts for such a separation.