Electron spin resonance g tensors for complexes of Ne and Ar with AlO: theoretical studies related to the large matrix effect observed for AlO

For Ne(n)-AlO (n=2, 4, 6, 8, 10) and Ar(n)-AlO clusters (n=2, 4, 6, 8), the perpendicular (relative to AlO) component of the g tensor was calculated by second-order perturbation theory, using multireference configuration-interaction wave functions. The rare-gas (Rg) atoms were placed axially and/or off axially (one or two rings of four Rg atoms each), and the distance of the Rg atoms from the Al and O atoms, or from the AlO axis, was varied from 4 to 12 bohrs. Rg atoms placed axially mostly increase g(perpendicular), whereas off-axially placed ones lower it below the gas-phase value of AlO. The largest deviations from g(perpendicular) of isolated AlO occur at Ne-Al,O distances of 5-6 bohrs, and Ar-Al,O distances of 6-9 bohrs, with maximal lowerings of about 1600 ppm for Ne and about 2200 ppm (estimated) for Ar in the case of two axial and eight off-axial Rg atoms. Electron spin resonance studies by Knight and Weltner found large matrix effects for AlO, with downshifts of g(perpendicular) observed to be about 450 and 1150 ppm in Ne and Ar matrices, respectively.     

Basis set limits of the second order Moller-Plesset correlation energies of water, methane, acetylene, ethylene, and benzene

We report second order Moller-Plesset (MP2) and MP2-F12 total energies on He, Ne, Ar, H(2)O, CH(4), C(2)H(2), C(2)H(4), and C(6)H(6), using the correlation consistent basis sets, aug-cc-pVXZ (X=D-7). Basis set extrapolation techniques are applied to the MP2 and MP2-F12/B methods. The performance of the methods is tested in the calculations of the atoms, He, Ne, and Ar. It is indicated that the two-point extrapolation of MP2-F12/B with the basis sets (X=5,6) is the most reliable. Similar accuracy is obtained using two-point extrapolated conventional MP2 with the basis sets (X=6,7). For the molecules investigated the valence MP2 correlation energy is estimated within 1 mE(h).     

Ab initio potential energy surface and predicted microwave spectra for Ar--OCS dimer and structures of Arn--OCS (n = 2-14) clusters

An ab initio potential energy surface for the Ar--OCS dimer was calculated using the coupled-cluster singles and doubles with noniterative inclusion of connected triples [CCSD(T)] with a large basis set containing bond functions. The interaction energies were obtained by the supermolecular approach with the full counterpoise correction for the basis set superposition error. The CCSD(T) potential was found to have two minima corresponding to the T-shaped and the collinear Ar--SCO structures. The two-dimensional discrete variable representation method was employed to calculate the rovibrational energy levels for five isotopomers Ar--OCS, Ar--OC34S, Ar--O13CS, Ar--18OCS, and Ar--17OCS. The calculated pure rotational transition frequencies for the vibrational ground state of the five isotopomers are in good agreement with the observed values. The corresponding microwave spectra show that the b-type transitions (Delta Ka = +/-1) are significantly stronger than the a-type transitions (Delta Ka = 0). Minimum-energy structures of the Ar2--OCS trimer were been determined with MP2 optimization, whereas the minimum-energy structures of the Arn--OCS clusters with n = 3-14 were obtained with the pairwise additive potentials. It was found that there are two minima corresponding to one distorted tetrahedral structure and one planar structure for the ternary complex. The 14 nearest neighbor Ar atoms form the first solvation shell around the OCS molecule.     

Cooperativity between S···π and Rg···π in the OCS···C6H6···Rg (Rg = He, Ne, Ar, and Kr) van der Waals complexes

The complexes OCS···C(6)H(6), C(6)H(6)···Rg, and OCS···C(6)H(6)···Rg (Rg = He, Ne, Ar, and Kr) have been studied by means of MP2 calculations and QTAIM analyses. The optimized geometries of the title complexes have C(6v) symmetry. The intermolecular interactions in the OCS···C(6)H(6)···Rg complexes are comparatively stronger than that in the OCS···C(6)H(6) complex, which prove that the He, Ne, Ar, and Kr atoms have the ability to form weak bonds with the benzene molecule. In QTAIM studies, the π-electron density of benzene was separated from the total electron density. The molecular graphs and topological parameters of the OCS···πC(6)H(6), πC(6)H(6)···Rg, and OCS···πC(6)H(6)···Rg complexes indicate that the interactions are mainly attributed to the electron density provided by the π-bonding electrons of benzene and the top regions of the S and Rg atoms. Charge transfer is observed from the benzene molecule to SCO/Rg in the formation of the OCS···C(6)H(6), C(6)H(6)···Rg, and OCS···C(6)H(6)···Rg complexes. Molecular electrostatic potential (MEP) analyses suggest that the electrostatic energy plays a pivotal role in these intermolecular interactions.     

Theoretical Study of Decomposition Pathways for Rare-gas-containing Compounds HRgX （Rg = He, Ne, Ar, Kr; X = Cl, Br）

1 INTRODUCTION transfer species with common formula H-Rgδ -Xδ, where X represents a strongly electrone-gative atom Since xenon hexafluoraplatinate, XePtF6 , the [1] or fragment, and Rg is a rare-gas atom. These species first rare gas-containing compound was discovered have linear equi- librium geometries and are mainly by Bartlett in 1962, rare gases are getting more and bound up by strong columbic attraction between (H- more attention and have been found to be possible to Rg) and…     

Rare gas effects on hyperfine coupling constants of BO, AlO, and GaO

Using density functional theory methods and large basis sets, we calculated hyperfine coupling constants (HFCCs) for the (11)B, (17)O, (27)Al, and (69)Ga nuclei of the radicals BO, AlO, and GaO (XO), embedded in 2-14 rare gas (Rg) Ne and Ar atoms. Kr atoms were included for AlO. The distance of the Rg atoms from XO was varied from 4 to 12 bohr. Matrix effects cause A(iso)(X) to increase, accompanied by decreases in A(dip)(X) and A(dip)(O), while A(iso)(O) remains close to zero. Changes are largest for AlO, slightly smaller for GaO, and very small for BO, in line with the molecular polarizabilities. Observed changes of A(iso)(X) and A(dip)(X) for BO in Ne matrixes and for AlO in Ne, Ar, and Kr matrixes are reproduced in complexes with 12 Rg atoms at distances of 5-6 bohr or 14 Rg atoms at distances of 6-7 bohr. For GaO, experimental data are available only in Ne matrixes. Theoretical results obtained for HFCCs of (17)O could not be verified due to insufficient experimental information. Estimates of HFCCs in matrixes not yet experimentally studied and for GaO in the gas phase have been made. Due to the interaction with rare gas atoms, p-spin density on the X and O atoms of XO is converted into s-spin density on X, thereby causing an increase (in magnitude) of A(iso)(X), accompanied by decreases in A(dip) of X and O. The higher polarizability of XO along the bond axis is reflected in complexes that have axial Rg atoms showing larger changes in HFCCs than comparable complexes without axial Rg atoms.     

Photoionization of helium nanodroplets doped with rare gas atoms

Photoionization of He droplets doped with rare gas atoms (Rg=Ne, Ar, Kr, and Xe) was studied by time-of-flight mass spectrometry, utilizing synchrotron radiation from the Advanced Light Source from 10 to 30 eV. High resolution mass spectra were obtained at selected photon energies, and photoion yield curves were measured for several ion masses (or ranges of ion masses) over a wide range of photon energies. Only indirect ionization of the dopant rare gas atoms was observed, either by excitation or charge transfer from the surrounding He atoms. Significant dopant ionization from excitation transfer was seen at 21.6 eV, the maximum of He 2p 1P absorption band for He droplets, and from charge transfer above 23 eV, the threshold for ionization of pure He droplets. No Ne+ or Ar+ signal from droplet photoionization was observed, but peaks from HenNe+ and HenAr+ were seen that clearly originated from droplets. For droplets doped with Rg=Kr or Xe, both Rg+ and HenRg+ ions were observed. For all rare gases, Rg2+ and HenRgm+ (n,m> or =1) were produced by droplet photoionization. Mechanisms of dopant ionization and subsequent dynamics are discussed.     

Atmospheric-pressure Penning ionization of aliphatic hydrocarbons

A study has been made of the atmospheric-pressure Penning ionization (APPeI) of aliphatic hydrocarbons (pentane, hexane, heptane, and octane) with long-lived rare gas atoms (Rg*). The metastable rare gas atoms (He*, Ne*, Ar* and Kr*) were generated by the negative-mode corona discharge of atmospheric-pressure rare gases. In the Rg*APPeI mass spectra for aliphatic hyrocarbons, the relative abundances of fragment ions were found to increase in the order of He* --> Ne* --> Ar* --> Kr*. The order is in the opposite direction to the internal energies of the Rg*. The less fragmentation observed for He* may be because the nascent molecular ions [M(+.)]* formed by Penning ionization have lifetimes long enough for them to be collisionally deactivated in the atmospheric-pressure ion source. It was found that the relative abundances of fragment ions in Ar*APPeI increased when the sample pressure in the ion source was reduced. This is attributed to the collision of Ar* with molecular ions followed by fragmentation.     

Theoretical Study of Rg·NO (Rg=He,Ne, Ar and Kr) Complexes

Rg·NO (Rg=He, Ne, Ar and Kr) complexes were studied using ab initio calculations. The neutral Rg·NO complex geometry and vibrational frequencies were calculated with the cc-pVDZ basis set at the CCSD(T) level of theory. The calculations show that the geometry of the Rg·NO complexes is a skewed T-shape with the Rg atom on the oxygen side of the NO molecule, and that the RgNO bond angle increases with mass. The dissociation energies (DE) and ionization energies (IE) of the neutral Rg·NO complexes, and the dissociation energies of Rg·NO+ ionic complexes were calculated using Gaussian-2 (G2) methods and a high accuracy energy model. The ionization energies of the neutral Rg·NO complexes range from 9.265 eV for He·NO to 9.132 eV for Kr·NO and the dissociation energies of Rg·NO+ range from 0.017 eV for He·NO+ to 0.156 eV for Kr·NO+, in line with the expectation based on the increasing polarizability of the Rg atom.     

Rotational excitations of N2O in small helium clusters and the role of Bose permutation symmetry

We present a detailed study of the energetics, structures, and Bose properties of small clusters of (4)He containing a single nitrous oxide (N(2)O) molecule, from N=1 (4)He up to sizes corresponding to completion of the first solvation shell around N(2)O (N=16 (4)He). Ground state properties are calculated using the importance-sampled rigid-body diffusion Monte Carlo method, rotational excited state calculations are made with the projection operator imaginary time spectral evolution method, and Bose permutation exchange and associated superfluid properties are calculated with the finite temperature path integral method. For N< or =5 the helium atoms are seen to form an equatorial ring around the molecular axis, at N=6 helium density starts to occupy the second (local) minimum of the N(2)O-He interaction at the oxygen side of the molecule, and N=9 is the critical size at which there is onset of helium solvation all along the molecular axis. For N> or =8 six (4)He atoms are distributed in a symmetric, quasirigid ring around N(2)O. Path integral calculations show essentially complete superfluid response to rotation about the molecular axis for N> or =5, and a rise of the perpendicular superfluid response from zero to appreciable values for N> or =8. Rotational excited states are computed for three values of the total angular momentum, J=1-3, and the energy levels fitted to obtain effective spectroscopic constants that show excellent agreement with the experimentally observed N dependence of the effective rotational constant B(eff). The non-monotonic behavior of the rotational constant is seen to be due to the onset of long (4)He permutation exchanges and associated perpendicular superfluid response of the clusters for N> or =8. We provide a detailed analysis of the role of the helium solvation structure and superfluid properties in determining the effective rotational constants.     

Replacement equivalence of H- and argon in small (Ar)nH- clusters from optimized structure calculations

The structural properties of some of the smaller ionic clusters of argon atoms containing the atomic impurity H-, ArnH- with n from 2 up to 7, are examined using different modeling for the interactions within each cluster and by employing different theoretical treatments, both classical and quantum, for the energetics. The same calculations are also carried out for the corresponding neutral homogeneous clusters Ar(n+1). The results of the calculations, the physical reliability of the interactions modeling, and the similarities and the difference between the anionic and the neutral complexes are discussed in some detail. The emerging picture shows that, due to specific features of the employed atom-atom potentials, the ArnH- and Ar(n+1) clusters present very similar structures, where the H- dopant substitutes for one of the outer Ar atoms but does not undergo as yet solvation within such small clusters.     

Noble‐Gas Complexes: Theoretical Investigation of Multicenter Polynuclear Species

Ab initio calculations at the MP2 level of theory disclose the conceivable existence of neutral complexes containing four or five distinct noble gases (Ng) each bound to a distinct Be‐atom. These multicenter polynuclear Ng molecules are formally obtained by replacing the H‐atoms of CH₄ and but‐2‐yne with ? NBeNg moieties, which behave as independent monovalent ‘functional groups’. Our investigated complexes include the five homotetranuclear [C(NBeNg)₄] complexes 1 – 5 (Ng=He? Xe), the five heterotetranuclear complexes [CN₄Be₄(He)(Ne)(Ar)(Kr)] ( 6 ), [CN₄Be₄(He)(Ne)(Ar)(Xe)] ( 7 ), [CN₄Be₄(He)(Ne)(Kr)(Xe)] ( 8 ), [CN₄Be₄(He)(Ar)(Kr)(Xe)] ( 9 ), and [CN₄Be₄(Ne)(Ar)(Kr)(Xe)] ( 10 ), and the heteropentanuclear complex [HC₄N₅Be₅(He)(Ne)(Ar)(Kr)(Xe)] ( 11 ). We also investigated the five model complexes [H₃CNBeNg] (Ng=He? Xe) containing a single ? NBeNg moiety. The geometries and vibrational frequencies of all these species, invariably characterized as minimum‐energy structures, were computed at the MP2(full)/6‐31G(d,p)/SDD level of theory, and their stability with respect to the loss of the various Ng‐atoms was evaluated by single‐point calculations at the MP2(full)/6‐311G(d)/SDD level of theory. The beryllium‐Ng binding energies range from ca. 17 (Ng=He) to ca. 63 (Ng=Xe) kJ/mol, and the results of natural‐bond‐orbital (NBO) and atoms‐in‐molecules (AIM) analysis reveal that the Be? Ng interaction is essentially electrostatic for helium, neon, argon, and krypton, and has probably a small covalent contribution for xenon.     

Energy transfer kinetics of the np5(n + 1)p excited states of Ne and Kr

Energy transfer rate constants for Ne(2p(5)3p) and Kr(4p(5)5p) atoms colliding with ground state rare gas atoms (Rg) have been measured. In part, this study is motivated by the possibility of using excited rare gas atoms as the active species in optically pumped laser systems. Rg(np(5)(n + 1)s) metastable states may be produced using low-power electrical discharges. The potential then exits for optical pumping and laser action on the np(5)(n + 1)p ? np(5)(n + 1)s transitions. Knowledge of the rate constants for collisional energy transfer and deactivation of the np(5)(n + 1)p states is required to evaluate the laser potential for various Rg + buffer gas combinations. In the present study we have characterized energy transfer processes for Ne (2p(5)3p) + He for the six lowest energy states of the multiplet. Rate constants for state-to-state transfer have been determined. Deactivation of the lowest energy level of Kr (4p(5)5p) by He, Ne, and Kr has also been characterized. Initial results suggest that Kr (4p(5)5p) + Ne mixtures may be the best suited for optically pumped laser applications.     

Structures and vibrational spectroscopy of partially reduced gas-phase cerium oxide clusters

This work demonstrates that the most stable structures of even small gas-phase aggregates of cerium oxide with 2-5 cerium atoms show structural motifs reminiscent of the bulk ceria. This is different from main group and transition metal oxide clusters, which often display structural features that are distinctly different from the bulk structure. The structures of Ce(2)O(2)(+), Ce(3)O(4)(+), and (CeO(2))(m)CeO(+) clusters (m = 0-4) are unambiguously determined by a combination of global structure optimizations at the density functional theory level and infrared vibrational predissociation spectroscopy of the cluster-rare gas atom complexes. The structures of Ce(2)O(2)(+) and Ce(2)O(3)(+) exhibit a Ce-O-Ce-O four-membered ring with characteristic absorptions between 430 and 680 cm(-1). Larger clusters have common structural features containing fused Ce-O-Ce-O four-membered rings which lead to intense absorption bands at around 500 and 650 cm(-1). Clusters containing a terminal Ce=O bond show a characteristic absorption band between 800 and 840 cm(-1). For some cluster sizes multiple isomers are observed. Their individual infrared signatures are identified by tuning their relative population through the choice of He, Ne or Ar messenger atoms. The present results allow us to benchmark different density functionals which yield different degrees of localization of unpaired electrons in Ce 4f states.     

Internal rotation in NH4(+)-Rg dimers (Rg = He, Ne, Ar): potential energy surfaces and IR spectra of the nu 3 band

The intermolecular potential energy surfaces for the electronic ground states of the ammonium ion-rare gas dimers NH4(+)-He and NH4(+)-Ne are calculated at the MP2 and CCSD(T)/aug-cc-pVXZ (X = D/T/Q) levels of theory. The global minima of both potentials correspond to proton (vertex)-bound structures, Re = 3.13 A, De = 171 cm-1 (He) and Re = 3.21 A, De = 302 cm-1 (Ne). The face- and edge-bound structures are local minima and transition states for the internal rotation dynamics, corresponding to barriers of approximately 20 (He) and 50 cm-1 (Ne). The ab initio potentials are employed in numerical solutions to the rotation-intermolecular vibration Hamiltonian to determine the term values and the rotational and distortion constants for the lowest bound levels in the intramolecular ground vibrational state of both complexes. The results are used to assess the accuracy of two-dimensional (fixed-R) representations of the potentials for determining the internal rotor levels in the ground and nu 3 vibrational states. This model is employed to produce simulations of the IR nu 3 transitions, which are compared to the experimental spectra recorded using photofragmentation spectroscopy. In the case of NH4(+)-Ne the potential parameters are least-squares fitted to the experimental spectrum. The trends within the NH4(+)-Rg series (Rg = He, Ne, Ar) revealed by both the IR spectra and theoretical calculations are discussed.     

Ab initio studies on the mechanism of the size-dependent hydrogen-loss reaction in Mg+(H2O)n

The mechanism of size-dependent intracluster hydrogen loss in the cluster ions Mg(+)(H(2)O)(n), which is switched on around n=6, and off around n=14, was studied by ab initio calculations at the MP2/6-31G* and MP2/6-31G** levels for n=1-6. The reaction proceeds by Mg(+)-assisted breaking of an H-O bond in one of the H(2)O molecules. The reaction barrier is dependent on both the cluster size and the solvation structure. As n increases from 1 to 6, there is a dramatic drop in the reaction barrier, from greater than 70 kcal mol(-1) for n=1 to less than 10 kcal mol(-1) for n=6. In the transition structures, the Mg atom is close to the oxidation state of +2, and H(2)O molecules in the first solvation shell are much more effective in stabilizing the transition structures and lowering the reaction barriers than H(2)O molecules in the other solvation shells. While the reaction barrier for trimer core structures with only three H(2)O molecules in the first shell is greater than 24 kcal mol(-1), even for Mg(+)(H(2)O)(6), it drops considerably for clusters with four-six H(2)O molecules in the first shell. The more highly coordinated complexes have comparable or slightly higher energy than the trimer core structures, and the presence of such high coordination number complexes is the underlying kinetic factor for the switching on of the hydrogen-loss reaction around n=6. For clusters with trimer core structures, the hydrogen loss reaction is much easier when it is preceded by an isomerization step that increases the coordination number around Mg(+). Delocalization of the electron on the singly occupied molecular orbital (SOMO) away from the Mg(+) ion is observed for the hexamer core structure, while at the same time this isomer is the most reactive for the hydrogen-loss reaction, with an energy barrier of only 2.7 kcal mol(-1) at the MP2/6-31G** level.     

Structure and photoabsorption properties of cationic alkali dimers solvated in neon clusters

We present a theoretical investigation of the structure and optical absorption of M(2)(+) alkali dimers (M=Li,Na,K) solvated in Ne(n) clusters for n=1 to a few tens Ne atoms. For all these alkali, the lowest-energy isomers are obtained by aggregation of the first Ne atoms at the extremity of the alkali molecule. This particular geometry, common to other M(2)(+)-rare gas clusters, is intimately related to the shape of the electronic density of the X (2)Σ(g)(+) ground state of the bare M(2)(+) molecules. The structure of the first solvation shell presents equilateral Ne(3) and capped pentagonal Ne(6) motifs, which are characteristic of pure rare gas clusters. The size and geometry of the complete solvation shell depend on the alkali and were obtained at n=22 with a D(4h) symmetry for Li and at n=27 with a D(5h) symmetry for Na. For K, our study suggests that the closure of the first solvation shell occurs well beyond n=36. We show that the atomic arrangement of these clusters has a profound influence on their optical absorption spectrum. In particular, the XΣ transition from the X (2)Σ(g)(+) ground state to the first excited (2)Σ(u)(+) state is strongly blueshifted in the Frank-Condon area.     

A theoretical study of 1:1 and 1:2 complexes of acetylene with nitrosyl hydride

Ab initio calculations at MP2 computational level using aug-cc-pVTZ basis set were used to analyze the interactions between 1:1 and 1:2 complexes of acetylene and nitrosyl hydride. The structures obtained have been analyzed with the atoms in molecules and the density functional theory–symmetry adapted perturbation theory methodologies. Four minima were located on the potential energy surface of the 1:1 complex. Twenty-four different structures have been obtained for the 1:2 complexes. Five types of interactions are observed, CH···O, CH···N, NH···π hydrogen bonds and orthogonal interactions between the π clouds of triple bond, or the lone pair of oxygen with the electron-deficient region of the nitrogen atom. Stabilization energies of the 1:1 and 1:2 clusters including basis set superposition error and ZPE are in the range 3–8 and 6–17 kJ mol⁻¹ at MP2/aug-cc-pVTZ computational level, respectively. Blue shift of NH bond upon complex formation in the ranges between 18–30 and 20–96 cm⁻¹ is predicted for 1:1 and 1:2 clusters, respectively. The total nonadditive energy in the 1:2 cluster, calculated as the sum of the supermolecular nonadditive MP2 energy and the three-body dispersion energy, presents values between −1.48 and 1.20 kJ mol⁻¹.     

Modeling and spectral simulation of matrix-isolated molecules by density functional calculations: a case study on formic acid dimer

The supermolecule approach has been used to model molecules embedded in solid argon matrix, wherein interaction between the guest and the host atoms in the first solvation shell is evaluated with the use of density functional calculations. Structural stability and simulated spectra have been obtained for formic acid dimer (FAD)-Ar(n) (n = 21-26) clusters. The calculations at the B971∕6-31++G(3df,3pd) level have shown that the tetrasubstitutional site on Ar(111) plane is likely to incorporate FAD most stably, in view of consistency with the matrix shifts available experimentally.