Mechanism and Kinetics of The Degradation of Nitazoxanide and Hydroxychloroquine Drugs by Hydroxyl Radicals: Theoretical Approach to Ecotoxicity

The efforts of contrasting the effects caused by the Covid-19 (coronavirus disease 2019) pandemic increased the disposal of active pharmaceutical ingredients. This paper reports the mechanisms and kinetics of the degradation in aqueous environments induced by  OH of two drugs, among those most widely probed at the outbreak of coronavirus, nitazoxanide and hydroxychloroquine. The investigation exploits quantum chemistry techniques and a reaction rate theory combined with diffusion-controlled processes and quantum mechanical tunneling. The reaction rate constants are obtained in an environmentally relevant temperature range. The results show that (i) the deacetylation of nitazoxanide with formation of tizoxanide is kinetically the most favorable channel, in agreement with experimental work; (ii) for hydroxychloroquine, the present theoretical calculations show that the most favorable channel is the addition of  OH at the aromatic ring. The half-life time degradation products are for both cases in the range between 12 to 138 days. Both drugs presented toxicities between harmful and toxic as obtained by computational toxicology calculations: The toxicity is also calculated for the degradation products: (i) in the nitazoxanide degradation process, tizoxanide was characterized as more toxic, while (ii) in the case of hydroxychloroquine, the major degradation product showed a decrease in the toxicity.

UV Photodissociation of Halothane in a Focused Molecular Beam: Space-Speed Slice Imaging of Competitive Bond Breaking into Spin-Orbit-Selected Chlorine and Bromine Atoms

A molecular beam of halothane (2-bromo-2-chloro-1,1,1-trifluoroethane) is focused by a hexapolar electrostatic field and photolyzed by UV laser radiation at 234 nm. Angular and speed distributions of chlorine and bromine photofragments emitted from halothane are measured for both spin-orbit states independently. Although the dissociation energy of the C-Cl bond is larger than that of C-Br, the relative yield of Cl to Br was found to be approximately 2. Measured speed and angular distributions of atomic fragments show distinct kinetic energy release and scattering characteristics: for bromine, observed fast and aligned fragments exhibit a signature of a direct mode of dissociation for the C-Br bond, via the electronically excited potential energy surface denoted nσ*(C-Br), of repulsive nature; for chlorine, a variation in the features is observed for the dissociation pathway through nσ*(C-Cl), from a modality similar to the bromine case, leading to fragments with appreciable kinetic energy release and pronounced directionality, to a modality involving slow products, nearly isotopically distributed. The origin of this behavior can be attributed to nonadiabatic interaction operating between the nσ*(C-Br) and nσ*(C-Cl) surfaces. These results are not only relevant for a detailed understanding of adiabatic versus diabatic coupling mechanisms in the manifold of excited states populated by photon absorption, but they also point out the possibility of selectively inducing specific dissociation pathways, even when involving energetically unfavorable outcomes, such as, in this case, the prevailing rupture of the stronger C-Cl bond against that of the weaker C-Br bond.

Benchmark Quantum Kinetics at Low Temperatures toward Absolute Zero and Role of Entrance Channel Wells on Tunneling, Virtual States, and Resonances: The F + HD Reaction

This paper reports a study of the quantum reaction dynamics and kinetics of the F + HD reaction at low and ultralow temperatures, focusing on the range from the Wigner limit up to 50 K. Close coupling time-independent quantum reactive scattering calculations for the production of HF and DF molecules have been carried out on two potential energy surfaces differing in the description of the reaction entrance channel. This case is computationally more demanding than the cases of F with H₂ and D₂ ( De Fazio et al. Frontiers in Chemistry 2019 , 7 , 328 ) but offers a wider phenomenology regarding the roles of quantum mechanical effects of tunneling, of virtual states, and of resonances. The results show that at the temperatures in the cold and ultracold regimes small changes in the entrance channel long-range interaction induce surprising near threshold features. The presence of a virtual state close to the reactive threshold gives rise to a marked anti-Arrhenius behavior of the rate constants below 100 mK. This effect enhances reaction rates by about 2 orders of magnitude, making them of the same order as those at room temperature and confining the onset of the Wigner regime in the microkelvin region.

Stereodynamic Imaging of Bromine Atomic Photofragments Eliminated from 1-Bromo-2-methylbutane Oriented via Hexapole State Selector

Both single-laser and two-laser experiments were conducted to look into the ion-imaging of Br*(²P_1/2) and Br(²P_3/2) photofragmented from 1-bromo-2-methylbutane in the range 232-240 nm via a detection scheme of (2+1) resonance-enhanced multiphoton ionization. The angular analysis of these photofragment distributions yields the anisotropy parameter β = 1.88 ± 0.06 for the Br* excited state which arises from a parallel transition, while β = 0.63 ± 0.09 for the Br ground state indicates the contribution from both a perpendicular transition and a non-adiabatic transition. When a hexapole coupled with an orienting field was implemented, the parent molecules are spatially oriented to yield an orientation efficiency |⟨cos θ⟩| of 0.15. Besides the χ angle between the recoil velocity v and the transition dipole moment μ, orienting molecules allows for the evaluation of the angle α between v and the permanent molecular dipole moment d. The angular analysis of Br* photofragment distribution yields χ = 11.5° and α in the range from 160° to 180° with weak dependency. In the two-laser experiments, the angular anisotropy of Br photofragment distribution was found to be smaller (0.38 ± 0.10) when the photolysis wavelength was red-shifted to 240 nm, suggesting the increasing contributions from perpendicular transitions.

Vectorial imaging of the photodissociation of 2-bromobutane oriented via hexapolar state selection

Molecular orientation techniques are becoming available in the study of elementary chemical processes, in order to highlight those structural and dynamical properties that would be concealed by random rotational motions. Recently successful orientation was achieved for asymmetric-top and chiral molecules of much larger complexity than hitherto. In this work, we report and discuss the correlation between the vectors' photofragment recoil velocity v, transition dipole moment μ, and permanent dipole moment d in a dissociation experiment on hexapole oriented 2-bromobutane, photoinitiated by a linearly polarized laser. The sliced ion images of the Br*(2P1/2) and Br(2P3/2) photofragments were acquired at 234.0 and 254.1 nm, respectively, by a (2 + 1) resonance-enhanced multiphoton ionization technique. A detailed analysis of the sliced ion images obtained at a tilting angle 45° of laser polarization provides information on the correlation of the three vectors, which are confined by two polar angles α and χ and one azimuthal angle φμd in the recoil frame. The sliced ion images of Br fragments eliminated individually from the enantiomers at 254.1 nm yield an asymmetric factor close to zero; for this reason the photofragment angular distributions do not show significant differences. The elimination of the Br* fragment at 234.0 nm is mainly correlated with a parallel transition, giving rise to a large anisotropy parameter of 1.85, and thus can be considered as a single state excitation. The resulting recoil frame angles are optimized to 163° ± 8° and 164° ± 1° for α and χ, respectively, whereas φμd is approaching 0° for the best fit. Since for the present molecule, the three vectors have an only slight spatial arrangement, the photofragment angular distributions of the two enantiomers do not show appreciable differences. Theoretical and computational simulations provide us the basis to state that oriented enantiomers can be discriminated on-the-fly in photodissociation processes even initiated by non-circularly polarized light, provided that the three vectors encountered above have specific three-dimensional arrangements. The fact that Br fragment elimination involves a multi-potential dissociation carries uncertainties in theoretical estimates of the vector direction. Therefore, this work represents a preliminary but significant step on the road to chiral discrimination on-the-fly, which is shown to be best propitiated in molecules where vectors are far from having degenerate mutual angular directions.

Quantum Dynamics and Kinetics of the F + H2 and F + D2 Reactions at Low and Ultra-Low Temperatures

Integral cross sections and rate constants for the prototypical chemical reactions of the fluorine atom with molecular hydrogen and deuterium have been calculated over a wide interval of collision energy and temperature ranging from the sub-thermal (50 K) down to the ultra-cold regimes (0.5 mK). Rigorous close coupling time-independent quantum reactive scattering calculations have been carried out on two potential energy surfaces, differing only at long-range in the reactants' channel. The results show that tunnel, resonance and virtual state effects enhance under-barrier reactivity giving rise to pronounced deviations from the Arrhenius law as temperature is lowered. Within the ultra-cold domain (below 1 mK), the reactivity is governed by virtual state effects and by tunneling through the reaction barrier; in the cold regime (1 mK–1 K), the shape resonances in the entrance channel of the potential energy surface make the quantum tunneling contribution larger so enhancing cross sections and rate constants by about one order of magnitude; at higher temperatures (above 10 K), the tunneling pathway enhanced by the constructive interference between two Feshbach resonances trapped in the reaction exit channel competes with the thermally activated mechanism, as the energy gets closer to the reaction barrier height. The results show that at low temperatures cross sections and rate constants are extremely sensitive to small changes in the long-range intermolecular interaction in the entrance channel of the potential energy surface, as well as to isotopic substitution.

Temperature Dependence of Rate Processes Beyond Arrhenius and Eyring: Activation and Transitivity

Advances in the understanding of the dependence of reaction rates from temperature, as motivated from progress in experiments and theoretical tools (e. g., molecular dynamics), are needed for the modeling of extreme environmental conditions (e.g., in astrochemistry and in the chemistry of plasmas). While investigating statistical mechanics perspectives (Aquilanti et al., 2017b, 2018), the concept of transitivity was introduced as a measure for the propensity for a reaction to occur. The Transitivity plot is here defined as the reciprocal of the apparent activation energy vs. reciprocal absolute temperature. Since the transitivity function regulates transit in physicochemical transformations, not necessarily involving reference to transition-state hypothesis of Eyring, an extended version is here proposed to cope with general types of transformations. The transitivity plot permits a representation where deviations from Arrhenius behavior are given a geometrical meaning and make explicit a positive or negative linear dependence of transitivity for sub- and super-Arrhenius cases, respectively. To first-order in reciprocal temperature, the transitivity function models deviations from linearity in Arrhenius plots as originally proposed by Aquilanti and Mundim: when deviations are increasingly larger, other phenomenological formulas, such as Vogel-Fulcher-Tammann, Nakamura-Takayanagi-Sato, and Aquilanti-Sanches-Coutinho-Carvalho are here rediscussed from the transitivity concept perspective and with in a general context. Emphasized is the interest of introducing into this context modifications to a very successful tool of theoretical kinetics, Eyring's Transition-State Theory: considering the behavior of the transitivity function at low temperatures, in order to describe deviation from Arrhenius behavior under the quantum tunneling regime, a "d-TST" formulation was previously introduced (Carvalho-Silva et al., 2017). In this paper, a special attention is dedicated to a derivation of the temperature dependence of viscosity, making explicit reference to feature of the transitivity function, which in this case generally exhibits a super-Arrhenius behavior. This is of relevance also for advantages of using the transitivity function for diffusion-controlled phenomena.

Kinetics of the OH+HCl→H2 O+Cl reaction: Rate determining roles of stereodynamics and roaming and of quantum tunneling

The OH + HCl → H₂ O + Cl reaction is one of the most studied four-body systems, extensively investigated by both experimental and theoretical approaches. Here, as a continuation of our previous work on the OH + HBr and OH + HI reactions, which manifest an anti-Arrhenius behavior that was explained by stereodynamic and roaming effects, we extend the strategy to understand the transition to the sub-Arrhenius behavior occurring for the HCl case. As previously, we perform first-principles on-the-fly Born-Oppenheimer molecular dynamics calculations, thermalized at four temperatures (50, 200, 350, and 500 K), but this time we also apply a high-level transition-state-theory, modified to account for tunneling conditions. We find that the theoretical rate constants calculated with Bell tunneling corrections are in good agreement with extensive experimental data available for this reaction in the ample temperature range: (i) simulations show that the roles of molecular orientation in promoting this reaction and of roaming in finding the favorable path are minor than in the HBr and HI cases, and (ii) dominating is the effect of quantum mechanical penetration through the energy barrier along the reaction path on the potential energy surface. The discussion of these results provides clarification of the origin on different non-Arrhenius mechanisms observed along this series of reactions. © 2018 Wiley Periodicals, Inc.

Double photoionization of propylene oxide: A coincidence study of the ejection of a pair of valence-shell electrons

Propylene oxide, a favorite target of experimental and theoretical studies of circular dichroism, was recently discovered in interstellar space, further amplifying the attention to its role in the current debate on protobiological homochirality. In the present work, a photoelectron-photoion-photoion coincidence technique, using an ion-imaging detector and tunable synchrotron radiation in the 18.0-37.0 eV energy range, permits us (i) to observe six double ionization fragmentation channels, their relative yields being accounted for about two-thirds by the couple (C₂H₄⁺, CH₂O⁺) and one-fifth by (C₂H₃⁺, CH₃O⁺); (ii) to measure thresholds for their openings as a function of photon energy; and (iii) to unravel a pronounced bimodality for a kinetic-energy-released distribution, fingerprint of competitive non-adiabatic mechanisms.

Kinetics of low-temperature transitions and a reaction rate theory from non-equilibrium distributions

This article surveys the empirical information which originated both by laboratory experiments and by computational simulations, and expands previous understanding of the rates of chemical processes in the low-temperature range, where deviations from linearity of Arrhenius plots were revealed. The phenomenological two-parameter Arrhenius equation requires improvement for applications where interpolation or extrapolations are demanded in various areas of modern science. Based on Tolman's theorem, the dependence of the reciprocal of the apparent activation energy as a function of reciprocal absolute temperature permits the introduction of a deviation parameter d covering uniformly a variety of rate processes, from those where quantum mechanical tunnelling is significant and d < 0, to those where d > 0, corresponding to the Pareto-Tsallis statistical weights: these generalize the Boltzmann-Gibbs weight, which is recovered for d = 0. It is shown here how the weights arise, relaxing the thermodynamic equilibrium limit, either for a binomial distribution if d > 0 or for a negative binomial distribution if d < 0, formally corresponding to Fermion-like or Boson-like statistics, respectively. The current status of the phenomenology is illustrated emphasizing case studies; specifically (i) the super-Arrhenius kinetics, where transport phenomena accelerate processes as the temperature increases; (ii) the sub-Arrhenius kinetics, where quantum mechanical tunnelling propitiates low-temperature reactivity; (iii) the anti-Arrhenius kinetics, where processes with no energetic obstacles are rate-limited by molecular reorientation requirements. Particular attention is given for case (i) to the treatment of diffusion and viscosity, for case (ii) to formulation of a transition rate theory for chemical kinetics including quantum mechanical tunnelling, and for case (iii) to the stereodirectional specificity of the dynamics of reactions strongly hindered by the increase of temperature.This article is part of the themed issue 'Theoretical and computational studies of non-equilibrium and non-statistical dynamics in the gas phase, in the condensed phase and at interfaces'.

Deformed transition-state theory: Deviation from Arrhenius behavior and application to bimolecular hydrogen transfer reaction rates in the tunneling regime

A formulation is presented for the application of tools from quantum chemistry and transition-state theory to phenomenologically cover cases where reaction rates deviate from Arrhenius law at low temperatures. A parameter d is introduced to describe the deviation for the systems from reaching the thermodynamic limit and is identified as the linearizing coefficient in the dependence of the inverse activation energy with inverse temperature. Its physical meaning is given and when deviation can be ascribed to quantum mechanical tunneling its value is calculated explicitly. Here, a new derivation is given of the previously established relationship of the parameter d with features of the barrier in the potential energy surface. The proposed variant of transition state theory permits comparison with experiments and tests against alternative formulations. Prescriptions are provided and implemented to three hydrogen transfer reactions: CH₄  + OH → CH₃  + H₂ O, CH₃ Cl + OH → CH₂ Cl + H₂ O and H₂  + CN → H + HCN, widely investigated both experimentally and theoretically. © 2016 Wiley Periodicals, Inc.

Hexapole-Oriented Asymmetric-Top Molecules and Their Stereodirectional Photodissociation Dynamics

Molecular orientation is a fundamental requisite in the study of stereodirected dynamics of collisional and photoinitiated processes. In this past decade, variable hexapolar electric filters have been developed and employed for the rotational-state selection and the alignment of molecules of increasing complexity, for which the main difficulties are their mass, their low symmetry, and the very dense rotational manifold. In this work, for the first time, a complex molecule such as 2-bromobutane, an asymmetric top containing a heavy atom (the bromine), was successfully oriented by a weak homogeneous field placed downstream from the hexapolar filter. Efficiency of the orientation was characterized experimentally, by combining time-of-flight measurements and a slice-ion-imaging detection technique. The application is described to the photodissociation dynamics of the oriented 2-bromobutane, which was carried out at a laser wavelength of 234 nm, corresponding to the breaking of the C-Br bond. The Br photofragment is produced in both the ground Br ((2)P3/2) and the excited Br ((2)P1/2) electronic states, and both channels are studied by the slice imaging technique, revealing new features in the velocity and angular distributions with respect to previous investigations on nonoriented molecules.

Stereodirectional Origin of anti-Arrhenius Kinetics for a Tetraatomic Hydrogen Exchange Reaction: Born-Oppenheimer Molecular Dynamics for OH + HBr

Among four-atom processes, the reaction OH + HBr → H2O + Br is one of the most studied experimentally: its kinetics has manifested an unusual anti-Arrhenius behavior, namely, a marked decrease of the rate constant as the temperature increases, which has intrigued theoreticians for a long time. Recently, salient features of the potential energy surface have been characterized and most kinetic aspects can be considered as satisfactorily reproduced by classical trajectory simulations. Motivation of the work reported in this paper is the investigation of the stereodirectional dynamics of this reaction as the prominent reason for the peculiar kinetics: we started in a previous Letter ( J. Phys. Chem. Lett. 2015 , 6 , 1553 - 1558 ) a first-principles Born-Oppenheimer "canonical" molecular dynamics approach. Trajectories are step-by-step generated on a potential energy surface quantum mechanically calculated on-the-fly and are thermostatically equilibrated to correspond to a specific temperature. Here, refinements of the method permitted a major increase of the number of trajectories and the consideration of four temperatures -50, +200, +350, and +500 K, for which the sampling of initial conditions allowed us to characterize the stereodynamical effect. The role is documented of the adjustment of the reactants' mutual orientation to encounter the entrance into the "cone of acceptance" for reactivity. The aperture angle of this cone is dictated by a range of directions of approach compatible with the formation of the specific HOH angle of the product water molecule; and consistently the adjustment is progressively less effective the higher the kinetic energy. Qualitatively, this emerging picture corroborates experiments on this reaction, involving collisions of aligned and oriented molecular beams, and covering a range of energies higher than the thermal ones. The extraction of thermal rate constants from this molecular dynamics approach is discussed and the systematic sampling of the canonical ensemble is indicated as needed for quantitative comparison with the kinetic experiments.

Interaction of O2 with CH4, CF4, and CCl4 by Molecular Beam Scattering Experiments and Theoretical Calculations

Gas phase collisions of O2 by CH4, CF4, and CCl4 have been investigated with the molecular beam technique by measuring both the integral cross section value, Q, and its dependence on the collision velocity, v. The adopted experimental conditions have been appropriate to resolve the oscillating "glory" pattern, a quantum interference effect controlled by the features of the intermolecular interaction, for all the three case studies. The analysis of the Q(v) data, performed by adopting a suitable representation of the intermolecular potential function, provided the basic features of the anisotropic potential energy surfaces at intermediate and large separation distances and information on the relative role of the physically relevant types of contributions to the global interaction. The present work demonstrates that while O2-CH4 and O2-CF4 are basically bound through the balance between size (Pauli) repulsion and dispersion attraction, an appreaciable intermolecular bond stabilization by charge transfer is operative in O2-CCl4. Ab initio calculations of the strength of the interaction, coupled with detailed analysis of electronic charge displacement promoted by the formation of the dimer, fully rationalizes the experimental findings. This investigation indicates that the interactions of O2, when averaged over its relative orientations, are similar to that of a noble gas (Ng), specifically Ar. We also show that the binding energy in the basic configurations of the prototypical Ng-CF4,CCl4 systems [ Cappelletti , D. ; Chem. Eur. J. 2015 , 21 , 6234 - 6240 ] can be reconstructed by using the interactions in Ng-F and Ng-Cl systems, previously characterized by molecular beam scattering experiments of state-selected halogen atom beams. This information is fundamental to approach the modeling of the weak intermolecular halogen bond. On the basis of the electronic polarizability, this also confirms [ Aquilanti , V. ; Angew. Chem., Int. Ed. 2005 , 44 , 2356 - 2360 ] that O2 can be taken as a proper reference partner for simulating the behavior of some basic noncovalent components of the interactions involving water. Present results are of fundamental relevance to build up the force field controlling the hydrophobic behavior of prototypical apolar CX4 (X = H, F, Cl) molecules.

Dynamical, spectroscopic and computational imaging of bond breaking in photodissociation: roaming and role of conical intersections

Recent experimental and theoretical advances in the study of the dissociation of excited molecules are revealing unexpected mechanisms, when their outcomes are tackled by combining (i) space-time ion imaging of translational features, with (ii) spectroscopic probing of rotational and vibrational distributions; crucial is the assistance of (iii) the quantum chemistry of structural investigations of rearrangements of chemical bonds, and of (iv) the simulations of molecular dynamics to follow the evolution of selective bond stretching and breaking. Here we present results of such an integrated approach to methyl formate, HCOOCH3, the simplest of esters; the main focus is on the rotovibrationally excited CO (v=1) product and in general on the energy distribution in the fragments. Previous laser studies of dissociation into CO and CH3OH at a sequence of various wavelengths discovered signatures of a roaming mechanism by the late arrival of CO (v=0) products in time-of-flight ion imaging. Subsequent detailed investigations as a function of excitation energy provided the assessment of the threshold, which opens for triple breakdown into CO and further fragments H and CH3O, as spectroscopically characterized by ion imaging and FTIR respectively. Accompanying quantum mechanical electronic structure calculations and classical molecular dynamics simulations clarify the origin of these fragments through "roaming" pathways involving incipient radical intermediates at energies below the triple fragmentation threshold: a specific role is played by nonadiabatic transitions at a conical intersection between ground and excited states; alternative pathways focalize our attention to regions of the potential energy surfaces other than those in the neighbourhoods of saddle points along minimum energy paths: eventually this leads us to look for avenues in reaction kinetics beyond those of venerable transition state theories.

Chirality of weakly bound complexes: the potential energy surfaces for the hydrogen-peroxide-noble-gas interactions

We consider the analytical representation of the potential energy surfaces of relevance for the intermolecular dynamics of weakly bound complexes of chiral molecules. In this paper we study the H2O2-Ng (Ng=He, Ne, Ar, Kr, and Xe) systems providing the radial and the angular dependence of the potential energy surface on the relative position of the Ng atom. We accomplish this by introducing an analytical representation which is able to fit the ab initio energies of these complexes in a wide range of geometries. Our analysis sheds light on the role that the enantiomeric forms and the symmetry of the H2O2 molecule play on the resulting barriers and equilibrium geometries. The proposed theoretical framework is useful to study the dynamics of the H2O2 molecule, or other systems involving O-O and S-S bonds, interacting by non-covalent forces with atoms or molecules and to understand how the relative orientation of the O-H bonds changes along collisional events that may lead to a hydrogen bond formation or even to selectivity in chemical reactions.

Stereodynamical Origin of Anti-Arrhenius Kinetics: Negative Activation Energy and Roaming for a Four-Atom Reaction

The OH + HBr → H2O + Br reaction, prototypical of halogen-atom liberating processes relevant to mechanisms for atmospheric ozone destruction, attracted frequent attention of experimental chemical kinetics: the nature of the unusual reactivity drop from low to high temperatures eluded a variety of theoretical efforts, ranking this one among the most studied four-atom reactions. Here, inspired by oriented molecular-beams experiments, we develop a first-principles stereodynamical approach. Thermalized sets of trajectories, evolving on a multidimensional potential energy surface quantum mechanically generated on-the-fly, provide a map of most visited regions at each temperature. Visualizations of rearrangements of bonds along trajectories and of the role of specific angles of reactants' mutual approach elucidate the mechanistic change from the low kinetic energy regime (where incident reactants reorient to find the propitious alignment leading to reaction) to high temperature (where speed hinders adjustment of directionality and roaming delays reactivity).

Theoretical reaction kinetics astride the transition between moderate and deep tunneling regimes: the F + HD case

For the reaction between F and HD, giving HF + D and DF + H, the rate constants, obtained from rigorous quantum scattering calculations at temperatures ranging from 350 K down to 100 K, show deviations from the Arrhenius behavior that have been interpreted in terms of tunneling of either H or D atoms through a potential energy barrier. The interval of temperature investigated extends from above to below a crossover value Tc, a transition temperature separating the moderate and deep quantum tunneling regimes. Below Tc, the rate of the H or D exchange reaction is controlled by the prevalence of tunneling over the thermal activation mechanism. In this temperature range, Bell's early treatment of quantum tunneling, based on a semiclassical approximation for the barrier permeability, provides a reliable tool to quantitatively account for the contribution of the tunneling effect. This treatment is here applied for extracting from rate constants properties of the effective tunneling path, such as the activation barrier height and width. We show that this is a way of parametrizing the dependence of the apparent activation energy on temperature useful for both calculated and experimental rate constants in an ample interval of temperature, from above to below Tc, relevant for modelization of astrophysical and in general very low-temperature environments.

Control of conformers combining cooling by supersonic expansion of seeded molecular beams with hexapole selection and alignment: experiment and theory on 2-butanol

Selection and alignment of rotamers and, more in general, of conformers in the gas phase is a challenge that we tackle experimentally by supersonic expansion of seeded molecular beams and hexapolar electrostatic fields with quadrupole mass detection. The studied system involves the nine conformers of the asymmetric-top molecule 2-butanol, which coexist because of nearly free rotations around a CC and a CO bond. From the measured time-of-flight of a 2-butanol supersonic molecular beam seeded in either He or Ar, the corresponding velocity distributions are obtained. The different nature and masses of the seeding gas decrease selectively the vibrational temperature and determine the population of the conformers, which is assessed on the basis of their statistical distribution, derived from high level accompanying quantum mechanical calculations. The use of a hexapolar electrostatic field permits us to induce a variation of the population distribution as a function of the applied voltage and of the selective focusing and alignment of the conformers. A technique, recently developed for treating asymmetric tops and involving extensive trajectory simulations, is applied to obtain the link between the focusing curves, i.e. the dependence of the beam intensity on the hexapole voltage, and the conformers' populations and alignment. Perspectives are provided for photo- and stereo-dynamics experiments, particularly appealing also on account that 2-butanol is the simplest chiral alcohol.

The He + H2+ → HeH+ + H reaction: ab initio studies of the potential energy surface, benchmark time-independent quantum dynamics in an extended energy range and comparison with experiments

In this work we critically revise several aspects of previous ab initio quantum chemistry studies [P. Palmieri et al., Mol. Phys. 98, 1835 (2000); C. N. Ramachandran et al., Chem. Phys. Lett. 469, 26 (2009)] of the HeH(2)(+) system. New diatomic curves for the H(2)(+) and HeH(+) molecular ions, which provide vibrational frequencies at a near spectroscopic level of accuracy, have been generated to test the quality of the diatomic terms employed in the previous analytical fittings. The reliability of the global potential energy surfaces has also been tested performing benchmark quantum scattering calculations within the time-independent approach in an extended interval of energies. In particular, the total integral cross sections have been calculated in the total collision energy range 0.955-2.400 eV for the scattering of the He atom by the ortho- and para-hydrogen molecular ion. The energy profiles of the total integral cross sections for selected vibro-rotational states of H(2)(+) (v = 0,...,5 and j = 1,...,7) show a strong rotational enhancement for the lower vibrational states which becomes weaker as the vibrational quantum number increases. Comparison with several available experimental data is presented and discussed.

Quantum chemistry of C(3)H(6)O molecules: structure and stability, isomerization pathways, and chirality changing mechanisms

Electronic structure calculations were carried out to study the various isomers of formula C(3)H(6)O, as a part of our current quantum chemical and dynamical approaches to intra- and intermolecular kinetics for the C(n)H(2n)O (n = 1, 2, 3) molecules. The usefulness of the GRRM (global reaction route mapping) program developed by Ohno and Maeda in predicting the structure of all isomers and of the transition states connecting them is fully exploited. All the isomers are identified as local minima on the MP2/CC-PVDZ potential energy surface. Acetone is the most stable isomer. In increasing order of stability the others are propanal, 2-propenol, 1-propenol, allyl alcohol, methyl vinyl ether, cyclopropanol, propylene oxide, and oxetane. Various isomerization paths connecting them are identified. All the transition states are fully characterized using intrinsic reaction coordinate calculations. The isomerization reactions may proceed through a single step or involve an intermediate species which is either a carbene or a diradical. Special attention is devoted to propylene oxide, a favorite molecule in current photochemical and stereodynamical studies because of its chiral nature. It is a rigid molecule, and chirality switching is found to be supported by its isomers. Two different chirality switching mechanisms which are assisted by propanal and allyl alcohol are presented.

Electrostatic hexapole state-selection of the asymmetric-top molecule propylene oxide

Rotational state-selection of the asymmetric-top molecule propylene oxide was carried out using an electrostatic hexapole field of 85-cm length. Molecular beam intensities were monitored by a quadrupole mass spectrometer. It was found that beam intensities of molecular beams for pure propylene oxide and those seeded in He and in Ar increased with increasing hexapole voltages. The hexapole voltage dependence of the beam intensity, which is called the focusing curve, was interpreted by computer simulation of the trajectories of molecules in the hexapolar field due to the Stark effect, as a function of rotational temperatures of molecular beams. The calculated best fit focusing curves, when compared with the experimental results, demonstrated that the rotational temperatures, associated with the distribution of states of a given rotational angular momentum J, are similar to the translational temperatures. It was found that the M = 0 states (where M is the projection of J along the direction of the electrostatic field) and negative values of the pseudoquantum number tau of propylene oxide can be selected using our experimental setup. These results suggest that the hexapole electric field is a tool even for the selection of rotational and orientation states of asymmetric-top molecules.