Solvation of fluoroform and fluoroform–dimethylether dimer in liquid krypton: A theoretical cryospectroscopic study

Abstract

A hybrid, sequential statistical physics-quantum mechanical electronic-quantum mechanical nuclei approach has been applied to study the C-H stretching frequencies of bare fluoroform dissolved in liquid krypton under cryogenic conditions (at ~130 K), as well as upon blue shifting hydrogen bonding interactions with dimethylether in the same solvent. The structure of the liquid at 130 K was generated by Monte Carlo simulations of cryogenic Kr solutions containing either fluoroform or fluoroform and dimethylether molecules. Statistically uncorrelated configurations were appropriately chosen from the equilibrated MC runs and supermolecular clusters containing solute and solvent molecules (either standalone or embedded in the "bulk" part of the solvent treated as a polarizable continuum) were subjected to quantum mechanical electronic (QMel) and subsequent quantum mechanical nuclei (QMnuc) calculations. QMel calculations were implemented to generate the in-liquid 1D intramolecular C-H stretching vibrational potential of the fluoroform moiety and subsequently in the QMnuc phase the corresponding anharmonic C-H stretching frequency was computed by diagonalization techniques. Finally, the constructed vibrational density of states histograms were compared to the experimental Raman bands. The calculated anharmonic vibrational frequency shifts of the fluoroform C-H stretching mode upon interaction with dimethylether in liquid Kr are in very good agreement with the experimental data (20.3 at MP2 level vs. 16.6 cm(-1) experimentally). Most of this relatively large frequency blue shift is governed by configurations characterized by a direct C-H···O contact between monomers. The second population detected during MC simulations, characterized by reversed orientation of the monomers, has a minor contribution to the spectral appearance. The experimentally observed trend in the corresponding bandwidths is also correctly reproduced by our theoretical approach. Solvation of the fluoroform monomer, according to experiment, results in small C-H stretching frequency red shift (~-2 cm(-1)), while our approach predicts a blue shift of about 10 cm(-1). By a detailed analysis of the anharmonic C-H stretching frequency dependence on the position of the nearest solvent krypton atom and also by analyzing the vibrational Stark effect induced by the local fluctuating field component parallel to the C-H axis, we have derived several conclusions related to these observations. The frequency vs. C···Kr distance dependence shows appreciable fluctuations and even changes in sign at R values close to the maximum of the C···Kr radial distribution function, so that most of the first-shell Kr atoms are located at positions at which the CH frequency shifts acquire either small negative or small positive values. It so happens, therefore, that even the actual sign of the frequency shift is strongly dependent on the correct description of the first solvation shell around CF3H by the Monte Carlo method, much more than the other in-liquid properties calculated by similar approaches.