Lipkin, N. ; Gerber, R. B. ; Moiseyev, N. ; Nathanson, G. M. ATOM SCATTERING STUDIOS OF LIQUID STRUCTURE AND DYNAMICS - COLLISIONS OF XE WITH A MODEL OF SQUALANE.
JOURNAL OF CHEMICAL PHYSICS 1994,
100, 8408-8417.
AbstractMolecular dynamics (MD) computer simulations are carried out for scattering of high-energy Xe atoms off liquid squalane, and the results are compared with those of molecular-beam scattering. experiments. A crude model for squalane is adopted, describing the hydrocarbon chain molecule as a sphere, and ignoring the role of internal modes. Good overall agreement is found between the results of the simulations and experiment, both for angular distributions and for trends in energy I transfer properties. In particular, excellent agreement is obtained for the dependence of the energy transfer on the deflection angle for in-plane scattering. Theory predicts less trapping events than found experimentally, probably due to the crude model adopted for the squalane molecules. The partial success of the model in predicting some properties and not others is discussed. The other main conclusions of the study are (1) The instantaneous local structure of the liquid surface is highly corrugated, giving rise to a broad annular distribution and to extensive out-of-plane scattering. (2) High-energy atoms undergo both a trapping desorption and also direct inelastic scattering, the latter. yielding information on liquid structure, (3) The angular distribution of atoms at a selected final velocity is sensitive to the local structure and dynamics of the surface. (4) The direct scattering can be conveniently interpreted in terms of contributions from single, double, and multiple collision events, these being roughly equal in relative weight. Forward scattering at grazing angle is dominated by single collisions, while double and multiple collisions have higher contribution at other directions. The double collision contribution in particular contains structural information. (5) There is a substantial yield per collision for sputtering of the squalane-like soft spheres. These results provide insight into the dynamics of gas-liquid collisions, and indicate the usefulness of beam scattering asa tool for studying liquid structure and dynamics.
Roitberg, A. E. ; Gerber, R. B. ; Ratner, M. A. A PERTURBED-MEAN-FIELD APPROACH TO THE DECAY-RATES OF EXCITED VIBRATIONAL-STATES IS EXTENDED SYSTEMS - AN APPLICATION TO I-2(NE)(N).
JOURNAL OF CHEMICAL PHYSICS 1994,
100, 4355-4366.
AbstractA method for calculating decay rates of vibrational modes in large polyatomic systems is proposed and tested. The high frequency excited vibration is treated quantum mechanically, and the soft modes are described classically. The initial state is described by the hybrid quantum/classical self-consistent-field (SCF) approximation. The formalism is based-on a golden-rule expression. The driving potential is the difference between the full Hamiltonian and the mean field Hamiltonian (SCF) causing the decay of the initial state to final mixed quantum/classical ;SCF states. These states are calculated using an extension of the usual static mean-field techniques to systems with mixed quantum and classical degrees of freedom. The formalism for obtaining the mean-field states and calculating the decay rates is presented, and the method is applied to a diatomic molecule treated quantum mechanically, embedded in a 1D model for a. rare gas cluster treated classically. The dependence of the eigenenergies of the quantum and the decay rates with temperature is studied. The influence on the system size is also presented and compared with the prediction of the isolated binary collision model. The effect of a change in the linear density of the cluster on the eigenenergies of the vibrational mode is presented.
Krylov, A. I. ; Gerber, R. B. PHOTODISSOCIATION OF ICN IN SOLID AND IN LIQUID AR - DYNAMICS OF THE CAGE EFFECT AND OF EXCITED-STATE ISOMERIZATION.
JOURNAL OF CHEMICAL PHYSICS 1994,
100, 4242-4252.
AbstractPhotodissociation of ICN by W excitation in solid and liquid Ar is studied by molecular dynamics simulations. The focus is on the differences between the cage effects on the CN photoproduct in the two phases, and on the excited state isomerization ICN* –> INC* dynamics in the solid matrix. Nonadiabatic transitions are neglected in this first study. The main results are: (1) No cage exit of the CN product is found in solid Ar, even in simulations at temperatures close to melting and for large excess energies. The result is in accord with recent experiments by Fraenkel and Haas. This should be contrasted with the large cage-exit probabilities found in many systems for atomic photofragments. The result is interpreted in terms of geometric and energy transfer considerations. It is predicted that complete caging of diatomic and larger photofragments will be typically the case for photodissociation in rare-gas matrices. (2) Almost 100% cage-exit probability for the CN product is found for ICN photolysis on the (II1)-I-1 potential surface in liquid Ar. On the other hand, photolysis on (II0+)-I-3 potential surface does not lead to cage exit on a time scale of 15 ps. The large differences between the reaction in the solid and in the liquid, and between the behavior of the process on the (II0+)-I-3 and the (II1)-I-1 potentials, respectively in the liquid, are interpreted. (3) CN rotational dynamics and subsequent relaxation leads to isomerization int he excited electronic states. On the (II0+)-I-3 potential surface one finds after t greater than or similar to 0.5 ps roughly equal amounts of the ICN and INC isomers. On the (II1)-I-1 surfaces only INC is found after t greater than or similar to 3.5 ps. This is explained in terms of the barriers for CN rotation in the two excited states, and in terms of time scales for rotational relaxation. The results throw light on the differences between cage effects for photochemical reactions in solid and in liquid solution, and on cage-induced isomerization dynamics in solid matrices.
Garciavela, A. ; Gerber, R. B. ; Buck, U. PHOTOLYSIS OF HCL IN THE AR2-HCL AND AR-HCL CLUSTERS - THE CLUSTER-SIZE EFFECT.
JOURNAL OF PHYSICAL CHEMISTRY 1994,
98, 3518-3526.
AbstractThe photodissociation dynamics of HCl in the clusters Ar2-HCl and Ar-HCl is studied in order to explore the cluster size effect. A quasi classical trajectory approach is employ to simulate the photofragmentation dynamics. interesting manifestations of the cage effect are found in the light and the heavy atom kinetic energy distributions, with important differences between the two clusters. The Ar and Cl distributions provide separate information on dynamical events which cannot be resolved in the hydrogen distribution. Only two different excitation wavelengths are used in the simulations. The cage effect is found to be strongly dependent on the wavelength employed to excite the HCl molecule. This is so to the extent that the trend of an increasing cage effect with the cluster size is reversed for the smaller wavelength. New types of resonance behavior are observed for the hydrogen motion in Ar2-HCl compared with the Ar-HCl cluster. An interesting difference between Ar-HCl and Ar2-HCl is that, in the latter case, Ar2 can form as product of the photodissociation, with a high yield for the two wavelengths. The dynamics of the process and the implications of the results are discussed.
McCoy, A. B. ; Gerber, R. B. ; Ratner, M. A. A QUANTITATIVE APPROXIMATION FOR THE QUANTUM DYNAMICS OF HYDROGEN-TRANSFER - TRANSITION-STATE DYNAMICS AND DECAY IN CLHCL-.
JOURNAL OF CHEMICAL PHYSICS 1994,
101, 1975-1987.
AbstractA general method for studying transition state spectroscopy and dynamics in hydrogen atom transfer reactions is presented. This approach is based on the time-dependent self-consistent field (TDSCF) approximation and is applied to a study of the CIHCl- photodetachment experiments of Metz et al. [Metz et al., J. Chem. Phys. 88, 1463 (1988)]. Comparison of results of exact time-dependent and TDSCF calculations are made for collinear and three-dimensional (J=0) approximations for the quantum dynamics. When ClHCl is constrained to be collinear, the TDSCF calculation overcorrelates the motions in the R atom displacement and ClCl extension coordinates. This results in relatively poor agreement with the exact result for many properties of the wavefunction. in contrast, when the system is propagated in the three vibrational coordinates of the system, the transition state dynamics are effectively overmuch more rapidly. Consequently, the TDSCF approximation yields results of very good quantitative accuracy over the time required for most of the wave function to decay off of the transition state. Comparison is also made between the wave function that results from the exact propagation and from TDSCF when the wave function in the ClCl stretch coordinate is approximated by a Gaussian wave packet. Here the magnitude of the overlap between the two TDSCF wave functions in the H atom coordinates, for quantum and semiclassical propagations of the wave function in the ClCl distance Coordinate, is greater than 0.98 over the time of the propagations. These TDSCF calculations are repeated fbi a wave function that is approximated by a product of a two-dimensional wave function in the hydrogen atom coordinates and a one-dimensional wave function in the ClCl extension coordinate and even better quantitative agreement with the exact propagation is achieved. The success of this method for studying ClHCl gives us confidence that TDSCF will provide a general powerful tool for studies of hydrogen and proton transfer reactions in large systems.
Vegiri, A. ; Alexander, M. H. ; Gregurick, S. ; McCoy, A. B. ; Gerber, R. B. QUANTUM MONTE-CARLO STUDIES OF SMALL B(H-2)(N) CLUSTERS.
JOURNAL OF CHEMICAL PHYSICS 1994,
101, 2577-2591.
AbstractThe structure and stability of clusters of a boron atom with one to eight H-2 molecules is investigated. For the simplest BH2 clusters, the lowest ab initio adiabatic potentials for o-H-2 and p-H-2 interacting with a boron atom are used. For the larger clusters (n=2-8), the p-Hz is treated as a sphere, and the total potential is taken to be the sum of pairwise additive B-H-2 and H-2-H-2 interactions which include, in the former case, an anisotropy due to the orientation of the unpaired B 2p electron. This electronic interaction is considerably more attractive when H-2 approaches the B atom in a plane perpendicular to the orientation of the 2p orbital. The local and global minima on these potential surfaces were located and diffusion quantum Monte Carlo simulations were used to determine the energies and properties of the ground state wave functions for these B-(H-2)(n) clusters. For the B(H-2) cluster, a comparison is made with the results of variational calculations.
Schmidt, B. ; Gerber, R. B. REACTIVE COLLISIONS AS A SIGNATURE FOR MELTINGLIKE TRANSITIONS IN CLUSTERS.
PHYSICAL REVIEW LETTERS 1994,
72, 2490-2493.
AbstractThe formation of HCl molecules in collisions of H atoms with Cl(Ar)12 Clusters is studied by classical trajectory calculations as a function of the cluster temperature. At sufficiently low impact energy (E = 10 kJ/mol) a pronounced threshold effect with temperature is observed: In the solidlike regime of the cluster (T less-than-or-equal-to 40 K) there is strong screening of the Cl atom by the icosahedral Ar solvation shell preventing H+Cl association. In the liquidlike phase of the aggregate (T almost-equal-to 45 K) a sharp onset of reactivity is found with a reactive cross section of 11.4 x 10(-20) m2.