Matrix isolation spectroscopy of H2O, D2O, and HDO in solid parahydrogen

Fast ortho-to-para conversion of molecular hydrogen in chemisorption and matrix-isolation systems

Molecular hydrogen has two nuclear-spin modifications called ortho and para. Because of the symmetry restriction with respect to permutation of the two protons, the ortho and para isomers take only odd and even values of the rotational quantum number, respectively. The ortho-to-para conversion is promoted in condensed systems, to which the excess rotational energy and spin angular momentum are transferred. We review recent studies on fast ortho-to-para conversion of hydrogen in molecular chemisorption and matrix isolation systems, discussing the conversion mechanism as well as rotational-relaxation pathways.

Ortho-para interconversion of nuclear states of H2O through replica transition state: prospect of quantum entanglement at homodromic Bjerrum defect site

CONTEXT

From a nuclear spin prospective, water exists as para and ortho nuclear spin isomers (isotopomers). Spin interconversions in isolated molecules of water are forbidden, but many recent reports have shown them to happen in bulk, through dynamic proton exchanges happening between interconnected networks of a large array of water molecules. In this contribution, a possible explanation for an unexpected slow or delayed interconversion of ortho-para water in ice observed in an earlier reported experiment is provided. Using the results of quantum mechanical investigations, we have discussed the roles played by Bjerrum defects in the dynamic proton exchanges and ortho-para spin state interconversions. We guess that at the sites of the Bjerrum defects, there are possibilities of quantum entanglements of states, through pairwise interactions. Based on the perfectly correlated exchange happening via a replica transition state, we speculate that it can have significant influences on ortho-para interconversions of water. We also conjecture that the overall ortho-para interconversion is not a continuous process, rather can be imagined to be happening serendipitously, but within the boundary of the rules of quantum mechanics.

METHODS

All computations were performed with Gaussian 09 program. B3LYP/6-31++G(d,p) methodology was used to compute all the stationary points. Further energy corrections were computed using CCSD(T)/aug-cc-pVTZ methodology. Intrinsic reaction coordinate (IRC) path computations were carried out for the transition states.

Infrared Spectroscopic Studies of Oxygen Atom Quantum Diffusion in Solid Parahydrogen

The thermally induced diffusion of atomic species in noble gas matrices was studied extensively in the 1990s to investigate low-temperature solid-state reactions and to synthesize reactive intermediates. In contrast, much less is known about the diffusion of atomic species in quantum solids such as solid parahydrogen (p-H₂). While hydrogen atoms were shown to diffuse in normal-hydrogen solids at 4.2 K as early as 1989, the diffusion of other atomic species in solid p-H₂ has not been reported in the literature. The in situ photogeneration of atomic oxygen, by ArF laser irradiation of an O₂-doped p-H₂ solid at 193 nm, is studied here to investigate the diffusion of O(³P) atoms in a quantum solid. The O(³P) atom mobility is detected by measuring the kinetics of the O(³P) + O₂ → O₃ reaction after photolysis via infrared spectroscopy of the O₃ reaction product. This reaction is barrierless and is thus assumed to be diffusion-controlled under these conditions such that the reaction rate constant can be used to estimate the oxygen atom diffusion coefficient. The O₃ growth curves are well fit by single exponential expressions allowing the pseudo-first-order rate constant for the O(³P) + O₂ → O₃ reaction to be extracted. The reaction rates are affected strongly by the p-H₂ crystal morphology and display a non-Arrhenius-type temperature dependence consistent with quantum diffusion of the O(³P) atom. The experimental results are compared to H(²S) atom reaction studies in p-H₂, analogous studies in noble gas matrices, and laboratory studies of atomic diffusion in astronomical ices and surfaces.

Formation of Para-H2O by Vacuum-UV Photolysis of O2 in Solid Hydrogen: Implication for Astrochemistry

The observation that the ortho to para ratio (OPR) of interstellar H₂O is smaller than 3 is an important yet unresolved subject in astronomy. We irradiated O₂ embedded in solid H₂ at 3 K with vacuum-ultraviolet (VUV) light and observed IR lines associated with para-H₂O (denoted as pH₂O) and nonrotating H₂O-(oH₂)_n (where oH₂ denotes ortho-H₂) but no lines associated with ortho-H₂O (denoted as oH₂O). After maintaining the matrix in darkness for ∼30 h, the amount of pH₂O decreased, accompanied by an increase in H₂O-(oH₂)_n via diffusion of oH₂. After that, the continuous nuclear-spin conversion from oH₂ to para-H₂ (denoted as pH₂) in solid H₂ over time resulted in the conversion of nonrotating H₂O-(oH₂)_n to rotating pH₂O in solid pH₂. The observation of the formation and conversion of pH₂O in our experiment suggests a plausible route in which VUV irradiation of O₂ and H₂ adsorbed on grain surfaces might be responsible for the smaller OPR of interstellar H₂O.

High resolution infrared spectroscopy of (HCl)2 and (DCl)2 isolated in solid parahydrogen: Interchange-tunneling in a quantum solid

Infrared spectroscopic studies of weakly bound clusters isolated in solid parahydrogen (pH₂) that exhibit large-amplitude tunneling motions are needed to probe how quantum solvation perturbs these types of coherent dynamics. We report high resolution Fourier transform infrared absorption spectra of (HCl)₂, HCl-DCl, and (DCl)₂ isolated in solid pH₂ in the 2.4-4.8 K temperature range. The (HCl)₂ spectra show a remarkable amount of fine structures that can be rigorously assigned to vibration-rotation-tunneling transitions of (HCl)₂ trapped in double substitution sites in the pH₂ matrix where end-over-end rotation of the cluster is quenched. The spectra are assigned using a combination of isotopically (H/D and ³⁵Cl/³⁷Cl) enriched samples, polarized IR absorption measurements, and four-line combination differences. The interchange-tunneling (IT) splitting in the ground vibrational state for in-plane and out-of-plane H³⁵Cl-H³⁷Cl dimers is 6.026(1) and 6.950(1) cm^-1, respectively, which are factors of 2.565 and 2.224 smaller than in the gas phase dimer. In contrast, the (DCl)₂ results show larger perturbations where the ground vibrational state IT splitting in D³⁵Cl-D³⁷Cl is 1.141(1) cm^-1, which is a factor of 5.223 smaller than in the gas phase, and the tunneling motion is quenched in excited intramolecular vibrational states. The results are compared to similar measurements on (HCl)₂ made in liquid helium nanodroplets to illustrate the similarities and differences in how both these quantum solvents interact with large amplitude tunneling motions of an embedded chromophore.

Infrared spectroscopy of an endohedral water in fullerene

An infrared absorption spectroscopy study of the endohedral water molecule in a solid mixture of H₂O@C₆₀ and C₆₀ was carried out at liquid helium temperature. From the evolution of the spectra during the ortho-para conversion process, the spectral lines were identified as para-H₂O and ortho-H₂O transitions. Eight vibrational transitions with rotational side peaks were observed in the mid-infrared: ω₁, ω₂, ω₃, 2ω₁, 2ω₂, ω₁ + ω₃, ω₂ + ω₃, and 2ω₂ + ω₃. The vibrational frequencies ω₂ and 2ω₂ are lower by 1.6% and the rest by 2.4%, as compared to those of free H₂O. A model consisting of a rovibrational Hamiltonian with the dipole and quadrupole moments of H₂O interacting with the crystal field was used to fit the infrared absorption spectra. The electric quadrupole interaction with the crystal field lifts the degeneracy of the rotational levels. The finite amplitudes of the pure v₁ and v₂ vibrational transitions are consistent with the interaction of the water molecule dipole moment with a lattice-induced electric field. The permanent dipole moment of encapsulated H₂O is found to be 0.50 ± 0.05 D as determined from the far-infrared rotational line intensities. The translational mode of the quantized center-of-mass motion of H₂O in the molecular cage of C₆₀ was observed at 110 cm^-1 (13.6 meV).

VUV photochemistry and nuclear spin conversion of water and water-orthohydrogen complexes in parahydrogen crystals at 4 K

Samples of H2O, HDO, and D2O were isolated in solid parahydrogen (pH2) matrices and irradiated by vacuum ultraviolet (VUV) radiation at 147 nm. Fourier-Transform Infrared (FTIR) spectra showed a clear depletion of D2O and an enrichment of both HDO and H2O by 147 nm irradiation. These irradiation-dependent changes are attributed to the production of OH and/or OD radicals through photodissociations of H2O, HDO, and D2O. The radicals subsequently react with the hydrogen matrix, leading to the observed enrichment of H2O. No trace of isolated OH or OD was detected in the FTIR spectra, indicating that the OH/OD radicals react with the surrounding matrix hydrogen molecules via quantum tunneling within our experimental timescale. The observed temporal changes in concentrations, especially the increase of HDO concentration during VUV irradiation, can be interpreted by a model with a rapid conversion from orthohydrogen (oH2) to pH2 in water-oH2 complexes upon VUV photodissociation, indicating either the acceleration of the nuclear spin conversion (NSC) of H2 due to the magnetic moment of the intermediate OH/OD radical, or the preferential reaction of the OH/OD radical with a nearby oH2 molecule over other pH2 molecules. We have also identified and quantified an anomalously slow NSC of H2O and D2O complexed with oH2 in solid pH2.

Hydrogen atom quantum diffusion in solid parahydrogen: The H + N2O → cis-HNNO → trans-HNNO reaction

The diffusion and reactivity of hydrogen atoms in solid parahydrogen at temperatures between 1.5 K and 4.3 K are investigated by high-resolution infrared spectroscopy. Hydrogen atoms are produced within solid parahydrogen as the by-products of the 193 nm in situ photolysis of N₂O, which induces a two-step tunneling reaction, H + N₂O → cis-HNNO → trans-HNNO. The second-order rate constant for the first step to form cis-HNNO is found to be inversely proportional to the N₂O concentration after photolysis, indicating that the hydrogen atoms move through solid parahydrogen via quantum diffusion. This reaction only readily occurs at temperatures below 2.8 K, not due to an increased rate constant for the first reaction step at low temperatures but rather due to an increased selectivity to the reaction. The rate constant for the second step of the reaction mechanism involving unimolecular isomerization is shown to be independent of the N₂O concentration as expected. The inverse concentration dependence of the rate constant for the reaction step that involves the hydrogen atom demonstrates clearly that quantum diffusion influences the reactivity of the hydrogen atoms in solid parahydrogen, which does not have an analogy in classical reaction kinetics.

Matrix Isolation Spectroscopy and Nuclear Spin Conversion of Propyne Suspended in Solid Parahydrogen

Parahydrogen (pH₂) quantum solids are excellent matrix isolation hosts for studying the rovibrational dynamics and nuclear spin conversion (NSC) kinetics of molecules containing indistinguishable nuclei with nonzero spin. The relatively slow NSC kinetics of propyne (CH₃CCH) isolated in solid pH₂ is employed as a tool to assign the rovibrational spectrum of propyne in the 600-7000 cm^-1 region. Detailed analyses of a variety of parallel (ΔK = 0) and perpendicular (ΔK=±1) bands of propyne indicate that the end-over-end rotation of propyne is quenched, but K rotation of the methyl group around the C₃ symmetry axis still persists. However, this single-axis K rotation is significantly hindered for propyne trapped in solid pH₂ such that the energies of the K rotational states do not obey simple energy-level expressions. The NSC kinetics of propyne follows first-order reversible kinetics with a 287(7) min effective time constant at 1.7 K. Intensity-intensity correlation plots are used to determine the relative line strengths of individual ortho- and para-propyne rovibrational transitions, enabling an independent estimation of the ground vibrational state effective A″ constant of propyne.

THz spectroscopy of weakly bound cluster molecules in solid para-hydrogen: a sensitive probe of van der Waals interactions

The present work demonstrates that 99.9% enriched solid para-H2 below 3 K provides an excellent inert and transparent medium for the exploration of large-amplitude intermolecular vibrational motion of weakly bound van der Waals cluster molecules in the THz spectral region. THz absorption spectra have been generated for CO2/H2O and CS2/H2O mixtures embedded in enriched solid para-H2 and numerous observed transitions associated with large-amplitude librational motion of the weakly bound binary CO2H2O and CS2H2O van der Waals cluster molecules have been assigned together with tentative assignments for the ternary CS2(H2O)2 system. The interaction strength, directionality and anharmonicity of the weak van der Waals "bonds" between the molecules can be characterized via these THz spectral signatures and yield rigorous benchmarks for high-level ab initio methodologies. It is suggested that even a less stable linear conformation of the ternary CS2(H2O)2 system, where one H2O molecule is linked to each S atom of the CS2 subunit, may be formed due to the kinetics associated with the mobility of free H2O molecules in the soft para-H2 medium. In addition, the spectroscopic observations confirm a linear and planar global intermolecular potential energy minimum for the binary CS2H2O system with C2v symmetry, where the O atom on the H2O molecule is linked to one of the S atoms on the CS2 subunit. A semi-experimental value for the vibrational zero-point energy contribution of 1.93 ± 0.10 kJ mol-1 from the class of large-amplitude intermolecular vibrational modes is proposed. The combination with CCSD(T)/CBS electronic energy predictions provides a semi-experimental estimate of 5.08 ± 0.15 kJ mol-1 for the binding energy D0 of the CS2H2O van der Waals system.

Large amplitude motions within molecules trapped in solid parahydrogen

Molecules of the β-diketone and β-dialdehyde families were trapped in solid parahydrogen (pH2) to investigate the vibrational behavior of systems containing an intramolecular hydrogen bond (IHB). In the simplest β-diketone, acetylacetone (AcAc), H transfer related to the IHB is coupled with methyl torsions. In pH2, the study of nuclear spin conversion (NSC) in methyl groups allows the characterisation of the influence of these large amplitude motions on the vibrational modes. The deuteration of the OH group involved in the IHB has important consequences on the vibrational spectrum of the molecule and evidence of NSC in methyl groups is difficult to obtain. In the chlorine derivative (3-chloroacetylacetone), the H-transfer is no longer coupled with methyl torsion, and NSC has undetectable effects on the IR spectrum. A search of H tunnelling splitting in the IR spectra of β-dialdehydes trapped in pH2 was performed. A few modes of 2-chloromalonaldehyde appear as doublets and were assigned to tunnelling levels. The spectroscopic results related to large amplitude motions are detailed and discussed, highlighting puzzling effects.

Diffusion of Water Molecules in Quantum Crystals

With the use of solid parahydrogen in matrix isolation spectroscopy becoming more commonplace over the past few decades, it is increasingly important to understand the behavior of molecules isolated in this solid. The mobility of molecules in solid parahydrogen can play an important role in the dynamics of the system. Water molecules embedded in solid parahydrogen as deposited were found to be mobile at 4.0 K on the time scale of a few days. The diffusion at this temperature must be due to quantum tunneling in solid parahydrogen. The diffusion dynamics were analyzed based on the theory of nucleation. The concentration dependence on the diffusion rate indicates that there might be correlated motion of water molecules, a signature of quantum diffusion. We find that both water monomers and water dimers migrate in solid parahydrogen and provide insight into the behavior of molecules embedded in this quantum crystal.

Signatures of a quantum diffusion limited hydrogen atom tunneling reaction

We are studying the details of hydrogen atom (H atom) quantum diffusion in parahydrogen quantum solids in an effort to better understand H atom transport and reactivity under these conditions.

Infrared Studies of the Symmetry Changes of the 28SiH4 Molecule in Low-Temperature Matrixes. Fundamental, Combination, and Overtone Transitions

Infrared spectra of ²⁸SiH₄ in argon and nitrogen matrixes at low temperature 6.5-20 K in the region of overtone and combination transitions were recorded for the first time. Additionally, the high-resolution spectra were obtained in the fundamental region. The frequencies and the relative intensities of all bands were determined. The set of experimental data suggests that the symmetry of molecules studied in the matrixes is different from the symmetry of the free molecules because of an interaction with the environment. The symmetry of ²⁸SiH₄ changes from T_d to C_3v on transition from the gas phase to a nitrogen matrix and to D_2d on transition to an argon matrix. A modeling of SiH₄ molecule force fields explains the experimental data as a change of a force constant of the selected SiH bond in the case of SiH₄ in the nitrogen matrix or force constants of two opposite angles in the case of SiH₄ in the argon matrix. In spite of small values of these changes, they result in noticeable spectroscopic effects: the band splitting and appearance of new bands in matrix spectra compared with spectra of free SiH₄. The interpretation of transitions in the fundamental and combination regions was performed.

Collisional cross-section of water molecules in vapour studied by means of 1H relaxation in NMR

In gas phase, collisions that affect the rotational angular momentum lead to the return of the magnetization to its equilibrium (relaxation) in Nuclear Magnetic Resonance (NMR). To the best of our knowledge, the longitudinal relaxation rates R₁ = 1/T₁ of protons in H₂O and HDO have never been measured in gas phase. We report R₁ in gas phase in a field of 18.8 T, i.e., at a proton Larmor frequency ν₀ = 800 MHz, at temperatures between 353 and 373 K and pressures between 9 and 101 kPa. By assuming that spin rotation is the dominant relaxation mechanism, we estimated the effective cross-section σ_J for the transfer of angular momentum due to H₂O-H₂O and HDO-D₂O collisions. Our results allow one to test theoretical predictions of the intermolecular potential of water in gas phase.

Challenges in preparing, preserving and detecting para-water in bulk: overcoming proton exchange and other hurdles

Para-water is an analogue of para-hydrogen, where the two proton spins are in a quantum state that is antisymmetric under permutation, also known as singlet state. The populations of the nuclear spin states in para-water are believed to have long lifetimes just like other Long-Lived States (LLSs). This hypothesis can be verified by measuring the relaxation of an excess or a deficiency of para-water, also known as a "Triplet-Singlet Imbalance" (TSI), i.e., a difference between the average population of the three triplet states T (that are symmetric under permutation) and the population of the singlet state S. In analogy with our recent findings on ethanol and fumarate, we propose to adapt the procedure for Dissolution Dynamic Nuclear Polarization (D-DNP) to prepare such a TSI in frozen water at very low temperatures in the vicinity of 1.2 K. After rapid heating and dissolution using an aprotic solvent, the TSI should be largely preserved. To assess this hypothesis, we studied the lifetime of water as a molecular entity when diluted in various solvents. In neat liquid H2O, proton exchange rates have been characterized by spin-echo experiments on oxygen-17 in natural abundance, with and without proton decoupling. One-dimensional exchange spectroscopy (EXSY) has been used to study proton exchange rates in H2O, HDO and D2O mixtures diluted in various aprotic solvents. In the case of 50 mM H2O in dioxane-d8, the proton exchange lifetime is about 20 s. After dissolving, one can observe this TSI by monitoring intensities in oxygen-17 spectra of H2O (if necessary using isotopically enriched samples) where the AX2 system comprising a "spy" oxygen A and two protons X2 gives rise to binomial multiplets only if the TSI vanishes. Alternatively, fast chemical addition to a suitable substrate (such as an activated aldehyde or ketone) can provide AX2 systems where a carbon-13 acts as a spy nucleus. Proton signals that relax to equilibrium with two distinct time constants can be considered as a hallmark of a TSI. We optimized several experimental procedures designed to preserve and reveal dilute para-water in bulk.

Heavy snow: IR spectroscopy of isotope mixed crystalline water ice

Temperature and isotopic dependence of simulated and experimental spectra shed light on the vibrational modes of crystalline water ice.

Quantum Diffusion-Controlled Chemistry: Reactions of Atomic Hydrogen with Nitric Oxide in Solid Parahydrogen

Our group has been working to develop parahydrogen (pH2) matrix isolation spectroscopy as a method to study low-temperature condensed-phase reactions of atomic hydrogen with various reaction partners. Guided by the well-defined studies of cold atom chemistry in rare-gas solids, the special properties of quantum hosts such as solid pH2 afford new opportunities to study the analogous chemical reactions under quantum diffusion conditions in hopes of discovering new types of chemical reaction mechanisms. In this study, we present Fourier transform infrared spectroscopic studies of the 193 nm photoinduced chemistry of nitric oxide (NO) isolated in solid pH2 over the 1.8 to 4.3 K temperature range. Upon short-term in situ irradiation the NO readily undergoes photolysis to yield HNO, NOH, NH, NH3, H2O, and H atoms. We map the postphotolysis reactions of mobile H atoms with NO and document first-order growth in HNO and NOH reaction products for up to 5 h after photolysis. We perform three experiments at 4.3 K and one at 1.8 K to permit the temperature dependence of the reaction kinetics to be quantified. We observe Arrhenius-type behavior with a pre-exponential factor of A = 0.036(2) min(-1) and Ea = 2.39(1) cm(-1). This is in sharp contrast to previous H atom reactions we have studied in solid pH2 that display definitively non-Arrhenius behavior. The contrasting temperature dependence measured for the H + NO reaction is likely related to the details of H atom quantum diffusion in solid pH2 and deserves further study.

Nuclear spin conversion of water inside fullerene cages detected by low-temperature nuclear magnetic resonance

The water-endofullerene H2O@C60 provides a unique chemical system in which freely rotating water molecules are confined inside homogeneous and symmetrical carbon cages. The spin conversion between the ortho and para species of the endohedral H2O was studied in the solid phase by low-temperature nuclear magnetic resonance. The experimental data are consistent with a second-order kinetics, indicating a bimolecular spin conversion process. Numerical simulations suggest the simultaneous presence of a spin diffusion process allowing neighbouring ortho and para molecules to exchange their angular momenta. Cross-polarization experiments found no evidence that the spin conversion of the endohedral H2O molecules is catalysed by (13)C nuclei present in the cages.

Matrix isolation spectroscopy and nuclear spin conversion of NH3 and ND3 in solid parahydrogen

We present matrix isolation infrared absorption spectra of NH3 and ND3 trapped in solid parahydrogen (pH2) at temperatures around 1.8 K. We used the relatively slow nuclear spin conversion (NSC) of NH3 and ND3 in freshly deposited pH2 samples as a tool to assign the sparse vibration-inversion-rotation (VIR) spectra of NH3 in the regions of the ν2, ν4, 2ν4, ν1, and ν3 bands and ND3 in the regions of the ν2, ν4, ν1, and ν3 fundamentals. Partial assignments are also presented for various combination bands of NH3. Detailed analysis of the ν2 bands of NH3 and ND3 indicates that both isotopomers are nearly free rotors; that the vibrational energy is blue-shifted by 1-2%; and that the rotational constants and inversion tunneling splitting are 91-94% and 67-75%, respectively, of the gas-phase values. The line shapes of the VIR absorptions are narrow (0.2-0.4 cm(-1)) for upper states that cannot rotationally relax and broad (>1 cm(-1)) for upper states that can rotationally relax. We report and assign a number of NH3-induced infrared absorption features of the pH2 host near 4150 cm(-1), along with a cooperative transition that involves simultaneous vibrational excitation of a pH2 molecule and rotation-inversion excitation of NH3. The NSCs of NH3 and ND3 were found to follow first-order kinetics with rate constants at 1.8 K of k = 1.88(16) × 10(-3) s(-1) and k = 1.08(8) × 10(-3) s(-1), respectively. These measured rate constants are compared to previous measurements for NH3 in an Ar matrix and with the rate constants measured for other dopant molecules isolated in solid pH2.

Phonon coupling of water monomers in a solid nitrogen matrix

The infrared absorption spectra of the H2O, HDO, and D2O monomers isolated in solid N2 have been recorded at various temperatures between 4 and 30 K. A study of the absorption features of the ν1, ν2, and ν3 vibrational modes for each monomer shows their optical line shapes to be strongly temperature dependent. For all three modes, a decrease in the absorption amplitude and a proportional broadening of the linewidth was observed with increasing temperature, while the integrated absorbance remained constant. These observations were explained in terms of phonon coupling, by which high frequency intramolecular modes decay by exciting matrix phonons. Fits of the linewidth for the lowest frequency ν2 vibrational mode to the predicted vibrational relaxation rate in a solid medium gave average phonon mode frequencies consistent with the Debye frequency for solid N2.

Transient H2O Infrared Satellite Peaks Produced in UV Irradiated Formic Acid Doped Solid Parahydrogen

We report newly identified satellite features of the R(0) rovibrational transition of all the fundamental modes of HDO and the ν3 mode of H2O measured via FTIR spectroscopy immediately after the 193 nm in situ photolysis of formic acid (HCOOH and DCOOD) in solid parahydrogen. The intensities of these satellite features decay slowly with a time constant of τ = 121(7) min after photolysis, even when the sample is maintained below 2 K. We propose that the van der Waals complex H···H2O (H···HDO) is the carrier of the satellite peaks and that these metastable complexes are produced after the low-temperature tunneling reaction of the OH (OD) photoproduct with the parahydrogen host.

Isotope effects of reactions in quantum solids initiated by IR + UV lasers: quantum model simulations for Cl((2)P(3/2)) + X(2)(ν) → XCl + X in X(2) matrices (X = H, D)

Six isotope effects (i)-(vi) are discovered for the reactions Cl + H(2)(ν) → HCl + H in solid para-H(2) ( 1 ) versus Cl + D(2)(ν) → DCl + D in ortho-D(2) ( 2 ), by means of quantum reaction dynamics simulations, within the frame of our simple model ( J. Phys. Chem. A 2009 , 113 , 7630 . ). Experimentally, the reactions may be initiated for ν = 0 and ν ≥ 1, by means of "UV only" photodissociation of the matrix-isolated precursor, Cl(2), or by "IR + UV" coirradiation ( Kettwich , S. C. , Raston , P. L. , and Anderson , D. T. J. Phys. Chem. A 2009 , 113 , 7621 . ), respectively. Specifically, (i) various shape and Feshbach reaction resonances correlate with vibrational thresholds of reactants and products, due to the near-thermoneutrality and low barrier of the system. The energetic density of resonances increases as the square root of mass, from M(X) = M(H) to M(D). (ii) The state selective reaction ( 1 ), ν = 1, is supported by a shape resonance, whereas this type of resonance is absent in ( 2 ), ν = 1. As a consequence, time-resolved measurements should monitor different three-step versus direct error-function type evolutions of the formation of the products. (iii) The effective barrier is lower for reaction 1 , ν = 0, enhancing the tunneling rate, as compared to that for reaction 2 , ν = 0. (iv) For reference, the reaction probabilities P versus total energy E(tot) in the gas exhibit sharp resonance peaks or zigzag behaviors of the reaction probability P versus total energy, near the levels of resonances ( Persky , A. and Baer , M. J. Chem. Phys . 1974 , 60 , 133 . ). These features tend to be washed out and broadened for reaction 1 , and even more so for reaction 2 . For comparison, they disappear for reactions in classical solids. (v) The slopes of P versus E(tot) below the potential barrier increase more steeply for reaction 1 , ν = 0, than for reaction 2 , ν = 0. This enhances the tunneling rate of the heavier isotopomer, reaction 2 , ν = 0, compared to that for reaction 1 . (vi) For a given value of the UV frequency, the translational energy E(trans) increases with mass M(X). Again, this effect supports tunneling of the heavier isotopomer. The isotope effects (i)-(iii), (iv)-(v), and (vi) may be classified as energetic, translational amplitude, and kinematic, respectively. Specifically, the effects (iv)-(v) are due to a systematic decrease of the amplitudes of translational motions of the reactant molecules, from quasi infinite in the gas via still rather large values of para-H(2)(ν) and smaller values for ortho-D(2)(ν) to very small values in classical solids. These isotope effects are special phenomena in quantum solids, which do not occur, neither in the gas phase nor in classical solids. Quantitative predictions, e.g., for the effects of increasing UV frequency on the ratio of reactions probabilities for the UV only versus IR + UV experiments, must account for the interplay of various isotope effects, e.g., (vi) combined with the antagonistic effects (iii) versus (iv) and (v).

Spectral line shape profile of rovibrational transitions of CO embedded in p-H2 crystals studied by high resolution IR diode laser spectroscopy

Line profiles of rovibrational transitions of CO embedded in p-H(2) crystals were studied by high resolution midinfrared diode laser spectroscopy. The line profile analysis for the R(0)(parallel), R(0)(perpendicular), P(1)(parallel), and P(1)(perpendicular) transitions shows that spectral line shapes are well reproduced by a convolution of Gaussian and Lorentzian functions. The temperature dependence of the Lorentzian Gamma(L)(T) and Gaussian widths Gamma(G)(T) shows that there is a nonzero linewidth contribution to each at the T = 0 K limit. The main part of the Lorentzian width Gamma(L)(T = 0) shows anisotropy in the hcp structure and is explained by spontaneous decay of the rotational excited state energy to phonon modes. A smaller part of Gamma(L)(T = 0) is attributed to inhomogeneous broadening due to the point defects of other CO molecules in the crystal. On the other hand, the Gaussian width Gamma(G)(T = 0) is explained by inhomogeneous broadening due to dislocations. In the T > 0 region, Gamma(L)(T) shows strong temperature dependence but Gamma(G)(T) does not. The center frequencies of the R(0)(perpendicular) and P(1)(parallel) transitions show blueshifts and those of the R(0)(parallel) and P(1)(perpendicular) transitions show redshifts with increasing temperature. This phenomenon is explained by a decrease in the anisotropy in the crystal field, which is caused by the averaging of thermal lattice fluctuations. Furthermore, the contribution of vibration and rotation to the linewidth is discussed.

Multiplet splittings and intensities of fine structure components of the Q1(0)H2 + S0(0)N2 transition in a solid parahydrogen matrix

A comprehensive analysis of theoretical multiplet splittings and intensities of the fine structure components of the Q(1)(0)H(2) + S(0)(0)N(2) transition in a solid parahydrogen crystal is presented. The consideration of higher order anisotropic term responsible for splittings is essential to explain the observed splitting of the three components. The pair interaction parameters DeltaB and DeltaC have been determined by comparing the theoretical splittings with the experimental values. The information about the small splittings (approximately 0.1 cm(-1)) due to crystal-field interaction is completely obscured due to fast hopping of v(') = 1, J(') = 0 H(2) vibron. Also, the theoretical expressions are derived for the intensities of the fine structure components of the Q(v(H(2)))(0) + S(v(N(2)))(0) transition and the theoretical results are compared with the experimental findings.

The Cl + H2 --> HCl + H reaction induced by IR + UV irradiation of Cl2 in solid para-H2: quantum model simulation

Recent experimental investigations by the group of D. T. Anderson (Kettwich, S. C.; Raston, P. L.; Anderson, D. T. J. Phys. Chem. A 2009, 113, DOI 10.1021/jp811206a) show that the reaction Cl + H(2) --> HCl + H in the para-H(2) crystal can be induced by infrared (IR) + ultraviolet (UV) coirradiations causing vibrational pre-excitation of the molecular reactant, H(2)(v=1), and generation of the atomic reactant, Cl((2)P(3/2)), by near-resonant photodissociation of a matrix-isolated Cl(2) molecule in the C (1)Pi(u) state, respectively. The corresponding reaction probability P(v=1) for the reactants Cl + H(2)(v=1) is approximately 0.15; this is approximately 25 times larger than P(v=0) for Cl + H(2)(v=0) (as initiated by pure UV irradiation). We present a simple three-step quantum model which accounts for some important parts of the experimental results and allows predictions for other scenarios, for example, UV photodissociation of the Cl(2) molecule by a laser pulse. The first step, vibrational pre-excitation of H(2), yields the molecular initial state which is described using the Einstein model of the para-H(2) crystal. The second step, photodissociation of Cl(2), generates the Cl((2)P(3/2)) atom approaching H(2)(v=1). In the third step, Cl reacts with H(2)(v=1) much more efficiently than with H(2)(v=0) close to threshold. The ultrashort time domains (approximately 100 fs) of steps 2 plus 3 support one- and then two-dimensional models of photodissociation of Cl(2) by short laser pulses and of the subsequent reaction of the system Cl-H-H embedded in frozen environments. The widths of the corresponding wave function describing the translational motion of the reactants is revealed as a significant parameter which is determined not only by the duration of the laser pulse but, even more importantly, by the width of the Gaussian-type distribution of the center of mass of the H(2) molecule in its Einstein cell. As a consequence, the resulting P(v) are quite robust versus variations of the UV pulse durations, allowing extrapolations to continuous wave irradiation. Quantum dynamics simulations of the reaction reveal that the experimental results are due to energetic and dynamical effects.