High-resolution spectroscopy of HCl and DCl isolated in solid parahydrogen: Direct, induced, and cooperative infrared transitions in a molecular quantum solid

Exploring HCl-HCl interactions: QZVPP calculations, improved Lennard-Jones potential, and second virial coefficient analysis for thermodynamics and industrial applications

In this paper, we present a comprehensive analysis of HCl-HCl interactions, including QZVPP calculations, energy fitting, conformation validation, and the determination of the second virial coefficient B using improved Lennard-Jones (ILJ) potential parameters. To acquire accurate interaction energies, initial QZVPP calculations are performed on approximately 1851 randomly generated HCl-HCl conformations. Then, these energies are used to fit an improved Lennard-Jones potential energy surface, allowing for a robust description of HCl-HCl interactions. The ILJ potential parameters are then used to validate particular HCl dimer conformations, ensuring their stability and consistency with experimental observations. The correlation between calculated and experimental conformations strengthens the validity of the ILJ potential parameters. In addition, the second viral coefficient B is calculated at various temperatures using the ILJ potential. The obtained B values are compared to experimental data, demonstrating close agreement, and validating the ILJ potential's ability to accurately capture the intermolecular interactions and gas-phase behavior of the HCl-HCl system. The results of this study demonstrate the effective implementation of QZVPP calculations, energy fitting, and ILJ potential parameters in validating HCl-HCl conformations and accurately determining the second virial coefficient B. The high degree of concordance between calculated B values and experimental data demonstrates the validity of the ILJ potential and its suitability for modeling HCl-HCl interactions. This research contributes to a greater comprehension of HCl-HCl interactions and their implications for numerous chemical and atmospheric processes. The validated conformations, energy fitting method, and calculated second virial coefficients provide valuable instruments for future research and pave the way for more accurate modeling and simulations of HCl-HCl systems.

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.

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.

Quantitative Absorption Spectroscopy of Aluminum Atoms in Cryogenic Parahydrogen Solids

We report results of quantitative ultraviolet (UV) and infrared (IR) absorption spectroscopy experiments on Al-atom-doped cryogenic parahydrogen (pH₂) solids produced by codeposition of Al vapor and pH₂ gas. For Al-atom concentrations [Al] ≲200 parts-per-million (ppm), the Al/pH₂ solids are optically transparent and primarily contain isolated Al atoms, with a small admixture of AlH, Al₂H₂, and Al₂H₄ molecules formed by UV irradiation and Al-atom recombination/reaction. We assign the Al/pH₂ UV absorption spectrum by invoking a large (≈0.6 eV) gas-to-matrix blue shift to accompany the increase in principal quantum number in the 4s ²S ← 3p ²P_1/2 transition, as previously discussed for boron-atom-doped pH₂ solids. We assign a series of sharp features observed in the 4140-4155 cm^-1 IR region to Al-atom-induced Q₁(0) and Q₁(1) absorptions of the pH₂ solid. We use the solid pH₂ Q₁(0) + S₀(0) absorption to determine the sample thickness and to establish a constant pH₂ deposition efficiency independent of the flow rate. Using all of these absorption features in concert, we show that the Al-atom flux delivered by the effusive source is well described by the Knudsen-Langmuir equation, calculate an absolute Al-atom deposition yield per mass of aluminum evaporated, and demonstrate both constant Al-atom deposition and isolation efficiencies for [Al] ≲200 ppm. We discuss this unexpected constant Al-atom isolation efficiency in detail and speculate that it indicates nonuniform Al-atom recombination/reaction on the surface of the accreting sample, perhaps dominated by processes occurring near pH₂ crystallite grain boundaries. We demonstrate "control" over the deposition process, which we define as the ability to set, achieve, and verify "targets" for the final Al-atom concentrations and pH₂ solid thicknesses. This ability is key to sorting out the very different phenomena observed in samples targeting [Al] ≳300 ppm, which are described in the immediately following companion manuscript.

High-resolution infrared spectroscopy of atomic bromine in solid parahydrogen and orthodeuterium

This work extends our earlier investigation of the near-infrared absorption spectroscopy of atomic bromine (Br) trapped in solid parahydrogen (pH2) and orthodeuterium (oD2) [S. C. Kettwich, L. O. Paulson, P. L. Raston, and D. T. Anderson, J. Phys. Chem. A 112, 11153 (2008)]. We report new spectroscopic observations on a series of double transitions involving excitation of the weak Br-atom spin-orbit (SO) transition ((2)P(1/2) ← (2)P(3/2)) in concert with phonon, rotational, vibrational, and rovibrational excitation of the solid molecular hydrogen host. Further, we utilize the rapid vapor deposition technique to produce pH2 crystals with a non-equilibrium mixture of face centered cubic (fcc) and hexagonal closed packed (hcp) crystal domains in the freshly deposited solid. Gentle annealing (T = 4.3 K) of the pH2 sample irreversibly converts the higher energy fcc crystal domains to the slightly more stable hcp structure. We follow the extent of this conversion process using the intensity of the U1(0) transition of solid pH2 and correlate crystal structure changes with changes in the integrated intensity of Br-atom absorption features. Annealing the pH2 solid causes the integrated intensity of the zero-phonon Br SO transition to increase approximately 45% to a value that is 8 times larger than the gas phase value. We show that the magnitude of the increase is strongly correlated to the fraction of hcp crystal domains within the solid. Theoretical calculations presented in Paper II show that these intensity differences are caused by the different symmetries of single substitution sites for these two crystal structures. For fully annealed Br-atom doped pH2 solids, where the crystal structure is nearly pure hcp, the Br-atom SO transition sharpens considerably and shows evidence for resolved hyperfine structure.

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.

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.

Nuclear spin conversion of methane in solid parahydrogen

The nuclear spin conversion of CH(4) and CD(4) isolated in solid parahydrogen was investigated by high resolution Fourier transform infrared spectroscopy. From the analysis of the temporal changes of rovibrational absorption spectra, the nuclear spin conversion rates associated with the rotational relaxation from the J=1 state to the J=0 state for both species were determined at temperatures between 1 and 6 K. The conversion rate of CD(4) was found to be 2-100 times faster than that of CH(4) in this temperature range. The faster conversion in CD(4) is attributed to the quadrupole interaction of D atoms in CD(4), while the conversion in CH(4) takes place mainly through the nuclear spin-nuclear spin interaction. The conversion rates depend on crystal temperature strongly above 3.5 K for CH(4) and above 2 K for CD(4), while the rates were almost constant below these temperatures. The temperature dependence indicates that the one-phonon process is dominant at low temperatures, while two-phonon processes become important at higher temperatures as a cause of the nuclear spin conversion.

Infrared Spectra and Theoretical Calculations of Lithium Hydride Clusters in Solid Hydrogen, Neon, and Argon

A matrix isolation IR study of laser-ablated lithium atom reactions with H2 has been performed in solid para-hydrogen, normal hydrogen, neon, and argon. The LiH molecule and (LiH)(2,3,4) clusters were identified by IR spectra with isotopic substitution (HD, D(2), and H(2) + D(2)) and comparison to frequencies calculated by density functional theory and the MP2 method. The LiH diatomic molecule is highly polarized and associates additional H(2) to form primary (H(2))(2)LiH chemical complexes surrounded by a physical cage of solid hydrogen where the ortho and para spin states form three different primary complexes and play a role in the identification of the bis-dihydrogen complex and in characterization of the matrix cage. The highly ionic rhombic (LiH)(2) dimer, which is trapped in solid matrices, is calculated to be 22 kcal/mol more stable than the inverse hydrogen bonded linear LiH-LiH dimer, which is not observed here. The cyclic lithium hydride trimer and tetramer clusters were also observed. Although the spontaneous reaction of 2 Li and H(2) to form (LiH)(2) occurs on annealing in solid H(2), the formation of higher clusters requires visible irradiation. We observed the simplest possible chemical reduction of dihydrogen using two lithium valence electrons to form the rhombic (LiH)(2) dimer.

Ultrafast dynamics of halogens in rare gas solids

We perform time resolved pump-probe spectroscopy on small halogen molecules ClF, Cl2, Br2, and I2 embedded in rare gas solids (RGS). We find that dissociation, angular depolarization, and the decoherence of the molecule is strongly influenced by the cage structure. The well ordered crystalline environment facilitates the modelling of the experimental angular distribution of the molecular axis after the collision with the rare gas cage. The observation of many subsequent vibrational wave packet oscillations allows the construction of anharmonic potentials and indicate a long vibrational coherence time. We control the vibrational wave packet revivals, thereby gaining information about the vibrational decoherence. The coherence times are remarkable larger when compared to the liquid or high pressure gas phase. This fact is attributed to the highly symmetric molecular environment of the RGS. The decoherence and energy relaxation data agree well with a perturbative model for moderate vibrational excitation and follow a classical model in the strong excitation limit. Furthermore, a wave packet interferometry scheme is applied to deduce electronic coherence times. The positions of those cage atoms, excited by the molecular electronic transitions are modulated by long living coherent phonons of the RGS, which we can probe via the molecular charge transfer states.

High-Resolution Infrared Spectroscopy of HCN−Znn (n = 1−4) Clusters: Structure Determination and Comparisons with Theory

High-resolution infrared laser spectroscopy has been used to obtain rotationally resolved spectra of HCN-Zn(n) (n = 1-4) complexes formed in helium nanodroplets. In the present study the droplets passed through a metal oven, where the zinc vapor pressure was adjusted until one or more atoms were captured by the droplets. A second pickup cell was then used to dope the droplets with a single HCN molecule. Rotationally resolved infrared spectra are obtained for all of these complexes, providing valuable information concerning their structures. Stark spectra are reported and used to determine the corresponding permanent electric dipole moments. Ab initio calculations are also reported for these complexes for comparison with the experimental results.

Infrared-induced reaction of Cl atoms trapped in solid parahydrogen

The 355 nm photodissociation of Cl(2) trapped in a solid parahydrogen matrix at 2 K leads to the formation of isolated Cl photofragments. At these low temperatures (k(B)T approximately 1.4 cm(-1)), the Cl atoms can not react with the parahydrogen matrix since the reaction Cl + H(2)(v = 0, j = 0) --> HCl(v = 0, j = 0) + H is endothermic by 360 cm(-1). Irradiation of the Cl atom doped parahydrogen solid with broadband infrared radiation from 4000 cm(-1) to 5000 cm(-1) induces reaction of atomic Cl with the parahydrogen matrix to form HCl. The infrared-induced chemistry is attributed to solid parahydrogen absorptions that lead to the creation of vibrationally excited H(2)(v = 1), which supply the necessary energy to induce reaction. The kinetics of this low temperature infrared-induced reaction is studied using Fourier Transform infrared spectroscopy of the HCl reaction product. The HCl formation kinetics is first-order and the magnitude of the effective rate constant for the infrared-induced reaction depends on the properties of the near infrared radiation.

Infrared spectra of seeded hydrogen clusters: (para-H2)N–N2O and (ortho-H2)N–N2O, N=2–13

High-resolution infrared spectra of clusters containing para-H2 and/or ortho-H2 and a single nitrous oxide molecule are studied in the 2225-cm(-1) region of the upsilon1 fundamental band of N2O. The clusters are formed in pulsed supersonic jet expansions from a cooled nozzle and probed using a tunable infrared diode laser spectrometer. The simple symmetric rotor-type spectra generally show no resolved K structure, with prominent Q-branch features for ortho-H2 but not para-H2 clusters. The observed vibrational shifts and rotational constants are reported. There is no obvious indication of superfluid effects for para-H2 clusters up to N=13. Sharp transitions due to even larger clusters are observed, but no definite assignments are possible. Mixed (para-H2)N-(ortho-H2)M-N2O cluster line positions can be well predicted by linear interpolation between the corresponding transitions of the pure clusters.

Small para-hydrogen clusters doped with carbon monoxide: Quantum Monte Carlo simulations and observed infrared spectra

The structures and rotational dynamics of clusters of a single carbon monoxide molecule solvated in para-hydrogen, (paraH(2))(N)-CO, have been simulated for sizes up to N=17 using the reptation Monte Carlo technique. The calculations indicate the presence of two series of R(0) rotational transitions with J=1<--0 for cold clusters, similar to those predicted and observed in the case of He(N)-CO. Infrared spectra of these clusters have been observed in the region of the C-O stretch ( approximately 2143 cm(-1)) in a pulsed supersonic jet expansion using a tunable diode laser probe. With the help of the calculations, the observed R(0) rotational transitions have been assigned up to N=9 for the b-type series and N=14 for the a-type series. Theory and experiment agree rather well, except that theory tends to overestimate the b-type energies. The (paraH(2))(12)-CO cluster is calculated to be particularly stable and (relatively) rigid, corresponding to completion of the first solvation shell, and it is observed to have the strongest a-type transition.

Infrared spectra of seeded hydrogen clusters: (paraH2)N–OCS, (orthoH2)N–OCS, and (HD)N–OCS, N=2–7

Infrared spectra of hydrogen-carbonyl sulfide clusters containing paraH2, orthoH2, or HD have been studied in the 2060 cm(-1) region of the C-O stretching vibration. The clusters were formed in pulsed supersonic jet expansions and probed using a tunable infrared diode laser spectrometer. Simple symmetric rotor type spectra were observed and assigned for clusters containing up to N = 7 hydrogen molecules. There was no resolved K structure, and Q-branch features were present for orthoH2 and HD but absent for paraH2. These characteristics can be rationalized in terms of near symmetric rotor structures, very low effective rotational temperatures (0.15 to 0.6 K), and nuclear spin statistics. The observed vibrational shifts were compared with those from recent observations on the same clusters embedded in helium nanodroplets. The observed rotational constants for the paraH2 clusters are in good agreement with a recent quantum Monte Carlo simulation. Some mixed clusters were also observed, such as HD-HD-He-OCS and paraH2 - orthoH2 - OCS.

Tunneling chemical reactions in solid parahydrogen: Direct measurement of the rate constants of R+H2→RH+H (R=CD3,CD2H,CDH2,CH3) at 5 K

Tunneling chemical reactions between deuterated methyl radicals and the hydrogen molecule in a parahydrogen crystal have been studied by Fourier transform infrared spectroscopy. The tunneling rates of the reactions R + H2 --> RH + H (R = CD3,CD2H,CDH2) in the vibrational ground state were determined directly from the temporal change in the intensity of the rovibrational absorption bands of the reactants and products in each reaction in solid parahydrogen observed at 5 K. The tunneling rate of each reaction was found to differ definitely depending upon the degree of deuteration in the methyl radicals. The tunneling rates were determined to be 3.3 x 10(-6) s(-1), 2.0 x 10(-6) s(-1), and 1.0 x 10(-6) s(-1) for the systems of CD3, CD2H, and CDH2, respectively. Conversely, the tunneling reaction between a CH3 radical and the hydrogen molecule did not proceed within a week's time. The upper limit of the tunneling rate of the reaction of the CH3 radical was estimated to be 8 x 10(-8) s(-1).

Infrared matrix-isolation spectroscopy using pulsed deposition of p-H2

We employed pulsed deposition of p-H2 onto a cold target to form a matrix sample suitable for measurements of infrared absorption. Unlike the method of rapid vapor deposition at approximately 2.5 K, developed by Fajardo et al., this method can be performed at a temperature as high as 5.5 K, achievable with a closed-cycle refrigerator; pumping on liquid helium in a cryostat is eliminated. Compared with the enclosed-cell method developed by Oka, Shida, Momose, and co-workers, this method is more versatile in sample preparation, especially for samples at a greater concentration or with high reactivity. Two experiments were tested: the pulse-deposited sample of CH4/p-H2 yields an infrared absorption spectrum nearly identical to that recorded with rapid vapor deposition, and a sample of vinyl chloride (C2H3Cl) in solid p-H2 irradiated with laser emission at 193 nm yields C2H5, in contrast to formation of HCl, C2H2, and a complex of HClC2H2 observed upon photolysis of C2H3Cl in an Ar matrix. These experiments are also compared with those with n-H2 or Ne as the matrix host.