Quantum Chemical Calculations of the First- and Second-Order Hyperpolarizabilities of Molecules in Solutions

A study of the effect of solute/solvent interactions on static and dynamic molecular hyperpolarizabilities (β and γ) of series of prototypical π-conjugated donor−acceptor chromophores is presented. The solvent effect was included via a recently proposed discrete quantum-mechanical Langevin dipoles/Monte Carlo method. The nonlinear optical properties (NLO) were computed using the INDO-like Hamiltonian based finite-field (FF) and sum-over-states (SOS) methods implemented in the GRINDOL code. The calculated β(γ) values are compared with experimental data determined in solution phase EFISH(THG) measurements and other published theoretical works. Generally, a reasonable agreement between calculated and experimental data was obtained. Moreover, our results indicate that the QM/LD/MC model gives a correct description of the solvent effect on the nonlinear optical response of molecules. It suggests that this level of theory can be used as an effective tool for investigation of NLO properties in condensed phases.