Polarizabilities fixed-charge method

Moments and polarizabilities can also be obtained by the fixed-charge method [76]. This technique allows for the single-step incorporation of the nonuniform electric field contributions due to gradients and higher order field derivatives. One or more charges are placed around the molecule in regions where the molecular wavefunctions are negligible. It is important that the basis set used for the field-free molecule be the same as that used in the presence of the field and that the molecule basis be adequate to describe any... 

Nevertheless, these methods are mostly applied with fixed charges (even if these are chosen in a sophisticated way) and with pairwise additivity approximation as well as with the neglect of nuclear quantum effects. Suggestions for polarizable models appeared in literature mainly for water [23], The quality of potential parameterization varies from system to system and from quantity to quantity, raising the question of transferability. Spontaneous events like reactions cannot appear in simulations unless the event is included in the parameterization. Despite these problems, it is possible to reproduce important quantities as structural, thermodynamic and transport properties with traditional MD (MC) mainly due to the condition of the nanosecond time scale and the large system size in which the simulation takes place [24],... 

Several approaches for calculating excited states in protein environments were proposed to improve the ordinary QM/MM description. The effect of polarization was included as a classical force field [27], and the excitation energy calculated for bacteriorhodopsin (bR) was 0.34 eV less than that from a fixed-charge non-polarizable QM/MM method [27]. Later, a triple-layer QM1/QM2/MM approach was proposed, and DFT(PBEO) calculations were performed for the QM2 layer, which consisted of the amino acids 4 A from the retinal PSB [28]. The calculated excitation energy of bR was only 0.08 eV smaller than that obtained using the ordinary QM/MM method [28]. In another study, an empirical polarization model combined with the QM/MM calculation produced a red shift of 0.14-0.17 eV [29]. However, these pioneering studies neglected the CT effects between the retinal and the protein environments. 

We developed expressions for EOM-CCSD-PCM in the SS formalism, and a series of approximations to reduce the considerable computational cost of this approach. Christiansen and Mikkelsen originally developed the LR formalism for the CCSD wave function in solution for a simple continuum solvation model. Later, they extended it to their flavour of explicit polarizable solvation model. Cammi has also presented several interesting developments in this research area, including a rederivation of the LR-CCSD expressions for PCM, and we presented the first implementation of the method. Other examples of CC methods combined with (non-)polarizable solvation models [e.g., fixed point charges) are also available in the literature. ... 

---