This research activity focuses on the nonlinear optical response properties of molecules, exploring the complex interactions between light and matter at the molecular level. The study aims to derive and implement quadratic response theory within the combined coupled cluster/dielectric continuum framework, providing a deeper understanding of the underlying mechanisms and their implications for various applications.