towards a spectroscopically flexible water dimer potential energy surface ross e. a. kelly, and...

Towards a Spectroscopically Towards a Spectroscopically Flexible Water Dimer Potential Flexible Water Dimer Potential

Energy SurfaceEnergy Surface

Ross E. A. KellyRoss E. A. Kelly, and Jonathan Tennyson, and Jonathan TennysonDepartment of Physics & AstronomyDepartment of Physics & Astronomy

University College LondonUniversity College London

Gerrit C. Groenenboom and Ad van der Avoird Gerrit C. Groenenboom and Ad van der Avoird Theoretical Chemistry, Institute for Molecules and Theoretical Chemistry, Institute for Molecules and

Materials, Materials, Radboud University, Nijmegen.Radboud University, Nijmegen.

Imperial College, December 2008

OutlineOutline

I. MotivationsI. Motivations II. The Water DimerII. The Water Dimer III. Vibration-Rotation TunnellingIII. Vibration-Rotation Tunnelling IV. Dimer PotentialsIV. Dimer Potentials V. Theoretical MethodV. Theoretical Method VI. VRT StatesVI. VRT States VII. Monomer Band OriginsVII. Monomer Band Origins VIII. Monomer Corrected SurfaceVIII. Monomer Corrected Surface IX. Conclusions & Further WorkIX. Conclusions & Further Work

I. MotivationsI. Motivations to understand to understand water dimer absorption water dimer absorption

throughout visible and IR region in the throughout visible and IR region in the atmosphere.atmosphere.

To create a To create a high accuracy water dimer high accuracy water dimer spectraspectra in agreement with in agreement with experimentsexperiments..

To create a To create a linelistlinelist of all possible of all possible water water dimer transitionsdimer transitions..

• To get all Vibration Rotation Energy Levels To get all Vibration Rotation Energy Levels Calculations require:Calculations require:

1. Some theoretical methodology for solving 1. Some theoretical methodology for solving the Hamiltonianthe Hamiltonian 2. High accuracy Potential Energy Surface2. High accuracy Potential Energy Surface

Preferably fully flexible (12D)Preferably fully flexible (12D)

•4 O-H distances4 O-H distances•2 H-O-H Angles2 H-O-H Angles•O-O distanceO-O distance•5 Euler Angles5 Euler Angles

II. Water DimerII. Water Dimer

Extremely Complex LandscapeExtremely Complex Landscape

X. Huang, B. J. Braams, J. M. Bowman, J. Phys. Chem. A 110, 445 (2006).

II. Water DimerII. Water Dimer

Complicated further by tunnelling effects!Complicated further by tunnelling effects! Tunnelling between eqiulavent states in the PES is Tunnelling between eqiulavent states in the PES is

feasible!feasible! Acceptor Tunnelling:Acceptor Tunnelling:

No bond breaking hereNo bond breaking hereLowest tunnelling barrier Lowest tunnelling barrier

Also, by breaking the Hydrogen bond, other tunnelling Also, by breaking the Hydrogen bond, other tunnelling paths possible: paths possible: Donor-Acceptor interchangeDonor-Acceptor interchangeDonor Bifurcation TunnellingDonor Bifurcation Tunnelling

III. Vibration-Rotation TunnellingIII. Vibration-Rotation Tunnelling

Can be represented by Permutation-Inversion Group GCan be represented by Permutation-Inversion Group G1616..

1 1

1 1

5 5

5 5

2 2

2 2

6 6

6 6

6 6

6 6

5

5 5

5

4

4

4

4

3

3

3

3

3 3

3 3

4

4

4

4

1 1

1 1

2

2 2

2

Isomorphic to D4h

with Irreducible Elements:

A1

+, A2

+, A1

-, A2

-, B1

+, B2

+, B1

-, B2

-, E+, E-

-> Water Dimer Spectroscopic Labels

III. Labelling Water Dimer StatesIII. Labelling Water Dimer States

[1] R. Bukowski, K. Szalewicz, G. C. Groenenboom, A. van der Avoird, Science 315, p1249-1252 (2007).

[2] X. Huang, B. J. Braams, J. M. Bowman, J. Phys. Chem. A 110, 445 (2006).

[3] X. Huang, B. J. Braams, J. M. Bowman, R. E. A. Kelly, J. Tennyson, G. C. Groenenboom, A. van der Avoird, J. Chem. Phys. 128,

034312 (2008)

Benchmark:Benchmark:●6D6D CC-polCC-pol Potential [1] - Potential [1] - Rigid MonomersRigid Monomers

New (2) New (2) 12D PESs12D PESs by Huang, Braams & Bowman by Huang, Braams & Bowman [2,3][2,3]

●HBB0HBB0 [2] [2] ●and and HBBHBB [3] – [3] – NEW 12D PESNEW 12D PES

- All completely - All completely ab initioab initio- Good agreement with experiment- Good agreement with experiment

Many PES’s available!Many PES’s available!

IV. Dimer Potentials AvailableIV. Dimer Potentials Available

• 30 000 configurations. 30 000 configurations.

• Calculated at coupled-cluster, single and double and Calculated at coupled-cluster, single and double and

perturbative treatment of triple excitations method.perturbative treatment of triple excitations method.

• augmented, correlation consistent, polarized triple zeta augmented, correlation consistent, polarized triple zeta

basis set.basis set.

• Polynomial fit with 5227 coefficients.Polynomial fit with 5227 coefficients.

• However, compared to However, compared to CCpolCCpol potential (benchmark): potential (benchmark):● Lower relative grid coverage of the surface than Lower relative grid coverage of the surface than

benchmark 6D PES.benchmark 6D PES.● Not extrapolated to the complete basis set limit Not extrapolated to the complete basis set limit

(CBS)(CBS)● No bond functions in analytical fitNo bond functions in analytical fit● Dissociation less accurately describedDissociation less accurately described

12DIV. IV. HBBHBB PES for H PES for H

44OO22

[1] X. Huang, B. J. Braams, J. M. Bowman, R. E. A. Kelly, J. Tennyson, G. C. Groenenboom, A. van der Avoird, J. Chem. Phys. 128, 034312

(2008).

• Dissociation Energy = 1665.82 cmDissociation Energy = 1665.82 cm-1-1 / 1757.90 cm / 1757.90 cm-1-1

• Other tests? Other tests?

• Compare with Low temperature high-resolution Compare with Low temperature high-resolution

Tetrahertz Spectroscopy (prepared in supersonic Tetrahertz Spectroscopy (prepared in supersonic

molecular beams), around 5 K.molecular beams), around 5 K.

• How can this be done theoretically?How can this be done theoretically?

IV. IV. HBBHBB PES for H PES for H44OO22

Rigid monomer Hamiltonian [1]:

Only for the Intermolecular modes

Used for water dimer previously, detailed account [2]

Coupled product of Symmetric rotor functions (Wigner-D functions) for the Angular coordinates

Radial basis: sinc Discrete Variable Representation (DVR)

[1] G. Brocks, A. van der Avoird, B. T. Sutcliffe, J. Tennyson, Mol. Phys. 50, 1025 (1983).

[2] G. C. Groenenboom, et al., JCP 113, 6702 (2000).

V. Theoretical Method for VRT levelsV. Theoretical Method for VRT levels

VI. Ground State VRT states for HVI. Ground State VRT states for H44OO22

[1] X. Huang, B. J. Braams, J. M. Bowman, R. E. A. Kelly, J. Tennyson, G. C. Groenenboom, A. van der Avoird, J. Chem. Phys. 128, 034312

(2008).

Very good agreement with:

Ground State Tunnelling splittings Rotational Constants

Not so good agreement with:

Acceptor Tunnelling

Excellent agreement with:

Ground State Tunnelling splittings Rotational Constants

Not so good agreement with:

Acceptor Tunnelling

[1] X. Huang, B. J. Braams, J. M. Bowman, R. E. A. Kelly, J. Tennyson, G. C. Groenenboom, A. van der Avoird, J. Chem. Phys. 128, 034312

(2008).

VI. Ground State VRT for DVI. Ground State VRT for D44OO22

In cm-1

Red – ab initio potential Black – experimental

GS – ground state

DT – donor torsion

AW – acceptor wag

AT – acceptor twist

DT2 – donor torsion overtone

VI. More low level VRT States for HVI. More low level VRT States for H44OO22

In cm-1

Red – ab initio potential Black – experimental

GS – ground state

DT – donor torsion

AW – acceptor wag

AT – acceptor twist

DT2 – donor torsion overtone

VI. More low level VRT States for DVI. More low level VRT States for D44OO22

VII. Monomer Band OriginsVII. Monomer Band Origins New 12D Huang New 12D Huang et al.et al. PES seems to work well: PES seems to work well:

– for low-level dimer VRT statesfor low-level dimer VRT states How about for H2O monomer energy levels?How about for H2O monomer energy levels? Use DVR3D [1] for Water monomer levels:Use DVR3D [1] for Water monomer levels:

ExciteExciteFixFix100 bohr

1. J. Tennyson et al., Comp. Phys. Comm. 2004, 163, 85-116.

Comparison is not so good.Comparison is not so good.

VII. Monomer Band OriginsVII. Monomer Band Origins

Blue HBB dimerBlue HBB dimerGreen HBB0Green HBB0Red Shirin 2008.Red Shirin 2008.(Monomer).(Monomer).

VIII. Adding Monomer CorrectionVIII. Adding Monomer Correction

)],,(),,([)](2)()([),( 21212121 eeecorr VVVVVV qqqqqqqqq

Correction for monomer modes:Correction for monomer modes:

New Potential Expression:New Potential Expression:

),(),,(),,(),,( 21212121 qqqqrqqrqq corrcorr VVVV

Tests for PotentialTests for Potential Evaluation of the saddle points.Evaluation of the saddle points. Evaluation of the monomer & dimer VRT states.Evaluation of the monomer & dimer VRT states.

VIII. Monomer Corrected SurfaceVIII. Monomer Corrected SurfaceCharacteristicCharacteristic HBBHBB HBB+MCHBB+MC

DDee -1657.12-1657.12 -1665.82-1665.82

Point 2Point 2 198.47198.47 193.46193.46

Point 3Point 3 243.98243.98 241.09241.09

Point 4Point 4 340.60340.60 358.06358.06

Point 5Point 5 1178.921178.92 1181.761181.76

Point 6Point 6 588.19588.19 590.19590.19

Point 7Point 7 894.83894.83 898.27898.27

Not Particularly worse than HBB potential.

Still a very good Still a very good agreement with agreement with exp.exp.

Nothing is Nothing is changed changed significantly.significantly.

VIII. Ground VIII. Ground State VRT State VRT

levels for Hlevels for H44OO22

VIII. More VIII. More VRT States VRT States for HBB+MCfor HBB+MC

surfacesurface

Again, Again, agreement is agreement is not not significantly significantly worse.worse.

IX. Conclusions & Further WorkIX. Conclusions & Further Work

Introduction of monomer correction makes PES Introduction of monomer correction makes PES more transferable for spectroscopic purposes.more transferable for spectroscopic purposes.

Little effect on dimer characteristics.Little effect on dimer characteristics. Monomer band origins significantly improved.Monomer band origins significantly improved. We are working towards a model which We are working towards a model which

incorporates monomer excitations into the dimer incorporates monomer excitations into the dimer states, so spectra can be produced.states, so spectra can be produced.