1 / 31

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 meeting Cambridge, Feb 2011

karah
Download Presentation

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

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. 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 meeting Cambridge, Feb 2011 CAVIAR - Continuum Absorption at Visible and Infrared wavelengths and its Atmospheric Relevance

  2. Water Dimer Method • (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. HBB – X. Huang et al. J. Chem. Phys. 128, 034312 (2008). HBB2 – X. Huang et al. J. Chem. Phys.130, 144314 (2009).

  3. (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).

  4. Solving the 6D intermolecular problem • Brocks et al. Hamiltonian* Solved using Dimer code of Groenboom and van derAvoird: Gives dimer “VRT” states * G. Brocks et al. Mol. Phys. 50, 1025 (1983).

  5. Acceptor Twist (AT) Acceptor Wag (AW) Donor Torsion (DT) Solving the 6D intermolecular problem Solving the 6D intermolecular problem In Plane Bend (IPB) Out-of-Plane Bend (OPB) Stretch 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

  6. 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

  7. Solving the 6D intermolecular problem • Now we can vibrationally average the potential • Input for 6D calculations • donor acceptor • State m State n • How well does it perform for |0 0> calculations

  8. Vibrational Averaging • 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).

  9. 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

  10. Full model for high frequency absorption Approximate separation between monomer and dimer modes Franck-Condon approximation for vibrational fine structure Rotational band model

  11. 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

  12. Model for high frequency absorption Approximate separation between monomer and dimer modes Franck-Condon approximation for vibrational fine structure Rotational band model

  13. Franck-Condon Approx for overtone spectra Assume monomer m1 excited, m2 frozen m2i = m2f I a (2) Franck-Condon factor (square of overlap integral): Gives dimer vibrational fine structure (1) Monomer vibrational band Intensity

  14. Allowed Transitions in our Model Assume excitation localised on one monomer 2. Acceptor 1. Donor All transitions from ground monomer vibrational states

  15. 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).

  16. 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)

  17. 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)

  18. 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)

  19. Outline of full problem • Calculate states for donor • Calculate states for acceptor • Vibrationally average potential for each state-state combination • Really only |0j> and |i0> • Need to ultimately solve (6D problem) • H=K+Veff • Veff sampled on a 6D grid

  20. (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. Averaging Techniques

  21. Averaging Techniques • Form of the wavefunction: • (I) Uncoupled free monomer • (II) Uncoupled perturbed (fixed) monomer • R. E. A. Kelly et al. To submit shortly.

  22. 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

  23. 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.

  24. 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.

  25. 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.

  26. 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 • .

  27. 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 6Dintermolecular problem • Obtain vibrational fine structure • For T=296 K requires allvibrational states

  28. Solving the 6D intermolecular problem:Allowed permutations for excited monomers 2 2 6 6 4 3 1 1 5 5 3 4 • 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

  29. Model for high frequency absorption Approximate separation between monomer and dimer modes Franck-Condon approximation for vibrational fine structure Rotational band model

  30. Solving the 6D intermolecular problem: For atmospheric temperatures • G4 Symmetry of Hamiltonian • Require allvibrational 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)

  31. Conclusions • 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

More Related