New developments in the symmetry-adapted algorithm of the Polarizable Continuum Model

We present recent developments in the symmetry implementation of the Polarizable Continuum Model (PCM). The structure of the matrix, which defines the PCM solvent response, is examined, and we demonstrate how this matrix can be transformed to a block diagonal form where each block belongs to different irreducible representations of the molecular point group. This development is especially important at the Multi-configurational Self-Consistent Field (MCSCF) level where symmetry is needed to avoid problems with symmetry breaking in the wave function and facilitate the optimization of electronic excited states. Moreover, although only the totally symmetric part of the solvent interaction is needed for energy calculations, in response or perturbation theory calculations of molecular properties, other irreps play an important role and the classification of solvent interaction terms by irrep is, therefore, desirable. In addition, the use of symmetry reduces the computational cost. The implementation presented here is illustrated with a series of calculations of absorption and emission processes in solution on the diazines pyrazine, pyrimidine, and pyridazine. These calculations allow us to illustrate both formal aspects of the implementation such as the choice of active spaces based on orbital symmetry as well as practical aspects such as the speed-up of the calculation.