Garcia-Borràs Marc, Solà Miquel, Lauvergnat David, Reis Heribert, Luis Josep M, Kirtman Bernard
Institut de Química Computacional and Departament de Química, Universitat de Girona, Campus Montilivi, 17071 Girona, Catalonia, Spain.
Laboratoire de Chimie Physique, CNRS, Bât. 349, UMR 8000, Orsay F-91405, France and Université Paris-Sud, Orsay F-91405, France.
J Chem Theory Comput. 2013 Jan 8;9(1):520-32. doi: 10.1021/ct300805p. Epub 2012 Nov 12.
Previously, a reduced dimensionality approach was used to determine the vibrational contribution to nonlinear optical properties for molecules with large amplitude anharmonic modes that takes into account tunneling between potential wells (Luis, J. M.; Reis, H.; Papadopoulos, M. G.; Kirtman, B. J. Chem. Phys.2009, 131, 034116). Here, the treatment is extended, again using ammonia as an example, to include the remaining modes at several approximate levels. It is shown that this extension is essential to obtaining the correct results. Our new approach fully accounts for tunneling and avoids possible convergence problems associated with the normal coordinate expansion of the potential energy surface in a single-well treatment. For accurate numerical values, a good treatment of electron correlation is required along with a flexible basis set including diffuse functions.