Computational study of the electrostatic potential and charges of multivalent ionic liquid molecules

Effect of structure and interaction on physicochemical properties of new [Emim][BF3X] complex anion ionic liquids studied by quantum chemistry

CONTEXT

One of the key challenges in the industrial application of ionic liquids (ILs) is their extreme characteristics, such as viscosity, glass transition temperatures, and conductivity. Understanding the relationship between ILs structure and physicochemical property is a crucial aspect of the directed design of ILs with good properties, which is a prerequisite for their successful implementation in industrial processes. In this work, high-level quantum-chemical research with for four pairs ionic liquids, [Emim][X] and [Emim][BF3X] (X = CH3SO3, EtSO4, HSO4, Tos), was performed, to analyze the stable structure, interionic interaction, and charge transfer and provide a new insight into the property variances at the molecular level. The result shows that the overall structural stability of ionic liquids is contributed with hydrogen bonding network between the protons in the C-H and N-H of the cation and oxygen atoms of the anion, as well as fluorine atoms. The nature and strength of the interionic interaction were measured via atoms in molecule analysis and sobEDAw method and results suggested that BF3 could waning interionic interaction of ion pairs. Moreover, a close relation between the binding energies of ion pairs and physicochemical properties was established: the weaker the interionic interaction, the lower is the viscosity and glass transition, and the higher is the conductivity.

METHODS

Quantum chemistry calculations were performed under B3LYP-D3/aug-cc-pVTZ level of DFT functional using the Gaussian 16 package (version C01). The Multiwfn 3.7 program was used to calculate the electrostatic potential, interaction region indicator, the information of bond critical points, core-valence bifurcation index, and ADCH charge. Visualization of structure and the region of interaction were achieved using VESTA and VMD.

Molecular Interactions of Perfluorinated and Branched Fluorine-Free Surfactants at Interfaces: Insights from a New Reliable Force Field

Per- and polyfluoroalkyl substances (PFAS) constitute a class of synthetic compounds with exceptional interfacial properties. Their widespread use in many industrial applications and consumer products, combined with their remarkable chemical and thermal stability, has led to their ubiquitous presence in environmental matrices, including surface water and groundwater. To replace PFAS with fluorine-free surfactants, it is necessary first to develop a deep molecular-level understanding of the mechanisms responsible for the exceptional properties of PFAS. For instance, it has been shown that fluorine-free surfactants with highly branched or methylated chains can achieve low surface tensions at air-water interfaces and can provide highly hydrophobic surface coatings. Although molecular simulations combined with experiments are promising for uncovering these mechanisms, the reliability of simulation results depends strongly on the accuracy of the force fields implemented. At the moment, atomistic force fields are not available to describe PFAS in a variety of environments. Ab initio methods could help fill this knowledge gap, but they are computationally demanding. As an alternative, ab initio calculations could be used to develop accurate force fields for atomistic simulations. In this work, a new algorithm is proposed, which, built from accurate ab initio calculations, yields force fields for perfluorinated sulfonic and perfluoroalkyl acids. The accuracy of the new force field was benchmarked against solvation free energy and interfacial tension data. The new force fields were then used to probe the interfacial behavior of the PFAS surfactants. The interfacial properties observed in our simulations were compared with those manifested by two branched fluorine-free surfactants. The good agreement achieved with experiments and ab initio calculations suggests that the proposed protocol could be implemented to study other perfluorinated substances and help in the design of fluorine-free surfactants for targeted applications.

Pro-inflammatory enzyme inhibition of lipoxygenases by flavonoid rich extract from Artemisia vulgaris

The peroxyl radicals generated by the activity of lipoxygenases (LOX) are mediators to trigger inflammatory diseases. Therefore, it is important to investigate potent LOX inhibitor for modulating the occurrence and resolving inflammatory processes. Artemisa vulgaris, is a herbal plant that is known for flavonoids, potentially inhibiting lipid peroxidation and scavenging radicals. The objectives of the present study were to obtain flavonoids rich extract from A. vulgaris, and determine the inhibitory mode of the extract against LOX. The flavonoids rich extract was optimized in an ultrasound assisted extraction using ionic liquids as extraction solvent. The results found that the optimum conditions; ratio of solid-to-liquid (1:10) and 30 min of extraction time could produce the high yield (10.14 %) and flavonoid content (5.30 mg QE/g). The LOX activity was demonstrated to follow a mixed mode of inhibition in the presence of the flavonoid rich extract as an inhibitor. The Michaelis-Menten constant (Km) was increased from 0.283 µM to 0.435 µM, whereas the maximum velocity was reduced from 0.22 µM/min to 0.058 µM/min in the inhibition. The flavonoids rich extract is likely to be a natural potent non-competitive inhibitor which may bind to free LOX or substrate-bound LOX.

Insights into the Effect of Charges on Hydrogen Bonds

Previous computational and experimental studies showed that charges located at the surroundings of hydrogen bonds can exert two opposite effects on them: rupture or strengthening of the hydrogen bond. This work aims to generalize the effect of charges in different hydrogen-bonded systems and to propose a coherent explanation of this effect. For these purposes, 19 systems with intra- and intermolecular hydrogen bonds were studied computationally with DFT. The FT-IR spectra of the systems were simulated, and two energy components of the hydrogen bond were studied separately to determine their variation upon the presence of a charge: charge transfer and molecular overlap. It was determined that either the breaking or strengthening of the hydrogen bond can be favored one over the other, for instance, depending on the heteroatom involved in the hydrogen bond. In addition, it is showed that the strengthening of the hydrogen bond by the presence of a charge is directly related to the decrease in charge transfer between the monomers, which is explained by an increase in molecular overlapping, suggesting a more covalent character of the interaction. The understanding of how hydrogen bonds are affected by charges is important, as it is a key towards a strategy to manipulate hydrogen bonds at convenience.

Molecular Modelling of Ionic Liquids: Situations When Charge Scaling Seems Insufficient

Charge scaling as an effective solution to the experiment-computation disagreement in molecular modelling of ionic liquids (ILs) could bring the computational results close to the experimental reference for various thermodynamic properties. According to the large-scale benchmark calculations of mass density, solvation, and water-ILs transfer-free energies in our series of papers, the charge-scaling factor of 0.8 serves as a near-optimal option generally applicable to most ILs, although a system-dependent parameter adjustment could be attempted for further improved performance. However, there are situations in which such a charge-scaling treatment would fail. Namely, charge scaling cannot really affect the simulation outcome, or minimally perturbs the results that are still far from the experimental value. In such situations, the vdW radius as an additional adjustable parameter is commonly tuned to minimize the experiment-calculation deviation. In the current work, considering two ILs from the quinuclidinium family, we investigate the impacts of this vdW-scaling treatment on the mass density and the solvation/partition thermodynamics in a fashion similar to our previous charge-scaling works, i.e., scanning the vdW-scaling factor and computing physical properties under these parameter sets. It is observed that the mass density exhibits a linear response to the vdW-scaling factor with slopes close to -1.8 g/mL. By further investigating a set of physiochemically relevant temperatures between 288 K and 348 K, we confirm the robustness of the vdW-scaling treatment in the estimation of bulk properties. The best vdW-scaling parameter for mass density would worsen the computation of solvation/partition thermodynamics, and a marginal decrease in the vdW-scaling factor is considered as an intermediate option balancing the reproductions of bulk properties and solvation thermodynamics. These observations could be understood in a way similar to the charge-scaling situation. i.e., overfitting some properties (e.g., mass density) would degrade the accuracy of the other properties (e.g., solvation free energies). Following this principle, the general guideline for applying this vdW-tuning protocol is by using values between the density-derived choice and the solvation/partition-derived solution. The charge and current vdW scaling treatments cover commonly encountered ILs, completing the protocol for accurate modelling of ILs with fixed-charge force fields.