The MBD method efficiently captures many-body polarization and dispersion effects to infinite order in perturbation theory using coupled quantum Drude oscillators (QDOs) to model collective electron density fluctuations. Each oscillator, effectively representing the response of valence electrons of an atom, is defined by three parameters: mass m, frequency ω, and charge q. To map atoms onto oscillators, we employ the optimized QDO parametrization (vdW-OQDO), designed to accurately describe atomic response to electric fields and vdW binding curves of atomic dimers . This allows us to determine {m, ω, q} solely from two response properties of atoms: the static polarizability (α 0 ) and the dipolar dispersion coefficient (C 6 ) -both obtained from reliable free-atom reference data and semi-local electron density .