A Monte Carlo-quantum mechanics study of the solvatochromic shifts of the lowest transition of benzene

Exploring the molecular solvatochromism, stability, reactivity, and non-linear optical response of resveratrol

CONTEXT

This work analyzes the isomerization effects and solvent contributions to the stability, electronic excitations, reactivity, and non-linear optical properties (NLO) of resveratrol molecules within the formalism of the Density Functional Theory. The findings suggest that resveratrol solvatochromism is significantly influenced by solvent polarization. The electronic and free energies (E and G) indicate that trans is the most stable conformer. The system is classified as a strong nucleophile. However, the analysis of the Fukui functions and the Mulliken charges indicate that cis-trans isomerization jointly affects the reactive indices of the carbon and hydrogen atoms. The results also suggest that solvent is relevant to solvatochromism and the NLO response. Both cis and trans conformers present strong π - π ∗ excitations that undergo a visible hypsochromic change when the polarity of the solvent increases. Once the absorption spectra are connected to the first hyperpolarization ( β ) by the Oudar and Chemla relation, the hypsochromism of resveratrol is the reason for the drop in the generation of the second harmonic when the ambient polarity decreases. The CAM-B3LYP DFT results suggest that resveratrol is interesting for NLO applications. Depending on the choice of solvent, values ∼ 50 times those observed for urea ( β = 0.34 × 10 - 34 esu), which is a standard NLO material.

METHODS

The optimized geometries of cis and trans isomers of resveratrol in vacuum were obtained using Density Functional Theory (DFT) with the hybrid exchange-correlation function (CAM-B3LYP) and Pople basis set functions, specifically 6-311++G(d,p). The solvent effect on the geometries of both isomers was included using the polarizable continuum model (PCM) with the same level of QM calculation. Vibrational analysis was conducted to confirm that all optimized geometries correspond to the minimum energy. Various electronic properties, including dipole moments, molecular orbitals, transition energy, dipole polarizabilities, and global reactivity parameters, were calculated using both continuum and discrete solvation models based on the sequential QM/MM methodology. All QM calculations were performed with the Gaussian 09 program and the MC simulations with the DICE program. All NLO analysis was carried out using the Multiwfn code.

Novel Approach for Predicting Vertical Electron Attachment Energies in Bulk-Solvated Molecules

When low-energy electrons interact with molecules, they can give rise to transient anion states commonly known as resonances. These states are formed through vertical electron attachment processes and have the potential to induce various forms of DNA lesions, including base damage, single- and double-strand breaks, cross-links, and clustered lesions that are challenging to repair. So far, most experimental and theoretical studies have investigated the formation of resonances of (bio)molecules in the gas phase or in microsolvated environments. Since cellular environments are mainly composed of water molecules, it is crucial to understand how bulk water affects the resonances of (bio)molecules. Given the existing gap in studies on resonances of bulk-solvated molecules, we propose a novel theoretical-computational approach to address this void. Our approach combines the multibasis-set (time-dependent-)density functional theory and self-consistent sequential quantum mechanics/molecular mechanics polarizable electrostatic embedding methods. We apply this combined methodology to predict the vertical electron attachment energies of 1-methyl-5-nitroimidazole (1M5NI), a well-known radiosensitizer model, in bulk water. In addition, we analyze the rapid mutual polarization between the resonances (both shape- and core-excited) of 1M5NI and the surrounding bulk water environment. For comparison, we also studied the isolated and microsolvated 1M5NI. Overall, while the polarization of the environment is clearly sensitive to the solute charge, causing a significant impact on the vertical electron affinity and consequently on the attachment electron energies, it does not have a significant impact on the excitation energies of the anion.

Structure and Ultrafast X-ray Diffraction of the Hydrated Metaphosphate

We study the pathway of metaphosphate hydration when a metaphosphate anion is dissolved in liquid water with an explicit water model. For this purpose, we propose a sequential Monte Carlo algorithm incorporated with the ab initio quantum mechanics/molecular mechanics (QM/MM) method, which can reduce the amount of ab initio QM/MM sampling while retaining the accuracy of the simulation. We demonstrate the numerical calculation of the standard enthalpy change for the successive addition reaction PO3-·2H2O + H2O ⇌ PO3-·3H2O in the liquid phase, which helps to clarify the hydration pathway of the metaphosphate. With the obtained hydrated structure of the metaphosphate anion, we perform ab initio calculations for its relaxation dynamics upon vibrational excitation and characterize the energy transfer process in solution with simulated ultrafast X-ray diffraction signals, which can be experimentally implemented with X-ray free-electron lasers.

Enhanced Sampling in Molecular Dynamics Simulations: How Many MD Snapshots can be Needed to Reproduce the Biological Behavior?

Since its early days in the 19th century, medicinal chemistry has concentrated its efforts on the treatment of diseases, using tools from areas such as chemistry, pharmacology, and molecular biology. The understanding of biological mechanisms and signaling pathways is crucial information for the development of potential agents for the treatment of diseases mainly because they are such complex processes. Given the limitations that the experimental approach presents, computational chemistry is a valuable alternative for the study of these systems and their behavior. Thus, classical molecular dynamics, based on Newton's laws, is considered a technique of great accuracy, when appropriated force fields are used, and provides satisfactory contributions to the scientific community. However, as many configurations are generated in a large MD simulation, methods such as Statistical Inefficiency and Optimal Wavelet Signal Compression Algorithm are great tools that can reduce the number of subsequent QM calculations. Accordingly, this review aims to briefly discuss the importance and relevance of medicinal chemistry allied to computational chemistry as well as to present a case study where, through a molecular dynamics simulation of AMPK protein (50 ns) and explicit solvent (TIP3P model), a minimum number of snapshots necessary to describe the oscillation profile of the protein behavior was proposed. For this purpose, the RMSD calculation, together with the sophisticated OWSCA method was used to propose the minimum number of snapshots.

Theoretical investigation of solvent and oxidation/deprotonation effects on the electronic structure of a mononuclear Ru-aqua-polypyridine complex in aqueous solution

Mononuclear polypyridine ruthenium (Ru) complexes can catalyze various reactions, including water splitting, and can also serve as photosensitizers in solar cells. Despite recent progress in their synthesis, accurately modeling their physicochemical properties, particularly in solution, remains challenging. Herein, we conduct a theoretical investigation of the structural and electronic properties of a mononuclear Ru-aqua polypyridine complex in aqueous solution, considering five of its possible oxidation/protonation states species: [RuII(H2O)(py)(bpy)2]2+, [RuII(OH)(py)(bpy)2]+, [RuIII(H2O)(py)(bpy)2]3+, [RuIII(OH)(py)(bpy)2]2+ and [RuIV(O)(py)(bpy)2]2+, where py = pyridine and bpy = 2,2'-bipyridine. At first, we investigate the impact of proton-coupled and non-coupled electron transfer reactions on the geometry and electronic structure of the complexes in vacuum and in solution, using an implicit solvent model. Then, using a sequential multiscale approach that combines quantum mechanics and molecular mechanics (S-QM/MM), we examine the explicit solvent effects on the electronic excitations of the complexes, and compare them with the experimental results. The complexes were synthesized, and their absorption spectra measured in aqueous solution. To accurately describe the QM interactions between the metal center and the aqueous ligand in the MM simulations, we developed new force field parameters for the Ru atom. We analyze the solvent structure around the complexes and account for its explicit influence on the polarization and electronic excitations of the complexes. Notably, accounting for the explicit solvent polarization effects of the first solvation shells is essential to correctly describe the energy of the electronic transitions, and the explicit treatment of the hydrogen bonds at the QM level in the excitation calculations improves the accuracy of the description of the metal-to-ligand charge-transfer bands. Transition density matrix analysis is used to characterize all electronic transitions in the visible and ultraviolet ranges according to their charge-transfer (CT) character. This study elucidates the electronic structure of those ruthenium polypyridyl complexes in aqueous solution and underscores the importance of precisely describing solvent effects, which can be achieved employing the S-QM/MM method.

Reactivity properties of the HOSO and HSO2 isomers in liquid medium: a sequential Monte Carlo/quantum mechanics study

CONTEXT

The rationalization of acid rain formation steps is fundamental for mitigating its effects. It is believed the hydroxysulfinyl radical is an intermediate species for the production of atmospheric sulfuric acid. Two stable configurations HOSO and HSO₂ have been reported for such a radical in the gas phase. This work aims at studying these isomers in the aqueous medium. The electrical and reactivity quantities - electronic chemical potential ([Formula: see text]), chemical hardness ([Formula: see text]), and electrophilicity ([Formula: see text]) - are here calculated and compared. Considering first solvation shells (15 H₂O for HSO₂ and 23 H₂O for HOSO), an increase of 25% in the dipole moment of HSO₂ was obtained, while the dipole moment of HOSO decreases in 11%. Both solvated isomers grow softer ([Formula: see text] decreases) contrasted to the gas phase.

METHODS

HOSO and HSO₂ are studied through a sequential Monte Carlo/quantum mechanics approach. Lennard-Jones plus the Coulomb potentials were used to represent intermolecular potential interaction in the frame of the DICE package. Molecular structure calculations were performed at CASPT2/aug - cc - pV(T + d)Z level of theory using the MOLPRO suite of programs.

Investigation of the Second Harmonic Generation at the Water-Vacuum Interface by Using Multi-Scale Modeling Methods

The Sequential Quantum Mechanics/Molecular Mechanics scheme has been enacted to perform a systematic investigation of the polarizability (α) and first hyperpolarizability (β) responses at the water-vacuum interface. After performing classical molecular dynamics simulations to provide snapshots of the structures, quantum chemistry calculations of the linear and nonlinear optical responses have been performed for clusters of five water molecules at the time-dependent DFT level in combination with different embedding schemes, ranging from point charges to polarizable point charges, with and without local field effects. When going from the bulk to the interface, the main observations of these calculations encompass i) a modest increase of the average polarizability but an increase by about a factor of two of its anisotropy, ii) an increase by about 20 % of the βHRS response, accompanied by a small increase of its depolarization ratio, and iii) a net increase of the component of the β tensor normal to the interface (βzzz ) as well as of β// . Globally, the interfacial effects on β are localized at the first molecular layer while they are observed up to the fourth molecular layer on α.

Unraveling the acid-base characterization and solvent effects on the structural and electronic properties of a bis-bidentate bridging ligand

Understanding the interactions and the solvent effects on the distribution of several species in equilibrium and how it can influence the ¹H-NMR properties, spectroscopy (UV-vis absorption), and the acid-base equilibria can be especially challenging. This is the case of a bis-bidentate bridging ligand bis(2-pyridyl)-benzo-bis(imidazole), where the two pyridyl and four imidazolyl nitrogen atoms can be protonated in different ways, depending on the solvent, generating many isomeric/tautomeric species. Herein, we report a combined theoretical-experimental approach based on a sequential quantum mechanics/molecular mechanics procedure that was successfully applied to describe in detail the acid-base characterization and its effects on the electronic properties of such a molecule in solution. The calculated free-energies allowed the identification of the main species present in solution as a function of the solvent polarity, and its effects on the magnetic shielding of protons (¹H-NMR chemical shifts), the UV-vis absorption spectra, and the acid-base equilibrium constants (pK_as) in aqueous solution. Three acid-base equilibrium constants were experimentally/theoretically determined (pK_a1 = 1.3/1.2, pK_a2 = 2.1/2.2 and pK_a5 = 10.1/11.3) involving mono-deprotonated and mono-protonated cis and trans species. Interestingly, other processes with pK_a3 = 3.7 and pK_a4 = 6.0 were also experimentally determined and assigned to the protonation and deprotonation of dimeric species. The dimerization of the most stable neutral species was investigated by Monte Carlo simulations and its electronic effects were considered for the elucidation of the UV-vis absorption bands, revealing transitions mainly with the charge-transfer characteristic and involving both the monomeric species and the dimeric species. The good matching of the theoretical and experimental results provides an atomistic insight into the solvent effects on the electronic properties of this bis-bidentate bridging ligand.

Extending NMR Quantum Computation Systems by Employing Compounds with Several Heavy Metals as Qubits

Nuclear magnetic resonance (NMR) is a spectroscopic method that can be applied to several areas. Currently, this technique is also being used as an experimental quantum simulator, where nuclear spins are employed as quantum bits or qubits. The present work is devoted to studying heavy metal complexes as possible candidates to act as qubit molecules. Nuclei such 113Cd, 199Hg, 125Te, and 77Se assembled with the most common employed nuclei in NMR-QIP implementations (1H, 13C, 19F, 29Si, and 31P) could potentially be used in heteronuclear systems for NMR-QIP implementations. Hence, aiming to contribute to the development of future scalable heteronuclear spin systems, we specially designed four complexes, based on the auspicious qubit systems proposed in our previous work, which will be explored by quantum chemical calculations of their NMR parameters and proposed as suitable qubit molecules. Chemical shifts and spin–spin coupling constants in four complexes were examined using the spin–orbit zeroth-order regular approximation (ZORA) at the density functional theory (DFT) level, as well as the relaxation parameters (T1 and T2). Examining the required spectral properties of NMR-QIP, all the designed complexes were found to be promising candidates for qubit molecules.

A QM/MM study of the conformation stability and electronic structure of the photochromic switches derivatives of DHA/VHF in acetonitrile solution

We present a detailed theoretical study of the electronic absorption spectra and thermochemistry of molecular photoswitches composed of one and two photochromic units of dihydroazulene (DHA)/vinylheptafulvene (VHF) molecules. Six different isomers are considered depending on the ring opening/closure forms of the DHA units. The solvent effect of acetonitrile is investigated using a sequential Molecular Mechanics/Quantum Mechanics approach. The thermochemical investigations of these photochromic molecules were performed using the Free Energy Perturbation method, and the simulations were performed using Configurational Bias Monte Carlo. We show that to open the 5-member ring of the DHA, there is no significant gain in thermal release of energy for the back reaction when a unit or two DHA units are considered. Overall, we found agreement between the solvation free energy based on Monte Carlo simulations and the continuum solvent model. However, the cavitation term in the continuum model is shown to be a source of disagreement when the non-electrostatic terms are compared. The electronic absorption spectra are calculated using TDDFT CAM-B3LYP/cc-pVDZ. Agreement with experiment is obtained within 0.1 eV, considering statistically uncorrelated configurations from the simulations. Inhomogeneous broadening is also considered and found to be well described in all cases.

Second Harmonic Generation Responses of Ion Pairs Forming Dimeric Aggregates

A sequential approach combining molecular dynamics and density functional theory calculations has been worked out to unravel the second harmonic generation responses of anion-cation (AC) pairs when they form dimeric aggregates, where the cation is a stilbazolium derivative and the anions range from small inorganic iodide to medium-size organic p-toluenesulfonate. These complexes showed a strong self-aggregation behavior in molecular dynamics simulations within high-concentration conditions and formed stable dimeric aggregates, (AC)₂, which can adopt different structural shapes from stacked, Λ, to head-to-head configurations. These various structures are associated with different symmetries, which are shown to modulate the second- and third-order nonlinear optical (NLO) responses. By consolidating the NLO results of this work with those previously obtained for single AC pairs [ J. Chem. Inf. Model. 2020, 60, 4817-4826], we have been able to explain the experimentally observed variations of the electrical-field-induced second harmonic generation (EFISHG) responses of these complexes as a function of concentration [ ChemPhysChem 2010, 11, 495-507]. Moreover, results have highlighted that (i) the second-order contribution, μβ_//, dominates the global EFISHG response; (ii) the μβ_// responses of dimers are about half of those computed for the parent AC pairs, while the third-order contributions, γ_//, are reduced by only 10%; (iii) these distinct trends are ascribed to the formation of dimers adopting mainly Λ and head-to-head shapes, increasing the centrosymmetric character, in comparison to the monomers, a situation in which the second-order response cancels out as well as influences the dipole moment on μβ_//; (iv) the presence of a strong amino donor group in the cation enhances the μβ_// response by 1 order of magnitude and γ_// by about a factor of 2; and finally, (v) dimeric aggregation has similar effects on the hyper-Rayleigh scattering response, β_HRS, as on μβ_//, while it reduces the one-dimensional character of β_HRS. This work constitutes a step forward for the modeling of the NLO responses of AC aggregates in solution.

A new interpretation of the absorption and the dual fluorescence of Prodan in solution

Remarkable interest is associated with the interpretation of the Prodan fluorescent spectrum. A sequential hybrid Quantum Mechanics/Molecular Mechanics method was used to establish that the fluorescent emission occurs from two different excited states, resulting in a broad asymmetric emission spectrum. The absorption spectra in several solvents were measured and calculated using different theoretical models presenting excellent agreement. All theoretical models [semiempirical, time dependent density functional theory and and second-order multiconfigurational perturbation theory] agree that the first observed band at the absorption spectrum in solution is composed of three electronic excitations very close in energy. Then, the electronic excitation around 340 nm-360 nm may populate the first three excited states (π-π^*L_b, n-π^*, and π-π^*L_a). The ground state S₀ and the first three excited states were analyzed using multi-configurational calculations. The corresponding equilibrium geometries are all planar in vacuum. Considering the solvent effects in the electronic structure of the solute and in the solvent relaxation around the solute, it was identified that these three excited states can change the relative order depending on the solvent polarity, and following the minimum path energy, internal conversions may occur. A consistent explanation of the experimental data is obtained with the conclusive interpretation that the two bands observed in the fluorescent spectrum of Prodan, in several solvents, are due to the emission from two independent states. Our results indicate that these are the n-π^* S₂ state with a small dipole moment at a lower emission energy and the π-π^*L_b S₁ state with large dipole moment at a higher emission energy.

X-ray Photoelectron Fingerprints of High-Valence Ruthenium-Oxo Complexes along the Oxidation Reaction Pathway in an Aqueous Environment

Recent advances in operando-synchrotron-based X-ray techniques are making it possible to address fundamental questions related to complex proton-coupled electron transfer reactions, for instance, the electrocatalytic water splitting process. However, it is still a grand challenge to assess the ability of the different techniques to characterize the relevant intermediates, with minimal interference on the reaction mechanism. To this end, we have developed a novel methodology employing X-ray photoelectron spectroscopy (XPS) in connection with the liquid-jet approach to probe the electrochemical properties of a model electrocatalyst, [Ru^II(bpy)₂(py)(OH₂)]²⁺, in an aqueous environment. There is a unique fingerprint of the extremely important higher-valence ruthenium-oxo species in the XPS spectra along the oxidation reaction pathway. Furthermore, a sequential method combining quantum mechanics and molecular mechanics is used to illuminate the underlying physical chemistry of such systems. This study provides the basis for the future development of in-operando XPS techniques for water oxidation reactions.

Investigating the first hyperpolarizability of liquid carbon tetrachloride

Sequential QMMM calculations have been carried out to investigate the first hyperpolarizability of liquid CCl₄.

Quenching of Charge Transfer in Nitrobenzene Induced by Vibrational Motion

Although nitrobenzene is the smallest nitro-aromatic molecule, the nature of its electronic structure is still unclear. Most notably, the lowest-energy absorption band was assessed in numerous studies providing conflicting results regarding its charge-transfer character. In this study, we employ a combination of molecular dynamics and quantum chemical methods to disentangle the nature of the lowest-energy absorption band of nitrobenzene. Surprisingly, the charge-transfer transition from the benzene moiety to the nitro group is found to be quenched by a flow of charge into the opposite direction induced by vibrational motion. Beyond clarifying the charge-transfer character of nitrobenzene, we show that the widely used approach of analyzing the ground-state minimum-energy geometry provides completely wrong conclusions.

The first hyperpolarizability of nitrobenzene in benzene solutions: investigation of the effects of electron correlation within the sequential QM/MM approach

Molecular polarization due to solvation amplifies the impact of electron correlation on the first hyperpolarizability of nitrobenzene.

Which charge definition for describing the crystal polarizing field and the χ(1) and χ(2) of organic crystals?

Crystal optical susceptibilities are probes to assess the performance of the charge definition employed to describe the crystal polarizing field.

Theoretical study of the solvent effect on the aromaticity of benzene: a NICS analysis

Nucleus-independent chemical shift (NICS) quantities for benzene-benzene and benzene-water species were obtained and are discussed in gas phase and in solution. Besides standard polarizable continuum model (PCM) calculations, sequential Monte Carlo/quantum mechanics (S-MC/QM) were also performed. Benzene was shown to be slightly more aromatic in condensate phase when we considered the average solvent configuration (ASEC) approach with explicit molecules.

Self-aggregation and optical absorption of stilbazolium merocyanine in chloroform

Dipolar aggregation is in many cases detrimental for the functioning of optical materials. In this study we investigate self-aggregation and optical absorption of stilbazolium merocyanine (SM) in chloroform solution by performing classical Molecular Dynamics (MD) simulations under ambient conditions. The reversal solvatochromic shift, the large bathochromic shift, and the structured absorption band presented by SM in chloroform solution are all aspects of its optical absorption behavior for which the existence of self-aggregation is yet not completely understood. Moreover, the spectroscopic properties of SM oligomers and their occurrence in solvent of low polarity remain a relevant topic that deserves to be investigated. Our analysis of the aggregation behavior of SM in chloroform verified that the majority of the chromophores are involved in the formation of oligomers in solution, where the whole dimer and part of the trimer populations present a stable π-stacking structure. The optical properties of the monomers and oligomers in solution were evaluated by means of a discrete polarizable embedding quantum mechanical/molecular mechanical (PE-QM/MM) response scheme where the quantum part is described at the level of density functional theory. The visible absorption spectrum of SM in chloroform is simulated using time average values obtained for the monomeric and oligomeric forms of SM from the PE-QM/MM calculations performed on uncorrelated configurations extracted from the classical MD simulations. This study shows that the self-aggregation of SM in chloroform may exist, but it is not essential for reproducing the reversal solvatochromic shift in chloroform and that the process does not contribute to enhance the bathochromic shift nor explain the structure observed in its absorption band. Moreover, it is verified that since the electronic transitions of the monomer and oligomers are close together, changes in the interplane separation between the monomeric units of the stacked oligomers substantially affect the spectral resolution of their contribution to the optical absorption spectrum.

One- and two-photon absorption of fluorescein dianion in water: a study using S-QM/MM methodology and ZINDO method

One- and two-photon absorption (1PA and 2PA) of fluorescein dianion (FSD) in water were studied using a combined and sequential Quantum Mechanics/Molecular Dynamics methodology. Different sets of 250 statistically relevant (uncorrelated) configurations composed by the solute and several solvent molecules were sampled from the classical simulation. On these configurations, the electronic properties were calculated a posteriori using the Zerner's intermediate neglect of differential overlap (ZINDO) method. The linear and nonlinear absorption of FSD in water were calculated using discrete and explicit solvent models. In the largest case, the relevant configurations are composed by FSD and 47 explicit water molecules embedded in the electrostatic field of all remaining water molecules. Both INDO/CIS and INDO/CISD calculations were performed to study the absorption processes of FSD and the Sum-Over-States (SOS) model was used to describe the 2PA process. A semi-classical method for spectrum simulations was employed to simulate the 1PA and 2PA cross-section spectra of FSD in water. For comparison purposes, in the case of the 2PA process two approaches, the "full expression" and "resonant expression" methods, were employed to simulate the nonlinear spectrum. The last method assumes resonant conditions and on the computation point of view it represents an interesting option to study the 2PA process. The INDO/CI calculations give a satisfactory description of the 1PA spectrum of FSD and properly describe the unusual blue-shift of its first π→π(*) transition in water. In the case of 2PA, the introduction of doubly excited configuration interactions (INDO/CISD) has proven to be essential for an appropriate description of the process at the higher energy spectral region. It was observed that the solvent effects do not drastically change the cross-sections of both processes. The simulated 2PA cross-section spectrum provided by the "full expression" method presents a better definition of the bands which appear along the experimental spectrum than the one provided by the "resonant expression" method. However, both approaches provide similar description for the effect of the solvent environment on the 2PA process of FSD in water.

Solvation of fluoroform and fluoroform-dimethylether dimer in liquid krypton: a theoretical cryospectroscopic study

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.