Journal of Chemical Physics
Feedback | Help | Exit
The Journal of Chemical Physics
Search:
   
 
 
 
Previous Article
Single-molecule fluorescence autocorrelation experiments on pentacene: The dependence of intersystem crossing on isotopic composition
Single pentacene molecules containing 13C or 1H in a pentacene-d14 doped p-terphenyl crystal have been studied by fluorescence autocorrelation. The triplet dynamics has been analyzed and a systematic ...
Next Article
Pressure-induced phase transition in mixed-valence gold complexes Cs2Au2X6 (X = Cl and Br)
Effects of hydrostatic pressure on Raman spectra have been investigated in mixed-valence (MV) gold complexes Cs2Au2X6 (X = Cl and Br). We have observed pressure-induced deactivation of the Raman-activ...

Regularity in chaotic reaction paths. I. Ar6

J. Chem. Phys. 110, 9160 (1999); doi:10.1063/1.478838

Issue Date: 8 May 1999 | See: Erratum

You are not logged in to this journal. Log in

Tamiki Komatsuzaki and R. Stephen Berry
Department of Chemistry, The University of Chicago, Chicago, Illinois 60637
We scrutinize the saddle crossings of a simple cluster of six atoms to show (a) that it is possible to choose a coordinate system in which the transmission coefficient for the classical reaction path is unity at all energies up to a moderately high energy, above which the transition state is chaotic; (b) that at energies just more than sufficient to allow passage across the saddle, all or almost all the degrees of freedom of the system are essentially regular in the region of the transition state; and (c) that the degree of freedom associated with the reaction coordinate remains essentially regular through the region of the transition state, even to moderately high energies. Microcanonical molecular dynamics simulation of Ar6 bound by pairwise Lennard-Jones potentials reveals the mechanics of passage. We use Lie canonical perturbation theory to construct the nonlinear transformation to a hyperbolic coordinate system which reveals these regularities. This transform "rotates away" the recrossings and nonregular behavior, especially of the motion along the reaction coordinate, leaving a coordinate and a corresponding dividing surface in phase space which minimize recrossings and mode–mode mixing in the transition state region. The action associated with the reactive mode tends to be an approximate invariant of motion through the saddle crossings throughout a relatively wide range of energy. Only at very low energies just above the saddle could any other approximate invariants of motion be found for the other, nonreactive modes. No such local invariants appeared at energies at which the modes are all chaotic and coupled to one another. ©1999 American Institute of Physics.
History: Received 27 October 1998; accepted 16 February 1999
Permalink: http://link.aip.org/link/?JCPSA6/110/9160/1
BUY THIS ARTICLE   (US$24)
Download HTML Download Sectioned HTML Download PDF (344 kB) View Cart

ERRATUM

  1. Erratum: "Regularity in chaotic reaction paths. I. Ar6" [J. Chem. Phys. 110, 9160 (1999)]
    Tamiki Komatsuzaki et al.
    J. Chem. Phys. 116, 862 (2002)

EPAPS

KEYWORDS and PACS

Keywords
PACS
  • 36.40.Jn
    Studies of special atoms, molecules, and their ions; clusters Atomic and molecular clusters Reactivity of clusters
  • 82.20.Db
    Physical chemistry Chemical kinetics Statistical theories (including transition state)
  • YEAR: 1999

RELATED DATABASES


To view database links for this article,
you need to log in.
To view database links for this article,
you need to log in.

PUBLICATION DATA

ISSN:
0021-9606 (print)   1089-7690 (online)
Publisher:
AIP is a member of CrossRef AIP

REFERENCES (55)
View in Separate Window/Tab
For access to fully linked references, you need to log in. For access to fully linked references, you need to Log in.
  1. C. Amitrano and R. S. Berry, Phys. Rev. Lett. 68, 729 (1992);
    2. Phys. Rev. E 47, 3158 (1993).
  2. T. L. Beck, D. M. Leitner, and R. S. Berry, J. Chem. Phys. 89, 1681 (1988).
  3. R. J. Hinde and R. S. Berry, J. Chem. Phys. 99, 2942 (1993).
  4. R. S. Berry, Chem. Rev. 93, 237 (1993).
  5. R. S. Berry, Int. J. Quantum Chem. 58, 657 (1996).
  6. E. R. Lovejoy, S. K. Kim, and C. B. Moore, Science 256, 1541 (1992).
  7. R. A. Marcus, Science 256, 1523 (1992).
  8. I. Ohmine, J. Phys. Chem. 99, 6767 (1995).
  9. R. B. Shirts and W. P. Reinhardt, J. Chem. Phys. 77, 5204 (1982);
    10. 79, 3173 (1983).
  10. M. J. Davis and S. K. Gray, J. Chem. Phys. 84, 5389 (1986).
  11. S. Wiggins, Normally Hyperbolic Invariant Manifolds in Dynamical Systems (Springer, New York, 1991).
  12. R. E. Gillian and G. S. Ezra, J. Chem. Phys. 94, 2648 (1991).
  13. M. Toda, Phys. Rev. Lett. 74, 2670 (1995).
  14. M. Toda, Phys. Lett. A 227, 232 (1997).
  15. T. Komatsuzaki and M. Nagaoka, J. Chem. Phys. 105, 10838 (1996).
  16. T. Komatsuzaki and M. Nagaoka, Chem. Phys. Lett. 265, 91 (1997).
  17. T. Komatsuzaki and R. S. Berry (unpublished).
  18. A. J. Lichtenberg and M. A. Lieberman, Regular and Chaotic Dynamics, 2nd ed. (Springer, New York, 1992).
  19. G. Hori, Publ. Astron. Soc. Jpn. 18, 287 (1966);
    20. 19, 229 (1967).
  20. A. Deprit, Celest. Mech. 1, 12 (1969).
  21. A. J. Dragt and J. M. Finn, J. Math. Phys. 17, 2215 (1976);
    22. 20, 2649 (1979).
  22. J. R. Cary, Phys. Rep. 79, 130 (1981).
  23. L. E. Fried and G. S. Ezra, J. Chem. Phys. 86, 6270 (1987).
  24. L. E. Fried and G. S. Ezra, Comput. Phys. Commun. 51, 103 (1988).
  25. L. E. Fried and G. S. Ezra, J. Phys. Chem. 92, 3144 (1988).
  26. M. Page and J. W. McIver, Jr., J. Chem. Phys. 88, 922 (1988).
  27. At a nonstationary point there exist only three zero-frequency modes corresponding to the total translational motions and the normal coordinates of the infinitesimal rotational motions appear in all the terms including the zeroth-order Hamiltonian.
  28. The canonical perturbation theory can be applied to some nonzero total angular momentum systems by Dragt and Finn's two-step transformation—Ref. 21 The discussions will be given elsewhere.
  29. W. H. Miller, Faraday Discuss. Chem. Soc. 62, 40 (1977).
  30. T. Seideman and W. H. Miller, J. Chem. Phys. 95, 1768 (1991).
  31. R. Hernandez and W. H. Miller, Chem. Phys. Lett. 214, 129 (1993).
  32. H. Goldstein, Classical Mechanics (Addison–Wesley, Reading, MA, 1950).
  33. See AIP Document No. E-PAPS: EJCPSA6-110-509918for the full expansions of the transformed coordinates and momenta, in first and second order and in both saddle regions. E-PAPS document files may be retrieved free of charge from AIP's FTP server (http://www.aip.org/pubservs/paps.html) or from ftp.aip.org in the directory /epaps/. For further information: e-mail: paps@aip.org or fax: 516-576-2223. [EPAPS]
  34. The deltaV was determined empirically by checking whether it gives indistinguishable PES and LCPT transformations with smaller deltaV. For example, even with deltaV = 1 × 10–8, the total numbers of the three-, and four-body nonlinear couplings terms were less changed than with deltaV = 2 × 10–5, i.e., 107 and 380 for saddle I, and 189 and 700 for saddle II, respectively, the number of total terms in H-bar2 (saddle I) is the same with deltaV = 2 × 10–5, without any significant change in coefficients with each term.
  35. J. C. Keck, Discuss. Faraday Soc. 33, 173 (1962).
  36. J. B. Anderson, J. Chem. Phys. 58, 4684 (1973).
  37. J. P. Bergsma, J. R. Reimers, K. R. Wilson, and J. T. Hynes, J. Chem. Phys. 85, 5625 (1986).
  38. J. P. Bergsma, B. J. Gertner, K. R. Wilson, and J. T. Hynes, J. Chem. Phys. 86, 1356 (1987).
  39. B. J. Gertner, K. R. Wilson, and J. T. Hynes, J. Chem. Phys. 90, 3537 (1989).
  40. W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, Numerical Recipes, 2nd ed. (Cambridge University Press, New York, 1992).
  41. T. Komatsuzaki and R. S. Berry, Phys. Chem. Chem. Phys.1, 1387 (1999).
  42. Such a characteristic time can be estimated by a plateau time taup; taup = 1/(kappaMD(0) – kappaMD(tau))[integral]<sub>0</sub><sup>tau</sup> dt[kappaMD(t) – kappaMD(tau)] where tau is the simulation cutoff time in evaluating kappaMD(t). taup would be the 1/e relaxation time if kappaMD(t) were exponential (Ref. 38). The taup are estimated as 1.359(I), 1.322(II) (at E=0.1 epsilon), 1.041(I), 0.935(II) (0.5 epsilon), and 1.048(I), 1.133(II) (1.0 epsilon), where the number in parentheses denotes either saddle I or saddle II.
  43. D. J. Wales and R. S. Berry, J. Phys. B 24, L351 (1991).
  44. R. J. Hinde, R. S. Berry, and D. J. Wales, J. Chem. Phys. 96, 1376 (1992).
  45. The quantity omega-barF(p,q) might not necessarily maintain its approximate invariance even while the associated action J-barF does, because omega-barF depends on not only J-barF but also the other J-bark(k [not-equal] F) [see Eq. (10)]. The equation of motion of the kth mode itself [Eqs. (12) and (13)] is rigorously satisfied if the J-bark is conserved despite the constancy or nonconstancy of omega-bark(J-bar ). As seen in Fig. 7, in almost cases in this study the frequencies omega-bar1(p,q) were fairly well conserved during saddle crossings as were the actions J-bar1. Even if it were not so, we might still expect any nonconstancy of omega-barF to leave the separability of q-barF as unaffected as J-barF in the transition state. The reason is that to pass through a dividing surface from one side to the other, the crossing trajectories typically require, at most, a half-period of the reactive hyperbolic orbit ~ pi/omega-barF; the fluctuation of omega-barF should have little influence during such short time intervals. We are now analyzing this in detail.
  46. C. C. Martens, M. J. Davis, and G. S. Ezra, Chem. Phys. Lett. 142, 519 (1987).
  47. M. J. Davis, J. Chem. Phys. 83, 1016 (1985).
  48. D. G. Truhlar, B. C. Garrett, and S. J. Kippenstein, J. Phys. Chem. 100, 12771 (1996).
  49. D. C. Chatfield, R. S. Friedman, S. L. Mielke, G. C. Lynch, T. C. Allison, D. G. Truhlar, and D. W. Schwenke, in Dynamics of Molecules and Chemical Reactions, edited by R. E. Wyatt and J. Z. H. Zhang (Dekker, New York, 1996), p. 323.
  50. J. Villa and D. G. Truhlar, Theor. Chem. Acc. 97, 317 (1997).
  51. In some trajectories at moderately high energies it was found that H-bar0([bold p]-bar2nd,[bold q]-bar2nd) = [summation]<sub>k = 1</sub><sup>12</sup> omegakJ-bar<sub>k</sub><sup>2nd</sup>(p,q) tend to maintain an approximate constant of motion during saddle crossings in addition to J-bar<sub>1</sub><sup>2nd</sup>(p,q). It implies that while there might not exist any approximate constant of action for each bath DOF, the system could possess an additional approximate constancy in some portions of the phase space in the transition state.
  52. G. A. Baker Jr., Essentials of Padé Approximants (Academic, New York, 1975).
  53. D. Farrelly and T. Uzer, J. Chem. Phys. 85, 308 (1986).
  54. L. E. Fried and G. S. Ezra, J. Chem. Phys. 90, 6378 (1989).
  55. As an another approach altering the zeroth-order integrable Hamiltonian, see I. D. Feranchuk and L. I. Komarov, Phys. Lett. A 88, 211 (1982).

CITING ARTICLES
For access to citing articles, you need to log in.
For access to citing articles, you need to Log in.


Copyright © American Institute of Physics