Anharmonic Vibrational Spectroscopy Calculations for Novel Rare-Gas-Containing Compounds: HXeH, HXeCl, HXeBr, and HXeOH

Complexes of HXeY with HX (Y, X = F, Cl, Br, I): Symmetry-Adapted Perturbation Theory Study and Anharmonic Vibrational Analysis

A comprehensive analysis of the intermolecular interaction energy and anharmonic vibrations of 41 structures of the HXeY⋯HX (X, Y = F, Cl, Br, I) family of noble-gas-compound complexes for all possible combinations of Y and X was conducted. New structures were identified, and their interaction energies were studied by means of symmetry-adapted perturbation theory, up to second-order corrections: this provided insight into the physical nature of the interaction in the complexes. The energy components were discussed, in connection to anharmonic frequency analysis. The results show that the induction and dispersion corrections were the main driving forces of the interaction, and that their relative contributions correlated with the complexation effects seen in the vibrational stretching modes of Xe-H and H-X. Reasonably clear patterns of interaction were found for different structures. Our findings corroborate previous findings with better methods, and provide new data. These results suggest that the entire group of the studied complexes can be labelled as "naturally blueshifting", except for the complexes with HI.

My Trajectory in Molecular Reaction Dynamics and Spectroscopy

This is the story of a career in theoretical chemistry during a time of dramatic changes in the field due to phenomenal growth in the availability of computational power. It is likewise the story of the highly gifted graduate students and postdoctoral fellows that I was fortunate to mentor throughout my career. It includes reminiscences of the great mentors that I had and of the exciting collaborations with both experimentalists and theorists on which I built much of my research. This is an account of the developments of exciting scientific disciplines in which I was involved: vibrational spectroscopy, molecular reaction mechanisms and dynamics, e.g., in atmospheric chemistry, and the prediction of new, exotic molecules, in particular noble gas molecules. From my very first project to my current work, my career in science has brought me the excitement and fascination of research. What a wonderful pursuit!

A theoretical study on the infrared signatures of proton-bound rare gas dimers (Rg-H+-Rg), Rg = {Ne, Ar, Kr, and Xe}

The infrared spectrum of proton-bound rare gas dimers has been extensively studied via matrix isolation spectroscopy. However, little attention has been paid on their spectrum in the gas phase. Most of the Rg₂H⁺ has not been detected outside the matrix environment. Recently, Ar_nH⁺ (n = 3-7) has been first detected in the gas-phase [D. C. McDonald et al., J. Chem. Phys. 145, 231101 (2016)]. In that work, anharmonic theory can reproduce the observed vibrational structure. In this paper, we extend the existing theory to examine the vibrational signatures of Rg₂H⁺, Rg = {Ne, Ar, Kr, and Xe}. The successive binding of Rg to H⁺ was investigated through the calculation of stepwise formation energies. It was found that this binding is anti-cooperative. High-level full-dimensional potential energy surfaces at the CCSD(T)/aug-cc-pVQZ//MP2/aug-cc-pVQZ were constructed and used in the anharmonic calculation via discrete variable representation. We found that the potential coupling between the symmetric and asymmetric Rg-H⁺ stretch (ν₁ and ν₃ respectively) causes a series of bright n₁ν₁ + ν₃ progressions. From Ne₂H⁺ to Xe₂H⁺, an enhancement of intensities for these bands was observed.

Computational Structures and SAPT Interaction Energies of HXeSH···H2Y (Y=O or S) Complexes

Ab initio calculations of the structures, vibrational spectra and supermolecular and symmetry-adapted perturbation theory (SAPT) interaction energies of the HXeOH and HXeSH complexes with H2O and H2S molecules are presented. Two minima already reported in the literature were reproduced and ten new ones were found together with some transition states. All complexes show blue shift in Xe–H stretching mode upon complexation. The computed spectra suggest that it should be possible to detect and distinguish the complexes experimentally. The structures where H2O or H2S is the proton-donor were found to be the most stable for all complex compositions. The SAPT analysis shows significant differences between the complexes with H2O and H2S indicating much larger dispersion and exchange contributions in the complexes with H2S.

A Noble-Gas Hydride in a Nitrogen Medium: Structure, Spectroscopy, and Intermolecular Vibrations of HXeBr@(N2)22

Noble-gas hydrides have been extensively studied in noble gas matrices. However, little is known on their stability and properties in molecular hosts. Here, HXeBr in the N2 environment is modeled at the B3LYP-D level of theory in a complete single shell of 22 N2 molecules. The system is compared to similar models of HXeBr in CO2 and Xe clusters. The optimized structure of (HXeBr)@(N2)22 is of low symmetry and is highly anisotropic. None of the N2 molecules are freely rotating, and the host molecules are not symmetrically positioned with respect to the HXeBr axis. The axes of the N2 molecules are nonuniformly distributed. The computed anharmonic H-Xe stretching frequency of HXeBr in the N2 cluster is in good accord with the experimental value. The soft-mode frequencies of the cluster including both intermolecular vibrations and librations, have a broad distribution that ranges from 8.7 to 107 cm(-1). It is expected that these findings and specifically, the single-shell model, may shed light also on the local structure and vibrations of other impurities in a molecular media.

Matrix site effects on vibrational frequencies of HXeCCH, HXeBr, and HXeI: a hybrid quantum-classical simulation

The matrix shifts of the H-Xe stretching frequency of noble-gas hydrides, HXeCCH, HXeBr, and HXeI in various noble-gas matrices (in Ne, Ar, Kr, and Xe matrices) are investigated via the hybrid quantum-classical simulations. The order of the H-Xe stretching frequencies is found to be ν(gas) < ν(Ne) < ν(Xe) < ν(Kr) < ν(Ar) for HXeCCH and HXeBr, while it is ν(gas) < ν(Ne) < ν(Xe) < ν(Ar) < ν(Kr) for HXeI. This order is anomalous with respect to the matrix dielectric constants, and the calculated results reproduce the experimentally observed shifts quite successfully. We also find that the matrix shifts from the gas-phase values are Δν(HXeCCH) ≈ Δν(HXeCl) < Δν(HXeBr) < Δν(HXeI) in the same noble-gas matrix environments, which implies that the weakly bound molecules exhibit large matrix shifts. The local trapping site is analyzed in detail, and it is shown that a realistic modeling of the surrounding matrix environments is essential to describe the unusual matrix shifts accurately.

On the vibrational linear and nonlinear optical properties of compounds involving noble gas atoms: HXeOXeH, HXeOXeF, and FXeOXeF

The vibrational (hyper)polarizabilities of some selected Xe derivatives are studied in the context of Bishop-Kirtman perturbation theory (BKPT) and numerical finite field methodology. It was found that for this set of rare gas compounds, the static vibrational properties are quite large, in comparison to the corresponding electronic ones, especially those of the second hyperpolarizability. This also holds for the dc-Pockels β(-ω;ω,0), Kerr γ(-ω;ω,0,0) and electric field second harmonic generation γ (-2ω;ω,ω,0) effects, although the computed nuclear relaxation (nr) vibrational contributions are smaller in magnitude than the static ones. HXeOXeH was used to study the effects of electron correlation, basis set, and geometry. Geometry effects were found to lead to noticeable changes of the vibrational and electronic second hyperpolarizability. A limited study of the effect of Xe insertion to the nr vibrational properties is also reported. Assessment of the results revealed that Xe insertion has a remarkable effect on the nr (hyper)polarizabilities. In terms of the BKPT, this is associated with a remarkable increase of the electrical and mechanical anharmonicity terms. The latter is consistent with the anharmonic character of several vibrational modes reported for rare gas compounds.

THEORETICAL DETERMINATION OF ROVIBRATIONAL ENERGIES AND ANOMALOUS ISOTOPIC EFFECT OF WEAKLY BOUND CLUSTER HXeOH

The rovibrational energy levels of HXeOH and its isotopomers have been calculated using a two-layer Lanczos algorithm in Radau-diatom-Jacobi coordinates based on a high level ab initio potential energy surface. The surface is obtained by fitting to 1229 RCCSD(T)/SDB-cc-pVQZ energy points. The equilibrium geometry of HXeOH is determined to have a trans configuration with a nearly collinear HXeO bond angle of 177.32 degrees. The well depth of this minimum is only 0.6123 eV with respect to the OH + Xe + H dissociation limit. The results show that the H–Xe stretching frequency of HXeOH and HXeOD has an anomalous isotopic shift upon 18 O isotope substitution, whereas the D–Xe stretch in DXeOH and DXeOD displays a normal isotopic shift. This trend is consistent with the experimental observations in a Xe solid matrix. The present results predict a lower frequency for the H–Xe stretch than was observed in the solid matrix experiment. Either these results are too low or there is a strong blue shift due to the matrix in the experimental values.

On the origin of the large electron correlation contribution to the hyperpolarizabilities of some diacetylene rare gas compounds

A comprehensive study of the linear and nonlinear molecular optical properties of HRgC(4)H, where Rg = Ar,Kr,Xe, has been performed. Dynamical electron correlation effects were computed by employing the coupled cluster methodology. A large electron correlation contribution to the nonlinear properties of HArC(4)H has been revealed. This contribution decreases by increasing the atomic number of the inserted rare gas atom. In order to interpret the origin of this noteworthy property, the complete active space self-consistent field method was employed. We have performed a systematic study of the linear and nonlinear electric properties by modifying the active space. The calculations have shown the significant contribution of the doubly excited sigma(*2) configuration and a negligible contribution of pi(*2). A quite remarkable discrepancy between numerically and analytically evaluated hyperpolarizabilities has also been observed for HArC(4)H. This was attributed to the contribution of near degenerate states.

Formation of novel rare-gas molecules in low-temperature matrices

Progress in the study of a new class of chemically bound compounds of noble-gas atoms is reviewed. The focus is on rare-gas molecules of the form HNgY, where Ng is a noble-gas atom and Y is an electronegative group, prepared by photolysis of HY in the rare-gas matrix. Other related types of new molecules of noble-gas atoms are discussed as well. Topics discussed in this review include: (a) The nature of bonding and the energetic stability of the compounds. (b) The vibrational spectroscopy of the molecules, and its role in identification of the species. (c) The mechanism and dynamics of photochemical formation of HNgY in the matrix, and the pathways for thermal and infrared (IR)-induced decomposition. Specifically, attention is given to the issue of "direct" formation following photolysis of HY versus "delayed" formation involving H atom diffusion. (d) Molecules of the lighter rare gases Ar, Ne, and He, focusing on the experimentally prepared HArF and on theoretical predictions suggesting the existence of other molecules. (e) The most-recently discovered photochemically induced insertion compounds of Ng into hydrocarbons, such as HXeCCH. (f) Clusters of HNgY with other molecules. The possible existence of neat aggregates and crystals of HNgY. The reviewed state-of-the-art suggests this field is at an early stage of development with major open questions bearing on the surprising properties of the molecules and on the formation mechanisms. These are part of the challenge for the future.

The Dipole Moment, Polarizabilities, and First Hyperpolarizabilities of HArF. A Computational and Comparative Study

The electronic and vibrational contributions to the dipole moment, polarizabilities, and first hyperpolarizabilities of HArF are reported. These have been computed by using a series of systematically built basis sets and a hierarchy of computational methods. HArF has a very large first hyperpolarizability along the z axis. This has been rationalized by invoking the difference in the electronic structure between the ground and the first excited state. The argon fluorohydride has been recently derived and characterized. The present study provides complementary data for the understanding of the electronic structure of this interesting argon derivative.

An ab initio potential energy surface and predissociative resonances of HArF

A three-dimensional potential energy surface of the ground electronic state HArF is constructed from more than 2000 ab initio points at the multireference averaged quadratic coupled-cluster level employing an augmented large basis set. The calculations indicate that the linear HArF molecule is metastable with a barrier of 0.643 eV in the atomization (HArF --> H + Ar + F) channel and a barrier of 1.017 eV in the dissociation (HArF --> Ar + HF) channel. Variational calculations of low-lying predissociative resonances of both HArF and DArF are performed on the three-dimensional potential energy surface using a complex-symmetric Lanczos propagation method, which yields both positions and widths of the resonance states. The resonance lifetime generally decreases with energy, but strong mode selectivity exists. Reasonably good agreement with experiment confirms the accuracy of our potential. These calculations provide valuable information on the stability and dynamics of HArF/DArF in its ground electronic state.

Organo-noble-gas hydride compounds HKrCCH, HXeCCH, HXeCC, and HXeCCXeH: Formation mechanisms and effect of [sup 13]C isotope substitution on the vibrational properties

We investigate the formation mechanism of HXeCCXeH in a Xe matrix. Our experimental results show that the HXeCCXeH molecules are formed in the secondary reactions involving HXeCC radicals. The experimental data on the formation of HXeCCXeH is fully explained based on the model involving the HXeCC+Xe+H-->HXeCCXeH reaction. This reaction is the first case when a noble-gas hydride molecule is formed from another noble-gas molecule. In addition, we investigate the (12)C/(13)C isotope effect on the vibrational properties of organo-noble-gas hydrides (HKrCCH, HXeCCH, HXeCC, and HXeCCXeH) in noble-gas matrixes. The present experimental results and ab initio calculations on carbon isotope shifts of the vibrational modes support the previous assignments of these molecules. Upon (12)C to (13)C isotope substitution, we observed a pronounced effect on the H-Kr stretching mode of HKrCCH (downshift of 1.0-3.6 cm(-1), depending on the matrix site) and a small anomalous shift (+0.1 cm(-1)) of the H-Xe stretching mode of HXeCCH and HXeCCXeH.

Chemical Compounds Formed from Diacetylene and Rare-Gas Atoms: HKrC4H and HXeC4H

New organic rare-gas compounds, HRgC4H (Rg = Kr or Xe), are identified in matrix-isolation experiments supported by ab initio calculations. These compounds are the largest molecules among the known rare-gas hydrides. They are prepared in low-temperature rare-gas matrixes via UV photolysis of diacetylene and subsequent thermal mobilization of H atoms at approximately 30 and 45 K for Kr and Xe, respectively. The strongest IR absorption bands of the HRgC4H molecules are the H-Rg stretches with the most intense components at 1290 cm(-1) for HKrC4H and at 1532 cm(-1) for HXeC4H, and a number of weaker absorptions (C-H stretching, C-C-C bending, and C-C-H bending modes) are also found in agreement with the theoretical predictions. As probably the most important result, the IR absorption spectra indicate some further stabilization of the HRgC4H molecules as compared with the corresponding HRgC2H species identified recently (Khriachtchev et al. J. Am. Chem. Soc. 2003, 125, 4696 and Khriachtchev et al. J. Am. Chem. Soc. 2003, 125, 6876). The computational energetic results support this trend. HXeC4H is predicted to be 2.5 eV lower in energy than H + Xe + C4H, which is approximately 1 eV larger than the corresponding value for HXeC2H. We expect that the larger molecules HRgC6H and HRgC8H are even more stable and the HRgC2nH species are good candidates for bulk organic rare-gas material.

Noble gas-actinide complexes of the CUO molecule with multiple Ar, Kr, and Xe atoms in noble-gas matrices

Laser-ablated U atoms react with CO in excess argon to produce CUO, which is trapped in a triplet state in solid argon at 7 K, based on agreement between observed and relativistic density functional theory (DFT) calculated isotopic frequencies ((12)C(16)O, (13)C(16)O, (12)C(18)O). This observation contrasts a recent neon matrix investigation, which trapped CUO in a linear singlet state calculated to be about 1 kcal/mol lower in energy. Experiments with krypton and xenon give results analogous to those with argon. Similar work with dilute Kr and Xe in argon finds small frequency shifts in new four-band progressions for CUO in the same triplet states trapped in solid argon and provides evidence for four distinct CUO(Ar)(4-n)(Ng)(n) (Ng = Kr, Xe, n = 1, 2, 3, 4) complexes for each Ng. DFT calculations show that successively higher Ng complexes are responsible for the observed frequency progressions. This work provides the first evidence for noble gas-actinide complexes, and the first example of neutral complexes with four noble gas atoms bonded to one metal center.

Noble gas-actinide compounds: evidence for the formation of distinct CUO(Ar)(4-n)(Xe)(n) and CUO(Ar)(4-n)(Kr)(n) (n = 1, 2, 3, 4) complexes

Laser-ablated U atoms react with CO in excess argon to produce CUO, which gives rise to 852.5 and 804.3 cm-1 infrared absorptions for the triplet state CUO(Ar)n complex in solid argon at 7 K. Relativistic density functional calculations show that the CUO(Ar) complex is stable and that up to four or five argon atoms can complex to CUO. When 1-3% Xe is added to the argon/CO reagent mixture, strong absorptions appear at 848.0 and 801.3 cm-1 and dominate new four-band progressions, which increase on annealing to 35-50 K as Xe replaces Ar in the intimate coordination sphere. Analogous spectra are obtained with 1-2% Kr added. This work provides evidence for eight distinct CUO(Ng)n(Ar)4-n (Ng = Kr, Xe, n = 1, 2, 3, 4) complexes and the first characterization of neutral complexes involving four noble-gas atoms on one metal center.

Noble gas-actinide compounds: complexation of the CUO molecule by Ar, Kr, and Xe atoms in noble gas matrices

The CUO molecule, formed from the reaction of laser-ablated U atoms with CO in a noble gas, exhibits very different stretching frequencies in a solid argon matrix [804.3 and 852.5 wave numbers (cm(-1))] than in a solid neon matrix (872.2 and 1047.3 cm(-1)). Related experiments in a matrix consisting of 1% argon in neon suggest that the argon atoms are interacting directly with the CUO molecule. Relativistic density functional calculations predict that CUO can bind directly to one argon atom (U-Ar = 3.16 angstroms; binding energy = 3.2 kilocalories per mole), accompanied by a change in the ground state from a singlet to a triplet. Our experimental and theoretical results also suggest that multiple argon atoms can bind to a single CUO molecule.