Rearrangements and tunneling splittings of protonated water trimer
|
|
- Rosalyn Bailey
- 6 years ago
- Views:
Transcription
1 JOURNAL OF CHEMICAL PHYSICS VOLUME, NUMBER 8 8 NOVEMBER 999 Rearrangements and tunneling splittings of protonated water trimer David J. Wales University Chemical Laboratories, Lensfield Road, Cambridge CB2 EW, United Kingdom Received 3 April 999; accepted August 999 Stationary points and rearrangement mechanisms are characterized for protonated water trimer using a variety of basis sets and density functional theory to describe electron correlation. For the largest basis sets there are three distinct low-lying minima separated in energy by only a few wave numbers. Ten distinct transition states were found with barriers spanning nearly three orders of magnitude. Several of these mechanisms should produce observable tunneling splittings. 999 American Institute of Physics. S I. INTRODUCTION The promise of new experimental results for protonated water clusters from far-infrared vibration rotation tunneling FIR VRT spectroscopy 2 6 seems likely to spark a new wave of research in this field, as it did for neutral water clusters. The present work, and a preceding paper on H 5 O 2, 7 are intended to assist in assignment of these new spectra by characterizing rearrangement pathways and the associated tunneling splitting patterns. The magnitude of the splittings presents a difficult problem in quantum nuclear dynamics, and in the present work we only infer the likely feasibility of different rearrangements from the calculated barriers and path lengths. Although simplistic, this approach served to explain the observed splittings of water trimer 8 and predict those recently found for water pentamer. 9,0 The theoretical framework is identical to the procedure described in previous accounts of rearrangements and tunneling splittings in neutral 9, 6 and protonated 7 water clusters. Several experiments have previously probed the infrared spectrum of H 7 O 3 in the OH stretching region In particular, Okumura et al. confirmed the assignment by Schwarz of a band at 2670 cm to H 7 O Comparisons of the present results with previous ab initio calculations 2 30 are made in Sec. III. To deduce the molecular symmetry MS group 3 for a nonrigid molecule we must characterize the reaction graph for permutational isomerization, i.e., the transition states and pathways which connect the minima in question. As before, we adopt the notation of Bone et al. 32 where a structure is understood to mean a particular molecular geometry and a version is a particular labeled permutational isomer of a given structure. Versions that are directly connected by a single transition state are said to be adjacent with respect to the corresponding mechanism, and rearrangements which produce observable tunneling splittings are termed feasible. 3 The largest tunneling splittings are expected for degenerate rearrangements 33 which link permutational isomers of the same structure via a single transition state, i.e., adjacent versions. For the two largest basis sets three minima were found for H 7 O 3 in the present work, the lowest pair separated by a very small energy difference on the order of a few wave numbers, with two nondegenerate rearrangements linking them. The remaining eight transition states all mediate degenerate rearrangements, seven of which are symmetric with the two sides of the path related by a symmetry operation. 34 Throughout this paper we follow Murrell and Laidler s definition of a transition state as a stationary point with a single negative Hessian eigenvalue. 35 The reaction pathways are then defined by steepest-descent paths from the transition states which are, of course, independent of mass, temperature and coordinate system within the Born Oppenheimer approximation when defined properly in terms of covariant derivatives. 36,37 II. GEOMETRY OPTIMIZATION AND PATHWAYS As in our previous studies of water clusters all the stationary points were located using eigenvector-following in Cartesian coordinates and the scheme described previously for water pentamer. 9 Analytic first and second derivatives of the energy were used at every step with no symmetry constraints and were generated with the CADPAC program. 44 Three basis sets were considered, as for the previous study of H 5 O 2. 7 The first, DZPdiff, is based upon a double- 45,46 plus polarization DZP basis, with polarization functions consisting of a single set of p orbitals on each hydrogen atom exponent.0 and a set of six d orbitals on each oxygen atom exponent 0.9. To these functions were added a diffuse s orbital on each hydrogen atom exponent and diffuse s and p orbitals on each oxygen atom exponents and for s and p, respectively, 47 giving a total of 02 functions for H 7 O 3. We also employed the standard aug-cc-pvdz and aug-cc-pvtz basis sets, 48,49 with a total of 38 and 340 basis functions, respectively. Correlation corrections were obtained through density functional theory DFT, since MP2 second derivative calculations were impossible due to disk space limitations. We employed the Becke nonlocal exchange functional 50 and the Lee Yang Parr correlation functional 5 together referred to as BLYP. Derivatives of the grid weights were not included and the core electrons were not frozen. Numerical integration of the BLYP functionals was performed using the CADPAC HIGH option. Calculations were deemed to be converged when the root-mean-square gradient fell below 20 6 atomic units. Since derivatives of the grid weights were not included the zero frequencies were sometimes as large as /99/(8)/8429/9/$ American Institute of Physics
2 8430 J. Chem. Phys., Vol., No. 8, 8 November 999 David J. Wales cm. Pathways were calculated at the DZPdiff/BLYP and aug-cc-pvdz/blyp levels of theory and all thirteen stationary points were also optimized at the aug-cc-pvtz/blyp level. One refinement to the optimization scheme used for H 5 O 2 was necessary. It was impossible to calculate analytic second derivatives for the largest basis set due to program limitations. The aug-cc-pvtz/blyp geometries were therefore relaxed from those obtained at the aug-cc-pvdz/blyp level using the Murtagh Sargent Hessian update scheme 52 and an initial Hessian from the aug-cc-pvtz/blyp calculation. The convergence criteria were met within four or five eigenvector-following steps using the approximate Hessian, and the Hessian index did not change. Three parameters are useful to describe the rearrangement mechanisms. The first is the integrated path length, S, calculated as a sum over eigenvector-following steps, m: S m Q m Q m, where Q m is the 3n-dimensional position vector for n nuclei in Cartesian coordinates at step m. The second is the distance between the two minima in nuclear configuration space, D: DQsQ f, 2 where Q(s) and Q( f ) are the 3n-dimensional position vectors of the minima at the start and finish of the path. The third is the moment ratio of displacement, 53, which gives a measure of the cooperativity of the rearrangement: n iq i sq i f 4 i Q i sq i f 2, 3 2 where Q i (s) is the position vector of atom i in starting minimum s, etc. If every atom undergoes the same displacement then, while if only one atom moves then n. III. REARRANGEMENTS OF PROTONATED WATER TRIMER Three minima and ten transition states were identified for the two largest basis sets Fig. and Table I. However, the higher-lying minimum, min3, and the transition state that links it to min were not located with the DZPdiff basis. The rotational constants of the three minima are given in Table II. The point group symmetries and Hessian indices are the same for equivalent stationary points with all three basis sets. Harmonic frequencies for the minima and transition states at the aug-cc-pvdz/blyp level are collected in Tables III and IV. Since the low frequency vibrations, in particular, are expected to be rather anharmonic, the harmonic frequencies and zero-point energies tabulated here should be interpreted with caution. Path lengths and cooperativity indices for each rearrangement are given in Table V. Pathways were first calculated by taking small displacements of order 0.05a 0 away from a transition state parallel or antiparallel to the transition vector, and then employing eigenvector-following energy minimization to find the associated minimum. 54 The pathways obtained by this procedure have been compared to steepest-descent paths and pathways FIG.. Side views of the three minima found at the aug-cc-pvdz/blyp and aug-cc-pvtz/blyp levels of theory. that incorporate a kinetic metric 37 in previous work the mechanism is usually found to be represented correctly. 55,56 To test this result for the present system all the pathways were recalculated using the quadratic steepest-descent algorithm of Page and McIver. 57 The mechanisms were unchanged, except for the bifurcation mediated by tsb, described below, where a small perturbation changes one of the minima from min to min2 for the DZPdiff basis, and from min to min3 for the aug-cc-pvdz basis. min3 does not appear to exist for the smaller basis set. Aside from this result the values of the distance and cooperativity indices, D and, are practically identical for the eigenvector-following and steepest-descent paths. The integrated path length, however, is always somewhat shorter for the steepest-descent paths, typically by 0% 20%, and in one case by 40%. In view of these results all the tabulated data is for the steepestdescent paths, and all the illustrations are for the larger augcc-pvdz basis. In several previous studies a C 2v symmetry chain structure has been classified as a minimum. 2,23,24,29 In the present work this structure was found to be a transition state tsf2, in agreement with the MP2 calculations of Lee and Dyke. 26 Lee and Dyke also characterize a C s minimum which corresponds to min2 or min3. Two other studies have reported single minima which appear to have C symmetry. 27,30 Geissler et al. 30 also report two transition states for the Stillinger David potential 58 which correspond to tsrc and tsrc2 in the present study, and sampled dynamical transition paths for the corresponding mechanisms. They also report rotational exchange of protons on the terminal water monomers but did not characterize these processes further. The facile monomer inversion process corresponding to tsinv is shown in Fig. 2 and connects min and min2. The
3 J. Chem. Phys., Vol., No. 8, 8 November 999 Protonated water trimer 843 TABLE I. Energies/hartree Ref. 54 and point groups of the protonated water trimer minima and transition states for three different basis sets using the BLYP functional. Zero-point energies of the minima are also given, along with the energy separation of each stationary point from min at the same level of theory, E the value in brackets includes harmonic zero-point terms; the aug-cc-pvdz frequencies were used for the unknown aug-cc-pvtz zero-point energies. Group DZPdiff aug-cc-pvdz aug-cc-pvtz min C ZPE/cm min2 C s E/cm 7(65) 4(39) 2(38) min3 C s E/cm 54(39) 64(29) tsinv: monomer inversion C E/cm (90) 23(46) 29(5) tsrot: internal rotation C E/cm 90(32) 48(59) 36(48) tsf: flip C s E/cm 202(242) 29(274) 279(333) tsf2: flip 2 C 2v E/cm 204(392) 23(83) 274(244) tsf3: flip 3 C 2v E/cm 268(55) 334(2) tsrc: ring closure C E/cm 2208(999) 223(2004) 265(956) tsrc2: ring closure 2 C s E/cm 25(2725) 2466(249) 2420(2445) tsb: bifurcation C E/cm 2662(249) 265(2393) 2623(2402) tsb2: double bifurcation C s E/cm 6409(6740) 625(6463) 699(6447) tsb3: bifurcation C s E/cm 739(757) 779(74) 7237(7200) barriers are even smaller than for H 5 O 2. 7 Only one internal rotation mechanism was found in the present work, in contrast to two for H 5 O 2. 7 It also connects min and min2 and involves low barriers Fig. 3. Profiles for these two paths are shown in Fig. 4. Symmetric degenerate rearrangements corresponding to flipping of the unbound proton on the central oxygen were found for all three minima Fig. 5. These processes have relatively low barriers. The five remaining pathways all involve barriers that are an order of magnitude larger. The two mechanisms involving ring closure first found by Geissler et al. 30 are symmetric degenerate rearrangements of min in the present calculations Fig. 6. The rearrangements mediated by tsb, tsb2, and tsb3 all involve bifurcations, two in the case of tsb2. The barriers involved for tsb are similar to those for ring TABLE II. Rotational constants/cm of the three minima for the three different basis sets. DZPdiff aug-cc-pvdz aug-cc-pvtz C min B A C min2 B A C min3 B A closure, while the others are significantly larger. All three processes are shown in Fig. 7. There are a total of 23!7! distinct versions of the H 7 O 3 or D 5 O 2 ) C global minimum, where the first TABLE III. Harmonic frequencies/cm of min, min2, and min3 at the aug-cc-pvdz/blyp level of theory. Intensities in km/mol are given in square brackets. min min2 min
4 8432 J. Chem. Phys., Vol., No. 8, 8 November 999 David J. Wales TABLE IV. Harmonic frequencies/cm of the H 7 O 3 transition states at the aug-cc-pvdz/blyp level of theory. tsinv tsrot tsf tsf2 tsf3 tsrc tsrc2 tsb tsb2 tsb3 44i 90i 253i 27i 238i 73i 88i 370i 557i 64i factor accounts for the inversion operation and there are 3! and 7! permutations of the oxygen and hydrogen or deuterium atoms, respectively. For min2 and min3 we must divide by two to account for the order of the point group, 32 and there are only distinct versions. The reaction graph for one set of labeled minima involving the most likely feasible mechanisms, i.e. those mediated by monomer inversion, internal rotation and the two flips tsinv, tsrot, tsf and tsf2, is shown in Fig. 8. min3 is not included to simplify the analysis and because it lies somewhat higher in energy than min and min2. The presence of two low energy minima complicates the molecular symmetry group analysis. One way to treat such cases is to use effective generators which relate permuta- TABLE V. Properties of the steepest-descent pathways found for H 7 O 3 with the two smaller basis sets. S, D, and are all defined in Sec. II. Parameter DZPdiff aug-cc-pvdz DZPdiff aug-cc-pvdz tsinv: monomer inversion tsrot: internal rotation S D tsf: flip tsf2: flip 2 S D tsf: flip 3 tsrc: ring closure S D tsrc2: ring closure 2 tsb: bifurcation S D tsb2: double bifurcation tsb3: bifurcation 3 S D
5 J. Chem. Phys., Vol., No. 8, 8 November 999 Protonated water trimer 8433 FIG. 4. Energy profiles of the aug-cc-pvdz/blyp steepest-descent paths for monomer inversion tsinv and internal rotation tsrot. terms of the effective tunneling matrix elements 2 and f2, and assuming a Hückel-type approximation, the splitting pattern is quite simple: FIG. 2. Two views of the monomer inversion pathway for H 7 O 3 calculated at the aug-cc-pvdz/blyp level. tional isomers of the same structure. 59 For min2 suitable generators are AC35627* reflection symmetry, 2 tsinv followed by tsrot, or vice versa and AC72635 the flip, tsf2, and the versions of min2 are then linked in sets of eight. The resulting MS group is the G(6) group first characterized by Dyke for neutral water dimer, 60 which is also appropriate for H 5 O 2. 7 The character table and nuclear spin weights are given in Table VI. In 2 2 f2 2 2 f2 f2 f2 2 2 f2 A B E E B f2 A 2. FIG. 3. Two views of the internal rotation pathway of H 7 O 3 calculated at the aug-cc-pvdz/blyp level. FIG. 5. Flip rearrangements of H 7 O 3 corresponding to tsf, tsf2, and tsf3 calculated at the aug-cc-pvdz/blyp level.
6 8434 J. Chem. Phys., Vol., No. 8, 8 November 999 David J. Wales FIG. 6. Ring closure rearrangements of H 7 O 3 corresponding to tsrc left and tsrc2 right calculated at the aug-cc-pvdz/blyp level. FIG. 7. Pathways of H 7 O 3 involving bifurcated transition states calculated at the aug-cc-pvdz/blyp level. If f2 is smaller than 2 the result is a triplet of doublets where the doublets are split by 2 f2 and the components of the triplet are split by 2 2. For min five generators are required to specify all the likely connections, namely 267* for the flip tsf and AC35672*, AC35726*, AC35726* and AC35627* for the indirect processes. The versions of min are then linked in sets of 6 and the appropriate MS group is again G(6). An accidental degeneracy arises in the tunneling spectrum which complicates the symmetry assignments. However, it can be removed by including the generator operation AC35627 as a small perturbation. The resulting spectrum includes a set of levels which match those found for min2, above, and eight more levels with wavefunctions that are not delocalized over the min2 isomers. Since min and min2 are so close in energy we will present a full analysis of the combined spectrum rather than focusing further on the results for the different minima. If min2 lies at an energy above min, and we adopt a Hückel-type approximation for the secular determinant, then an analytic form can be found for the energy levels: 2 8 inv rot 2 f f2 2 2 f f inv rot 2 f f2 2 2 f f inv 2 rot f f2 2 2 f f inv 2 rot f f2 2 2 f f inv rot 2 f f2 2 2 f f inv rot 2 f f2 2 2 f f2 2 f f A B E E B 2 A 2 A 2,B,E A,B 2,E
7 J. Chem. Phys., Vol., No. 8, 8 November 999 Protonated water trimer inv rot 2 f f2 2 2 f f2 2 B inv rot 2 f f2 2 2 f f2 2 A inv 2 rot f f2 2 2 f f2 2 E inv 2 rot f f2 2 2 f f2 2 E 2 8 inv rot 2 f f2 2 2 f f2 2 A 2 8 inv rot 2 f f2 2 2 f f2 2 B. FIG. 8. Reaction graph for H 7 O 3 min and min2 including pathways mediated by monomer inversion, internal rotation and the two flips.
8 8436 J. Chem. Phys., Vol., No. 8, 8 November 999 David J. Wales TABLE VI. Character table for the MS group G(6) and nuclear spin statistical weights for H 7 O 3 and D 7 O 3. G(6) E (2)(67) (2) 67 (AC)(35)(726) AC35627 (AC)(35)(7)(26) AC35627 E* (2)(67)* (2)* 67* (AC)(35)(726)* AC35627* (AC)(35)(7)(26)* AC35627* H 7 O 3 D 7 O 3 A A B B E A A B B E From the calculated barrier heights and path lengths a reasonable guess for the magnitudes of the tunneling matrix elements in wavenumbers for H 7 O 3 corresponding to tsinv-tsf2 is probably inv : rot : f : f2 0::0.:0.. The energy levels are then arranged in the order given above. Whatever the value of the levels always appear in doublets split by f f2. This doubling may be the smallest observable splitting if none of the other mechanisms corresponding to tsrc-tsb3 are feasible. In the limit 0 the energy levels obey a mirror relation about the energy zero with two triplets of doublets at around 2 inv, assuming that inv is the largest matrix element Fig. 9. The triplet spacing is roughly 2 rot if we also neglect f and f2. When is large compared to the tunneling matrix elements the two triplets of doublets move to energies around zero and Fig. 9 as the wavefunctions become more strongly localized on either min or min2 isomers. The triplet splitting is around 4 inv rot /. The expected intensities of the lines in the triplet of doublets are 9:9:30:30:24:26 or 8:60:50 if the splitting due to the flip is not resolved. For D 7 O 3 the corresponding ratios are 339:336:540:540:26:26 and 675:080:432. If any of the other mechanisms corresponding to the remaining transition states are feasible then the MS group will be enlarged and further splittings of the tunneling states will occur. However, the present results suggest that the effects of ring closure and bifurcation are probably beyond the limit of current experimental resolution, especially for D 7 O 3. Although bifurcation tunneling is resolved in (H 2 O) 3,(D 2 O) 3,2,6 63 and probably (H 2 O) 5, 9,0 the barriers in H 7 O 3 are higher. IV. CONCLUSIONS The present calculations of rearrangement pathways in protonated water trimer reveal a more complicated situation than for the protonated dimer, with two low energy minima interconverting via a variety of rearrangements and a third minimum lying only a little higher in energy before zeropoint effects are included. The nondegenerate monomer inversion and internal rotation processes seem likely to produce observable splittings via indirect tunneling 4 in which both min and min2 are involved. The degenerate flip rearrangements involve higher barriers and should result in an additional doublet splitting if they are feasible. The other mechanisms, including bifurcation processes, involve significantly larger barriers and do not appear to have been observed in previous quantum nuclear dynamics simulations of either H 5 O 2 64 or a hydrated proton in water. 65 Calculations on protonated water clusters are notoriously sensitive to changes in basis set and treatment of electron correlation, 66 and the additional stationary points found with the two largest basis sets provide further evidence of this problem for H 7 O 3. Extensive searches for lower energy pathways connecting min3 to either min or min2 were not conducted. However, the patterns predicted in the present work would probably not change much if such additional paths exist or if the details of the low energy pathways were altered, so long as the connectivity of the reaction graph is unaffected. New results of high resolution spectroscopy in the far infrared should soon test this theory. ACKNOWLEDGMENTS The author gratefully acknowledges the support of the Royal Society of London and the EPSRC. FIG. 9. Splitting pattern for H 7 O 3 calculated using a Hückel-type approximation for 0, 20, and 40 cm. The other parameters were fixed at inv 0, rot.5 and f f2 0.4, all in wave numbers. R. J. Saykally personal communication. 2 R. C. Cohen and R. J. Saykally, J. Phys. Chem. 94, N. Pugliano and R. J. Saykally, J. Chem. Phys. 96, R. J. Saykally and G. A. Blake, Science 259, K. Liu, J. D. Cruzan, and R. J. Saykally, Science 27,
9 J. Chem. Phys., Vol., No. 8, 8 November 999 Protonated water trimer J. B. Paul, C. P. Collier, R. J. Saykally, J. J. Scherer, and A. O Keefe, J. Phys. Chem. A 0, D. J. Wales, J. Chem. Phys. 0, D. J. Wales, J. Am. Chem. Soc. 5, D. J. Wales and T. R. Walsh, J. Chem. Phys. 05, M. G. Brown, F. N. Keutsch, and R. J. Saykally, J. Chem. Phys. 09, D. J. Wales, J. Am. Chem. Soc. 5, T. R. Walsh and D. J. Wales, J. Chem. Soc., Faraday Trans. 92, D. J. Wales and T. R. Walsh, J. Chem. Phys. 06, D. J. Wales, in Advances in Molecular Vibrations and Collision Dynamics, edited by J. M. Bowman and Z. Bačić JAI Press, Stamford, 998, Vol. 3, pp D. J. Wales, in Theory of Atomic and Molecular Clusters, edited by J. Jellinek Springer-Verlag, Heidelberg, 999, pp D. J. Wales, in Recent Theoretical and Experimental Advances in Hydrogen-Bonded Clusters, edited by S. Xantheas Kluwer, Dordrecht, in press available from the Los Alamos preprint server at URL xxx.lanl.gov/abs/physics/ H. A. Schwarz, J. Chem. Phys. 67, M. Okumura, L. I. Yeh, J. D. Myers, and Y. T. Lee, J. Chem. Phys. 85, L. I. Yeh, M. Okumura, J. D. Myers, J. M. Price, and Y. T. Lee, J. Chem. Phys. 9, M. Okumura, L. I. Yeh, J. D. Myers, and Y. T. Lee, J. Phys. Chem. 94, M. D. Newton and S. Ehrenson, J. Am. Chem. Soc. 93, M. D. Newton, J. Chem. Phys. 67, C. A. Deakyne, M. Meot-Ner, C. L. Campbell, M. G. Hughes, and S. P. Murphy, J. Chem. Phys. 84, G. R. J. Williams, J. Mol. Struct.: THEOCHEM 38, E. P. F. Lee, J. M. Dyke, A. E. Wilders, and P. Watts, Mol. Phys. 7, E. P. F. Lee and J. M. Dyke, Mol. Phys. 73, D. Wei and D. R. Salahub, J. Chem. Phys. 0, H.-P. Cheng, J. Phys. Chem. A 02, U. W. Schmitt and G. A. Voth, J. Phys. Chem. B 02, P. L. Geissler, C. Dellago, and D. Chandler, Phys. Chem. Chem. Phys., H. C. Longuet-Higgins, Mol. Phys. 6, R. G. A. Bone, T. W. Rowlands, N. C. Handy, and A. J. Stone, Mol. Phys. 72, R. E. Leone and P. v. R. Schleyer, Angew. Chem. Int. Ed. Engl. 9, J. G. Nourse, J. Am. Chem. Soc. 02, J. N. Murrell and K. J. Laidler, Trans. Faraday Soc. 64, B. Friedrich, Z. Herman, R. Zahradnik, and Z. Havlas, Adv. Quantum Chem. 9, A. Banerjee and N. P. Adams, Int. J. Quantum Chem. 43, C. J. Cerjan and W. H. Miller, J. Chem. Phys. 75, J. Simons, P. Jorgensen, H. Taylor, and J. Ozment, J. Phys. Chem. 87, D. O Neal, H. Taylor, and J. Simons, J. Phys. Chem. 88, A. Banerjee, N. Adams, J. Simons, and R. Shepard, J. Phys. Chem. 89, J. Baker, J. Comput. Chem. 7, J. Baker, J. Comput. Chem. 8, R. D. Amos, I. L. Alberts, J. S. Andrews, S. M. Colwell, N. C. Handy, D. Jayatilaka, P. J. Knowles, R. Kobayashi, K. E. Laidig, G. Laming, A. M. Lee, P. E. Maslen, C. W. Murray, J. E. Rice, E. D. Simandiras, A. J. Stone, M.-D. Su, and D. J. Tozer, The Cambridge analytic derivatives package Issue T. H. Dunning, J. Chem. Phys. 53, S. J. Huzinaga, J. Chem. Phys. 42, J. E. Fowler and H. F. Schaefer, J. Am. Chem. Soc. 7, T. H. Dunning, J. Chem. Phys. 90, R. A. Kendall, T. H. Dunning, and R. J. Harrison, J. Chem. Phys. 96, A. D. Becke, Phys. Rev. A 38, C. Lee, W. Yang, and R. G. Parr, Phys. Rev. B 37, B. A. Murtagh and R. W. H. Sargent, Comput. J. UK 3, F. H. Stillinger and T. A. Weber, Phys. Rev. A 28, The atomic units of distance and energy, a 0 and h, are m and J, respectively. 55 T. R. Walsh and D. J. Wales, Z. Phys. D 40, G. N. Merrill and M. S. Gordon, J. Phys. Chem. A 02, M. Page and J. W. McIver, J. Chem. Phys. 88, F. H. Stillinger and C. W. David, J. Chem. Phys. 69, B. J. Dalton and P. D. Nicholson, Int. J. Quantum Chem. 9, T. R. Dyke, J. Chem. Phys. 66, N. Pugliano and R. J. Saykally, Science 257, E. H. T. Olthof, A. Van der Avoird, P. E. S. Wormer, K. Liu, and R. J. Saykally, J. Chem. Phys. 05, A. Van der Avoird, E. H. T. Olthof, and P. E. S. Wormer, J. Chem. Phys. 05, M. E. Tuckerman, D. Marx, M. L. Klein, and M. Parrinello, Science 275, D. Marx, M. E. Tuckerman, J. Hutter, and M. Parrinello, Nature 397, E. F. Valeev and H. F. Schaefer, J. Chem. Phys. 08,
Rearrangements and tunneling splittings of protonated water dimer
JOURNAL OF CHEMICAL PHYSICS VOLUME 110, NUMBER 21 1 JUNE 1999 Rearrangements and tunneling splittings of protonated water dimer David J. Wales University Chemical Laboratories, Lensfield Road, Cambridge
More informationRearrangements of Water Dimer and Hexamer
. Rearrangements of Water Dimer and Hexamer David J. Wales University Chemical Laboratories, Lensfield Road, Cambridge, CB2 1EW, United Kingdom Abstract. Rearrangement mechanisms of the water dimer and
More informationTheoretical study of the water tetramer
Theoretical study of the water tetramer David J. Wales and Tiffany R. Walsh University Chemical Laboratories, Lensfield Road, Cambridge CB2 1EW, United Kingdom Received 18 December 1996; accepted 22 January
More informationThe far-infrared vibration rotation tunneling spectrum of the water tetramer-d8
The far-infrared vibration rotation tunneling spectrum of the water tetramer-d8 J. D. Cruzan, M. G. Brown, K. Liu, L. B. Braly, and R. J. Saykally Department of Chemistry, University of California at Berkeley,
More informationDFT calculations of NMR indirect spin spin coupling constants
DFT calculations of NMR indirect spin spin coupling constants Dalton program system Program capabilities Density functional theory Kohn Sham theory LDA, GGA and hybrid theories Indirect NMR spin spin coupling
More informationExample questions for Molecular modelling (Level 4) Dr. Adrian Mulholland
Example questions for Molecular modelling (Level 4) Dr. Adrian Mulholland 1) Question. Two methods which are widely used for the optimization of molecular geometies are the Steepest descents and Newton-Raphson
More informationNon-covalent force fields computed ab initio
Non-covalent force fields computed ab initio Supermolecule calculations Symmetry-adapted perturbation theory (SAPT) Supermolecule calculations Requirements: E = E AB E A E B. Include electron correlation,
More informationOslo node. Highly accurate calculations benchmarking and extrapolations
Oslo node Highly accurate calculations benchmarking and extrapolations Torgeir Ruden, with A. Halkier, P. Jørgensen, J. Olsen, W. Klopper, J. Gauss, P. Taylor Explicitly correlated methods Pål Dahle, collaboration
More informationChemistry 483 Lecture Topics Fall 2009
Chemistry 483 Lecture Topics Fall 2009 Text PHYSICAL CHEMISTRY A Molecular Approach McQuarrie and Simon A. Background (M&S,Chapter 1) Blackbody Radiation Photoelectric effect DeBroglie Wavelength Atomic
More informationProfessor N. C. Handy and Dr. R. D. Amos University Chemical Laboratory Lensfield Road Cambridge CB2 lew UK. 1,P r, '.*-, -..
/A THE DEVELOPMENT OF QUANTUM CHEMISTRY CODES AD-A277 149 * IEIIIllil Professor N. C. Handy and Dr. R. D. Amos University Chemical Laboratory Lensfield Road Cambridge CB2 lew UK 1,P r, '.*-, -.. *. Grant
More informationAssignment 1: Molecular Mechanics (PART 2 25 points)
Chemistry 380.37 Fall 2015 Dr. Jean M. Standard September 2, 2015 Assignment 1: Molecular Mechanics (PART 2 25 points) In this assignment, you will perform some additional molecular mechanics calculations
More informationQuantitative characterization of the water trimer torsional manifold by terahertz laser spectroscopy and theoretical analysis. II.
JOURNAL OF CHEMICAL PHYSICS VOLUME 111, NUMBER 17 1 NOVEMBER 1999 Quantitative characterization of the water trimer torsional manifold by terahertz laser spectroscopy and theoretical analysis. II. H 2
More informationTerahertz vibration-rotation-tunneling spectroscopy of the water tetramer-d8: Combined analysis of vibrational bands at 4.1 and 2.
THE JOURNAL OF CHEMICAL PHYSICS 128, 094302 2008 Terahertz vibration-rotation-tunneling spectroscopy of the water tetramer-d8: Combined analysis of vibrational bands at 4.1 and 2.0 THz Wei Lin, a Jia-Xiang
More informationExploring the energy landscape
Exploring the energy landscape ChE210D Today's lecture: what are general features of the potential energy surface and how can we locate and characterize minima on it Derivatives of the potential energy
More informationMO Calculation for a Diatomic Molecule. /4 0 ) i=1 j>i (1/r ij )
MO Calculation for a Diatomic Molecule Introduction The properties of any molecular system can in principle be found by looking at the solutions to the corresponding time independent Schrodinger equation
More informationWater Pentamer: Characterization of the Torsional-Puckering Manifold by Terahertz VRT Spectroscopy
J. Phys. Chem. A 2005, 109, 6483-6497 6483 Water Pentamer: Characterization of the Torsional-Puckering Manifold by Terahertz VRT Spectroscopy Heather A. Harker, Mark R. Viant, Frank N. Keutsch, Ernest
More informationA contribution to the understanding of the structure of xenon hexafluoride
A contribution to the understanding of the structure of xenon hexafluoride T. Daniel Crawford, a) Kristen W. Springer, b) and Henry F. Schaefer III Center for Computational Quantum Chemistry, University
More informationWater dimer hydrogen bond stretch, donor torsion overtone, and in-plane bend vibrations
JOURNAL OF CHEMICAL PHYSICS VOLUM19, NUMBER 17 1 NOVEMBER 2003 Water dimer hydrogen bond stretch, donor torsion overtone, and in-plane bend vibrations Frank N. Keutsch Department of Chemistry and Chemical
More informationChemistry 3502/4502. Final Exam Part I. May 14, 2005
Advocacy chit Chemistry 350/450 Final Exam Part I May 4, 005. For which of the below systems is = where H is the Hamiltonian operator and T is the kinetic-energy operator? (a) The free particle
More informationUptake of OH radical to aqueous aerosol: a computational study
Uptake of OH radical to aqueous aerosol: a computational study Grigory Andreev Karpov Institute of Physical Chemistry 10 Vorontsovo pole, Moscow, 105064, Russia Institute of Physical Chemistry and Electrochemistry
More informationJack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah
1. Born-Oppenheimer approx.- energy surfaces 2. Mean-field (Hartree-Fock) theory- orbitals 3. Pros and cons of HF- RHF, UHF 4. Beyond HF- why? 5. First, one usually does HF-how? 6. Basis sets and notations
More informationSTRUCTURAL DEFECTS IN IMIDATES : AN AB INITIO STUDY
Int. J. Chem. Sci.: 9(4), 2011, 1763-1767 ISSN 0972-768X www.sadgurupublications.com STRUCTURAL DEFECTS IN IMIDATES : AN AB INITIO STUDY M. FATHIMA BEGUM, HEMA TRESA VARGHESE a, Y. SHEENA MARY a, C. YOHANNAN
More informationLecture 9 Electronic Spectroscopy
Lecture 9 Electronic Spectroscopy Molecular Orbital Theory: A Review - LCAO approximaton & AO overlap - Variation Principle & Secular Determinant - Homonuclear Diatomic MOs - Energy Levels, Bond Order
More informationAB INITIO MODELING OF THE STRUCTURAL DEFECTS IN AMIDES
Int. J. Chem. Sci.: 9(4), 2011, 1564-1568 ISSN 0972-768X www.sadgurupublications.com AB INITIO MODELING OF THE STRUCTURAL DEFECTS IN AMIDES M. FATHIMA BEGUM, HEMA TRESA VARGHESE a, Y. SHEENA MARY a, C.
More informationAtom-molecule molecule collisions in spin-polarized polarized alkalis: potential energy surfaces and quantum dynamics
Atom-molecule molecule collisions in spin-polarized polarized alkalis: potential energy surfaces and quantum dynamics Pavel Soldán, Marko T. Cvitaš and Jeremy M. Hutson University of Durham with Jean-Michel
More information2m dx 2. The particle in a one dimensional box (of size L) energy levels are
Name: Chem 3322 test #1 solutions, out of 40 marks I want complete, detailed answers to the questions. Show all your work to get full credit. indefinite integral : sin 2 (ax)dx = x 2 sin(2ax) 4a (1) with
More informationChemistry 334 Part 2: Computational Quantum Chemistry
Chemistry 334 Part 2: Computational Quantum Chemistry 1. Definition Louis Scudiero, Ben Shepler and Kirk Peterson Washington State University January 2006 Computational chemistry is an area of theoretical
More informationRotations and vibrations of polyatomic molecules
Rotations and vibrations of polyatomic molecules When the potential energy surface V( R 1, R 2,..., R N ) is known we can compute the energy levels of the molecule. These levels can be an effect of: Rotation
More information( )! rv,nj ( R N )! ns,t
Chapter 8. Nuclear Spin Statistics Notes: Most of the material presented in this chapter is taken from Bunker and Jensen (2005) Chap. 9 and Bunker and Jensen (1998) Chap. 8. 8.1 The Complete Internal Wave
More informationPhysical Chemistry Chemical Physics
Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2018 Physical Chemistry Chemical Physics Supporting Information Thermochemistry of the
More informationSymmetry: Translation and Rotation
Symmetry: Translation and Rotation The sixth column of the C 2v character table indicates the symmetry species for translation along (T) and rotation about (R) the Cartesian axes. y y y C 2 F v (x) T x
More informationQuantum mechanics can be used to calculate any property of a molecule. The energy E of a wavefunction Ψ evaluated for the Hamiltonian H is,
Chapter : Molecules Quantum mechanics can be used to calculate any property of a molecule The energy E of a wavefunction Ψ evaluated for the Hamiltonian H is, E = Ψ H Ψ Ψ Ψ 1) At first this seems like
More informationConcerted Hydrogen-Bond Breaking by Quantum Tunneling in the Water Hexamer Prism
Concerted Hydrogen-Bond Breaking by Quantum Tunneling in the Water Hexamer Prism Jeremy O. Richardson, 1,2 Cristóbal Pérez, 3 Simon Lobsiger, 3 Adam A. Reid, 1 Berhane Temelso, 4 George C. Shields, 4 Zbigniew
More informationChemistry 881 Lecture Topics Fall 2001
Chemistry 881 Lecture Topics Fall 2001 Texts PHYSICAL CHEMISTRY A Molecular Approach McQuarrie and Simon MATHEMATICS for PHYSICAL CHEMISTRY, Mortimer i. Mathematics Review (M, Chapters 1,2,3 & 4; M&S,
More informationChemistry 3502/4502. Final Exam Part I. May 14, 2005
Chemistry 3502/4502 Final Exam Part I May 14, 2005 1. For which of the below systems is = where H is the Hamiltonian operator and T is the kinetic-energy operator? (a) The free particle (e) The
More informationRenner-Teller Effect in Tetra-Atomic Molecules
Groupe de Chimie Théorique du MSME Renner-Teller Effect in Tetra-Atomic Molecules Laurent Jutier, G. Dhont, H. Khalil and C. Léonard jutier@univ-mlv.fr (non linear) Outline General Presentation Structure
More informationStructures and infrared spectra of fluoride hydrogen sulfide clusters from ab initio calculations: F -(H 2 S) n, n = 1 5w
RESEARCH PAPER Structures and infrared spectra of fluoride hydrogen sulfide clusters from ab initio calculations: F -(H 2 S) n, n = 1 5w D. A. Wild* and T. Lenzer PCCP www.rsc.org/pccp MPI fu r biophysikalische
More informationARTICLES. Normal-mode analysis without the Hessian: A driven molecular-dynamics approach
JOURNAL OF CHEMICAL PHYSICS VOLUME 119, NUMBER 2 8 JULY 2003 ARTICLES Normal-mode analysis without the Hessian: A driven molecular-dynamics approach Joel M. Bowman, a) Xiubin Zhang, and Alex Brown Cherry
More informationTools for QM studies of large systems
Tools for QM studies of large systems Automated, hessian-free saddle point search & characterization QM/MM implementation for zeolites Shaama Mallikarjun Sharada Advisors: Prof. Alexis T Bell, Prof. Martin
More informationJack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah
1. Born-Oppenheimer approx.- energy surfaces 2. Mean-field (Hartree-Fock) theory- orbitals 3. Pros and cons of HF- RHF, UHF 4. Beyond HF- why? 5. First, one usually does HF-how? 6. Basis sets and notations
More informationFigure 1: Transition State, Saddle Point, Reaction Pathway
Computational Chemistry Workshops West Ridge Research Building-UAF Campus 9:00am-4:00pm, Room 009 Electronic Structure - July 19-21, 2016 Molecular Dynamics - July 26-28, 2016 Potential Energy Surfaces
More informationChem 673, Problem Set 5 Due Thursday, November 29, 2007
Chem 673, Problem Set 5 Due Thursday, November 29, 2007 (1) Trigonal prismatic coordination is fairly common in solid-state inorganic chemistry. In most cases the geometry of the trigonal prism is such
More informationComputational Methods. Chem 561
Computational Methods Chem 561 Lecture Outline 1. Ab initio methods a) HF SCF b) Post-HF methods 2. Density Functional Theory 3. Semiempirical methods 4. Molecular Mechanics Computational Chemistry " Computational
More informationQUANTUM CHEMISTRY PROJECT 3: PARTS B AND C
Chemistry 460 Fall 2017 Dr. Jean M. Standard November 6, 2017 QUANTUM CHEMISTRY PROJECT 3: PARTS B AND C PART B: POTENTIAL CURVE, SPECTROSCOPIC CONSTANTS, AND DISSOCIATION ENERGY OF DIATOMIC HYDROGEN (20
More informationTerahertz Laser Spectroscopy of the Water Pentamer: Structure and Hydrogen Bond Rearrangement Dynamics
J. Phys. Chem. A 1997, 101, 9011-9021 9011 Terahertz Laser Spectroscopy of the Water Pentamer: Structure and Hydrogen Bond Rearrangement Dynamics K. Liu, M. G. Brown, J. D. Cruzan, and R. J. Saykally*
More informationQuantum chemistry and vibrational spectra
Chapter 3 Quantum chemistry and vibrational spectra This chapter presents the quantum chemical results for the systems studied in this work, FHF (Section 3.) and OHF (Section 3.3). These triatomic anions
More informationFinite Element Analysis of Molecular Rydberg States
Excerpt from the Proceedings of the COMSOL Conference 2009 Boston Finite Element Analysis of Molecular Rydberg States M. G. Levy, R. M. Stratt, and P. M. Weber* Brown University Department of Chemistry
More informationChem120a : Exam 3 (Chem Bio) Solutions
Chem10a : Exam 3 (Chem Bio) Solutions November 7, 006 Problem 1 This problem will basically involve us doing two Hückel calculations: one for the linear geometry, and one for the triangular geometry. We
More informationAn Introduction to Quantum Chemistry and Potential Energy Surfaces. Benjamin G. Levine
An Introduction to Quantum Chemistry and Potential Energy Surfaces Benjamin G. Levine This Week s Lecture Potential energy surfaces What are they? What are they good for? How do we use them to solve chemical
More informationChem 442 Review of Spectroscopy
Chem 44 Review of Spectroscopy General spectroscopy Wavelength (nm), frequency (s -1 ), wavenumber (cm -1 ) Frequency (s -1 ): n= c l Wavenumbers (cm -1 ): n =1 l Chart of photon energies and spectroscopies
More informationAb initio study of rearrangements between C 60 fullerenes
Chemical Physics Letters 374 (2003) 125 131 www.elsevier.com/locate/cplett Ab initio study of rearrangements between C 60 fullerenes Yuko Kumeda, David J. Wales * The University Chemical Laboratory, Lensfield
More informationBeyond the Hartree-Fock Approximation: Configuration Interaction
Beyond the Hartree-Fock Approximation: Configuration Interaction The Hartree-Fock (HF) method uses a single determinant (single electronic configuration) description of the electronic wavefunction. For
More informationNPTEL/IITM. Molecular Spectroscopy Lectures 1 & 2. Prof.K. Mangala Sunder Page 1 of 15. Topics. Part I : Introductory concepts Topics
Molecular Spectroscopy Lectures 1 & 2 Part I : Introductory concepts Topics Why spectroscopy? Introduction to electromagnetic radiation Interaction of radiation with matter What are spectra? Beer-Lambert
More informationMaking Electronic Structure Theory work
Making Electronic Structure Theory work Ralf Gehrke Fritz-Haber-Institut der Max-Planck-Gesellschaft Berlin, 23 th June 2009 Hands-on Tutorial on Ab Initio Molecular Simulations: Toward a First-Principles
More informationPractical Advice for Quantum Chemistry Computations. C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology
Practical Advice for Quantum Chemistry Computations C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology Choice of Basis Set STO-3G is too small 6-31G* or 6-31G** 6 probably
More informationChem 673, Problem Set 5 Due Thursday, December 1, 2005
otton, Problem 9.3 (assume D 4h symmetry) Additional Problems: hem 673, Problem Set 5 Due Thursday, December 1, 2005 (1) Infrared and Raman spectra of Benzene (a) Determine the symmetries (irreducible
More informationThe Potential Energy Surface
The Potential Energy Surface In this section we will explore the information that can be obtained by solving the Schrödinger equation for a molecule, or series of molecules. Of course, the accuracy of
More informationClose coupling results for inelastic collisions of NH3 and Ar. A stringent test of a spectroscopic potential
12 August 1994 ELSEVIER Chemical Physics Letters 226 ( 1994) 22-26 CHEMICAL PHYSICS LETTERS Close coupling results for inelastic collisions of NH3 and Ar. A stringent test of a spectroscopic potential
More informationDensity Functional Theory
Density Functional Theory March 26, 2009 ? DENSITY FUNCTIONAL THEORY is a method to successfully describe the behavior of atomic and molecular systems and is used for instance for: structural prediction
More informationTheoretical determination of the heat of formation of methylene
Theoretical determination of the heat of formation of methylene Nikos L. Doltsinis and Peter J. Knowles School of Chemistry, University of Birmingham, Edgbaston, Birmingham B5 2TT, United Kingdom The heat
More informationExperiment 11: NUCLEAR MAGNETIC RESONANCE SPECTROSCOPY
Experiment 11: NUCLEAR MAGNETIC RESONANCE SPECTROSCOPY Purpose: This is an exercise to introduce the use of nuclear magnetic resonance spectroscopy, in conjunction with infrared spectroscopy, to determine
More informationHONOUR SCHOOL OF NATURAL SCIENCE. Final Examination GENERAL PHYSICAL CHEMISTRY I. Answer FIVE out of nine questions
HONOUR SCHOOL OF NATURAL SCIENCE Final Examination GENERAL PHYSICAL CHEMISTRY I Monday, 12 th June 2000, 9.30 a.m. - 12.30 p.m. Answer FIVE out of nine questions The numbers in square brackets indicate
More informationSpectroscopy Of Hydrogen-bonded Formanilide Clusters In A Supersonic Jet: Solvation Of A Model Trans Amide
Bowling Green State University ScholarWorks@BGSU Chemistry Faculty Publications Chemistry 6-2000 Spectroscopy Of Hydrogen-bonded Formanilide Clusters In A Supersonic Jet: Solvation Of A Model Trans Amide
More information6.2 Polyatomic Molecules
6.2 Polyatomic Molecules 6.2.1 Group Vibrations An N-atom molecule has 3N - 5 normal modes of vibrations if it is linear and 3N 6 if it is non-linear. Lissajous motion A polyatomic molecule undergoes a
More informationChapter IV: Electronic Spectroscopy of diatomic molecules
Chapter IV: Electronic Spectroscopy of diatomic molecules IV.2.1 Molecular orbitals IV.2.1.1. Homonuclear diatomic molecules The molecular orbital (MO) approach to the electronic structure of diatomic
More informationMolecular Orbital Theory This means that the coefficients in the MO will not be the same!
Diatomic molecules: Heteronuclear molecules In heteronuclear diatomic molecules, the relative contribution of atomic orbitals to each MO is not equal. Some MO s will have more contribution from AO s on
More informationElectronic structure simulations of water solid interfaces
Electronic structure simulations of water solid interfaces Angelos Michaelides London Centre for Nanotechnology & Department of Chemistry, University College London www.chem.ucl.ac.uk/ice Main co-workers:
More informationInfrared spectra of small biomolecules from first-principle molecular dynamics simulations and effective normal mode analysis
Infrared spectra of small biomolecules from first-principle molecular dynamics simulations and effective normal mode analysis R. Vuilleumier, M.-P. Gaigeot and D. Borgis Département de chimie, Ecole Normale
More informationAN INTRODUCTION TO QUANTUM CHEMISTRY. Mark S. Gordon Iowa State University
AN INTRODUCTION TO QUANTUM CHEMISTRY Mark S. Gordon Iowa State University 1 OUTLINE Theoretical Background in Quantum Chemistry Overview of GAMESS Program Applications 2 QUANTUM CHEMISTRY In principle,
More informationIntroduction to Theories of Chemical Reactions. Graduate Course Seminar Beate Flemmig FHI
Introduction to Theories of Chemical Reactions Graduate Course Seminar Beate Flemmig FHI I. Overview What kind of reactions? gas phase / surface unimolecular / bimolecular thermal / photochemical What
More informationZwitterions in Neutral Arginine Clusters
Chapter 3 43 Cooperative Salt Bridge Stabilization of Gas Phase Zwitterions in Neutral Arginine Clusters Portions published previously in: Julian R. R.; Beauchamp J. L.; Goddard W. A. J. Phys. Chem. A
More informationStructure of diatomic molecules
Structure of diatomic molecules January 8, 00 1 Nature of molecules; energies of molecular motions Molecules are of course atoms that are held together by shared valence electrons. That is, most of each
More informationAb initio interpolated potential energy surface and classical reaction dynamics for HCO + +H, HOC + +H, and deuterated analogues
THE JOURNAL OF CHEMICAL PHYSICS 124, 124318 2006 Ab initio interpolated potential energy surface and classical reaction dynamics for HCO + +H, HOC + +H, and deuterated analogues Gloria E. Moyano, Seth
More informationSupplementary information
Matthias Heger, Tina Scharge, and Martin A. Suhm Institute of Physical Chemistry, Georg-August-Universität, Tammannstraße 6, 37077 Göttingen, Germany. E-mail: msuhm@gwdg.de Current address: Gesellschaft
More informationRethinking Hybridization
Rethinking Hybridization For more than 60 years, one of the most used concepts to come out of the valence bond model developed by Pauling was that of hybrid orbitals. The ideas of hybridization seemed
More informationNew Models for Aqueous Systems: Construction of Vibrational Wave Functions for use in Monte Carlo Simulations.
New Models for Aqueous Systems: Construction of Vibrational Wave Functions for use in Monte Carlo Simulations. Maria A. Gomez and Lawrence R. Pratt T-12 and CNLS Theoretical Division Los Alamos National
More informationFull configuration interaction potential energy curves for breaking bonds to hydrogen: An assessment of single-reference correlation methods
JOURNAL OF CHEMICAL PHYSICS VOLUME 118, NUMBER 4 22 JANUARY 2003 Full configuration interaction potential energy curves for breaking bonds to hydrogen: An assessment of single-reference correlation methods
More informationBe H. Delocalized Bonding. Localized Bonding. σ 2. σ 1. Two (sp-1s) Be-H σ bonds. The two σ bonding MO s in BeH 2. MO diagram for BeH 2
The Delocalized Approach to Bonding: The localized models for bonding we have examined (Lewis and VBT) assume that all electrons are restricted to specific bonds between atoms or in lone pairs. In contrast,
More informationVIBRATION-ROTATION SPECTRUM OF CO
Rice University Physics 332 VIBRATION-ROTATION SPECTRUM OF CO I. INTRODUCTION...2 II. THEORETICAL CONSIDERATIONS...3 III. MEASUREMENTS...8 IV. ANALYSIS...9 April 2011 I. Introduction Optical spectroscopy
More informationIntroduction to Hartree-Fock Molecular Orbital Theory
Introduction to Hartree-Fock Molecular Orbital Theory C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology Origins of Mathematical Modeling in Chemistry Plato (ca. 428-347
More informationCharacteristics of the interaction in azulene (H 2 X) n=1,2 (X=O,S) clusters.
Characteristics of the interaction in azulene (H 2 X) n=1,2 (X=O,S) clusters. Enrique M. Cabaleiro-Lago (a), Ángeles Peña-Gallego (b), Jesús Rodríguez-Otero (b), M. Merced Montero-Campillo (b) (a) Departamento
More information26 Group Theory Basics
26 Group Theory Basics 1. Reference: Group Theory and Quantum Mechanics by Michael Tinkham. 2. We said earlier that we will go looking for the set of operators that commute with the molecular Hamiltonian.
More informationV( x) = V( 0) + dv. V( x) = 1 2
Spectroscopy 1: rotational and vibrational spectra The vibrations of diatomic molecules Molecular vibrations Consider a typical potential energy curve for a diatomic molecule. In regions close to R e (at
More informationA dimer method for finding saddle points on high dimensional potential surfaces using only first derivatives
JOURNAL OF CHEMICAL PHYSICS VOLUME 111, NUMBER 15 15 OCTOBER 1999 A dimer method for finding saddle points on high dimensional potential surfaces using only first derivatives Graeme Henkelman and Hannes
More informationElectronic quantum effect on hydrogen bond geometry in. water dimer
Electronic quantum effect on hydrogen bond geometry in water dimer Danhui Li 1,2, Zhiyuan Zhang 1,2 Wanrun Jiang 1,2 Depeng Zhang 1,2 Yu Zhu 1,2 and Zhigang Wang 1,2* 1 Institute of Atomic and Molecular
More information5.61 Physical Chemistry Final Exam 12/16/09. MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Chemistry Chemistry Physical Chemistry
5.6 Physical Chemistry Final Exam 2/6/09 MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Chemistry Chemistry - 5.6 Physical Chemistry Final Examination () PRINT your name on the cover page. (2) It
More informationCalculate a rate given a species concentration change.
Kinetics Define a rate for a given process. Change in concentration of a reagent with time. A rate is always positive, and is usually referred to with only magnitude (i.e. no sign) Reaction rates can be
More informationComputational Material Science Part II. Ito Chao ( ) Institute of Chemistry Academia Sinica
Computational Material Science Part II Ito Chao ( ) Institute of Chemistry Academia Sinica Ab Initio Implementations of Hartree-Fock Molecular Orbital Theory Fundamental assumption of HF theory: each electron
More informationPublicat ion I. Publication I
Publicat ion I I Publication I K. J. Jalkanen, V. W. Jurgensen and I. M. Degtyarenko, Linear response properties required to simulate vibrational spectra of biomolecules in various media: (R)-Phenyloxirane
More informationVSCF/VCI vibrational spectroscopy of H 7 O + 3 and H 9O + 4 using high-level, many-body potential energy surface and dipole moment surfaces.
VSCF/VCI vibrational spectroscopy of H 7 O + 3 and H 9O + 4 using high-level, many-body potential energy surface and dipole moment surfaces Qi Yu 1 and Joel M. Bowman 1, 1 Department of Chemistry and Cherry
More informationWolfgang Demtroder. Molecular Physics. Theoretical Principles and Experimental Methods WILEY- VCH. WILEY-VCH Verlag GmbH & Co.
Wolfgang Demtroder Molecular Physics Theoretical Principles and Experimental Methods WILEY- VCH WILEY-VCH Verlag GmbH & Co. KGaA v Preface xiii 1 Introduction 1 1.1 Short Historical Overview 2 1.2 Molecular
More informationExcited States Calculations for Protonated PAHs
52 Chapter 3 Excited States Calculations for Protonated PAHs 3.1 Introduction Protonated PAHs are closed shell ions. Their electronic structure should therefore be similar to that of neutral PAHs, but
More informationForbidden Electric Dipole Transitions in the Hydrogen Molecular Ion First Estimates
Bulg. J. Phys. 44 (2017) 76 83 Forbidden Electric Dipole Transitions in the Hydrogen Molecular Ion First Estimates P. Danev Institute for Nuclear Research and Nuclear Energy, Bulgarian Academy of Sciences,
More informationINPUT DESCRIPTION FOR SQM version 2.0
INPUT DESCRIPTION FOR SQM version 2.0 INTRODUCTION SQM is an add-on module for the PQS program which scales force constants to produce a Scaled Quantum Mechanical (SQM) Force Field. This can correct for
More informationTHE JOURNAL OF CHEMICAL PHYSICS 123,
THE JOURNAL OF CHEMICAL PHYSICS 23, 2240 2005 Efficient methods for finding transition states in chemical reactions: Comparison of improved dimer method and partitioned rational function optimization method
More informationSUPPLEMENTARY INFORMATION
DOI: 10.1038/NCHEM.1680 On the nature and origin of dicationic, charge-separated species formed in liquid water on X-ray irradiation Stephan Thürmer, 1 Milan Ončák, 2 Niklas Ottosson, 3 Robert Seidel,
More informationIsomerization of the Newly Discovered Interstellar Molecule SiCN to SiNC Through Two Transition States
Liberty University DigitalCommons@Liberty University Faculty Publications and Presentations Department of Biology and Chemistry 2003 Isomerization of the Newly Discovered Interstellar Molecule SiCN to
More informationINPUT DESCRIPTION FOR SQM version 1.0
INPUT DESCRIPTION FOR SQM version 1.0 INTRODUCTION SQM is an add-on module for the PQS program which scales force constants to produce a Scaled Quantum Mechanical (SQM) Force Field. This can correct for
More informationAb initio calculations for potential energy surfaces. D. Talbi GRAAL- Montpellier
Ab initio calculations for potential energy surfaces D. Talbi GRAAL- Montpellier A theoretical study of a reaction is a two step process I-Electronic calculations : techniques of quantum chemistry potential
More informationMP2 Basis Set Limit Binding Energy Estimates of Hydrogen-bonded Complexes from Extrapolation-oriented Basis Sets
386 Bull. Korean Chem. Soc. 2007, Vol. 28, No. 3 Young Choon Park and Jae Shin Lee MP2 Basis Set Limit Binding Energy Estimates of Hydrogen-bonded Complexes from Extrapolation-oriented Basis Sets Young
More information