  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.
