slide1
Download
Skip this Video
Download Presentation
A global ab initio -based potential energy surface for of H 5 + , vibrational zero-point, and reaction dynamics of H 3 + + HD

Loading in 2 Seconds...

play fullscreen
1 / 34

A global ab initio-based potential energy surface for of H5, vibrational zero-point, and reaction dynamics of H3 HD - PowerPoint PPT Presentation


  • 113 Views
  • Uploaded on

A global ab initio -based potential energy surface for of H 5 + , vibrational zero-point, and reaction dynamics of H 3 + + HD. JMB, Bas Braams, Zhen Xie, Emory University, . National Science Foundation, Office of Naval Research. Issues in Astrophysics. Deuterium fraction Spectroscopy

loader
I am the owner, or an agent authorized to act on behalf of the owner, of the copyrighted work described.
capcha
Download Presentation

PowerPoint Slideshow about 'A global ab initio-based potential energy surface for of H5, vibrational zero-point, and reaction dynamics of H3 HD' - Audrey


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.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.


- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -
Presentation Transcript
slide1

A global ab initio-based potential energy

surface for of H5+, vibrational zero-point,

and reaction dynamics of H3+ + HD

JMB, Bas Braams, Zhen Xie,

Emory University,

NationalScience Foundation, Office of Naval Research

slide2

Issues in Astrophysics

  • Deuterium fraction
  • Spectroscopy
  • Thermochemistry
  • Reaction rates

D. Gerlich, E. Herbst, and E. Roueff, Planet Space Sci. 50, 1275 (2002)

D. Gerlich, S. Schlemmer, , Planet. Space Science 50, 1287(2002)

slide3

Challenges for Theoretical Chemistry

  • Potential energy surfaces
  • Spectroscopy
  • Thermochemistry/ZPE
  • Dynamics
slide4

Previous Theoretical Work

Ab initio electronic structure

R. Prosmiti, P. Villarreal, G. Delgado-barrio, (2001,2003) - DIM-fit PES

H. Müller and W. Kutzelnigg, PCCP (2000) - CC-R12, stationary pts, De

Y. Yamaguchi, J. F. Gaw, R. B. Remington, H. F. Schaefer III, (1987) -

Dynamics

I.Stich, D.Marx, M.Parinello, T.Terakura (1997) - PI/DFT

W.P.Kraemer, V.Spirko, and O.Bludsk (1994) - Red Dim vibrations

G. Moyano and M. Collins, (2003) - “GROW” QCT k(T)

slide5

H5+ Potential Landscape

Stationary Points

slide7

H5+ Potential Energy Surface (2005)

  • Ca 105ab initio energies [CCSD(T)/aug-cc-pVTZ]
  • Fit with a polynomial basis the is invariant wrt
  • any permutation of the 5 H atoms.
  • Many-body representation that gives the
  • fragments H3+ and H2
  • Use a switching function to describe
  • long-range electrostatic interaction.
  • Z. Xie, B. J. Braams, and J. M. Bowman J. Chem. Phys. 122, 224307 (2005).
slide8

3

4

6

2

1

5

Permutational Invariance: H2CO Example

H

4

H’

3

6

5

1

2

H

H’

slide11

H5+ Potential Energy Surface (2005)

a Muller &Kutzelnigg, b Prosmiti et al. eMoyano & Collins

slide12

Spectroscopy and Thermochemistry

HD + H3+  [H4D+]H2+ + H2D+

Thermochemistry - get the ZPEs accurately. HO?

Spectroscopy - detect reactant and products and maybe stabilizedH4D+

M. Okumura, L. I. Yeh, and Y. T. Lee (1988) - H5+

slide13

Thermochemistry

H5+ H3+ + H2 , D0= 6.9±0.3 kcal/mole

Recall D0 = De+ DZPE

ZPE H3++H2

ZPE H5+

De

slide14

Thermochemistry

H5+ H3+ + H2 , D0= 6.8±0.3 kcal/mol

  • The most accurate De = 8.58 kcal/mol
  • The present PES De = 8.30 kcal/mol
  • Using HO ZPEs PES gets D0 = 5.57 kcal/mol
  • (Not in good agreement with experiment)
slide15

Quantum Diffusion Monte Carlo ZPEs

Solution to the TDSE

E0 is the exact ZPE

Let t = it

slide16

H5+ H3+ + H2

H3+ + H2

H5+

ZPE H5+ = 20.6 kcal/mol(= 2.48De)

(Application of DMC requires a global PES)

slide17

H5+ H3+ + H2

From the DMC ZPEs we get D0 = 6.33±.03 kcal/mol

Exp* = 6.8±0.3 (HO = 5.57)

*Exp1 = 6.6 ±0.3,Exp2 = 7.0 ±0.1

Present Deis low by 0.28 compared

to most accurate ab intio value

So, our estimate would be 6.6 kcal/mol

slide18

More DMC-based Energetics

H4D+ H3+ + HD D0 = 6.63 kcal/mol

 H2D+ + H2 D0 = 6.37 kcal/mol

H3+ + HD  H2D+ + H2 , E = -90 cm-1 (-0.26 kcal/mol)

slide19

Dynamics and Reaction Rates

HD + H3+ H2+ + H2D+

slide20

HD + H3+ H4D+H2+ + H2D+

Structure of the Global Minimum

Three distinct localized locations for D

slide21

HD + H3+ H4D+H2+ + H2D+

Structure of the Global Minimum

Three distinct localized locations for D

slide23

HD + H3+  H4D+H2 + H2D+

( H’D + H3+)

Statistical expectations

Always form a “collision complex” - Langevin xsection

Classically 60% to H2 + H2D+, P(b) = 0.6

and40% to H’D+ H3+ withP(b) = 0.3 for H’≠H

P(b) = 0.1 H’=H

slide24

HD + H3+  H4D+H2 + H2D+

( H’D + H3+)

Virtually no exchange, P(b) 0.05

slide26

A sample direct trajectory

H3+(J=1)+HD, Ecoll = 100 cm-1

Ecoll= 95 K

slide29

Summary So Far

  • Global PES done - PI, dissociates, VLR ad hoc
  • DMC calcs of ZPEs and Do
  • QCT calculations of xsection and rate constant
  • (Not statistical, mostly proton hopping)

To Do

  • Quantum calculations of vibrational energies
  • and IR spectrum. (This is feasible with Multimode.)
  • Quantum calculation of the rate constant. (This is
  • very hard to do in full dimensionality.)
slide30

CH5+ CH3+ + H2

CH5+

Like H5+ there are 5! equivalent

Minima, saddle points.

Fit to 36173 CCSD(T)/aug-cc-pVTZ

Does CH5+ have a structure?

  • Oka’s unassigned spectrum
  • Marx and Parinello
  • Schaefer, Scheiner, Schlyer
  • Klopper and Kutzelnigg
slide32

A. Brown, B. J. Braams, K. Christoffel, Z. Jin,

and J. M. Bowman, J. Chem. Phys., 2003, 119, 8790. - fit and MD spectrum

A. Brown, A. B. McCoy, B. J. Braams, Z. Jin,

and J. M. Bowman, J. Chem. Phys., 2004, 121, 4105. - fit and DMC

A.B. McCoy, B .J. Braams, A. Brown

et al., J. Phys. Chem. A, 2004, 108, 4991. - DMC isotopmers

slide34

CH4D+ CH3+ + HD, CH2D+ + H2

Close to statistical and independent of initial geometry

ad