Computer-aided simulation of infrared spectra of ethanol conformations in gas, liquid and in CCl4 solution

Structure and intermolecular interactions in ionic liquid 1-ethyl-3-methylimidazolium bromide and its aqueous solutions investigated by vibrational spectroscopy and quantum chemical computations

The recently developed efficient protocol to explicit quantum mechanical modeling of structure and IR spectra of liquids and solutions (S. A. Katsyuba, S. Spicher, T. P. Gerasimova, S. Grimme, J. Phys. Chem. B 2020, 124, 6664) is applied to ionic liquid (IL) 1-ethyl-3-methylimidazolium bromide (EmimBr), its C2-deuterated analog [Emim-d]Br and its aqueous solutions. It is shown that the solvation strongly modifies frequencies and IR intensities of the CH/CD stretching vibrations (νCH/νCD) of the imidazolium ring. The main vibrational spectroscopic features of the neat IL are reproduced by the simulations for a cluster (EmimBr)9, in which all three imidazolium CH moieties of the solvated cation form short contacts with three Br- anions, and another two Br- anions are located on top and bottom of imidazolium ring. Cluster models of aqueous solutions reproduce the experimental vibrational frequencies of actual solutions, provided that the Br- anion of solvated contact ion pair (CIP) is situated on top of imidazolium ring, and CH/CD moieties of the latter participate in short contacts with surrounding water molecules. Both structural and spectroscopic analysis allow to interpret the short contacts CH/CD⋯Br- and CH/CD⋯OH2 as hydrogen bonds of approximately equal strength. Enthalpies of bonding of these liquid-state H-bonds, estimated with the use of empirical correlations, amount to ca. 1.4 kcal⋅mol-1, while the analogous estimates obtained for the gas-phase charged species [Emim]2Br+ increase to 5.6 kcal⋅mol-1. It is shown that formation of solvent-shared ion pair (SIP) in aqueous solution, where the counterions of IL are separated by two water molecules H-bonded to a Br- anion, produces frequency shifts ΔνCH/CD, strongly different from the case of CIP formation. This difference can be used for IR/Raman spectroscopic differentiation of the type of solvated ion pairs of EmimBr or other related ILs.

Structural flexibility of favipiravir and its structural analogues in solutions: experimental and computational insight

Combined UV-vis and quantum chemical studies of the structural flexibility and tautomerism of 6-R-3-hydroxy-2-pyrazine carboxamides in solutions revealed that their keto-enol transformations are accompanied by the deprotonation of enol tautomers and the formation of the corresponding anionic species. Both the solvent and the 6-R substituent strongly influence the relative abundance of the above forms in solutions. Anions are not formed in 1,2-dichloroethane (DCE), but the probability of deprotonation in neutral water and N,N-dimethylformamide (DMF) increases in the order R = H < F < NO2. Only enol tautomers of all solutes are found in DCE. DMF stabilizes keto forms only moderately and assists much strongly in the deprotonation of all three compounds. Water tends to stabilize both keto tautomers and deprotonated anions: the keto form dominates in the case of R = H (antiviral drug T-1105), the anions are found exclusively for R = NO2, and the aqueous solution of another antiviral drug, favipiravir (R = F), contains both the keto tautomer and the anionic form. The results of quantum chemical free energy calculations are in agreement with the experimental observations.

Computationally assisted vibrational spectroscopy of nucleic acid bases. 2. Thymine

As in the case of cytosine [Phys. Chem. Chem. Phys. 2023, 25, 24121-24128], Raman and infrared (IR) spectra of aqueous thymine and its N-deuterated derivative, thymine-d2 have been computationally reproduced and interpreted with the use of the recently developed efficient protocol to explicit quantum mechanical modeling of structure and IR spectra of liquids and solutions [J. Phys. Chem. B, 2020, 124, 6664-6670]. A cluster model of a solute surrounded by 30 water molecules is shown to be sufficient to reproduce experimental vibrational frequencies and relative Raman intensities with the use of B3LYP-D3/def2-TZVP or B3LYP-D3/aug-cc-pVDZ simulations. Analogous PBE-D3 computations provided a less good, but still reasonably accurate, modeling of Raman spectra. It is shown that strong changes of frequencies and relative intensities of the Raman bands of thymine, caused by its hydration, can be interpreted mainly as a result of hydrogen bonding with 6 nearest water molecules. Non-negligible improvement of the quality of simulations for larger clusters comprising water molecules that do not have direct contacts with the solute, suggests that spectroscopic effects of hydration should be ascribed to the joined action of solute-solvent and solvent-solvent interactions. Nevertheless, the moderate number of water molecules required for successful simulations of the Raman spectra of aqueous thymine, suggests that the vibrational modes and derivatives of the polarizability of the solute are mainly locally influenced, while the effect of bulk water is rather modest.

Structure of Tris[2-(4-pyridyl)ethyl]phosphine, Tris[2-(2-pyridyl)ethyl]phosphine, and Their Chalcogenides in Solution: Dipole Moments, IR Spectroscopy, and DFT Study

Tris(hetaryl)substituted phosphines and their chalcogenides are promising polydentate ligands for the design of metal complexes. An experimental and theoretical conformational analysis of tris[2-(4-pyridyl)ethyl]phosphine, tris[2-(2-pyridyl)ethyl]phosphine, and their chalcogenides was carried out by the methods of dipole moments, IR spectroscopy and DFT B3PW91/6-311++G(df,p) calculations. In solution, these compounds exist as an equilibrium of mainly non-eclipsed (synclinal or antiperiplanar) forms with a predominance of a symmetrical conformer having a gauche-orientation of the Csp3-Csp3 bonds of pyridylethyl substituents relative to the P=X bond (X = lone pair, O, S, Se) and a gauche-orientation of the pyridyl rings relative to the zigzag ethylene bridges. Regardless of the presence and nature of the chalcogen atom (oxygen, sulfur, or selenium) in the studied molecules with many axes of internal rotation, steric factors-the different position of the nitrogen atoms in the pyridyl rings and the configuration of ethylene bridges-determine the realization and spatial structure of preferred conformers.

Efficient electrocatalytic CO2 reduction to ethanol through the proton coupled electron transfer process of PVnMo(12-n) (n = 1, 2, 3) over indium electrode

The multistep proton-coupled electron transfer (PCET) processes are beneficial for products distribution and selectivity of the electrocatalytic CO2 reduction reaction (CO2RR), which are affected by the nature of the catalyst and electrolyte at electrode-electrolyte interface. Polyoxometalates (POMs) are electron regulators of PCET processes, which can catalyze CO2RR effectively. Accordingly, the commercial indium electrodes are combined in this work with a series of Keggin-type POMs (PVnMo(12-n)O40)(n+3)-, n = 1, 2, 3) to process CO2RR with Faradaic efficiency toward ethanol reaching 93.4% at -0.3 V (vs. RHE). The cyclic voltammetry and X-ray photoelectron spectroscopy results reveal the activation of CO2 molecules by the first PCET process of the VⅤ/Ⅳ in POM. Subsequently, the PCET process of MoⅥ/Ⅴ results the oxidation of the electrode, causing the loss of In0 active sites. Electrochemical in-situ infrared spectroscopy confirms the weak adsorption of *CO at the later stage of electrolysis due to the oxidation of the In0 active sites. The indium electrode in PV3Mo9 system retains more In0 active sites owing to the highest V-substitution ratio, thereby ensuring a high adsorption ratio of *CO and CC coupling. In sum, the regulation of the interface microenvironment by POM electrolyte additives can be used to boost the performance of CO2RR.

Computational analysis of the vibrational spectra and structure of aqueous cytosine

The recently developed efficient protocol for the explicit quantum mechanical modeling of the structure and IR spectra of liquids and solutions [Katsyuba et al., J. Phys. Chem. B, 2020, 124, 6664-6670] is used to describe aqueous solutions of cytosine. The same cluster model of a solute surrounded by the first solvation shell of solvent molecules was shown to be sufficient to reproduce experimental vibrational frequencies and relative IR and Raman intensities. An equally good quality of Raman spectra was provided by B3LYP-D3/def2-TZVP and B3LYP-D3/aug-cc-pVDZ simulations. Computations using the PBE functional were sufficient for modeling of the IR spectra but failed in the simulations of Raman scattering. It is shown that strong changes of frequencies and relative intensities of Raman and IR bands of cytosine, caused by its hydration, cannot be completely assigned to the influence of hydrogen bonds (HBs) with 7 or 8 closest water molecules. They are rather ascribed to the combined effect of solute-solute and solute-solvent HBs with the participation of at least 30 water molecules separating cytosine from the bulk solvent. This suggests that the vibrational modes and derivatives of the polarizability and dipole moment of the solute are mainly locally influenced by its first hydration shell, while the influence of bulk water is rather modest.

Semiclassical and VSCF/VCI Calculations of the Vibrational Energies of trans- and gauche-Ethanol Using a CCSD(T) Potential Energy Surface

A recent full-dimensional Δ-Machine learning potential energy surface (PES) for ethanol is employed in semiclassical and vibrational self-consistent field (VSCF) and virtual-state configuration interaction (VCI) calculations, using MULTIMODE, to determine the anharmonic vibrational frequencies of vibration for both the trans and gauche conformers of ethanol. Both semiclassical and VSCF/VCI energies agree well with the experimental data. We find significant mixing between the VSCF basis states due to Fermi resonances between bending and stretching modes. The same effects are also accurately described by the full-dimensional semiclassical calculations. These are the first high-level anharmonic calculations using a PES, in particular a "gold-standard" CCSD(T) one.

Quantum Calculations on a New CCSD(T) Machine-Learned Potential Energy Surface Reveal the Leaky Nature of Gas-Phase Trans and Gauche Ethanol Conformers

Ethanol is a molecule of fundamental interest in combustion, astrochemistry, and condensed phase as a solvent. It is characterized by two methyl rotors and trans (anti) and gauche conformers, which are known to be very close in energy. Here we show that based on rigorous quantum calculations of the vibrational zero-point state, using a new ab initio potential energy surface (PES), the ground state resembles the trans conformer, but substantial delocalization to the gauche conformer is present. This explains experimental issues about identification and isolation of the two conformers. This “leak” effect is partially quenched when deuterating the OH group, which further demonstrates the need for a quantum mechanical approach. Diffusion Monte Carlo and full-dimensional semiclassical dynamics calculations are employed. The new PES is obtained by means of a Δ-machine learning approach starting from a pre-existing low level density functional theory surface. This surface is brought to the CCSD(T) level of theory using a relatively small number of ab initio CCSD(T) energies. Agreement between the corrected PES and direct ab initio results for standard tests is excellent. One- and two-dimensional discrete variable representation calculations focusing on the trans–gauche torsional motion are also reported, in reasonable agreement with experiment.

Automated Molecular Cluster Growing for Explicit Solvation by Efficient Force Field and Tight Binding Methods

An automated and broadly applicable workflow for the description of solvation effects in an explicit manner is introduced. This method, termed quantum cluster growth (QCG), is based on the semiempirical GFN2-xTB/GFN-FF methods, enabling efficient geometry optimizations and MD simulations. Fast structure generation is provided using the intermolecular force field xTB-IFF. Additionally, the approach uses an efficient implicit solvation model for the electrostatic embedding of the growing clusters. The novel QCG procedure presents a robust cluster generation tool for subsequent application of higher-level (e.g., DFT) methods to study solvation effects on molecular geometries explicitly or to average spectroscopic properties over cluster ensembles. Furthermore, the computation of the solvation free energy with a supermolecular approach can be carried out with QCG. The underlying growing process is physically motivated by computing the leading-order solute-solvent interactions first and can account for conformational and chemical changes due to solvation for low-energy barrier processes. The conformational space is explored with the NCI-MTD algorithm as implemented in the CREST program, using a combination of metadynamics and MD simulations. QCG with GFN2-xTB yields realistic solution geometries and reasonable solvation free energies for various systems without introducing many empirical parameters. Computed IR spectra of some solutes with QCG show a better match to the experimental data compared to well-established implicit solvation models.

What quantum chemical simulations tell us about the infrared spectra, structure and interionic interactions of a bulk ionic liquid

The recently developed efficient protocol to explicit quantum mechanical modeling of the structure and IR spectra of liquids and solutions [Katsyuba et al., J. Phys. Chem. B, 2020, 124, 6664-6670] is applied to ionic liquid 1-ethyl-3-methyl-imidazolium tetrafluoroborate [Emim][BF₄], and its C2-deuterated analog [Emim-d][BF₄]. It is shown that the solvation strongly modifies the frequencies and IR intensities of both cationic and anionic components of the ionic liquids. The main features of the bulk spectra are reproduced by the simulations for cluster ([Emim][BF₄])₈, representing an ion pair solvated by the first solvation shell. The geometry of the cluster closely resembles the solid-state structure of the actual ionic liquid and is characterized by short contacts of all CH moieties of the imidazolium ring with [BF₄]^- anions. Both structural and spectroscopic analyses allow the contacts to be interpreted as hydrogen bonds of approximately equal strength. The enthalpies of these liquid-state H-bonds, estimated with the use of empirical correlations, amount to 1.2-1.5 kcal mol^-1, while the analogous estimates obtained for the gas-phase charged species [Emim][BF₄]₂^- and [Emim]₂[BF₄]⁺ increase to 3.6-3.9 kcal mol^-1.

Synthesis and Computationally Assisted Spectroscopic Study of Tautomerism in 3-(Phenyl(2-arylhydrazineylidene)methyl)quinoxalin-2(1H)-ones

The recently developed efficient protocol combining implicit and explicit, accurate quantum mechanical modeling of the condensed state [Katsyuba et al., J. Chem. Phys. 155, 024507 (2021)] is used to describe...