Theory of ozone isotopic effects and various electron transfer reactions

Three separate topics, each stimulated by experiments, are treated theoretically in this dessertation: isotopic effects of ozone, electron transfer at interfaces, and intramolecular directional electron transfer in a supramolecular system.

The strange mass-independent isotope effect for the enrichment of ozone, which has been a puzzle in the literature for some 20 years, and the equally puzzling unconventional strong mass-dependent effect of individual reaction rate constants are studied as different aspects of a symmetry-driven behavior. A statistical (RRKM-based) theory with a hindered-rotor transition state is used. The individual rate constant ratios of recombination reactions at low pressures are calculated using the theory involving (1) small deviation from the statistical density of states for symmetric isotopomers, and (2) weak collisions for deactivation of the vibrationally excited ozone molecules. The weak collision and partitioning among exit channels play major roles in producing the large unconventional isotope effect in "unscrambled" systems. The enrichment studies reflect instead the non-statistical effect in "scrambled" systems. The theoretical results of low-pressure ozone enrichments and individual rate constant ratios obtained from these calculations are consistent with the corresponding experimental results. The isotopic exchange rate constant for the reaction ^(16)O + ^(18)O ^(18)O→+ ^(16)O ^(18)O + ^(18)O provides information on the nature of a variationally determined hindered-rotor transition state using experimental data at 130 K and 300 K. Pressure effects on the recombination rate constant, on the individual rate constant ratios and on the enrichments are also investigated. The theoretical results are consistent with the experimental data. The temperature dependence of the enrichment and rate constant ratios is also discussed, and experimental tests are suggested. The desirability of a more accurate potential energy surface for ozone in the transition state region is also noted.

Electron transfer reactions at semiconductor /liquid interfaces are studied using a tight-binding model for the semiconductors. The slab method and a z-transform method are employed in obtaining the tight-binding electronic structures of semiconductors having surfaces. The maximum electron transfer rate constants at Si/viologen^(2-/+) and InP /Me_(2)Fc^(+/O) interfaces are computed using the tight-binding type calculations for the solid and the extended-Huckel for the coupling to the redox agent at the interface. These electron transfer reactions are also studied using a free electron model for the semiconductor and the redox molecule, where Bardeen's method is adapted to calculate the coupling matrix element between the molecular and semiconductor electronic states. The calculated results for maximum rate constant of the electron transfer from the semiconductor bulk states are compared with the experimentally measured values of Lewis and coworkers, and are in reasonable agreement, without adjusting parameters. In the case of InP /liquid interface, the unusual current vs applied potential behavior is additionally interpreted, in part, by the presence of surface states.

Photoinduced electron transfer reactions in small supramolecular systems, such as 4-aminonaphthalimide compounds, are interesting in that there are, in principle, two alternative pathways (directions) for the electron transfer. The electron transfer, however, is unidirectional, as deduced from pH-dependent fluorescence quenching studies on different compounds. The role of electronic coupling matrix element and the charges in protonation are considered to explain the directionality of the electron transfer and other various results. A related mechanism is proposed to interpret the fluorescence behavior of similar molecules as fluorescent sensors of metal ions.

Efficient sampling of atomic configurational spaces.

We describe a method to explore the configurational phase space of chemical systems. It is based on the nested sampling algorithm recently proposed by Skilling (AIP Conf. Proc. 2004, 395; J. Bayesian Anal. 2006, 1, 833) and allows us to explore the entire potential energy surface (PES) efficiently in an unbiased way. The algorithm has two parameters which directly control the trade-off between the resolution with which the space is explored and the computational cost. We demonstrate the use of nested sampling on Lennard-Jones (LJ) clusters. Nested sampling provides a straightforward approximation for the partition function; thus, evaluating expectation values of arbitrary smooth operators at arbitrary temperatures becomes a simple postprocessing step. Access to absolute free energies allows us to determine the temperature-density phase diagram for LJ cluster stability. Even for relatively small clusters, the efficiency gain over parallel tempering in calculating the heat capacity is an order of magnitude or more. Furthermore, by analyzing the topology of the resulting samples, we are able to visualize the PES in a new and illuminating way. We identify a discretely valued order parameter with basins and suprabasins of the PES, allowing a straightforward and unambiguous definition of macroscopic states of an atomistic system and the evaluation of the associated free energies.

On representing chemical environments

We review some recently published methods to represent atomic neighbourhood environments, and analyse their relative merits in terms of their faithfulness and suitability for fitting potential energy surfaces. The crucial properties that such representations (sometimes called descriptors) must have are differentiability with respect to moving the atoms, and invariance to the basic symmetries of physics: rotation, reflection, translation, and permutation of atoms of the same species. We demonstrate that certain widely used descriptors that initially look quite different are specific cases of a general approach, in which a finite set of basis functions with increasing angular wave numbers are used to expand the atomic neighbourhood density function. Using the example system of small clusters, we quantitatively show that this expansion needs to be carried to higher and higher wave numbers as the number of neighbours increases in order to obtain a faithful representation, and that variants of the descriptors converge at very different rates. We also propose an altogether new approach, called Smooth Overlap of Atomic Positions (SOAP), that sidesteps these difficulties by directly defining the similarity between any two neighbourhood environments, and show that it is still closely connected to the invariant descriptors. We test the performance of the various representations by fitting models to the potential energy surface of small silicon clusters and the bulk crystal.

High-spin isomeric structures in exotic odd-odd nuclei: Exploration of the proton drip line and beyond

We investigate the role of two-quasiparticle isomeric states along the proton drip line, using configuration-constrained potential-energy-surface calculations. In contrast to even-even nuclei, odd-odd nuclei can have coexisting low-lying two-quasiparticle states. The low excitation energy and high angular momentum can lead to long-lived isomers. Also, because of the hindrance by spin selection, the probabilities of beta and proton decays from high-spin isomers can be reduced significantly. The present calculations reproduce reasonably well the available data for observed isomers in such nuclei. Unobserved high-spin isomers are predicted, which could be useful for future experimental studies of exotic nuclei at and beyond the proton drip line.

beta-delayed proton decays and spin assignments for Sm-133 and Yb-149

The proton-rich isotope Sm-133 was produced via the fusion evaporation reaction Ca-40 + Ru-96. Its beta-delayed proton decay was studied by p-gamma coincidence in combination with a He-jet tape transport system, and half-lives, proton energy spectra, gamma-transitions following the proton emission, as well as beta-delayed proton branching ratios to the low-lying states in the grand-daughter nucleus were determined. Comparing the observed beta-delayed proton branching ratios with statistical model calculations, the best agreement is found assuming that only one level with the spin of 3/2 in Sm-133 decays or two levels with the spins of 1/2 and 5/2 decay with similar half-lives. The configuration-constrained nuclear potential energy surfaces of Sm-133 were calculated using the Woods-Saxon-Strutinsky method, which suggests a 1/2-ground state and a 5/2(+) isomer with an excitation energy of 120 keV. Therefore, the simple(EC+beta(+)) decay scheme of Sm-133 in Eur. Phys. J.A 11,277(2001) has been revised. In addition, our previous experimental data on the beta-delayed proton decay of Yb-149 reported in Eur. Phys. J. A 12,1 ( 2 0 0 1) was also analyzed using the same method. The spin-parity of Yb-149 is suggested to be 1/2(-).

beta-delayed proton decays and spin assignments for Tb-140, Dy-141 and Dy-143

The proton-rich isotopes Tb-140 and Dy-141 were produced via the fusion evaporation reaction Ca-40+ Cd-106. Their beta-delayed proton decays were studied by p-gamma coincidence in combination with a He-jet tape transport system, and half-lives, proton energy spectra, gamma-transitions following the proton emission, as well as beta-delayed proton branching ratios to the low-lying states in the grand-daughter nuclei were determined. Comparing the experimental data with statistical model calculations, the ground-state spins of Tb-140 and Dy-141 were found to be consistent with 7 and 9/2, respectively. The configuration-constrained nuclear potential energy surfaces (NPES) of Tb-140 and Dy-141 were calculated using the Woods-Saxon-Strutinsky method, which suggest the ground-state spins and parities of Tb-140 and Dy-141 to be 7(+) and 9/2(-), respectively. In addition, the configuration-constrained NPES of Dy-143 were calculated, which predict a 1/2(+) ground state and a 11/2(-) isomer with excitation energy of 198 keV. These findings are consistent with our previous experimental data on Dy-143 reported in Eur. Phys. J. A 16, 347 (2003).

Molecular dynamics and control following excitation with an ultra-short intense laser pulse

Vibronic excitations of the tri-atomic molecule OClO (A(2)A(2)(nu(1), nu(2), nu(3)) <-- (XB1)-B-2 (0, 0, 0)) with weak and strong ultra-short laser fields are studied within full quantum wavepacket dynamics in hyperspherical coordinates. Different dynamics is observed following excitation with laser pulses of different intensities. With a strong laser pulse, many vibrational states are excited and a spatially more localised wavepacket arises. The numerical results show that the population of different vibrational states of the wavepacket on the excited potential energy surface is altered by the intensity of the laser pulse. The numerical results also suggest a related effect on the phase of the wavepacket. These interesting phenomena can be understood by an analysis of the corresponding results for two model diatomic molecules. The possible physical mechanisms of control of chemical processes using strong laser fields are discussed. (C) 2004 Elsevier B.V. All rights reserved.

Theoretical study on the formation mechanism of Iso-CH2I-Cl

The dissociation and isomerization reaction mechanism on the ground-state potential energy surface for CH2ClI are investigated by ab initio calculations. It is found that the isomer iso-CH2I-Cl can be produced from either the recombination of the photodissociation. fragments or the isomerization reaction of CH2ClI, rather than from isomerization reaction of iso-CH2Cl-I. Further explanations of experimental results are also presented. (C) 2003 Wiley Periodicals, Inc.

Close-coupling time-dependent quantum dynamics study of the H+HCl reaction

The paper presents a theoretical study of the dynamics of the H + HCl system on the potential energy surface (PES) of Bian and Werner (Bian, W.; Werner, H. -J., J. Chem. Phys. 2000, 112, 220). A time-dependent wave packet approach was employed to calculate state-to-state reaction probabilities for the exchanged and abstraction channels. The most recent PES for the system has been used in the calculations. Reaction probabilities have also been calculated for several values of the total angular momentum J > 0. Those have then been used to estimate cross sections and rate constants for both channels. The calculated cross sections can be compared with the results of previous quasiclassical trajectory calculations and reaction dynamics experimental on the abstraction channel. In addition, the calculated rate constants are in the reasonably good agreement with experimental measurement.

Theoretical study on the mechanism of the reaction of CH4+MgO

The reactions of (1) CH4 + MgO --> MgOH. + CH3. and (2) CH4 + MgO --> Mg + CH3OH have been studied on the singlet spin state potential energy surface at the MP2/6-311+G(2d,2p) level. These two reaction channels, both involving intermediates and transition states, have been rationalized by the structures of the species involved, natural bond orbital (NBO), and vibrational frequency analysis. We have considered two initial interacting models between CH4 and MgO: a collinear C-H approach to the O end of the MgO forming the MgOCH4 complex with C-3nu symmetry and three hydrogen atoms of the methane point to the Mg end of the MgO forming the OMgCH4 complex with C-1 symmetry. The calculations predict that reactions 1 and 2 are exothermic by 39.8 and 86.5 kJ mol(-1), respectively. Also, the former reaction proceeds more easily than the latter, and the complex HOMgCH3 is energetically preferred in the reaction of MgO + CH4.

Accurate quantum mechanical calculation for the N+OH reaction

Accurate three-dimensional time-dependent quantum wave packet calculations for the N+OH reaction on the (3)A' potential energy surface [Guadagnini, Schatz, and Walch, J. Chem. Phys. 102, 774 (1995)] have been carried out. The calculations show for the first time that the initial state-selected reaction probabilities are dominated by resonance structures, and the lifetime of the resonance is generally in the subpicosecond time scale. The calculated reaction cross sections indicate that they are a decreasing function of the translational energy, which is in agreement qualitatively with the quasiclassical trajectory calculations. The rate constants obtained from the quantum mechanical calculations are consistent with the quasiclassical trajectory results and the experimental measurements. (C) 2003 American Institute of Physics.

Effects of the transition dipole moment function on the dynamics of ozone photodissociation: an exact 3D quantum mechanical study

The effects of the transition dipole moment function (TDMF) on the dynamics Of O-3 photodissociation in the Hartley band have been exploited by means of exact 3D time-dependent wavepacket method using the SW potential energy surface [J. Chem. Phys. 78 (1983) 7191]. The calculations show that the explicit inclusion of the TDMF results in slight uniform reductions for the intensities of recurrence peaks of the autocorrelation function and a slight broadening of the absorption spectrum, in comparison with the result where the TDMF is assumed to be constant. The pattern of recurrence structures of the autocorrelation function is essentially unaffected. (C) 2001 Elsevier Science B.V. All rights reserved.

A density function theory study on the properties and decomposition of CF3OF

The density function theory was used to calculate the potential energy surface for the decomposition of CF3OF. The geometries, vibrational frequencies and energies of all stationary points were obtained. The calculated harmonic frequencies agreed well with the experimental ones. Three decomposition channels of CF3OF were studied. The calculated reaction enthalpy (29.85 kcal/mol) of the elimination reaction CF3OF --> CF2O + F-2 was in good agreement with the experimental value (27.7 kcal/mol). The O-F bond of CF3OF is broken easily by comparing the energies, while the decomposition channel to yield the CF30 and F radicals is the main reaction path. (C) 2002 Published by Elsevier Science B.V.

Theoretical studies of product polarization and state distributions of the H+HCl reaction

The angular momentum polarization and rotational state distributions of the H-2 and HCl products from the H + HCl reaction are calculated at a relative translational energy of 1.6 eV by using quasiclassical trajectories on two potential energy surfaces, one from G3 surface [T.C. Allison et al., J. Phys. Chem. 100 (1996) 13575], and the other from BW2 surface [W. Bian, H.-J. Werner, J. Chem. Phys. 112 (2000) 220]. Product rotational distributions obtained on the G3 potential energy surface (PES) are much closer to the experimental results (P.M. Aker et al., J. Chem. Phys. 90 (1989) 4795; J. Chem. Phys. 90 (1989) 4809) than the distributions calculated on the BW2 PES. The distributions of P(phi(r)) for the H-2 and HCl products obtained on the G3 PES are similar, whereas the rotational alignment effect of the H-2 product is stronger than that of the HCl product. In contrast to the polarization distributions obtained on the G3 PES, the rotational alignment effect of the two products calculated on the BW2 PES is similar. However, the abstraction reaction is dominated by out-of-plane mechanisms, while the exchange reaction is dominated by in-plane mechanisms. The significant difference of the product rotational polarization obtained on the G3 and BW2 PESs implies that the studies of the dynamical stereochemistry can provide a sensitive test for the accuracy of the PES. (C) 2002 Elsevier Science B.V. All rights reserved.