A reliable new three-dimensional potential energy surface for H2–Kr

Full quantum calculation of the rovibrational states and intensities for a symmetric top-linear molecule dimer: Hamiltonian, basis set, and matrix elements

The rovibrational energy levels and intensities of the CH₃F-H₂ dimer have been obtained using our recent global intermolecular potential energy surface [X.-L. Zhang et al., J. Chem. Phys. 148, 124302 (2018)]. The Hamiltonian, basis set, and matrix elements are derived and given for a symmetric top-linear molecule complex. This approach to the generation of energy levels and wavefunctions can readily be utilized for studying the rovibrational spectra of other van der Waals complexes composed of a symmetric top molecule and a linear molecule, and may readily be extended to other complexes of nonlinear molecules and linear molecules. To confirm our method, the rovibrational levels of the H₂O-H₂ dimer have been computed and shown to be in good agreement with experiment and with previous theoretical results. The rovibrational Schrödinger equation has been solved using a Lanczos algorithm together with an uncoupled product basis set. As expected, dimers containing ortho-H₂ are more strongly bound than dimers containing para-H₂. Energies and wavefunctions of the discrete rovibrational levels of CH₃F-paraH₂ complexes obtained from the direct vibrationally averaged 5-dimensional potentials are in good agreement with the results of the reduced 3-dimensional adiabatic-hindered-rotor (AHR) approximation. Accurate calculations of the transition line strengths for the orthoCH₃F-paraH₂ complex are also carried out, and are consistent with results obtained using the AHR approximation. The microwave spectrum associated with the orthoCH₃F-orthoH₂ dimer has been predicted for the first time.

Noncovalent Complexes of the Noble-Gas Atoms: Analyzing the Transition from Physical to Chemical Interactions

The bonding character of the noncovalent complexes of the noble-gas (Ng) atoms ranges from nearly purely dispersive contacts to interactions featuring appreciable contributions of induction and charge transfer. In this study, we discuss a new quantitative index that seems peculiarly informative about these diverse bonding situations. This index was termed as the degree of polarization (DoP) of Ng, as it measures, in essence, the Ng polarization promoted by the binding partner. The definition of the DoP(Ng) relies on the analysis of the local electron energy density H(r), and its physical meaning was best appreciated by studying also the charge-displacement function and the molecular electrostatic potential of the investigated benchmark species, that include nearly 60 Ngs complexes of different bonding character. The DoP(Ng) appears of general applicability, and is also positively correlated with other bonding character indices. © 2019 Wiley Periodicals, Inc.

Insight into the halogen-bond nature of noble gas-chlorine systems by molecular beam scattering experiments, ab initio calculations and charge displacement analysis

We have carried out molecular-beam scattering experiments and high-level ab initio investigations on the potential energy surfaces of a series of noble-gas-Cl2 adducts. This effort has permitted the construction of a simple, reliable and easily generalizable analytical model potential formulation, which is based on a few physically meaningful parameters of the interacting partners and transparently shows the origin, strength, and stereospecificity of the various interaction components. The results demonstrate quantitatively beyond doubt that the interaction between a noble-gas (Ng) atom - even He - and Cl2 in a collinear configuration is characterized by weak halogen bond (XB) formation, accompanied by charge transfer (CT) from the Ng to chlorine. This characteristic, which stabilizes the adduct, rapidly disappears on going towards the T-shaped configuration, dominated by pure van der Waals (vdW) forces. Similarly, a pure vdW interaction takes place - with no CT component in any configuration - if Cl2 is present in the lowest πg* → σu* excited state, because the change in electron density that accompanies the excitation eliminates the Cl2 polar flattening and σ hole, making the XB interaction inaccessible.

Selective Emergence of the Halogen Bond in Ground and Excited States of Noble-Gas-Chlorine Systems

Molecular-beam scattering experiments and theoretical calculations prove the nature, strength, and selectivity of the halogen bonds (XB) in the interaction of halogen molecules with the series of noble gas (Ng) atoms. The XB, accompanied by charge transfer from the Ng to the halogen, is shown to take place in, and measurably stabilize, the collinear conformation of the adducts, which thus becomes (in contrast to what happens for other Ng-molecule systems) approximately as bound as the T-shaped form. It is also shown how and why XB is inhibited when the halogen molecule is in the ³ Π_u excited state. A general potential formulation fitting the experimental observables, based on few physically essential parameters, is proposed to describe the interaction accurately and is validated by ab initio computations.

Predicted infrared spectra in the HF stretching band of the H2-HF complex

The infrared spectra with hydrogen fluoride (HF) and deuterium fluoride (DF) (v₂ = 1 ← 0) for eight isotropic species of H₂-HF complex are predicted, based on our newly constructed high-accuracy ab initio potential energy surface [D. Yang et al., J. Chem. Phys. 148, 184301 (2018)]. The radial discrete variable representation/angular finite basis representation method and Lanczos algorithm were used to determine the ro-vibrational energy levels and wave functions for eight species of H₂-HF complex (para-H₂-HF, ortho-H₂-HF, para-D₂-HF, ortho-D₂-HF, para-H₂-DF, ortho-H₂-DF, para-D₂-DF, and ortho-D₂-DF) with separating the inter- and intra-molecular vibrations. Bound states properties including their dissociation energies and rotational constants were presented. The calculated band origins are all red shifted to the isolated HF molecule and in good agreement with available experimental values. The frequencies and line intensities of ro-vibrational transitions in the HF stretching band were further calculated, and the predicted infrared spectra are consistent with available observed spectra. Among them, the spectra for three isotopic species of H₂-HF (para-H₂-DF, para-D₂-DF, and ortho-D₂-DF) were predicted for the first time.

Interaction-induced electric (hyper)polarizability in the dihydrogen-neon pair: basis set and electron correlation effects

We investigated the interaction (hyper)polarizability of neon-dihydrogen pairs by performing high-level ab initio calculations with atom/molecule-specific, purpose-oriented Gaussian basis sets. We obtained interaction-induced electric properties at the SCF, MP2, and CCSD levels of theory. At the CCSD level, for the T-shaped configuration, around the respective potential minimum of 6.437 a₀, the interaction-induced mean first hyperpolarizability varies for 5 <  R/a₀ < 10 as[Formula: see text]Again, at the CCSD level, but for the L-shaped configuration around the respective potential minimum of 6.572 a₀, this property varies for 5 <  R/a₀ < 10 as[Formula: see text] Graphical Abstract Interaction-induced mean dipole polarizability ([Formula: see text]) for the T-shaped configuration of H₂-Ne calculated at the SCF, MP2, and CCSD levels of theory.

Ab initio potential energy surface and microwave spectra for the H2-HCCCN complex

We present a four-dimensional ab initio potential energy surface of the H₂-HCCCN complex at the coupled-cluster singles and doubles with noniterative inclusion of connected triples [CCSD(T)]-F12 level with a large basis set including an additional set of bond functions. The artificial neural networks method was extended to fit the intermolecular potential energy surface. The complex has a planar linear global minimum with the well depth of 199.366 cm^-1 located at R = 5.09 Å, φ = 0°, θ₁ = 0°, and θ₂ = 180°. An additional planar local minimum is also found with a depth of 175.579 cm^-1 that is located at R = 3.37 Å, φ = 0°, θ₁ = 110°, and θ₂ = 104°. The radial discrete variable representation/angular finite basis representation and the Lanczos algorithm were employed to calculate the rovibrational energy levels for four species of H₂-HCCCN (pH₂-HCCCN, oH₂-HCCCN, pD₂-HCCCN, and oD₂-HCCCN). The rotational frequencies and spectroscopic parameters were also determined for four complexes, which agree well with the experimental values.

An exchange-Coulomb model potential energy surface for the Ne-CO interaction. II. Molecular beam scattering and bulk gas phenomena in Ne-CO mixtures

Four potential energy surfaces are of current interest for the Ne-CO interaction. Two are high-level fully ab initio surfaces obtained a decade ago using symmetry-adapted perturbation theory and supermolecule coupled-cluster methods. The other two are very recent exchange-Coulomb (XC) model potential energy surfaces constructed by using ab initio Heitler-London interaction energies and literature long range dispersion and induction energies, followed by the determination of a small number of adjustable parameters to reproduce a selected subset of pure rotational transition frequencies for the (20)Ne-(12)C(16)O van der Waals cluster. Testing of the four potential energy surfaces against a wide range of available experimental microwave, millimeter-wave, and mid-infrared Ne-CO transition frequencies indicated that the XC potential energy surfaces gave results that were generally far superior to the earlier fully ab initio surfaces. In this paper, two XC model surfaces and the two fully ab initio surfaces are tested for their abilities to reproduce experiment for a wide range of nonspectroscopic Ne-CO gas mixture properties. The properties considered here are relative integral cross sections and the angle dependence of rotational state-to-state differential cross sections, rotational relaxation rate constants for CO(v=2) in Ne-CO mixtures at T=296 K, pressure broadening of two pure rotational lines and of the rovibrational lines in the CO fundamental and first overtone transitions at 300 K, and the temperature and, where appropriate, mole fraction dependencies of the interaction second virial coefficient, the binary diffusion coefficient, the interaction viscosity, the mixture shear viscosity and thermal conductivity coefficients, and the thermal diffusion factor. The XC model potential energy surfaces give results that lie within or very nearly within the experimental uncertainties for all properties considered, while the coupled-cluster ab initio surface gives results that agree similarly well for all but one of the properties considered. When the present comparisons are combined with the ability to give accurate spectroscopic transition frequencies for the Ne-CO van der Waals complex, only the XC potential energy surfaces give results that agree well with all extant experimental data for the Ne-CO interaction.

An exchange-Coulomb model potential energy surface for the Ne-CO interaction. I. Calculation of Ne-CO van der Waals spectra

Exchange-Coulomb model potential energy surfaces have been developed for the Ne-CO interaction. The initial model is a three-dimensional potential energy surface based upon computed Heitler-London interaction energies and literature results for the long-range induction and dispersion energies, all as functions of interspecies distance, the orientation of CO relative to the interspecies axis, and the bond length of the CO molecule. Both a rigid-rotor model potential energy surface, obtained by setting the CO bond length equal to its experimental spectroscopic equilibrium value, and a vibrationally averaged model potential energy surface, obtained by averaging the stretching dependence over the ground vibrational motion of the CO molecule, have been constructed from the full data set. Adjustable parameters in each model potential energy surface have been determined through fitting a selected subset of pure rotational transition frequencies calculated for the (20)Ne-(12)C(12)O isotopolog to precisely known experimental values. Both potential energy surfaces provide calculated results for a wide range of available experimental microwave, millimeter-wave, and midinfrared Ne-CO transition frequencies that are generally far superior to those obtained using the best current literature potential energy surfaces. The vibrationally averaged CO ground state potential energy surface, employed together with a potential energy surface obtained from it by replacing the ground vibrational state average of the CO stretching dependence of the potential energy surface by an average over the first excited CO vibrational state, has been found to be particularly useful for computing and/or interpreting mid-IR transition frequencies in the Ne-CO dimer.

Dipole oscillator strength distributions, properties, and dispersion energies for ethylene, propene, and 1-butene

A recommended isotropic dipole oscillator strength distribution (DOSD) has been constructed for the ethylene molecule through the use of quantum mechanical constraint techniques and experimental dipole oscillator strength (DOS) data; the DOS data employed are recent experimental results not available at the time of the original constrained DOSD analysis of this molecule. The constraints are furnished by molar refractivity data and the Thomas–Reiche–Kuhn sum rule. The DOSD is used to evaluate a variety of isotropic dipole oscillator strength sums, logarithmic dipole oscillator strength sums, and mean excitation energies for ethylene. Pseudo-DOSDs for this molecule, and for propene and 1–butene, which are based on an earlier constrained DOSD analysis for these molecules, are developed. They are used to obtain reliable results for the isotropic dipole–dipole dispersion-energy coefficients C6, for the interactions of the alkenes with each other and with 47 other species, and the triple-dipole dispersion-energy coefficients C9 for interactions involving any triple of molecules taken from ethylene, propene, and 1–butene.Key words: alkenes, dipole properties, pseudo-states, dipole–dipole and triple-dipole dispersion energies, long-range additive, non-additive interaction energies.

Accuracy of recent potential energy surfaces for the He–N2 interaction. I. Virial and bulk transport coefficients

A new exchange-Coulomb semiempirical model potential energy surface for the He-N2 interaction has been developed. Together with two recent high-level ab initio potential energy surfaces, it has been tested for the reliability of its predictions of second-virial coefficients and bulk transport phenomena in binary mixtures of He and N2. The agreement with the relevant available measurements is generally within experimental uncertainty for the exchange-Coulomb surface and the ab initio surface of Patel et al. [J. Chem. Phys. 119, 909 (2003)], but with slightly poorer agreement for the earlier ab initio surface of Hu and Thakkar [J. Chem. Phys. 104, 2541 (1996)].

Time-dependent coupled-cluster calculations of polarizabilities and dispersion energy coefficients

Time-dependent coupled cluster theory, with unrestricted electron spins and full treatment of orbital rotation, is used to calculate polarizabilities at imaginary frequencies for Li, Ar, HCl, CO, N(2), O(2), and H(2)O, and to obtain dispersion energy coefficients for their pair interactions. Results obtained with augmented quadruple-zeta basis sets agree well with the best literature values of the C(6) dispersion energy coefficients. Time-dependent coupled cluster with single and double excitations theory will be useful as a benchmark for evaluating more approximate theories. (c) 2007 Wiley Periodicals, Inc. J Comput Chem, 2008.

A five-dimensional potential energy surface and predicted infrared spectra for the N2O-hydrogen complexes

We present a five-dimensional potential energy surface for the N(2)O-hydrogen complex using supermolecular approach with the full counterpoise correction at the coupled-cluster singles and doubles with noniterative inclusion of connected triple level. The normal mode Q(3) for the nu(3) antisymmetric stretching vibration of the N(2)O molecule was included in the calculations of the potential energies. The radial discrete variable representation/angular finite basis representation method and Lanczos algorithm were employed to calculate the rovibrational energy levels for four species of N(2)O-hydrogen complexes (N(2)O-para-H(2), -ortho-H(2), -ortho-D(2), and -para-D(2)) without separating the inter- and intramolecular vibrations. The calculated band origins are all blueshifted relative to the isolated N(2)O molecule and in good agreement with the experimental values. The calculated rotational spectroscopic constants and molecular structures agree well with the available experimental results. The frequencies and line intensities of the rovibrational transitions in the nu(3) region of N(2)O for the van der Waals ground vibrational state were calculated and compared with the observed spectra. The predicted infrared spectra are consistent with the observed spectra and show that the N(2)O-H(2) complexes are mostly a-type transitions while both a-type and b-type transitions are significant for the N(2)O-D(2) complexes.

Three-dimensional ab initio potential-energy surface and rovibrational spectra of the H2–Kr complex

We report a three-dimensional ab initio potential-energy surface for the H2-Kr complex calculated using a supermolecular method. The electronic calculations were performed at the coupled-cluster singles and doubles level with noniterative inclusion of connected triples levels with a large basis set including midbond functions and the full counterpoise correction for the basis-set superposition error. The intermolecular potential energy between the H2 molecule and the Kr atom were evaluated at five potential-optimized discrete variable representation (DVR) grid points generated from the potential-energy curve of H2. The potential for other bond lengths of H2 could be deduced using polynomial interpolations. The complex is found to have a linear preferred structure with a rather flat energy barrier. The three-dimensional DVR method and the Lanczos propagation algorithm were employed to calculate the rovibrational states without separating the inter- and intramolecular nuclear motions. In addition, the rovibrational spectra from the H2 fundamental vibrational band were calculated. The calculated shift for the band origin is -1.50 cm-1, which is in good agreement with the experimental value of -1.706 cm-1, and the calculated transition frequencies in Q1(0) and S1(0) bands are within 3% of the observed values.

Spectroscopy of Ar–SH and Ar–SD. II. Determination of the three-dimensional intermolecular potential-energy surface

All the pure rotational transitions reported in the previous studies [J. Chem. Phys. 113, 10121 (2000); J. Mol. Spectrosc. 222, 22 (2003)] and newly observed rotation-vibration transitions, P = 1/2 <-- 3/2, for Ar-SH and Ar-SD [J. Chem. Phys. (2005), the preceding paper] have been simultaneously analyzed to determine a new intermolecular potential-energy surface of Ar-SH in the ground state. A Schrodinger equation considering the three-dimensional freedom of motion for an atom-diatom complex in the Jacobi coordinate, R, theta, and r, was numerically solved to obtain energies of the rovibrational levels using the discrete variable representation method. A three-dimensional potential-energy surface is determined by a least-squares fitting with initial values of the parameters for the potential obtained by ab initio calculations at the RCCSD(T)/aug-cc-pVTZ level of theory. The potential well reproduces all the observed data in the microwave and millimeter wave regions with parity doublings and hyperfine splittings. Several low-lying rovibrational energies are calculated using the new potential-energy surface. The dependence of the interaction energy between Ar and SH(2pi(i)) on the bond length of the SH monomer is discussed.