The minimal quantum model for spontaneous emission is given by the interaction Hamiltonian
$$ H = \int d\omega [g(\omega)a(\omega)\sigma^+ + g^*(\omega)a^\dagger(\omega)\sigma^-]$$
where $a(\omega), a^\dagger(\omega)$ are bosonic modes associated with the light field and $\sigma^\pm$ are the ladder operators of a two-level system. The coupling function $g(\omega)$ depends on the electromagnetic environment and can e.g. be evaluated for free space or for a cavity.
The sponteneous emission problem is then specifically given by the initial state where the atom is fully excited and the light field is in the vacuum state (or a thermal state at a given temperature, but I will only consider the vacuum case for simplicity).
With regards to solution methods, there are multiple levels of complexity that one can consider. The standard approach, which is probably covered in every modern quantum optics book, is Wigner-Weisskopf theory. The latter makes some approximations such as the Markov approximation, which make it applicable at weak light-matter coupling.
For strong light-matter coupling, the standard approach is to use an alternative starting Hamiltonian, namely the Jaynes-Cummings model as suggested by @NorbertSchuch in the comments. The reason is that strong coupling is often achieved by coupling the atom to a resonator where a single mode dominates.
As an advanced topic, I would like to note that the spontaneous emission problem, specifically, can also be described at strong coupling using the above Hamiltonian. There are a bunch of resources on this topic, I think one of the Cohen-Tannoudji quantum optics textbooks also contains a chapter. My personal favourite is Krimer et al., PRA 89, 033820 (2014). It essentially presents a generalization of the Wigner-Weisskopf method, where the Markov approximation is eliminated by solving certain integrals exactly using Laplace transforms. With regards to the comments, the paper also shows the transition from the Wigner-Weisskopf type emission to the Jaynes-Cummings type dynamics as you turn up the coupling by making the electromagnetic environment more confined.
Note that the central reason why this model is so useful and can be solved so accurately is the rotating wave approximation. It results in the Hamiltonian being excitation number conserving, such that only the states $|1_\mathrm{atom}\rangle|0_\omega\rangle$ and $|0_\mathrm{atom}\rangle|1_\omega\rangle$ contribute.
Even this treatment breaks down eventually. In the ultra-strong coupling regime, the minimal model given above does not work anymore, since counter-rotating terms and stuff like gauge invariance start to matter. There is a lot of literature on this topic, see e.g. What are the "strong", "ultrastrong" and "deep strong" coupling regimes of the Rabi model? and references therein.
