Linear response properties of solvated systems: a computational study

We present a computational study of static and dynamic linear polarizabilities in solution. We use different theoretical approaches to describe solvent effects, ranging from quantum mechanics/molecular mechanics (QM/MM) to quantum embedding approaches. In particular, we consider non-polarizable and polarizable QM/MM methods, the latter based on the fluctuating charge (FQ) force field. In addition, we use a quantum embedding method defined in the context of multilevel Hartree-Fock (MLHF), where the system is divided into active and inactive regions, and combine it with a third layer described by means of the FQ model. The multiscale approaches are then used as reference wave functions for equation-of-motion coupled cluster (EOM-CC) response properties, allowing for the account of electron correlation. The developed models are applied to the calculation of linear response properties of two organic moieties-namely, para-nitroaniline and benzonitrile-in non-aqueous solvents-1,4-dioxane, acetonitrile, and tetrahydrofuran. The computed polarizabilities are then discussed in terms of the physico-chemical solute-solvent interactions described by each method (electrostatic, polarization and Pauli repulsion), and finally compared with the available experimental references.