In this work we introduce a novel subsystem-based electronic structure embedding method that combines the projection-based block-orthogonalized Manby-Miller embedding (BOMME) with the density-based Frozen Density Embedding (FDE) methods. Our approach is effective for systems in which the building blocks interact at varying strengths while still maintaining a lower computational cost compared to a quantum simulation of the entire system. To evaluate the performance of our method, we assess its ability to reproduce the X-ray absorption spectra (XAS) of chloride and fluoride anions in aqueous solutions (based on a 50-water droplet model) via real-time time-dependent density functional theory (rt-TDDFT) calculations. We employ an ensemble approach to compute XAS for the K- and L-edges, utilizing multiple snapshots of configuration space obtained from classical molecular dynamics simulations with a polarizable force field. Configurational averaging influences both the broadening of spectral features and their intensities, with contributions to the final intensities originating from different geometry configurations. We found that embedding models that are too approximate for halide-water specific interactions, as in the case of FDE, fail to reproduce the experimental spectrum for chloride. Meanwhile, BOMME tends to overestimate intensities, particularly for higher energy features because of finite-size effects. Combining FDE for the second solvation shell and retaining BOMME for the first solvation shell mitigates this effect, resulting in an overall improved agreement within the energy range of the experimental spectrum. Additionally, we compute the transition densities of the relevant transitions, confirming that these transitions occur within the halide systems. Thus, our real-time QM/QM/QM embedding method proves to be a promising approach for modeling XAS of solvated systems.
