vibrational averaging techniques to calculate the role of water dimers in atmospheric absorption...

Vibrational averaging techniques to calculate the role of water dimers in atmospheric absorption

Jonathan Tennyson, Matt J. Barber, Ross E. A. Kelly, Lorenzo Lodi

Department of Physics and Astronomy, University College London

CAVIAR meetingCambridge, Feb 2011

CAVIAR - Continuum Absorption at Visible and Infrared wavelengths and its Atmospheric Relevance

• (1) Need to solve the nuclear motion Hamiltonian– 12D problem! Approximations required.

• (2) Fully dimensional potential energy surface required– Huang, Braams and Bowman (HBB) potentials

• 30-40,000 configurations sampled.• Calculated at coupled-cluster, single and double and perturbative

treatment of triple excitations method.• Augmented, correlation consistent, polarized triple zeta basis set.• Polynomial fit with 5227 coefficients.

Water Dimer Method

HBB – X. Huang et al. J. Chem. Phys. 128, 034312 (2008).HBB2 – X. Huang et al. J. Chem. Phys. 130, 144314 (2009).

(2) Fully Dimensional Water Dimer Potential

Monomer corrected* HBB potential • Corrects for monomer excitation

– Accurate modes for the monomer

* S. V. Shirin et al., J. Chem. Phys. 128, 224306 (2008).R.E.A. Kelly, J. Tennyson, G C. Groenenboom, A. Van der Avoird, JQRST, 111, 1043 (2010).

Solving the 6D intermolecular problem

• Brocks Brocks et alet al. Hamiltonian*. Hamiltonian*

Solved using Dimer code ofSolved using Dimer code of

Groenboom and van der Avoird:Groenboom and van der Avoird:

Gives dimer “VRT” statesGives dimer “VRT” states

* G. Brocks et al. Mol. Phys. 50, 1025 (1983).

Solving the 6D intermolecular problemSolving the 6D intermolecular problem

AcceptorTwist (AT)

AcceptorWag (AW)

Donor Torsion (DT)

In PlaneBend (IPB) Stretch

Out-of-PlaneBend (OPB)

Generated by Matt Hodges and Anthony Stone. C. Millot et al. J. Phys. Chem. A 1998,102, 754. http://www-stone.ch.cam.ac.uk/research/water.dimer/modes.html

Adiabatic Separation

• Approximate separation between monomer and dimer modes– Separate intermolecular and intramolecular modes.

mD – water donor vibrational wavefunction

mA – water acceptor vibrational wavefunction

d – dimer VRT wavefunction

dmm AD

)()(|);,(|)()()( BBAABABBAAmm

eff mmVmmV BA QQRQQQQR

• Now we can vibrationally average the potentialNow we can vibrationally average the potential

• Input for 6D calculationsInput for 6D calculations

donordonor acceptoracceptor

State mState m State nState n

• How well does it perform for |0 0> calculationsHow well does it perform for |0 0> calculations

Solving the 6D intermolecular problem

• 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 R.E.A. Kelly, J. Tennyson, G C. Groenenboom, A. Van der Avoird, JQRST, 111, 1043 (2010).

Vibrational Averaging

Vibrational Averaging: 6D Costs!

• Computation:

– typical number of DVR points with different Morse Parameters:

– {9,9,24} gives 1,080 points for monomer

– 1,0802 = 1,166,400 points for both monomers

– 1,166,400 x 2,894,301 intermolecular points

= 3,374,862,926,400 points• Same monomer wavefunctions for all grid points• Distributed computing: Condor 1000 computers, 10 days

But we have a way to probe high frequency dimer spectra

Full model for high frequency absorption

• Approximate separation between monomer and dimer modes

• Franck-Condon approximation for vibrational fine structure

• Rotational band model

Adiabatic Separation

• Approximate separation between monomer and dimer modes– Separate intermolecular and intramolecular modes.

mD – water donor vibrational wavefunction

mA – water acceptor vibrational wavefunction

d – dimer VRT wavefunction

dmm AD

Model for high frequency absorption

• Approximate separation between monomer and dimer modes

• Franck-Condon approximation for vibrational fine structure

• Rotational band model

2

2121

2fffiii

fi dmmdmmI

22

1122

fifi

mmddmmfi

Franck-Condon Approx for overtone spectra

Assume monomer m1 excited, m2 frozen

m2i = m2

f

I

(2) Franck-Condon factor

(square of overlap integral):

Gives dimer vibrational fine structure

(1) Monomer vibrational band Intensity

Allowed Transitions in our Model

1. Donor2. Acceptor

All transitions from ground monomer vibrational states

Assume excitation localised on one monomer

Franck-Condon factors

– Overlap between dimer states on adiabatic potential energy surfaces for water monomer initial and final states

– Need the dimer states (based on this model).

Transitions: Example

Donor – Vibrational ground state (VGS)Acceptor – VGS

Donor –VGS Acceptor – bend

Acceptor Twist (AT)

Acceptor Wag (AW)

Donor Torsion (DT)

Ground State (GS)

Acceptor Twist (AT)

Acceptor Wag (AW)

Donor Torsion (DT)

Ground State (GS)

Transitions: Example

Acceptor Twist (AT)

Acceptor Wag (AW)

Donor Torsion (DT)

Ground State (GS)

Acceptor Twist (AT)

Acceptor Wag (AW)

Donor Torsion (DT)

Ground State (GS)

Donor – Vibrational ground state (VGS)Acceptor – VGS

Donor –VGS Acceptor – bend

Transitions: Example

Acceptor Twist (AT)

Acceptor Wag (AW)

Donor Torsion (DT)

Ground State (GS)

Acceptor Twist (AT)

Acceptor Wag (AW)

Donor Torsion (DT)

Ground State (GS)

Donor – Vibrational ground state (VGS)Acceptor – VGS

Donor –VGS Acceptor – bend

Outline of full problem

• Need to ultimately solve (6D problem)

• H=K+Veff

• Veff sampled on a 6D grid

dd EH • Calculate states for donor

• Calculate states for acceptor

• Vibrationally average potential for each state-state combination– Really only |0j> and |i0>

(a) 6D averaging:

(b) 3D+3D averaging:

1 C Leforestier et al, J Chem Phys, 117, 8710 (2002)2 R. E. A. Kelly et al. To submit shortly.

);,()()(

)()(|);,(|)()()(22 RQQQQ

QQRQQQQR

BABBAqq

A

BBAABABBAAeff

Vmm

mmVmmV

BA

);,()(|);,(|)(

)(|);,(|)()(000

0

RQQQRQQQ

QRQQQR

BABBBABB

AABAAAeff

VmVm

mVmV

Averaging Techniques

Averaging Techniques

• Form of the wavefunction:– (I) Uncoupled free monomer

– (II) Uncoupled perturbed (fixed) monomer

R. E. A. Kelly et al. To submit shortly.

Problems with Fixed Wavefunction approach (uncoupled methods)

• Donor bend • (Donor) Free OH stretch • (Donor) Bound OH stretch

• (Donor) Free OH stretch • (Donor) Bound OH stretch

Averaging Techniques

• Form of the wavefunction:– (I) Uncoupled free monomer

– (II) Uncoupled perturbed (fixed) monomer

– (III) Coupled Adiabatic

R. E. A. Kelly et al. To submit shortly.

Averaging Techniques

• Form of the wavefunction:– (I) Uncoupled free monomer

– (II) Uncoupled perturbed (fixed) monomer

– (III) Coupled Adiabatic

R. E. A. Kelly et al. To submit shortly.

Averaging Techniques

• Form of the wavefunction:– (I) Uncoupled free monomer

– (II) Uncoupled perturbed (fixed) monomer

– (III) Coupled Adiabatic

R. E. A. Kelly et al. To submit shortly.

Averaging Techniques

• Form of the wavefunction:– (I) Uncoupled free monomer– (II) Uncoupled perturbed monomer – (III) Coupled Adiabatic

• Coupled Adiabatic methods are the most suitable– Requires wavefunction calculations at each intermolecular grid

point! 2,893,401 * 2 DVR3D calculations!– So we use cheaper (3+3)D averaging technique.

– Still costs! 500-700 CPUs for 3-4 weeks.

This part is complete

.

Calculating dimer spectra with FC approach

• Solved for monomers • Coupled adiabatic appoach

• Vibrationally averaged potential for donor-acceptor state-state combinations |0j> and |i0>• Input for 6D intermolecular problem

• Now we can solve 6D intermolecular problem• Obtain vibrational fine structure

For T=296 K requires all vibrational states

1 15 5

2 26 6

4

4

3

3• G16 Symmetry of Hamiltonian for GS monomers– > replaced with G4

• Greatly increases computational requirements

So far• Reduced angular basis• Small radial basis• 320 diagonalizations for 0-10,000 cm-1

• Each at 16 GB• 8 states per symmetry block

• Gives 20,480 transitions: results presented by Matt

Solving the 6D intermolecular problem:Allowed permutations for excited monomers

Model for high frequency absorption

• Approximate separation between monomer and dimer modes

• Franck-Condon approximation for vibrational fine structure

• Rotational band model

• G4 Symmetry of Hamiltonian• Require all vibrational bound states • Greatly increases computational requirements

So far• Reduced angular basis• Small radial basis• 320 diagonalizations• Each at matrix 320,000 x 320,000 means ~ 1 TB RAM• ~400 states per symmetry block• 3 – 5 days CPU• Ongoing (problems with UCL computers)

Solving the 6D intermolecular problem: For atmospheric temperatures

• New Model to probe near IR and visible regions of the water dimer spectra.– vibrational fine structure

• Aim for spectra for up to 15,000 cm-1 produced.

• And all states up to dissociation to be calculated.

-- Computer resources a big issue

Conclusions