Correlated electronic-structure methods, such as those based on coupled-cluster (CC) theory, provide a robust platform for computing ground and excited states in a variety of closedand open-shell systems. The scope of molecules that can be treated by these high-level methods is expanding, thanks to advances in algorithmic strategies and computing hardware. In the past three decades, CC methods such as the equation-of-motion coupled-cluster (EOM-CC) and the coupled-cluster response theory (CC-RSP) have been vigorously developed to afford computations of an ever-growing list of molecular properties, greatly enhancing our capabilities for reliable first-principles modeling of spectroscopic experiments. These methods are now routinely employed for computing ground-and excited-state properties in linear and non-linear regimes. Examples of nonlinear properties include multiphoton absorption cross sections, inelastic scattering cross sections, first and second hyperpolarizabilities, excited state polarizabilities and multiphoton transition moments between excited states, multiphoton and magnetic circular dichroism strengths, g-tensors, to mention just a few. In addition, the highlevel CC methods provide the gold standard for benchmarking approximate and computationally less expensive methods, such as those based on density functional theory (DFT), which can treat much larger systems. Concomitant developments of embedding schemes afford robust CC condensed-phase modeling of molecular processes and properties in solutions, protein environments, and periodic systems. Further, advances in orbital concepts and wave-function analysis tools facilitate robust molecular-level understanding of the molecular processes within the framework of correlated methods.