I am trying to calculate the electrostatic activation of an aromatic molecule containing H, C, N, S, and P. This is a cycloaddition reaction to an unsaturated structure. I am using the M06-2x/6-31G+(d,p) level of theory. How can I do it? SMD solvent system is used. Can anyone help me by giving suggestions for appropriate keywords for input?
Regards
Keshab K. Adhikary