Spectroscopy anharmonic vibrational

In resonant infrared multidimensional spectroscopies the excitation pulses couple directly to the transition dipoles. The lowest order possible technique in noncentrosymmetrical media involves three-pulses, and is, in general, three dimensional (Fig. 1A). Simulating the signal requires calculation of the third-order response function. In a small molecule this can be done by applying the sum-over-states expressions (see Appendix A), taking into account all possible Liouville space pathways described by the Feynman diagrams shown in Fig. IB. The third-order response of coupled anharmonic vibrations depends on the complete set of one- and two-exciton states coupled to thermal bath (18), and the sum-over-states approach rapidly becomes computationally more expensive as the molecule size is increased. 

Chaban GM, Gerber RB. Anharmonic vibrational spectroscopy of the glycine-water complex calculations for ab initio, empirical and hybrid quantum mechanics/molecular mechanics potentials. J Chem Phys 2001 115 1340-1348. 

First-principles calculations of anharmonic vibrational spectroscopy of large molecules... 

The structure of this review is as follows. In Section 9.2, we briefly discuss methods for computing vibrational states of systems having several coupled vibrational degrees of freedom. This will also cover methods that were not yet adapted for direct use with ab initio potentials, since in our view, such extensions may be possible in the future, at least for some of the algorithms. The focus will be on methods that seem potentially applicable to large polyatomics, rather than those of great accuracy for small systems. Section 9.3 also deals with computational methods for anharmonic vibrational spectroscopy that are applicable to potential surfaces from electronic structure calculations. Our main focus will be on the Vibrational Self-Consistent Field (VSCF) approach in several variants and extensions. The performance of the available method in the present state of the art is discussed in Section 9.4. Future directions are outlined in Section 9.5. 

A variety of methods for performing anharmonic vibrational spectroscopy computations were developed to address these and related systems. At the early stages, essentially all the methods were developed for potential surfaces available as explicit analytic functions of the coordinates. 

In briefly surveying several of the many methods proposed for anharmonic vibrational spectroscopy calculations, we mention also methods used so far only for analytic potential surfaces. For many of these methods, adaptations to algorithms may be feasible, and the methods seem promising in this respect. 

Matsunaga et al. [110] introduced VSCF-DPT2, a method that includes the effects of degeneracies in the anharmonic vibrational spectra. The essential extension is to use Degenerate Perturbation Theory (as opposed to Non-degenerate Perturbation Theory) in introducing correlation effects. Also this method was interfaced with electronic strucmre codes, and is incorporated in gamess. There have been several applications of ab initio spectroscopy calculations with this method. 

Major progress was made in recent years in ab initio calculations of anharmonic vibrational spectroscopy. One of the important indicators is the good agreement with experiment found in calculations for relatively large molecules, having more than 10 atoms (24 vibrational modes). Treatment of such large systems at a good anharmonic... 

In the spirit of the above observation, anharmonic vibrational spectroscopy calculations provide, by comparison with experiment, an evaluation of the quality of the potential surface used. It seems important to use spectroscopy as a way to compare the relative accuracies of different force fields. The results may, of course, differ from case to case. A crude general picture that seems to emerge from a limited set of small molecules (that includes HjO, HCOOH, CH3COOH) [104] indicates that MP2, B97 and B3LYP are very roughly comparable in the agreement with experimental spectroscopy, while HCTH and BLYP functionals do significantly less well (in this order). Much more can be learned about the quality of potential surfaces from different types of ab initio and DFT methods, but this will require systematic studies for various types of molecules. 

First-principles calculations of anharmonic vibrational spectroscopy of large molecules 187 Table 9.4 OH-stretching overtone excitation frequencies for HNO3... 

The field of Ab initio Anharmonic Vibrational Spectroscopy has seen rapid recent advancements. Already, several methods are at hand, and a variety of codes are available in several of the leading electronic structure suites of programs. It is tempting to speculate on the possible directions of progress in the next few years. 

With regard to the electronic structure methodology, major obstacles must be surmounted before improvements can be made. Calculations with Coupled-Cluster methods, an obvious next step, are far more computationally costly than the presently used MP2, or B3LYP methods. In fact, there are extremely few direct ab initio calculations of anharmonic vibrational spectroscopy at higher than MP2 or DPT levels, even for small polyatomics. From the point of view of ab initio anharmonic spectroscopy, the leap from MP2 to the Coupled-Cluster method seems a bottleneck. One can draw encouragement from faster Coupled-Cluster implementations, so far employed with the perturbation theory anharmonic analysis [116,117]. 

---