Search

This computational study of line-broadening effects is based on an accurate, analytical representation of the intermonomer potential energy surface (PES) of the CO 2  ⋅ Ar van der Waals (vdW) complex. The PES is employed to compute collisional broadening coefficients for rovibrational lines of CO 2 perturbed by Ar. The semiclassical computations are performed using the modified Robert-Bonamy approach, including real and imaginary terms, and the exact trajectory model. The lines investigated are in the 10001←00011, 01101←00001, 00011←00001, and 00031←00001 vibrational bands and the computations are repeated at multiple temperatures. The computed results are in good agreement with the available experimental values, validating both the intermonomer PES developed and the methodology used. For lines in the 01101←00001 band of CO 2 , temperature-dependent Ar-broadening coefficients are reported for the first time. The parameters presented should prove useful, among other applications, for the accurate experimental determination of CO 2 and Ar abundances in planetary atmospheres., (© 2023 The Authors. ChemPhysChem published by Wiley-VCH GmbH.)

Currently, the explicitly correlated coupled cluster method is used routinely to generate the multi-dimensional potential energy surfaces (mD-PESs) of van der Waals complexes of small molecular systems relevant for atmospheric, astrophysical and industrial applications. Although very accurate, this method is computationally prohibitive for medium and large molecules containing clusters. For instance, the recent detections of complex organic molecules (COMs) in the interstellar medium, such as benzonitrile, revealed the need to establish an accurate enough electronic structure approach to map the mD-PESs of these species interacting with the surrounding gases. As a benchmark, we have treated the case of the polar molecule benzonitrile interacting with helium, where we use post-Hartree-Fock and symmetry-adapted perturbation theory (SAPT) techniques. Accordingly, we show that MP2 and distinguishable-cluster approximation (DCSD) cannot be used for this purpose, whereas accurate enough PESs may be obtained using the corresponding explicitly correlated versions (MP2-F12 or DCSD-F12) with a reduction in computational costs. Alternatively, computations revealed that SAPT(DFT) is as performant as CCSD(T)-F12/aug-cc-pVTZ, making it the method of choice for mapping the mD-PESs of COMs containing clusters. Therefore, we have used this approach to generate the 3D-PES of the benzonitrile-He complex along the intermonomer Jacobi coordinates. As an application, we have incorporated the analytic form of this PES into quantum dynamical computations to determine the cross sections of the rotational (de-)excitation of benzonitrile colliding with helium at a collision energy of 10 cm -1 .

The four-dimensional-potential energy surface (4D-PES) of the CO2-CO van der Waals complex is generated using the explicitly correlated coupled cluster with single, double, and perturbative triple excitation (CCSD(T)-F12) method in conjunction with the augmented correlation-consistent triple zeta (aug-cc-pVTZ) basis set. This 4D-PES is developed over the set of inter-molecular coordinates and where the CO2 and CO monomers are treated as rigid rotors. Afterwards, analytic fits of this 4D-PES are carried out. In addition to the already known C-bound and O-bound stable structures of CO2-CO, we characterise a new isomer: it has a T-shaped structure where the O atom of the CO2 moiety points into the centre of mass of CO. We also find the saddle points connecting these minimal structures. This new isomer may play a role during the intramolecular isomerization processes at low energies. Then, the 4D-PES expansion is incorporated into bound vibrational state computations of C-bound and O-bound complexes. We also computed the temperature dependence of the second virial coefficient for CO2-CO. A good agreement with experiments is found.

The hydronium ion, H3O+, presents a crucial key to understanding the chemistry of the interstellar clouds where it has been frequently observed. The present paper is devoted to studying the inelastic scattering of both forms of the hydronium ion, o-H3O+ and p-H3O+, by helium atoms. The interaction potential between H3O+ and He was mapped in Jacobi coordinates leading to a new three dimensional potential energy surface (3D-PES) at the CCSD(T)/aug-cc-pVQZ+BF (CCSD(T)/AVQZ+BF) level of theory. Close coupling treatment was used to obtain rotational cross-sections for both o-H3O+ and p-H3O+ involving 9 rotational states in each case. The cross-sections were computed from energy thresholds up to 500 cm-1 for o-H3O+ and up to 300 cm-1 for p-H3O+, in order to generate rotational rate coefficients for the kinetic temperature range 5-50 K for both forms. The use of the present rates are viewed to be a good tool to estimate hydronium abundance.

In a Bernal graphene bilayer, carbon atoms belong to two inequivalent sublattices A and B, with atoms that are coupled to the other layer by [Formula: see text] bonds belonging to sublattice A and the other atoms belonging to sublattice B. We analyze the density of states and the conductivity of Bernal graphene bilayers when atoms of sublattice A or B only are randomly functionalized. We find that for a selective functionalization on sublattice B only, a mobility gap of the order of 0.5 eV is formed close to the Dirac energy at concentration of adatoms [Formula: see text]. In addition, at some other energies conductivity presents anomalous behaviors. We show that these properties are related to the bipartite structure of the graphene layer.

Dioxygen difluoride is a tough molecule that has defied accurate theoretical description for many decades. In the present work we have identified the reason for this resistance: the flatness of the OO, and more important OF, stretching potential energy curves, which make it difficult to localise the global minimum. It is not related to the weak multi-reference character. Using high-level CCSD(T)-F12/VTZ-F12 ab initio theory, the global minimum has been properly located and vibrationally averaged bond lengths obtained. These vibrationally averaged parameters agree with experimental data to within 0.01 Å. Averaging was found essential to achieve this unprecedented accuracy. We have then simulated the IR and UV spectra, which compare well with experimental data and permit identification of the observed transitions.

There has been a disagreement amongst experimentalists and between experimentalists and theoreticians as to the gas-phase structure of dimethyl peroxide. We have investigated this problem with high-level CCSD(T)-F12 and MRCI procedures. There can be no doubt anymore that, at the minimum of the potential energy surface, the COOC fragment has a trans-structure. The dynamical structure of the molecule can, however, be different and be explained by the very slow torsional motion. We have analysed the dynamical structure using numerical wavefunctions of the torsional motion and a fully optimized potential curve of MP2/aug-cc-pVTZ quality. Computational and all experimental results are shown to be in complete agreement. The problem that has persisted for more than thirty years, highlighted in a recent review article by Oberhammer titled "Gas phase structures of peroxides: experiments and computational problems", has been resolved.

Ab initio and Perdew, Burke, and Ernzerhof (PBE) density functional theory with dispersion correction (PBE-D3) calculations are performed to study N 2 -Ar n (n ≤ 3) complexes and N 2 trapped in Ar matrix (i.e., N 2 @Ar). For cluster computations, we used both Møller-Plesset (MP2) and PBE-D3 methods. For N 2 @Ar, we used a periodic-dispersion corrected model for Ar matrix, which consists on a slab of four layers of Ar atoms. We determined the equilibrium structures and binding energies of N 2 interacting with these entities. We also deduced the N 2 vibrational frequency shifts caused by clustering or embedding compared to an isolated N 2 molecule. Upon complexation or embedding, the vibrational frequency of N 2 is slightly shifted, while its equilibrium distance remains unchanged. This is due to the weak interactions between N 2 and Ar within these compounds. Our calculations show the importance of inclusion of dispersion effects for the accurate description of geometrical and spectroscopic parameters of N 2 isolated, in interaction with Ar surfaces, or trapped in Ar matrices.

We provide in this work, the absolute solvation enthalpies and the absolute solvation free energies of the proton in ammonia explicitly versus temperature. As a result, the absolute solvation free energy of the proton remains quite constant for temperatures below 200 K. Above this temperature, it increases as a linear function of the temperature: ΔG am (H + ,T)=-1265.832+0.210 T. This indicates that a temperature change of 100 K would induce a solvation free energy change of 21 kJ mol -1 . Thus, ignoring this free energy change would lead to a bad description of hydrogen bonds and an unacceptable error higher than 3.7 pK a units. However, the absolute solvation enthalpy of the proton in ammonia is not significantly affected by a temperature change and, the room temperature value is -1217 kJ mol -1 . The change of the solvation enthalpy is only within 3 kJ mol -1 for a temperature change up to 200 K.

Structures of protonated ammonia clusters (H + (NH 3 ) n ) are very important for the determination of pKa's and solvation energies of the proton in ammonia. In this work, their structures were investigated at M06-2X/6-31++g(d,p) level of theory, for n=10-16 and for temperatures ranging from 0 to 400 K. In the cluster community, this is the first theoretical study on the protonated ammonia clusters larger than the nonamer. We noted that the population of the investigated clusters is reproduced by branched cage or cage like structures at low temperatures, while branched linear and branched cyclic or branched double cyclic isomers are the only isomers responsible for the population at higher temperatures. In these isomers, the proton is highly and entirely solvated at the center of the cluster. In addition, protonated ammonia clusters are all Eigen structures and the first solvation shell of the related ammonium ion core is saturated by four ammonia molecules. Moreover, infrared (IR) spectra of all isomers have been investigated and these spectra show good agreement with the experiment. This allowed us to assign experimental peaks and to provide the constitution of the populations of the various clusters.

Ab initio studies on CH x (OH) 4-x (x = 0-3) polyols are carried out to derive their structural and spectroscopic features. Several stereoisomers (both equilibrium structures and transition states) are found. Some are predicted here for the first time. We determined hence their geometrical parameters, vibrational frequencies, electronic excitation energies for the singlet manifold, and IR spectra. While the IR spectra for all polyols present similar shapes, their UV spectra exhibit however distinct band origin that are specific to each polyol and more interestingly to each diasteroisomer. Stereoelectronic effects are also noticed and discussed. It is suggested that UV spectroscopy is an efficient probe to experimentally identify polyols in mixtures involving polyols.

The accurate determination of the solvation energies of a proton in ammonia is based on the precise knowledge of the structures of neutral and protonated ammonia clusters. In this work, we have investigated all the possible and stable structures of protonated ammonia clusters H + (NH 3 ) n=2-9 , along with their isomeric distribution at a specific temperature. New significant isomers are reported here for the first time and show that the structures of protonated ammonia clusters are not only branched linear as assumed by all previous authors. Branched linear structures are the only ones responsible for the population of protonated ammonia clusters for n = 4-6 at any temperature. However, for larger cluster sizes, these types of structures compete with branched cyclic, double cyclic, branched double cyclic and triple cyclic structures depending on the temperature. In addition, we have shown that protonated ammonia clusters are all Eigen structures and the first solvation shell of the related ammonium ion core is saturated by four ammonia molecules. We have also carried out a study of the hydrogen bond network of protonated ammonia clusters establishing the stability rule governing the various isomers of each cluster from estimated energies of the hydrogen bond types in H + (NH 3 ) n=2-9 . With all these results, a route for the accurate determination of the solvation energies of a proton in ammonia at a given temperature could be conceivable.

Benchmark ab initio calculations are performed to investigate the stable isomers of [O,N,S,S]. These computations are carried out using coupled cluster (RCCSD(T)) and explicitly correlated coupled cluster methods (RCCSD(T)-F12). In addition to the already known cis isomer of SSNO, nine other stable forms are predicted. The most stable isomer is cis-OSNS. Nine structures are chain bent-bent with relatively large dipole moments which make them detectable, as cis-SSNO, by infrared, far-infrared, and microwave spectroscopies. We found also a C2v isomer (NS2O). Since these species are strongly suggested to play an important role as intermediates during the bioactive reaction products of the NO/H2S interaction, the rotational and vibrational spectroscopic parameters are presented to help aid the in vivo identification and assignment of these spectra. Results from this work show that [O,N,S,S] may play key roles during nitric oxide transport and deliver in biological media, as well as, provide an explanation for the weak characteristic of disulfide bridges within proteins.

Small methane clusters (CH4)n are irradiated using intense femtosecond laser excitation at 624 nm. The ionized species and those resulting from their fragmentation are detected via time-of-flight mass spectrometry (TOF MS). We find evidence of bound, multiply charged methane molecules and clusters resulting from Coulomb explosion upon exposure to highly energetic, ultrafast radiation. The assignment of the mass spectra is done after first-principles calculations (at the (R)MP2/aug-cc-pVXZ (X = D,T) level) on the charged (CH4)n(q+) clusters (n = 1-4, q = 1-4). We also considered the cluster stabilities and fragments that may result from intracluster molecular reactivity. Complex intracluster ion-molecule reactions induced by photoionization are expected to occur. Interestingly, we show that multi charged small methane clusters undergo intracluster reactions and fragmentations which are different from those observed for isolated methane ions or for large ionized methane clusters.

A ro-vibrationally inelastic close coupling study of the rotational excitation of OH(+)(X(3)Σ(-)) by collisions with H((2)S) is presented. The two lowest potential energy surfaces of doublet and quadruplet spin multiplicity are involved. The former is the one we developed recently, and the latter is a modified version of the quadruplet surface of Martinez et al. to include the long-range charge-induced-dipole potential. The details of the modification of this surface are presented as well as the comparison of the rotational excitation resulting from collisions with hydrogen on these two surfaces. The effect of the coupling between vibration and rotation on the rotational excitation rate is also discussed, as the potential well depth of the doublet surface is quite large and allows the coupling between many vibrational channels of OH(+). As the hydrogen exchange reaction can occur for both potential energy surfaces, we discuss the reliability of the approximation made by the calculation of the cross sections with a quantum dynamics limited to the inelastic process. The relative importance of the collisions on the doublet or quadruplet surface within a given rotational transition is also discussed.

A hydrogen bond network in ammonia clusters plays a key role in understanding the properties of species embedded in ammonia. This network is dictated by the structures of neutral ammonia clusters. In this work, structures of neutral ammonia clusters (NH3)n(=2-10) have been studied at M06-2X/6-31++G(d,p) and MP2/6-31++g(d,p) levels of theory. The analysis of the relative stabilities of various hydrogen bond types has also been studied and vibrational spectroscopy of the ammonia pentamer and decamer is investigated. We noted that M06-2X provides lower electronic energies, greater binding energies and higher structural resolution than MP2. We also noted that at the M06-2X level of theory, the binding energy converges to the experimental vaporization enthalpy faster than that at the MP2 level of theory. As a result, it is found that the M06-2X functional could be more suitable than the MP2 ab initio method in the description of structures and energies of ammonia clusters. However, we found that the electronic energy differences obtained at both levels of computation follow a linear relation with n (number of ammonia molecules in a cluster). As far as the structures of ammonia clusters are concerned, we proposed new "significant" isomers that have not been reported previously. The most remarkable is the global minimum electronic energy structure of the ammonia hexamer, which has an inversion centre and confirms experimental observation. Moreover, we reported the relative stabilities of neutral ammonia clusters for temperatures ranging from 25 to 400 K. The stability of isomers changes with the increase of the temperature. As a result, the branched and less bonded isomers are the most favored at high temperatures and disfavored at low temperatures, while compact and symmetric isomers dominate the population of clusters at low temperatures. In fine, from this work, the global minimum energy structures of ammonia clusters are known for the first time at a given temperature (T ∼ 0-400 K) and at a reliable computational level of theory.

Results are presented that suggest that thiazyl hydride (HSN)/thionitrosyl hydride (sulfimide, HNS) can be used as light-sensitive compounds for NO-delivery in biological media, as well as markers for the possible detection of intermediates in nitrites + H2S reactions at the cellular level. They are expected to be more efficient than the HNO/HON isovalent species and hence they should be considered instead. A set of characteristic spectroscopic features are identified that could aid in the possible detection of these species in the gas phase or in biological environments. The possibility of intramolecular dynamical processes involving excited states that are capable of interconverting HNS and its isomeric form HSN is examined.

We investigate the lowest electronic states of doublet and quartet spin multiplicity states of HNS(-) and HSN(-) together with their parent neutral triatomic molecules. Computations were performed using highly accurate ab initio methods with a large basis set. One-dimensional cuts of the full-dimensional potential energy surfaces (PESs) along the interatomic distances and bending angle are presented for each isomer. Results show that the ground anionic states are stable with respect to the electron detachment process and that the long range parts of the PESs correlating to the SH(-) + N, SN(-) + H, SN + H(-), NH + S(-), and NH(-) + S are bound. In addition, we predict the existence of long-lived weakly bound anionic complexes that can be formed after cold collisions between SN(-) and H or SH(-) and N. The implications for the reactivity of these species are discussed; specifically, it is shown that the reactions involving SH(-), SN(-), and NH(-) lead either to the formation of HNS(-) or HSN(-) in their electronic ground states or to autodetachment processes. Thus, providing an explanation for why the anions, SH(-), SN(-), and NH(-), have limiting detectability in astrophysical media despite the observation of their corresponding neutral species. In a biological context, we suggest that HSN(-) and HNS(-) should be incorporated into H2S-assisted heme-catalyzed reduction mechanism of nitrites in vivo.

Through benchmark studies, we explore the performance of PBE density functional theory, with and without Grimme's dispersion correction (DFT-D3), in predicting spectroscopic properties for molecules interacting with rare gas matrices. Here, a periodic-dispersion corrected model of matrix embedding is used for the first time. We use PBE-D3 to determine the equilibrium structures and harmonic vibrational frequencies of carbon monoxide in interaction with small Ar clusters (CO-Arn, n = 1, 2, 3), with an Ar surface and embedded in an Ar matrix. Our results show a converging trend for both the vibrational frequencies and binding energies when going from the gas-phase to a fully periodic approach describing CO embedding in Ar. This trend is explained in terms of solvation effects, as CO is expected to alter the structure of the Ar matrix. Due to a competition between CO-Ar interactions and Ar-Ar interactions, perturbations caused by the presence of CO are found to extend over several Å in the matrix. Accordingly, it is mandatory to fully relax rare gas matrices when studying their interaction with embedded molecules. Moreover, we show that the binding energy per Ar is almost constant (∼-130 cm(-1) atom(-1)) regardless of the environment of the CO molecule. Finally, we show that the concentration of the solute into the cold matrix influences the spectroscopic parameters of molecules embedded into cold matrices. We suggest hence that several cautions should be taken before comparing these parameters to gas phase measurements and to theoretical data of isolated species.

Four-dimensional potential energy surface (4D-PES) of the atmospherically relevant CO2-N2 van der Waals complex is generated using the explicitly correlated coupled cluster with single, double, and perturbative triple excitation (CCSD(T)-F12) method in conjunction with the augmented correlation consistent triple zeta (aug-cc-pVTZ) basis set. This 4D-PES is mapped along the intermonomer coordinates. An analytic fit of this 4D-PES is performed. Our extensive computations confirm that the most stable form corresponds to a T-shape structure where the nitrogen molecule points towards the carbon atom of CO2. In addition, we located a second isomer and two transition states in the ground state PES of CO2-N2. All of them lay below the CO2 + N2 dissociation limit. This 4D-PES is flat and strongly anisotropic along the intermonomer coordinates. This results in the possibility of the occurrence of large amplitude motions within the complex, such as the inversion of N2, as suggested in the recent spectroscopic experiments. Finally, we show that the experimentally established deviations from the C2v structure at equilibrium for the most stable isomer are due to the zero-point out-of-plane vibration correction.

Using the recently developed explicitly correlated coupled cluster method in connection with the aug-cc-pVTZ basis set, we generated the three-dimensional potential energy surface (3D-PES) of the ground state of the Ar-BeO complex. This PES covers the regions of the global and local minima, the saddle point, and the dissociation of the complex. The PES is also used for the calculation of the rovibrational spectrum up to the dissociation limit. The high density of levels which is observed favors the mixing of the states and hence the occurrence of anharmonic resonances. The wavefunctions of the high rovibrational levels exhibit large amplitude motions in addition to strong anharmonic resonances. Our theoretical spectrum should be helpful in identifying the van der Waals modes of this complex in laboratory.

We performed accurate ab initio investigations of the geometric parameters and the vibrational structure of neutral HNS/HSN triatomics and their singly charged anions and cations. We used standard and explicitly correlated coupled cluster approaches in connection with large basis sets. At the highest levels of description, we show that results nicely approach those obtained at the complete basis set limit. Moreover, we generated the three-dimensional potential energy surfaces (3D PESs) for these molecular entities at the coupled cluster level with singles and doubles and a perturbative treatment of triple excitations, along with a basis set of augmented quintuple-zeta quality (aug-cc-pV5Z). A full set of spectroscopic constants are deduced from these potentials by applying perturbation theory. In addition, these 3D PESs are incorporated into variational treatment of the nuclear motions. The pattern of the lowest vibrational levels and corresponding wavefunctions, up to around 4000 cm(-1) above the corresponding potential energy minimum, is presented for the first time.

We present an extensive study of the four-dimensional potential energy surface (4D-PES) of the carbon dioxide dimer, (CO2)2. This PES is developed over the set of intermolecular coordinates. The electronic computations are carried out at the explicitly correlated coupled cluster method with single, double, and perturbative triple excitations [CCSD(T)-F12] level of theory in connection with the augmented correlation-consistent aug-cc-pVTZ basis set. An analytic representation of the 4D-PES is derived. Our extensive calculations confirm that "Slipped Parallel" is the most stable form and that the T-shaped structure corresponds to a transition state. Later on, this PES is employed for the calculations of the vibrational energy levels of the dimer. Moreover, the temperature dependence of the dimer second virial coefficient and of the first spectral moment of rototranslational collision-induced absorption spectrum is derived. For both quantities, a good agreement is found between our values and the experimental data for a wide range of temperatures. This attests to the high quality of our PES. Generally, our PES and results can be used for modeling CO2 supercritical fluidity and examination of its role in planetary atmospheres. It can be also incorporated into dynamical computations of CO2 capture and sequestration. This allows deep understanding, at the microscopic level, of these processes.

We mapped the global three-dimensional potential energy surface (3D-PES) of the water cation at the MRCI/aug-cc-pV5Z including the basis set superposition (BSSE) correction. This PES covers the molecular region and the long ranges close to the H + OH(+)(X(3)Σ(-)), the O + H2(+)(X(2)Σg(+)), and the hydrogen exchange channels. The quality of the PES is checked after comparison to previous experimental and theoretical results of the spectroscopic constants of H2O(+)(X(2)B1) and of the diatomic fragments, the vibronic spectrum, the dissociation energy, and the barrier to linearity for H2O(+)(X(2)B1). Our data nicely approach those measured and computed previously. The long range parts reproduce quite well the diatomic potentials. In whole, a good agreement is found, which validates our 3D-PES.

A theoretical study of rotational collision of LiH(X(1)Σ(+),v = 0, J) with Ar has been carried out. The ab initio potential energy surface (PES) describing the interaction between the Ar atom and the rotating LiH molecule has been calculated very accurately and already discussed in our previous work [Computational and Theoretical Chemistry 993 (2012) 20-25]. This PES is employed to evaluate the de-excitation cross sections. The ab initio PES for the LiH(X(1)Σ(+))-Ar((1)S) Van der waals system is calculated at the coupled-cluster [CCSD(T)] approximation for a LiH length fixed to an experimental value of 3.0139 bohrs. The basis set superposition error (BSSE) is corrected and the bond functions are placed at mid-distance between the center of mass of LiH and the Ar atom. The cross sections are then derived in the close coupling (CC) approach and rate coefficients are inferred by averaging these cross sections over a Maxwell-Boltzmann distribution of kinetic energies. The 11 first rotational levels of rate coefficients are evaluated for temperatures ranging from 10 to 300 K. We notice that the de-excitation rate coefficients appear large in the order 10(-10) cm(-3) s(-1) and show very low temperature dependence. The rate coefficients magnify significantly the propensity toward ∆ J = -1 transitions. These results confirm the same propensity already noted for the cross sections.

We closely compare the accuracy of multidimensional potential energy surfaces (PESs) generated by the recently developed explicitly correlated coupled cluster (CCSD(T)-F12) methods in connection with the cc-pVXZ-F12 (X = D, T) and aug-cc-pVTZ basis sets and those deduced using the well-established orbital-based coupled cluster techniques employing correlation consistent atomic basis sets (aug-cc-pVXZ, X = T, Q, 5) and extrapolated to the complete basis set (CBS) limit. This work is performed on the benchmark rare gas-hydrogen halide interaction (HCl-He) system. These PESs are then incorporated into quantum close-coupling scattering dynamical calculations in order to check the impact of the accuracy of the PES on the scattering calculations. For this system, we deduced inelastic collisional data including (de-)excitation collisional and pressure broadening cross sections. Our work shows that the CCSD(T)-F12/aug-cc-pVTZ PES describes correctly the repulsive wall, the van der Waals minimum and long range internuclear distances whereas cc-pVXZ-F12 (X = D,T) basis sets are not diffuse enough for that purposes. Interestingly, the collision cross sections deduced from the CCSD(T)-F12/aug-cc-pVTZ PES are in excellent agreement with those obtained with CCSD(T)/CBS methodology. The position of the resonances and the general shape of these cross sections almost coincide. Since the cost of the electronic structure computations is reduced by several orders of magnitude when using CCSD(T)-F12/aug-cc-pVTZ compared to CCSD(T)/CBS methodology, this approach can be recommended as an alternative for generation of PESs of molecular clusters and for the interpretation of accurate scattering experiments as well as for a wide production of collisional data to be included in astrophysical and atmospherical models.

The accurate evaluation of pKa's, or solvation energies of the proton in methanol at a given temperature is subject to the determination of the most favored structures of various isomers of protonated (H(+)(MeOH)n) and neutral ((MeOH)n) methanol clusters in the gas phase and in methanol at that temperature. Solvation energies of the proton in a given medium, at a given temperature may help in the determination of proton affinities and proton dissociation energies related to the deprotonation process in that medium and at that temperature. pKa's are related to numerous properties of drugs. In this work, we were interested in the determination of the most favored structures of various isomers of protonated methanol clusters in the gas phase and in methanol, at a given temperature. For this aim, the M062X/6-31++G(d,p) and B3LYP/6-31++G(d,p) levels of theory were used to perform geometries optimizations and frequency calculations on various isomers of (H(+)(MeOH)n) in both phases. Thermal effects were retrieved using our homemade FORTRAN code. Thus, we accessed the relative populations of various isomers of protonated methanol clusters, in both phases for temperatures ranging from 0 to 400 K. As results, in the gas phase, linear structures are entropically more favorable at high temperatures, while more compact ones are energetically more favorable at lower temperatures. The trend is somewhat different when bulk effects are taken into account. At high temperatures, the linear structure only dominates the population for n ≤ 6, while it is dominated by the cyclic structure for larger cluster sizes. At lower temperatures, compact structures still dominate the population, but with an order different from the one established in the gas phase. Hence, temperature effects dominate solvent effects in small cluster sizes (n ≤ 6), while the reverse trend is noted for larger cluster sizes.

pKa's, proton affinities, and proton dissociation free energies characterize numerous properties of drugs and the antioxidant activity of some chemical compounds. Even with a higher computational level of theory, the uncertainty in the proton solvation free energy limits the accuracy of these parameters. We investigated the thermochemistry of the solvation of the proton in methanol within the cluster-continuum model. The scheme used involves up to nine explicit methanol molecules, using the IEF-PCM and the strategy based on thermodynamic cycles. All computations were performed at B3LYP/6-31++G(dp) and M062X/6-31++G(dp) levels of theory. It comes out from our calculations that the functional M062X is better than B3LYP, on the evaluation of gas phase clustering energies of protonated methanol clusters, per methanol stabilization of neutral methanol clusters and solvation energies of the proton in methanol. The solvation free energy and enthalpy of the proton in methanol were obtained after converging the partial solvation free energy of the proton in methanol and the clustering free energy of protonated methanol clusters, as the cluster size increases. Finally, the recommended values for the solvation free energy and enthalpy of the proton in methanol are -257 and -252 kcal/mol, respectively.