HF Monomer Research Articles

HF trimer: 12D fully coupled quantum calculations of HF-stretch excited intramolecular and intermolecular vibrational states using contracted bases of intramolecular and intermolecular eigenstates.

We present the computational methodology, which for the first time allows rigorous twelve-dimensional (12D) quantum calculations of the coupled intramolecular and intermolecular vibrational states of hydrogen-bonded trimers of flexible diatomic molecules. Its starting point is the approach that we introduced recently for fully coupled 9D quantum calculations of the intermolecular vibrational states of noncovalently bound trimers comprised of diatomics treated as rigid. In this paper, it is extended to include the intramolecular stretching coordinates of the three diatomic monomers. The cornerstone of our 12D methodology is the partitioning of the full vibrational Hamiltonian of the trimer into two reduced-dimension Hamiltonians, one in 9D for the intermolecular degrees of freedom (DOFs) and another in 3D for the intramolecular vibrations of the trimer, and a remainder term. These two Hamiltonians are diagonalized separately, and a fraction of their respective 9D and 3D eigenstates is included in the 12D product contracted basis for both the intra- and intermolecular DOFs, in which the matrix of the full 12D vibrational Hamiltonian of the trimer is diagonalized. This methodology is implemented in the 12D quantum calculations of the coupled intra- and intermolecular vibrational states of the hydrogen-bonded HF trimer on an abinitio calculated potential energy surface (PES). The calculations encompass the one- and two-quanta intramolecular HF-stretch excited vibrational states of the trimer and low-energy intermolecular vibrational states in the intramolecular vibrational manifolds of interest. They reveal several interesting manifestations of significant coupling between the intra- and intermolecular vibrational modes of (HF)3. The 12D calculations also show that the frequencies of the v = 1, 2 HF stretching states of the HF trimer are strongly redshifted in comparison to those of the isolated HF monomer. Moreover, the magnitudes of these trimer redshifts are much larger than that of the redshift for the stretching fundamental of the donor-HF moiety in (HF)2, most likely due to the cooperative hydrogen bonding in (HF)3. The agreement between the 12D results and the limited spectroscopic data for the HF trimer, while satisfactory, leaves room for improvement and points to the need for a more accurate PES.

Theoretical and matrix isolation studies of infrared spectra of the H2CO∙∙∙HF hydrogen-bonded complex

The equilibrium geometry, the binding energy, and harmonic and anharmonic spectral parameters of H2CO and HF monomers and the H2CO∙∙∙HF complex are calculated in the MP2/6-311++G(3df,3pd) approximation with the basis set superposition error taken into account. Anharmonic calculations are carried out using the second-order vibrational perturbation theory. The changes in transition frequencies and intensities of monomeric IR bands upon formation of the 1:1 complex are analyzed. A high intensity value (53 km mol–1) is predicted for the overtone band of the out-of-plane HF libration. The variational method was used to solve anharmonic vibrational problems in 1D–4D subspaces involving the H–F stretch. The influence of anharmonicity on spectral and structural parameters was studied. The variational calculation confirmed the high intensity value for the overtone of the HF librational mode from the perturbative calculation. Matrix-isolation experiments were performed in N2 at T = 8 K to record the absorption spectra of pure H2CO and HF and H2CO/HF mixtures in the IR region. The spectral features related to the H–F, C=O, and C–H stretches and the in-plane HF librational motion were reliably identified. The predictions of our perturbative calculations for isolated compounds, especially for the band shifts upon complexation, are in satisfactory agreement with the experimental matrix-isolation findings. The only exception is the H–F stretching mode in the complex because of the high anharmonicity of HF and strong interaction of free HF molecules with a matrix environment. Of particular interest is the observation of a rather strong band at 1069 cm–1, which is close to 1072 cm–1 predicted for the overtone of HF out-of-plane librational mode.

Additivity and non-additivity of dissociation energies in intermolecular interactions. Theoretical studies on (H2)n, n = 2-8, (CO2)n, n = 2-6 and (HF)n, n = 2-8

CCSD(T) and MP2 results using the aug-cc-pV5Z basis set are reported for chain, cyclic and other structures of the clusters (H2)n, n = 2-8, (CO2)n, n = 2-6 and (HF)n, n = 2-8. In chain-like structures of (H2)n and (CO2)n, with the bonding type of the dimer maintained, the dissociation energy De of the dimer doubles for the trimer, triples for the tetramer, and so on. Due to these systems being dominated by short-range forces, interactions are essentially restricted to neighbouring monomers. For other types of (H2)n and (CO2)n structures, the multipliers relative to the dimerisation energy can be much higher. Dissociation energies for the hexamers in S6 symmetry of both H2 (379 cm−1) and CO2 (4925 cm−1) are over ten times the respective dimerisation energies. For the chain-like trimer of HF, however, De is in excess of double the dimer value. Mainly due to longer-range dipolar forces, the interactions reach beyond the neighbouring monomers. The interaction energy between HF monomers in chains follows an approximate R−2 (R being the F–F distance) relationship, The calculated dissociation energies of the HF octamer are 15,985 cm−1 (factor of 10.4) for the chain, and 21,003 cm−1 (factor of 13.7) for the C6h cyclic structure.

Structures, Stabilities, and Spectral Properties of Endohedral Borospherenes M@B40 0/- (M = H2, HF, and H2O).

The discovery of borospherene B40 leads to a new beginning for the study of boron chemistry and may lead to new boron-based nanomaterials. Based on density functional theory, the structures, electronic properties, infrared and Raman spectra, photoelectron spectra, and electronic absorption spectra of endohedral borospherenes M@B400/– (M = H2, HF, and H2O) are investigated. It is found that H2, HF, and H2O monomers can form stable endohedral borospherenes M@B400/– (M = H2, HF, and H2O). In addition, the calculated results indicate that the doped molecule at the off-center location can relax to the center location within the cage and the symcenter of the doped molecule is almost located in the center of the cage. Unlike endohedral metalloborospherene Ca@B40, which is a charge-transfer complex between Ca2+ and B402–, natural population analyses and chemical bonding analyses reveal that there is no significant charge transfer of the doped molecule. The calculated spectra indicate that doping of a molecule (H2, HF, or H2O) in borospherene B40 can change the photoelectron spectra and doping of a polar molecule (HF or H2O) in borospherene B40 can change the spectral properties. For instance, the addition of a molecule can increase infrared and Raman-active modes and cause a red shift or blue shift of electronic spectra. These spectral features can be compared with future experimental values of endohedral borospherenes M@B400/– (M = H2, HF, and H2O).

Intermolecular interactions and proton transfer in the hydrogen halide-superoxide anion complexes.

The superoxide radical anion O2(-) is involved in many important chemical processes spanning different scientific disciplines (e.g., environmental and biological sciences). Characterizing its interaction with various substrates to help elucidate its rich chemistry may have far reaching implications. Herein, we investigate the interaction between O2(-) (X[combining tilde] (2)Πg) and the hydrogen halides (X[combining tilde] (1)Σ) with coupled-cluster theory. In contrast to the short (1.324 Å) hydrogen bond formed between the HF and O2(-) monomers, a barrierless proton transfer occurs for the heavier hydrogen halides with the resulting complexes characterized as long (>1.89 Å) hydrogen bonds between halide anions and the HO2 radical. The dissociation energy with harmonic zero-point vibrational energy (ZPVE) for FHO2(-) (X[combining tilde] (2)A'') → HF (X[combining tilde] (1)Σ) + O2(-) (X[combining tilde] (2)Πg) is 31.2 kcal mol(-1). The other dissociation energies with ZPVE for X(-)HO2 (X[combining tilde] (2)A'') → X(-) (X[combining tilde] (1)Σ) + HO2 (X[combining tilde] (2)A'') are 25.7 kcal mol(-1) for X = Cl, 21.9 kcal mol(-1) for X = Br, and 17.9 kcal mol(-1) for X = I. Additionally, the heavier hydrogen halides can form weak halogen bonds H-XO2(-) (X[combining tilde] (2)A'') with interaction energies including ZPVE of -2.3 kcal mol(-1) for HCl, -8.3 kcal mol(-1) for HBr, and -16.7 kcal mol(-1) for HI.

Complexation of Zr and Hf monomers in supercritical aqueous solutions: Insights from ab initio molecular dynamics simulations

The enhanced mobilization of high field strength elements (HFSE) by certain geological fluids at high pressures, P, and high temperatures, T, as indicated by field observations and recent solubility measurements is most likely related to the type of complex formation in the fluid phase. However, only a few in situ experimental studies have been performed so far to constrain HFSE speciation in supercritical aqueous fluids. Here, we complement these investigations by ab initio molecular dynamics simulations to explore the complexation of Zr4+ and Hf4+ monomers in aqueous solutions at T=1000K and P∼1GPa. Regardless of the type of fluid (basic, neutral, acidic) both Zr and Hf seem to prefer an octahedral coordination. By systematically changing fluid composition and coordinating anions, the structure of various [(Zr,Hf)Clk(OH)m(H2O)n]4−k−m(k=0–6,m=2–6,n=0–4) complexes is investigated. For each complex that is stable on the time scale of the simulation, theoretical X-ray absorption spectra are calculated and compared to experimental data. From this comparison, the most likely complexes in HCl solutions are [(Zr,Hf)Cl3(Cl,OH)2(H2O)]−. The differences between experimental and theoretical XANES spectra present in the case of NaOH solutions indicate HFSE speciation beyond monomers in this system. Charge-neutral [(Zr,Hf)(OH)4(H2O)1–2] complexes, which may be the dominant species in neutral fluids at very low HFSE concentrations, show the lowest coordination (5.2(1)) of all investigated species. Finally, differences in the complexation of Zr and Hf in aqueous solutions at ambient and supercritical conditions and the possible formation of more complex oligomeric species are discussed.

Interatomic Coulombic decay in (n = 2–3) clusters using CAP/EOM-CCSD method

The equation-of-motion coupled-cluster method along with the complex absorbing potential has been applied to study the interatomic Coulombic decay mechanism in hydrogen-bonded clusters. We have applied this method to calculate the lifetime of the F 2s inner-valence ionised state of (n = 2–3) clusters. The lifetime is found to be very short and decreases substantially with increasing the number of HF monomer.

Correlated ab initio investigations on the intermolecular and intramolecular potential energy surfaces in the ground electronic state of the ${\rm O}_2^ - ({\rm X}{}^2\Pi _{\rm g}) - {\rm HF}({\rm X}{}^1\Sigma^ +)$O2−(XΠg2)− HF (XΣ+1) complex

This work reports the first highly correlated ab initio study of the intermolecular and intramolecular potential energy surfaces in the ground electronic state of the O(2)(-)(X(2)Π(g))-HF(X(1)Σ(+)) complex. Accurate electronic structure calculations were performed using the coupled cluster method including single and double excitations with addition of the perturbative triples correction [CCSD(T)] with the Dunning's correlation consistent basis sets aug-cc-pVnZ, n = 2-5. Also, the explicitly correlated CCSD(T)-F12a level of theory was employed with the AVnZ basis as well as the Peterson and co-workers VnZ-F12 basis sets with n = 2 and 3. Results of all levels of calculations predicted two equivalent minimum energy structures of planar geometry and C(s) symmetry along the A" surface of the complex, whereas the A' surface is repulsive. Values of the geometrical parameters and the counterpoise corrected dissociation energies (Cp-D(e)) that were calculated using the CCSD(T)-F12a/VnZ-F12 level of theory are in excellent agreement with those obtained from the CCSD(T)/aug-cc-pV5Z calculations. The minimum energy structure is characterized by a very short hydrogen bond of length of 1.328 Å, with elongation of the HF bond distance in the complex by 0.133 Å, and D(e) value of 32.313 Kcal/mol. Mulliken atomic charges showed that 65% of the negative charge is localized on the hydrogen bonded end of the superoxide radical and the HF unit becomes considerably polarized in the complex. These results suggest that the hydrogen bond is an incipient ionic bond. Exploration of the potential energy surface confirmed the identified minimum and provided support for vibrationally induced intramolecular proton transfer within the complex. The T-shaped geometry that possesses C(2v) symmetry presents a saddle point on the top of the barrier to the in-plane bending of the hydrogen above and below the axis that connects centers of masses of the monomers. The height of this barrier is 7.257 Kcal/mol, which is higher in energy than the hydrogen bending frequency by 909.2 cm(-1). The calculated harmonic oscillator vibrational frequencies showed that the H-F stretch vibrational transition in the complex is redshifted by 2564 cm(-1) and gained significant intensity (by at least a factor of 30) with respect to the transition in the HF monomer. These results make the O(2)(-)-HF complex an excellent prototype for infrared spectroscopic investigations on open-shell complexes with vibrationally induced proton transfer.

Calculation of interaction of the stretching and bending vibrations of HF in the hydrogen-bonded complex [F(HF)2

Vibrations of the [F(HF)2]− complex are calculated with allowance for the anharmonic interactions of the stretching vibrations of HF monomers and their rotations about the centers of gravity of HF in the plane of the complex. A four-dimensional vibrational Schrodinger equation is solved using a potential energy surface calculated in the MP2/6-311++G(3df,3pd) approximation with the superposition of atomic functions of the monomers taken into account. The equilibrium and vibrationally averaged structures of the complex are determined. The frequencies and intensities for spectral transitions from the ground state to a number of excited vibrational states are calculated. It is shown that, due to resonances between the excited states of the stretching modes and doubly excited states of the bending modes, the overtone transitions associated with the bending modes borrow a significant part of the intensity of fundamental stretching transitions.

Spectroscopy of HF and HF-Containing Clusters in Solid Parahydrogen

We report measurements of FT-IR absorption spectroscopy of HF, DF, and their clusters in solid parahydrogen (pH(2)). The observed spectra contain many absorption lines which were assigned to HF monomers, HF polymers, and clusters with other species, such as N(2), O(2), orthohydrogen (oH(2)), etc. The rotational constants of HF and DF monomers were determined from the cooperative transitions of the vibration of solid pH(2) and the rotation of HF and DF. Small reduction of the rotational constants indicates that HF and DF are nearly free rotors in solid pH(2). Time dependence of the spectra suggests that HF and DF monomers migrate in solid pH(2) and form larger polymers, probably via tunneling reactions through high energy barriers on inserting another monomer to the polymers. The line width of HF monomers in solid pH(2) was found to be 4 cm(-1), which is larger than that of other hydrogen halides in solid pH(2). This broad line width is explained by rapid rotational relaxation due to the accidental coincidence between the rotational energy of HF and the phonon energy with maximum density of states of solid pH(2) and the rotational-translational coupling in a trapping site.

Axial and equatorial hydrogen-bond conformers between (CH2)3S and H(D)F: Fourier transform infrared spectroscopy and ab initio calculations

The coexistence of axial and equatorial hydrogen-bonded conformers of 1 : 1 (CH(2))(3)S-HF (and -DF) has been observed in the same adiabatic expansion of a supersonic jet seeded with argon and in a static absorption cell at room temperature. High level calculations computed the axial conformer to be the most stable one with a small energy difference with respect to the equatorial one, in full agreement with previous microwave experiments. On the grounds of band contour simulations of FTIR spectra and ab initio energetic and anharmonic vibrational calculations, two pairs of ν(s) HF donor stretching bands, observed in a series of jet-FTIR spectra at 3457.9 and 3480.5 cm(-1) have been respectively assigned to the axial and equatorial forms of the 1 : 1 complex. In the jet-FTIR spectra series with HF, the assignment of an additional broad band (about 200 cm(-1) higher in frequency with respect to ν(s)) to a 1 : 2 complex has been supported by theoretical investigations. Experimental detection of both axial and equatorial forms of a cyclic trimer has been confirmed by calculated energetic and vibrational properties. The nature of hydrogen bonding has been examined within topological frameworks. The energetic partitioning within the 1 : 1 dimers has been elucidated with SAPT techniques. Interestingly, the interconversion pathway between two 1 : 1 structures has been explored and it was seen that the formation of the 1 : 1 complex affects the interconversion barrier on the ring puckering motion. The band contour analysis of gas phase FTIR experiments provided a consistent set of vibrational frequencies and anharmonic coupling constants, in good agreement with ab initio anharmonic vibrational calculations. Finally, from a series of cell-FTIR spectra recorded at different partial pressures of (CH(2))(3)S and HF monomers, the absorption signal of the 1 : 1 complex could be isolated which enabled to estimate the equilibrium constant K(p) = 0.023 at 298 K for the dimerization.

Infrared laser spectroscopy of CH3⋯HF in helium nanodroplets: The exit-channel complex of the F+CH4 reaction

High-resolution infrared laser spectroscopy is used to study the CH3...HF and CD3...HF radical complexes, corresponding to the exit-channel complex in the F + CH4 --> HF + CH3 reaction. The complexes are formed in helium nanodroplets by sequential pickup of a methyl radical and a HF molecule. The rotationally resolved spectra presented here correspond to the fundamental v = 1 <-- 0 H-F vibrational band, the analysis of which reveals a complex with C(3v) symmetry. The vibrational band origin for the CH3...HF complex (3797.00 cm(-1)) is significantly redshifted from that of the HF monomer (3959.19 cm(-1)), consistent with the hydrogen-bonded structure predicted by theory [E. Ya. Misochko et al., J. Am. Chem. Soc. 117, 11997 (1995)] and suggested by previous matrix isolation experiments [M. E. Jacox, Chem. Phys. 42, 133 (1979)]. The permanent electric dipole moment of this complex is experimentally determined by Stark spectroscopy to be 2.4+/-0.3 D. The wide amplitude zero-point bending motion of this complex is revealed by the vibrational dependence of the A rotational constant. A sixfold reduction in the line broadening associated with the H-F vibrational mode is observed in going from CH3...HF to CD3...HF. The results suggest that fast relaxation in the former case results from near-resonant intermolecular vibration-vibration (V-V) energy transfer. Ab initio calculations are also reported (at the MP2 level) for the various stationary points on the F + CH4 surface, including geometry optimizations and vibrational frequency calculations for CH3...HF.

High-rank quantum topological electrostatic potential: Molecular dynamics simulation of liquid hydrogen fluoride

Five bulk properties of liquid hydrogen fluoride were obtained by a molecular dynamics simulation at five temperatures between 203 and 273 K. The rigid-body interaction potential was designed with emphasis on the electrostatic contribution. The potential was expanded in spherical tensor multipole moments, and its corresponding Ewald summations, containing all possible interaction up to L=5 (i.e., monopole–hexadecapole, dipole–octopole, quadrupole–quadrupole). The moments were generated by a quantum chemical topological analysis of the electron density of the geometry-optimized monomer at the MP2/aug-cc-pVTZ level. Only two parameters in the Lennard-Jones part of the potential were adjusted to the experimental density and radial distribution functions. Agreement with experiment is excellent for the total energy and the density, and reasonable (with even overall performance) for the diffusion coefficient, the isobaric heat capacity and the thermal expansion coefficient. Chains of hydrogen-bonded HF monomers up to 10 molecules are quite common in the liquid and they more or less retain the relative orientation of the HF monomers in the gas-phase dimer.

Spectrum and vibrational predissociation of the HF dimer. I. Bound and quasibound states

We present full six-dimensional calculations of the bound states of the HF dimer for total angular momentum J=0,1 and of the quasibound states for J=0 that correspond with vibrational excitation of one of the HF monomers, either the donor or the acceptor in the hydrogen bond. Transition frequencies and rotational constants were calculated for all four molecular symmetry blocks. A contracted discrete variable representation basis was used for the dimer and monomer stretch coordinates R,rA,rB; the generation of the monomer basis in the dimer potential leads to significantly better convergence of the energies. We employed two different potential energy surfaces: the SQSBDE potential of Quack and Suhm and the SO-3 potential of Klopper, Quack, and Suhm. The frequencies calculated with the SO-3 potential agree very well with experimental data and are significantly better than those from the SQSBDE potential.

Calculating energy levels of isomerizing tetra-atomic molecules. I. The rovibrational bound states of Ar2HF

A general, six-dimensional computational method for the accurate calculation of rotationally and vibrationally excited states of tetra-atomic molecules is developed. The resulting program is particularly appropriate for molecules executing wide-amplitude motions and isomerizations. An application to the Ar2HF van der Waals trimer is presented in which the HF intramolecular stretching coordinate is separated out adiabatically and is not treated explicitly. Vibrational term values up to about 100 cm−1 with absolute convergence to better than 0.1 cm−1 are reported. These calculations employ more extensive vibrational basis sets and hence consider a much higher density of states than hitherto. States that sample Ar–Ar–HF linear configurations and approach Ar–HF–Ar linear configurations are characterized for the first time. Results for total angular momentum J=0 and 1 provide the first accurate calculations of rotational constants for this system. The rotational constants for the HF bending states of Ar2HF in the ground and first vibrationally excited states of the HF monomer are in good agreement with experiment, confirming the accuracy of the potential used in this work.

Solvated electrons in very small clusters of polar molecules: (HF)(3)(-).

Photoelectron spectra of (HF)(3)(-) reveal coexistence of two anionic isomers with vertical electron detachment energies (VDE) of 0.24 and 0.43 eV. The results of electronic-structure calculations, performed at the coupled cluster level of theory with single, double, and noniterative triple excitations, suggest that the two isomers observed experimentally are an open, zigzag, dipole-bound anion and an asymmetric solvated electron, in which the dipole-bound anion of (HF)(2) is solvated by one HF monomer at the side of the excess electron. The theoretical VDE of 0.21 and 0.44 eV, respectively, are in excellent agreement with the experimental data.

Extended Coupled Cluster Approach for Molecular Properties: Study of H2O and HF Complexes

In this paper, we study stationary variant of extended coupled-cluster response approach for properties. This has been studied at the singles and doubles approximation using cubic-truncated functional. This approximation has been studied earlier around equilibrium for small molecules. In this paper, efficacy of this approximation has been shown using perturbative arguments. Further we have calculated dipole moments and polarizabilities of weakly interacting dimers of HF, H2O and H2O - HF complex. Results of HF and H2O monomers have been presented at the same level for comparison. The results have been compared with experimental results, wherever available and other theoretical results.

Convergence of variational calculations of the ground-state energy of HF dimer

We report large-scale quantum mechanical calculations for the ground-state energy E 0 of the HF dimer on the SQSBDE potential energy surface of Quack and Suhm to compare with the recent calculations of Mladenovic and Lewerenz, who recently reported significantly lower energies than all previously published results. We confirm their results within 0.01 cm −1 with bases of 21 824–29 541 functions built on asymptotic HF monomer vibrational eigenfunctions.

Photodissociation of cyclic HF complexes: Pentamer through heptamer

Photofragment angular and final rotational state distributions have been measured, corresponding to the vibrational predissociation of HF polymers [(HF)N N=5–7]. The present experiments confirm the theoretical prediction that, for clusters in this size range, the only open dissociation channel corresponds to the “evaporation” of a single HF monomer. The pump–probe experiments reported here indicate that the resulting HF monomer is produced with only modest rotational excitation. Rotational temperatures of 180 and 140 K have been determined for the evaporated HF monomer following dissociation of the HF pentamer and hexamer, respectively. Energy conservation in these photodissociation studies provides upper limits for the dissociation energies (D0), namely, 2941 cm−1 and 2854 cm−1 for the pentamer and hexamer, respectively.

The infrared spectrum and internal rotation barrier in HF–BF 3

The rotationally resolved infrared spectrum of the HF stretch fundamental in combination with one quantum of internal rotation excitation has been investigated for the HF–BF 3 complex using a color-center-laser optothermal molecular-beam spectrometer. The infrared band origin is red-shifted by 3.9 cm −1 from the HF monomer origin. The magnitude of the Coriolis interaction in the m=±1 internal rotor state observed in combination with the HF stretch gives a barrier to internal rotation of 54 cm −1 in the HF stretch excited state.