Many-body perturbation theory calculations using the GW approximation starting from DFT wavefunctions (GW@DFT) are thus appealing due to their fully ab initio approach and well-known accuracy in predicting oneparticle properties, including ionization potential (IP), electron affinity (EA), and energy gaps in a wide variety of systems. Calculations using the GW approximation for periodic materials have typically been applied towards the study of bulk solids or molecules/surfaces in vacuum calculations. This is primarily due to the complexity in sampling solvent molecules configurations in explicit solvation approaches, and the lack of a complete theoretical formulation and computational tools for GW calculations with implicit solvation approaches.