A stochastic cumulant GW method is presented, allowing us to map the evolution of photoemission spectra, quasiparticle energies, lifetimes and emergence of collective excitations from molecules to bulk-like systems with up to thousands of valence electrons, including Si nanocrystals and nanoplatelet. The quasiparticle energies rise due to their coupling with collective shake-up (plasmon) excitations, and this coupling leads to significant spectral weight loss (up to 50% for the low energy states), shortening the lifetimes and shifting the spectral features to lower energy by as much as 0.6 eV. Such features are common to all the systems studied irrespective of their size and shape. For small and low dimensional systems the surface plasmon resonances affect the frequency of the collective excitation and position of the satellites.Recent developments in Green's function (GF) techniques have allowed for the description of charge excitations, i.e., quasiparticles (QPs) [1,2], in the bulk, over a wide range of QP energies. Band-edge excitations are well-described by the so called G 0 W 0 approximation [3-5], while at higher QP energies corrections are required to account for charge-density fluctuations and hole-plasmon coupling [5][6][7][8]. Photoemission experiments on solids reveal significant QP lifetime shortening and coupling to other collective excitations, manifested by satellite structures in the photoemission spectra [6,9,10]. The satellite structure and the QP lifetime shortening is often captured by the cumulant expansion (CE) ansatz to G 0 W 0 [6][7][8][10][11][12][13][14] .In confined systems, the QP spectrum near the bandedge is governed by the quantum confinement of electrons and holes. Higher energy, satellite-excitations are attributed to simultaneous ionization and excitation of the valence electrons ("shake-up" excitations) [13,[15][16][17][18]. Transition and differences between the satellite spectral features of molecules and nanostructures with "shake-up" signatures and bulk with collective plasmon resonances have been difficult to assess, as they require many-body treatment of systems with hundreds and thousands of electrons. In fact, the quantum confinement effect on the satellite transitions has received little attention if any.In this letter, we address this challenge by combining the well-known cumulant expansion (CE) ansatz [5, 7-9, 11, 19, 20] with the recent stochastic GW approach (sGW [21,22]), to obtain a nearly linear-scaling algorithm that reveals the changes of the QP spectra from a single molecule to covalently bonded nanocrystals (NCs) of unprecedented size (here up to 5288 valence electrons). The formalism is presented and assessed for the two size extremes (molecule and bulk), followed by the study of the effects of quantum confinement on the satellite structure in silicon NCs of different size and shape. In small NCs the satellite features are affected by the changes in the plasmon energy. For large NCs, we find observable quantum confinement effects on the satellite features b...