Radiation-induced embrittlement of bainitic steels is the lifetime limiting factor of reactor pressure vessels in existing nuclear light water reactors. The primary mechanism of embrittlement is the obstruction of dislocation motion produced by nanometric defect structures that develop in the bulk of the material due to irradiation. In view of improving the predictive capability of existing models it is necessary to understand better the mechanisms leading to the formation of these defects, amongst which the so-called “late blooming phases”. In this work we study the stability of the latter by means of density functional theory (DFT) calculations and Monte Carlo simulations based on a here developed quaternary FeCuNiMn interatomic potential. The potential is based on extensive DFT and experimental data. The reference DFT data on solute–solute interaction reveal that, while Mn–Ni pairs and triplets are unstable, larger clusters are kept together by attractive binding energy. The NiMnCu synergy is found to increase the temperature range of stability of solute atom precipitates in Fe significantly as compared to binary FeNi and FeMn alloys. This allows for thermodynamically stable phases close to reactor temperature, the range of stability being, however, very sensitive to composition.