Near-infrared and multifunctional fluorescent probe enabled by cyanopyridine cyanine dye for bisulfite recognition and biological imaging

SO2 derivatives, sulfite/bisulfite, are widely employed in both the food processing and drug synthesis industries. Despite their widespread application, excessive levels of sulfite/bisulfite can negatively impact human health. Most probes for detecting sulfite/bisulfite are restricted by their fluorescence within the visible spectrum range and poor solubility in aqueous solution, which limit their use in food testing and biological imaging. Herein, a near-infrared probe comprising of the cyanopyridine cyanine skeleton, 4-((Z)-2-((E)-2-chloro-3-(2-cyano-2-(1-methylpyridine-4(1H)-ylidene)ethylidene)cyclohex-1-en-1-yl)-1-cyanovinyl)-1-methylpyridin-1-ium (abbreviated as CCP), was developed. This probe enables precise quantification of bisulfite (HSO3-) in almost pure buffered solutions, showing a near-infrared fluorescence emission at 784 nm with an impressively low detection limit of 0.32 μM. The probe stands out for its exceptional selectivity, minimal susceptibility to interference, and strong adaptability. The probe CCP utilizes the CC bond to trigger a near-infrared fluorescence quenching reaction with HSO3- via nucleophilic addition, which effectively disrupts the large delocalization within the molecule for accurate HSO3- identification. Moreover, the probe has been successfully applied in detecting HSO3- in various food products and living cells, simplifying the measurement of HSO3- content in water samples. This advancement not only enhances the analytical capabilities but also contributes to ensuring food safety and environmental protection. ENVIRONMENTAL IMPLICATION: SO2 derivatives including sulfite/bisulfite, serving dual roles as preservatives and antioxidants, have widespread application across various sectors including food preservation, water sanitation, and the pharmaceutical industry. Despite their widespread application, excessive levels of sulfite/bisulfite can affect human health. Developing methods for precisely and sensitively detecting sulfite/bisulfite in food products and biological samples is important for ensuring food safety and environmental protection. Here, a sensitive near-infrared and multifunctional fluorescent probe in a 99.9 % buffered solution, along with water gel encapsulation, has been successfully applied for the detection of bisulfite in food, authentic water samples, and biological cells.

Active species in chloroaluminate ionic liquids catalyzing low-temperature polyolefin deconstruction

Chloroaluminate ionic liquids selectively transform (waste) polyolefins into gasoline-range alkanes through tandem cracking-alkylation at temperatures below 100 °C. Further improvement of this process necessitates a deep understanding of the nature of the catalytically active species and the correlated performance in the catalyzing critical reactions for the tandem polyolefin deconstruction with isoalkanes at low temperatures. Here, we address this requirement by determining the nuclearity of the chloroaluminate ions and their interactions with reaction intermediates, combining in situ 27Al magic-angle spinning nuclear magnetic resonance spectroscopy, in situ Raman spectroscopy, Al K-edge X-ray absorption near edge structure spectroscopy, and catalytic activity measurement. Cracking and alkylation are facilitated by carbenium ions initiated by AlCl3-tert-butyl chloride (TBC) adducts, which are formed by the dissociation of Al2Cl7- in the presence of TBC. The carbenium ions activate the alkane polymer strands and advance the alkylation cycle through multiple hydride transfer reactions. In situ 1H NMR and operando infrared spectroscopy demonstrate that the cracking and alkylation processes occur synchronously; alkenes formed during cracking are rapidly incorporated into the carbenium ion-mediated alkylation cycle. The conclusions are further supported by ab initio molecular dynamics simulations coupled with an enhanced sampling method, and model experiments using n-hexadecane as a feed.

Chloride and Hydride Transfer as Keys to Catalytic Upcycling of Polyethylene into Liquid Alkanes

Transforming polyolefin waste into liquid alkanes through tandem cracking-alkylation reactions catalyzed by Lewis-acid chlorides offers an efficient route for single-step plastic upcycling. Lewis acids in dichloromethane establish a polar environment that stabilizes carbenium ion intermediates and catalyzes hydride transfer, enabling breaking of polyethylene C-C bonds and forming C-C bonds in alkylation. Here, we show that efficient and selective deconstruction of low-density polyethylene (LDPE) to liquid alkanes is achieved with anhydrous aluminum chloride (AlCl3) and gallium chloride (GaCl3). Already at 60 °C, complete LDPE conversion was achieved, while maintaining the selectivity for gasoline-range liquid alkanes over 70 %. AlCl3 showed an exceptional conversion rate of 5000g L D P E m o l c a t - 1 h - 1 ${{{\rm g}}_{{\rm L}{\rm D}{\rm P}{\rm E}}{{\rm \ }{\rm m}{\rm o}{\rm l}}_{{\rm c}{\rm a}{\rm t}}^{-1}{{\rm \ }{\rm h}}^{-1}}$ , surpassing other Lewis acid catalysts by two orders of magnitude. Through kinetic and mechanistic studies, we show that the rates of LDPE conversion do not correlate directly with the intrinsic strength of the Lewis acids or steric constraints that may limit the polymer to access the Lewis acid sites. Instead, the rates for the tandem processes of cracking and alkylation are primarily governed by the rates of initiation of carbenium ions and the subsequent intermolecular hydride transfer. Both jointly control the relative rates of cracking and alkylation, thereby determining the overall conversion and selectivity.

Magnetic Ionic Liquids: Current Achievements and Future Perspectives with a Focus on Computational Approaches

Magnetic ionic liquids (MILs) stand out as a remarkable subclass of ionic liquids (ILs), combining the desirable features of traditional ILs with the unique ability to respond to external magnetic fields. The incorporation of paramagnetic species into their structures endows them with additional attractive features, including thermochromic behavior and luminescence. These exceptional properties position MILs as highly promising materials for diverse applications, such as gas capture, DNA extractions, and sensing technologies. The present Review synthesizes key experimental findings, offering insights into the structural, thermal, magnetic, and optical properties across various MIL families. Special emphasis is placed on unraveling the influence of different paramagnetic species on MILs' behavior and functionality. Additionally, the Review highlights recent advancements in computational approaches applied to MIL research. By leveraging molecular dynamics (MD) simulations and density functional theory (DFT) calculations, these computational techniques have provided invaluable insights into the underlying mechanisms governing MILs' behavior, facilitating accurate property predictions. In conclusion, this Review provides a comprehensive overview of the current state of research on MILs, showcasing their special properties and potential applications while highlighting the indispensable role of computational methods in unraveling the complexities of these intriguing materials. The Review concludes with a forward-looking perspective on the future directions of research in the field of magnetic ionic liquids.

Combined Spectroscopic, Thermodynamic, and Theoretical Approach for Detecting and Quantifying Hydrogen Bonding and Dispersion Interaction in Ionic Liquids

ConspectusIonic liquids (ILs) are attracting increasing interest in science and engineering due to their unique properties that can be tailored for specific applications. Clearly, a better understanding of their behavior on the microscopic scale will help to elucidate macroscopic fluid phenomena and thereby promote potential applications. The advantageous properties of these innovative fluids arise from the delicate balance of Coulomb interactions, hydrogen bonding, and dispersion forces. The development of these properties requires a fundamental understanding of the strength, location, and direction of the different types of interactions and their contribution to the overall phase behavior. Contrary to expectations, hydrogen bonding and dispersion interactions have a significant influence on the structure, dynamics, and phase behavior of ILs.The synergy between experimental and theoretical methods has now advanced to a stage where hydrogen bonds and dispersion effects as well as the competition between the two can be studied in detail. In this account, we demonstrate that a suitable combination of spectroscopic, thermodynamic, and theoretical methods enables the detection, dissection, and quantification of noncovalent interactions, even in complex systems such as ionic liquids. This approach encompasses far-infrared vibrational spectroscopy (FIR), various thermodynamic methods for determining enthalpies of vaporization, and quantum chemical techniques that allow us to switch dispersion contributions on or off when calculating the energies and spectroscopic properties of clusters.We briefly discuss these experimental and theoretical methods, before providing various examples illustrating how the mélange of Coulomb interaction, hydrogen bonds, and dispersion forces can be analyzed, and their individual contributions quantified. First, we demonstrated that both hydrogen bonding and dispersion interactions are manifested in the FIR spectra and can be quantified by observed shifts of characteristic spectral signatures. Through the selection of suitable protic ionic liquids (PILs) featuring anions with varying interaction strengths and alkyl chain lengths, we were able to demonstrate that dispersion interactions can compete with hydrogen bonding. The resultant transition enthalpy serves as a measure of the dispersion interaction. Contrary to expectations, PILs possess lower enthalpies of vaporization compared with aprotic ILs (AILs). The reason for this is simple: In protic ILs, ion pairs carry both the hydrogen bond and attractive dispersion between the cation and anion into the gas phase. By utilizing a well-curated set of protic ILs and molecular analogues, we successfully disentangled Coulomb interaction, hydrogen bonding, and dispersion interaction through purely thermodynamic methods.

Effect of isopropyl side chain branching and different anions on electronic structure, vibrational spectra, and hydrogen bonding of isopropyl-imidazolium-based ionic liquids: Experimental and theoretical investigations

In the present work, four branched methylated, 1,2-dimethyl-3-isopropyl-imidazolium (i-[C₃Dmim⁺]) and protonated,1-methyl-3-isopropyl-imidazolium (i-[C₃mim⁺])-based ionic liquids (ILs) with varying anion (Br^-, BF₄^-, PF₆^-, and NTf₂^-) were synthesized and investigated by NMR, infrared (IR) and Raman spectroscopy. Based on infrared and Raman spectroscopy, complete vibrational assignments have been performed. The IR and Raman analysis revealed that the vibrational spectra are virtually unaffected upon methylation, while significant frequency changes were observed by changing anion. Furthermore, to determine the electronic structure, energetic stability, and vibrational properties of these i-[C₃Dmim]Y, i-[C₃mim]Y (Y = Br, BF₄, PF₆, and NTf₂) ion pairs, quantum chemical calculations including the dispersion correction method are performed both on single ions and on ionic couples. The calculated electron density was analyzed to expose non-covalent intra- and interionic interactions by the quantum theory of atoms in molecules (AIM) and interpreted in terms of both anion dependence and type of interaction. Computational results suggest that for all ionic couples the most energetically stable configuration is obtained with the anions located close to the C2 position of the imidazolium cation. However, in the case of i-[C₃mim]NTf₂ and i-[C₃Dmim]BF₄, similar energies were obtained in configurations 2 and 3 where the anion is located above the imidazolium ring. For i-[C₃mim]Br a stronger hydrogen bond is predicted than for other studied ILs. Calculations indicate that a red shift of the CH stretching bands should occur due to hydrogen bonding; indeed, such displacement of bands is experimentally observed.

Microheterogeneity-Induced Vibrational Spectral Dynamics of Aqueous 1-Alkyl-3-methylimidazolium Tetrafluoroborate Ionic Liquids of Different Cationic Chain Lengths

We have monitored the impacts of an increment in the alkyl chain length of the imidazolium-based tetrafluoroborate ionic liquids on the local deuteroxyl probe modes of interest. For this study, we have taken 1-ethyl-3-methylimidazolium tetrafluoroborate [EMIm][BF₄], 1-butyl-3-methylimidazolium tetrafluoroborate [BMIm][BF₄], 1-octyl-3-methylimidazolium tetrafluoroborate [OMIm][BF₄], and 1-decyl-3-methylimidazolium tetrafluoroborate [DMIm][BF₄] ionic liquid solutions with 5% HOD in H₂O as the vibrational reporter of the associated ultrafast system dynamics. Classical molecular dynamics (MD) simulations were employed to determine molecular structure and dynamic properties, while the spectral profiles were derived by applying the wavelet analysis of classical trajectories. Spatial distribution functions reveal the heterogeneity within the molecular structures of the ionic liquids (ILs) with varying alkyl chain lengths. The intense position of the spectral peak, the frequency corresponding to the shoulder peak, and the spectral linewidth of the O-D stretch distribution are not influenced by the increment in the cationic chain length. In addition, the ionic liquid (IL) [BMIm][BF₄] exhibits a notable trend; the dynamic timescales are longer than the other studied systems. Therefore, we have performed the Voronoi decomposition analysis of the ionic and the polar-apolar domains, symmetrically increasing the length of alkyl chains on the IL cations. Domain analysis reveals structural microheterogeneity; the anions form discrete domains, and the ionic liquid constituting cations form continuous domains irrespective of the alkyl chain length on the imidazolium cations. Therefore, this computational ultrafast spectroscopy study aids in forming a molecular-level picture of the ionic liquid cations and anions in the liquid phase, providing a detailed interpretation of the spectral properties of the probe stretching vibrations.

Comparison between Ab Initio Molecular Dynamics and OPLS-Based Force Fields for Ionic Liquid Solvent Organization

OPLS-based force fields (FFs) have been shown to provide accurate bulk-phase properties for a wide variety of imidazolium-based ionic liquids (ILs). However, the ability of OPLS to reproduce an IL solvent structure is not as well-validated given the relative lack of high-level theoretical or experimental data available for comparison. In this study, ab initio molecular dynamics (AIMD) simulations were performed for three widely used ILs: the 1-butyl-3-methylimidazolium cation with chloride, tetrafluoroborate, or hexafluorophosphate anions, that is, [BMIM][Cl], [BMIM][BF₄], and [BMIM][PF₆], respectively, as a basis for further assessment of two unique IL FFs: the ±0.8 charge-scaled OPLS-2009IL FF and the OPLS-VSIL FF. The OPLS-2009IL FF employs a traditional all-atom functional form, whereas the OPLS-VSIL FF was developed using a virtual site that offloads negative charge to inside the plane of the ring with careful attention given to reproducing hydrogen bonding. Detailed comparisons between AIMD and the OPLS FFs were made based on radial distribution functions (RDFs), combined distribution functions (CDFs), and spatial distribution functions (SDFs) to examine cation-anion interactions and π⁺-π⁺ stacking between the imidazolium rings. While both FFs were able to correctly capture the general solvent structure of these popular ILs, the OPLS-VSIL FF quantitatively reproduced interaction distances more accurately. In addition, this work provides further insights into the different short- and long-range structure patterns of these popular ILs.

Thermally activated delayed fluorescence in luminescent cationic copper(i) complexes

In this work, the photophysical characteristics of [Cu(N^N)2]+ and [Cu(N^N)(P^P)]+ complexes were described. The concept of thermally activated delayed fluorescence (TADF) and its development throughout the years was also explained. The importance of ΔE (S1-T1) and spin-orbital coupling (SOC) values on the TADF behavior of [Cu(N^N)2]+ and [Cu(N^N)(P^P)]+ complexes is discussed. Examples of ΔE (S1-T1) values reported in the literature were collected and some trends were proposed (e.g. the effect of the substituents at the 2,9 positions of the phenanthroline ligand). Besides, the techniques (or calculation methods) used for determining ΔE (S1-T1) values were described. The effect of SOC in TADF was also discussed, and examples of the determination of SOC values by DFT and TD-DFT calculations are provided. The last chapter covers the applications of [Cu(N^N)2]+ and [Cu(N^N)(P^P)]+ TADF complexes and the challenges that are still needed to be addressed to ensure the industrial applications of these compounds.

Structures and properties of ionic crystals and condensed phase ionic liquids predicted with the generalized energy-based fragmentation method

The generalized energy-based fragmentation (GEBF) approach is extended to facilitate ab initio investigations of structures, lattice energies, vibrational spectra and ¹ H NMR chemical shifts of ionic crystals and condensed-phase ionic liquids (ILs) with the periodic boundary conditions (PBC). For selected periodic systems, our results demonstrate that the so-called PBC-GEBF approach can provide satisfactory descriptions on ground-state energies, structures, and vibrational spectra of ionic crystals and IL crystals. The PBC-GEBF approach is then applied to three realistic condensed phase systems. For three ionic crystals (LiCl, NaCl, and KCl), we apply the PBC-GEBF approach with MP2 theory as well as some popular DFT methods to investigate their crystal structures and lattice energies. Our calculations indicate that the crystal structures obtained with PBC-GEBF-MP2/6-311 + G** are very close to the corresponding X-ray structures, while PBC-GEBF-ωB97X-D/6-311 + G** provides satisfactory prediction for crystal structures and lattice energies. For two polymorphs of [n-C₄ mim][Cl] crystals, we find that the PBC-GEBF approach at the M06-2X/6-311 + G** level can give a satisfactory descriptions on structures and Raman spectra of these two crystals. Furthermore, for [C₂ mim][BF₄ ] ILs, we demonstrate that their ¹ H NMR chemical shifts can be estimated from averaging over 5 typical snapshots (extracted from MD simulations) with the PBC-GEBF approach at the B97-2/pcSseg-2 level. The calculated results account for the observed experimental data quite well. Therefore, we expect that the PBC-GEBF approach, combined with various quantum chemistry methods, will become an effective tool in predicting structures and properties of ionic crystals and condensed-phase ILs.

Inclusion of More Physics Leads to Less Data: Learning the Interaction Energy as a Function of Electron Deformation Density with Limited Training Data

Machine learning (ML) approaches to predicting quantum mechanical (QM) properties have made great strides toward achieving the computational chemist's holy grail of structure-based property prediction. In contrast to direct ML methods, which encode a molecule with only structural information, in this work, we show that QM descriptors improve ML predictions of dimer interaction energy, both in terms of accuracy and data efficiency, by incorporating electronic information into the descriptor. We present the electron deformation density interaction energy machine learning (EDDIE-ML) model, which predicts the interaction energy as a function of Hartree-Fock electron deformation density. We compare its performance with leading direct ML schemes and modern DFT methods for the prediction of interaction energies for dimers of varying charge type, size, and intermolecular separation. Under a low-data regime, EDDIE-ML outperforms other direct ML schemes and is the only model readily transferrable to larger, more complex systems including base pair trimers and porous cages. The underlying physical connection between the density and interaction energy enables EDDIE-ML to reach an accuracy comparable to modern DFT functionals in fewer training data points compared to other ML methods.

Carbon dioxide adsorption and activation on ionic liquid decorated Au(111) surface: A DFT study

We use first principle approaches to study the adsorption and catalytic activation mechanism of CO₂ on ionic liquids (ILs, [C_nMIm]⁺[Cl]^- (n = 0-6)) attached to a Au(111) surface. The adsorption of CO₂ at this liquid-solid model interface occurs via either (i) parallel π-stacking mode or (ii) CO₂ oxygen lone pair (lp)···π interaction. These CO₂ physisorption modes, which depend on the CO₂ landing angle at this interface, are identified as an efficient way to activate CO₂ and its further conversion into value-added products. For illustration, we discuss the conversion of CO₂ into formic acid where the ILs@Au(111) decorated interface allows reduction of the activation energy for the CO₂ + H₂ → HCOOH reaction. In sum, our electrode/electrolyte based interface model provides valuable information to design novel heterogeneous catalysts for CO₂ conversion. Indeed, our work establishes that a suitable interface material is enough to activate CO₂.

Spectroscopic Signatures of Hydrogen-Bonding Motifs in Protonic Ionic Liquid Systems: Insights from Diethylammonium Nitrate in the Solid State

Diethylammonium nitrate, [N0 0 2 2][NO3], and its perdeuterated analogue, [N D D 2 2] [NO3], were structurally characterized and studied by infrared, Raman, and inelastic neutron scattering (INS) spectroscopy. Using these experimental data along with state-of-the-art computational materials modeling, we report unambiguous spectroscopic signatures of hydrogen-bonding interactions between the two counterions. An exhaustive assignment of the spectral features observed with each technique has been provided, and a number of distinct modes related to NH···O dynamics have been identified. We put a particular emphasis on a detailed interpretation of the high-resolution, broadband INS experiments. In particular, the INS data highlight the importance of conformational degrees of freedom within the alkyl chains, a ubiquitous feature of ionic liquid (IL) systems. These findings also enable an in-depth physicochemical understanding of protonic IL systems, a first and necessary step to the tailoring of hydrogen-bonding networks in this important class of materials.

Unravelling the Interactions of Magnetic Ionic Liquids by Energy Decomposition Schemes: Towards a Transferable Polarizable Force Field

This work aims at unravelling the interactions in magnetic ionic liquids (MILs) by applying Symmetry-Adapted Perturbation Theory (SAPT) calculations, as well as based on those to set-up a polarisable force field model for these liquids. The targeted MILs comprise two different cations, namely: 1-butyl-3-methylimidazolium ([Bmim]⁺) and 1-ethyl-3-methylimidazolium ([Emim]⁺), along with several metal halides anions such as [FeCl₄]⁻, [FeBr₄]⁻, [ZnCl₃]⁻ and [SnCl₄]²⁻ To begin with, DFT geometry optimisations of such MILs were performed, which in turn revealed that the metallic anions prefer to stay close to the region of the carbon atom between the nitrogen atoms in the imidazolium fragment. Then, a SAPT study was carried out to find the optimal separation of the monomers and the different contributions for their interaction energy. It was found that the main contribution to the interaction energy is the electrostatic interaction component, followed by the dispersion one in most of the cases. The SAPT results were compared with those obtained by employing the local energy decomposition scheme based on the DLPNO-CCSD(T) method, the latter showing slightly lower values for the interaction energy as well as an increase of the distance between the minima centres of mass. Finally, the calculated SAPT interaction energies were found to correlate well with the melting points experimentally measured for these MILs.

Bi(III) halometallate ionic liquids: Interactions and speciation

Bismuth containing compounds are of particular interest for optical or photo-luminescent applications in sensing, bio-imaging, telecommunications, and opto-electronics and as components in non-toxic extremely dense liquids. Bismuth(III) halometallates form highly colored novel ionic liquid based solvents for which experimental characterization and fundamental understanding are limited. In this work, Bismuth(III) halometallates incorporating chloride, bromide, and iodide have been studied via density functional theory employing B3LYP-D3BJ/aug-cc-pVDZ. Lone anions, and anions in clusters with sufficient 1-ethyl-3-methyl-imidazolium [C₂C₁Im]⁺ counter-cations to balance the charge, have been investigated in the gas- phase, and with polarizable continuum solvation. Evaluation of speciation profiles indicates that dimeric or trimeric anions are prevalent. In contrast to analogous Al systems, anions of higher charge (-2, -3) are present. Speciation profiles are similar, but not identical with respect to the halide. The Bi based anions [Bi_mX_n]^x- in the gas phase and generalized solvation environment produce multiple low energy conformers; moreover, key structural interaction patterns emerge from an analysis of ion-pair and neutral-cluster structures (Bi_mX_n)^x-(C₂C₁Im)_x ⁺ for x = 1, 2, and 3. Cation-anion interactions are weak; with Coulombic and dispersion forces predominating, anion-π structures are favored, while significant hydrogen bonding does not occur. Anion to cation charge transfer is minimal, but mutual polarization is significant, leading to local positive regions in the anion electrostatic potential surface. The key features of experimental x-ray photoelectron, UV-Vis spectra, and Raman spectra are reproduced, validating the computational results and facilitating rationalization of key features.

Adsorption of the drug bempedoic acid over different 2D/3D nanosurfaces and enhancement of Raman activity enabling ultrasensitive detection: First principle analysis

The nanocluster-based drug delivery system is of much importance, now days. This manuscript studies the interaction of pristine/substituted/doped GQDs, fullerene, helicene and CNT with bempedoic acid, which is an effective alternative of statins in the treatment of hypercholesteremia. The adsorption energies are calculated at B3LYP-D3/6-311G+(2d,p) level in order to study the adsorption of bempedoic acid over the surfaces of the nanoclusters incorporating Grimme's dispersion correction. Surface enhanced Raman scattering (SERS), which is a sound approach to vibrational spectroscopy, is used in order to detect bempedoic acid. All the studies signify that bempedoic acid can be detected with these nanoclusters and the negative adsorption energies advocate for the possible use of these nanoclusters as effective drug delivery system in case of bempedoic acid. Adsorption energy of bempedoic acid over helicene was found to be the most negative among the mentioned nanocluster systems, while adsorption on the surface of CNT was found to be the least negative.

Benchmarking London dispersion corrected density functional theory for noncovalent ion-π interactions

The strongly attractive noncovalent interactions of charged atoms or molecules with π-systems are important binding motifs in many chemical and biological systems. These so-called ion-π interactions play a major role in enzymes, molecular recognition, and for the structure of proteins. In this work, a molecular test set termed IONPI19 is compiled for inter- and intramolecular ion-π interactions, which is well balanced between anionic and cationic systems. The IONPI19 set includes interaction energies of significantly larger molecules (up to 133 atoms) than in other ion-π test sets and covers a broad range of binding motifs. Accurate (local) coupled cluster values are provided as reference. Overall, 19 density functional approximations, including seven (meta-)GGAs, eight hybrid functionals, and four double-hybrid functionals combined with three different London dispersion corrections, are benchmarked for interaction energies. DFT results are further compared to wave function based methods such as MP2 and dispersion corrected Hartree-Fock. Also, the performance of semiempirical QM methods such as the GFNn-xTB and PMx family of methods is tested. It is shown that dispersion-uncorrected DFT underestimates ion-π interactions significantly, even though electrostatic interactions dominate the overall binding. Accordingly, the new charge dependent D4 dispersion model is found to be consistently better than the standard D3 correction. Furthermore, the functional performance trend along Jacob's ladder is generally obeyed and the reduction of the self-interaction error leads to an improvement of (double) hybrid functionals over (meta-)GGAs, even though the effect of the SIE is smaller than expected. Overall, the double-hybrids PWPB95-D4/QZ and revDSD-PBEP86-D4/QZ turned out to be the most reliable among all assessed methods for the description of ion-π interactions, which opens up new perspectives for systems where coupled cluster calculations are no longer computationally feasible.

Extension of an Atom-Atom Dispersion Function to Halogen Bonds and Its Use for Rational Design of Drugs and Biocatalysts

A dispersion function Das in the form of a damped atom-atom asymptotic expansion fitted to ab initio dispersion energies from symmetry-adapted perturbation theory was improved and extended to systems containing heavier halogen atoms. To illustrate its performance, the revised Das function was implemented in the multipole first-order electrostatic and second-order dispersion (MED) scoring model. The extension has allowed applications to a much larger set of biocomplexes than it was possible with the original Das. A reasonable correlation between MED and experimentally determined inhibitory activities was achieved in a number of test cases, including structures featuring nonphysically shortened intermonomer distances, which constitute a particular challenge for binding strength predictions. Since the MED model is also computationally efficient, it can be used for reliable and rapid assessment of the ligand affinity or multidimensional scanning of amino acid side-chain conformations in the process of rational design of novel drugs or biocatalysts.

Ab initio and force field molecular dynamics study of bulk organophosphorus and organochlorine liquid structures

We performed ab initio molecular dynamics (AIMD) simulations to benchmark bulk liquid structures and to evaluate results from all-atom force field molecular dynamics (FFMD) simulations with the generalized Amber force field (GAFF) for organophosphorus (OP) and organochlorine (OC) compounds. Our work also addresses the current and important topic of force field validation, applied here to a set of nonaqueous organic liquids. Our approach differs from standard treatments, which validate force fields based on thermodynamic data. Utilizing radial distribution functions (RDFs), our results show that GAFF reproduces the AIMD-predicted asymmetric liquid structures moderately well for OP compounds that contain bulky alkyl groups. Among the OCs, RDFs obtained from FFMD overlap well with AIMD results, with some offsets in position and peak structuring. However, re-parameterization of GAFF for some OCs is needed to reproduce fully the liquid structures predicted by AIMD. The offsets between AIMD and FFMD peak positions suggest inconsistencies in the developed force fields, but, in general, GAFF is able to capture short-ranged and long-ranged interactions of OPs and OCs observed in AIMD. Along with the local coordination structure, we also compared enthalpies of vaporization. Overall, calculated bulk properties from FFMD compared reasonably well with experimental values, suggesting that small improvements within the FF should focus on parameters that adjust the bulk liquid structures of these compounds.

r2SCAN-3c: A "Swiss army knife" composite electronic-structure method

The recently proposed r²SCAN meta-generalized-gradient approximation (mGGA) of Furness and co-workers is used to construct an efficient composite electronic-structure method termed r²SCAN-3c. To this end, the unaltered r²SCAN functional is combined with a tailor-made triple-ζ Gaussian atomic orbital basis set as well as with refitted D4 and geometrical counter-poise corrections for London-dispersion and basis set superposition error. The performance of the new method is evaluated for the GMTKN55 database covering large parts of chemical space with about 1500 data points, as well as additional benchmarks for non-covalent interactions, organometallic reactions, and lattice energies of organic molecules and ices, as well as for the adsorption on polar salt and non-polar coinage-metal surfaces. These comprehensive tests reveal a spectacular performance and robustness of r²SCAN-3c: It by far surpasses its predecessor B97-3c at only twice the cost and provides one of the best results of all semi-local density-functional theory (DFT)/QZ methods ever tested for the GMTKN55 database at one-tenth of the cost. Specifically, for reaction and conformational energies as well as non-covalent interactions, it outperforms prominent hybrid-DFT/QZ approaches at two to three orders of magnitude lower cost. Perhaps, the most relevant remaining issue of r²SCAN-3c is self-interaction error (SIE), owing to its mGGA nature. However, SIE is slightly reduced compared to other (m)GGAs, as is demonstrated in two examples. After all, this remarkably efficient and robust method is chosen as our new group default, replacing previous composite DFT and partially even expensive high-level methods in most standard applications for systems with up to several hundreds of atoms.

Model Studies on the Formation of the Solid Electrolyte Interphase: Reaction of Li with Ultrathin Adsorbed Ionic-Liquid Films and Co3 O4 (111) Thin Films

In this work we aim towards the molecular understanding of the solid electrolyte interphase (SEI) formation at the electrode electrolyte interface (EEI). Herein, we investigated the interaction between the battery‐relevant ionic liquid (IL) 1‐butyl‐1‐methylpyrrolidinium bis(trifluoromethylsulfonyl)imide (BMP‐TFSI), Li and a Co₃O₄(111) thin film model anode grown on Ir(100) as a model study of the SEI formation in Li‐ion batteries (LIBs). We employed mostly X‐ray photoelectron spectroscopy (XPS) in combination with dispersion‐corrected density functional theory calculations (DFT‐D3). If the surface is pre‐covered by BMP‐TFSI species (model electrolyte), post‐deposition of Li (Li⁺ ion shuttle) reveals thermodynamically favorable TFSI decomposition products such as LiCN, Li₂NSO₂CF₃, LiF, Li₂S, Li₂O₂, Li₂O, but also kinetic products like Li₂NCH₃C₄H₉ or LiNCH₃C₄H₉ of BMP. Simultaneously, Li adsorption and/or lithiation of Co₃O₄(111) to Li_nCo₃O₄ takes place due to insertion via step edges or defects; a partial transformation to CoO cannot be excluded. Formation of Co⁰ could not be observed in the experiment indicating that surface reaction products and inserted/adsorbed Li at the step edges may inhibit or slow down further Li diffusion into the bulk. This study provides detailed insights of the SEI formation at the EEI, which might be crucial for the improvement of future batteries.

Quantification of Noncovalent Interactions in Azide-Pnictogen, -Chalcogen, and -Halogen Contacts

The noncovalent interactions between azides and oxygen‐containing moieties are investigated through a computational study based on experimental findings. The targeted synthesis of organic compounds with close intramolecular azide–oxygen contacts yielded six new representatives, for which X‐ray structures were determined. Two of those compounds were investigated with respect to their potential conformations in the gas phase and a possible significantly shorter azide–oxygen contact. Furthermore, a set of 44 high‐quality, gas‐phase computational model systems with intermolecular azide–pnictogen (N, P, As, Sb), –chalcogen (O, S, Se, Te), and –halogen (F, Cl, Br, I) contacts are compiled and investigated through semiempirical quantum mechanical methods, density functional approximations, and wave function theory. A local energy decomposition (LED) analysis is applied to study the nature of the noncovalent interaction. The special role of electrostatic and London dispersion interactions is discussed in detail. London dispersion is identified as a dominant factor of the azide–donor interaction with mean London dispersion energy‐interaction energy ratios of 1.3. Electrostatic contributions enhance the azide–donor coordination motif. The association energies range from −1.00 to −5.5 kcal mol⁻¹.

Macroscopic quantum electrodynamics and density functional theory approaches to dispersion interactions between fullerenes

The processing and material properties of commercial organic semiconductors, for e.g. fullerenes is largely controlled by their precise arrangements, specially intermolecular symmetries, distances and orientations, more specifically, molecular polarisabilities. These supramolecular parameters heavily influence their electronic structure, thereby determining molecular photophysics and therefore dictating their usability as n-type semiconductors. In this article we evaluate van der Waals potentials of a fullerene dimer model system using two approaches: (a) Density Functional Theory and, (b) Macroscopic Quantum Electrodynamics, which is particularly suited for describing long-range van der Waals interactions. Essentially, we determine and explain the model symmetry, distance and rotational dependencies on binding energies and spectral changes. The resultant spectral tuning is compared using both methods showing correspondence within the constraints placed by the different model assumptions. We envision that the application of macroscopic methods and structure/property relationships laid forward in this article will find use in fundamental supramolecular electronics.

A Systematic Study of DFT Performance for Geometry Optimizations of Ionic Liquid Clusters

Clusters of two ion pairs of imidazolium-based ionic liquids were optimized with 43 different levels of theory, including DFT functionals and MP2-based methods combined with varying Dunning's basis sets, and added dispersion corrections. Better preforming DFT functionals were then applied to clusters consisting of four ion pairs. Excellent performance of some DFT functionals for the two ion pair clusters did not always match that of the four ion-paired clusters despite interionic distances remaining constant between the optimized two and four ion-paired clusters of the same ionic liquid. Combinations of DFT functional and basis set such as ωB97X-D/cc-pVDZ, M06-2X/aug-cc-pVDZ, B3LYP-D3/cc-pVTZ, and TPSS-D3/cc-pVTZ gave excellent results for geometry optimization of two ion-paired clusters of imidazolium ionic liquids but gave larger deviations when applied to the four ion-paired clusters of varying ionic liquids. Empirical dispersion corrections were seen to be crucial in correctly capturing correlation effects in the studied ionic liquid clusters, becoming more important in larger clusters. Dunning's double-ζ basis set, cc-pVDZ, is associated with the smallest root mean squared deviations for geometries; however, it also produces the largest deviations in total electronic energies. ωB97X-D and M06-2X produced the best performance with the augmented version of this basis set. The triple-ζ basis set, cc-pVTZ, leads to the best performance of most of the DFT functionals (especially the dispersion-corrected ones) used, whereas its augmented version, aug-cc-pVTZ, was not seen to improve results. The combinations of functional and basis set that gave the best geometry and energetics in both two and four ion-paired clusters were PBE-D3/cc-pVTZ, ωB97X-D/aug-cc-pVDZ, and BLYP-D3/cc-pVTZ. All three combinations are recommended for geometry optimizations of larger clusters of ionic liquids. PBE-D3/cc-pVTZ performed the best with an average deviation of 2.3 kJ mol^-1 and a standard deviation of 3.4 kJ mol^-1 for total electronic energy when applied to four ion-paired clusters. Geometries optimized with FMO2-SRS-MP2/cc-pVTZ produced total energy within 2.0 kJ mol^-1 off the benchmark in two ion-paired clusters, with the cc-pVDZ basis set performing unsurprisingly poorly with the same method. The error increased to 4.8 kJ mol^-1 on average in four ion-paired clusters, with the smallest RMSD deviations in geometries when compared to the benchmark ones. This study is the first report that investigated the performance of DFT functionals for two and four ion-paired clusters of a wide range of ionic liquids consisting of commonly used cations such as pyrrolidinium, imidazolium, pyridinium, and ammonium. It also identified the importance of assessing the performance of quantum chemical methods for ionic liquids on a variety of cation-anion combinations.

Modeling the Physicochemical Properties of Natural Deep Eutectic Solvents

Natural deep eutectic solvents (NADES) are mixtures of naturally derived compounds with a significantly decreased melting point owing to specific interactions among the constituents. NADES have benign properties (low volatility, flammability, toxicity, cost) and tailorable physicochemical properties (by altering the type and molar ratio of constituents); hence, they are often considered to be a green alternative to common organic solvents. Modeling the relation between their composition and properties is crucial though, both for understanding and predicting their behavior. Several efforts have been made to this end. This Review aims at structuring the present knowledge as an outline for future research. First, the key properties of NADES are reviewed and related to their structure on the basis of the available experimental data. Second, available modeling methods applicable to NADES are reviewed. At the molecular level, DFT and molecular dynamics allow density differences and vibrational spectra to be interpreted, and interaction energies to be computed. Additionally, properties at the level of the bulk medium can be explained and predicted by semi-empirical methods based on ab initio methods (COSMO-RS) and equation of state models (PC-SAFT). Finally, methods based on large datasets are discussed: models based on group-contribution methods and machine learning. A combination of bulk-medium and dataset modeling allows qualitative prediction and interpretation of phase equilibria properties on the one hand, and quantitative prediction of melting point, density, viscosity, surface tension, and refractive index on the other. Multiscale modeling, combining molecular and macroscale methods, is expected to strongly enhance the predictability of NADES properties and their interaction with solutes, and thus yield truly tailorable solvents to accommodate (bio)chemical reactions.

The CO2 Absorption in Flue Gas Using Mixed Ionic Liquids

Because of the appealing properties, ionic liquids (ILs) are believed to be promising alternatives for the CO₂ absorption in the flue gas. Several ILs, such as [NH₂emim][BF₄], [C₄mim][OAc], and [NH₂emim[OAc], have been used to capture CO₂ of the simulated flue gas in this work. The structural changes of the ILs before and after absorption were also investigated by quantum chemical methods, FTIR, and NMR technologies. However, the experimental results and theoretical calculation showed that the flue gas component SO₂ would significantly weaken the CO₂ absorption performance of the ILs. SO₂ was more likely to react with the active sites of the ILs than CO₂. To improve the absorption capacity, the ionic liquid (IL) mixture [C₄mim][OAc]/ [NH₂emim][BF₄] were employed for the CO₂ absorption of the flue gas. It is found that the CO₂ absorption capacity would be increased by about 25%, even in the presence of SO₂. The calculation results suggested that CO₂ could not compete with SO₂ for reacting with the IL during the absorption process. Nevertheless, SO₂ might be first captured by the [NH₂emim][BF₄] of the IL mixture, and then the [C₄mim][OAc] ionic liquid could absorb more CO₂ without the interference of SO₂.

Heterogeneity in the microstructure and dynamics of tetraalkylammonium hydroxide ionic liquids: insight from classical molecular dynamics simulations and Voronoi tessellation analysis

Microscopic structural and dynamic heterogeneities were investigated for three ionic liquids (ILs), tetraethylammonium hydroxide, tetrapropylammonium hydroxide, and tetrabutylammonium hydroxide employing classical molecular dynamics (MD) simulations.

Counting cations involved in cationic clusters of hydroxy-functionalized ionic liquids by means of infrared and solid-state NMR spectroscopy

Sensitive probe of like-charge attraction: analyzing infrared spectra allows counting the number of cations involved in clusters of opposite (c–a) and like-charged (c–c) ions in ionic liquids. This approach is also applicable to molecular liquids.

Hydrogen Bonding Between Ions of Like Charge in Ionic Liquids Characterized by NMR Deuteron Quadrupole Coupling Constants-Comparison with Salt Bridges and Molecular Systems

We present deuteron quadrupole coupling constants (DQCC) for hydroxyl-functionalized ionic liquids (ILs) in the crystalline or glassy states characterizing two types of hydrogen bonding: The regular Coulomb-enhanced hydrogen bonds between cation and anion (c-a), and the unusual hydrogen bonds between cation and cation (c-c), which are present despite repulsive Coulomb forces. We measure these sensitive probes of hydrogen bonding by means of solid-state NMR spectroscopy. The DQCCs of (c-a) ion pairs and (c-c) H-bonds are compared to those of salt bridges in supramolecular complexes and those present in molecular liquids. At low temperatures, the (c-c) species successfully compete with the (c-a) ion pairs and dominate the cluster populations. Equilibrium constants obtained from molecular-dynamics (MD) simulations show van't Hoff behavior with small transition enthalpies between the differently H-bonded species. We show that cationic-cluster formation prevents these ILs from crystallizing. With cooling, the (c-c) hydrogen bonds persist, resulting in supercooling and glass formation.

Molecular Simulation of the Catalytic Regeneration of nBuLi through a Hydrometalation Route

Efficient regeneration of organolithium compounds is a challenging aspect in the process of novel organometathetical catalytic cycles. One of these catalytic cycles is a newly suggested method for Mg production from seawater that capitalizes on the rich chemistry of Grignard reagents. The proposed three-step catalytic cycle with Cp₂ MCl L catalyst ( M = Ti, Zr; L = select organic ligands) requires the regeneration of ⁿBuLi from Li(s), butene, and H₂. The potential of this approach is evaluated with density functional theory-based molecular simulations. The results reveal that the high affinity of Li toward Cl and N results in the formation of alkanes, and the strong coupling between the catalyst and BuLi leads to catalyst deactivation. To improve its catalytic performance, we proposed the use of a diamine cocatalyst and a modified catalyst with a ligand that does not contain N, which would help release BuLi from the vicinity of the catalytic center. Ab initio molecular dynamics simulations at 298 K in explicit solvent (THF) were used to estimate the Gibbs free energetics and equilibrium constants obtained from the vibrational density of states using velocity autocorrelation functions. The results show a marked improvement in the free energetics with lower barriers toward the completion of the catalytic cycle and suppression of deactivation channels.

Cation influence on heterocyclic ammonium ionic liquids: a molecular dynamics study

Four different ionic liquids (ILs) consisting of the bis(trifluoromethanesulfonyl)imide ([NTf₂]⁻) anion, with structurally similar systematically varying cations, are investigated herein through classical molecular dynamics.

Solvation Properties of Silver and Copper Ions in a Room Temperature Ionic Liquid: A First-Principles Study

Understanding the behavior of metal ions in room temperature ionic liquids (ILs) is essential for predicting and optimizing performance for technologies like metal electrodeposition; however, many mechanistic details remain enigmatic, including the solvation properties of the ions in ILs and how they are governed by the intrinsic interaction between the ions and the liquid species. Here, we utilize first-principles molecular dynamics simulations to unravel and compare the key structural properties of Ag⁺ and Cu⁺ ions in a common room temperature IL, 1-ethyl-3-methylimidazolium trifluoromethanesulfonate. We find that, when compared to Cu⁺, the larger Ag⁺ shows a more disordered and flexible solvation structure with a more frequent exchange of the IL species between its solvation shells. In addition, our simulations reveal an interesting analog in the solvation behavior of the ions in the IL and aqueous environments, particularly in the effect of the ion electronic structures on their solvation properties. This work provides fundamental understanding of the intrinsic properties of the metal ions in the IL, while offering mechanistic understanding and strategy for future selection of ILs for metal electrodeposition processes.

Multiscale Simulation of the Interaction and Adsorption of Ions on a Hydrophobic Graphene Surface

The adsorption of ions on a graphene surface is very important to control relevant graphene-based processes. In this work, a multiscale simulation was carried out to study the adsorption of Na⁺ /Cl^- ions on graphene by combining quantum mechanics calculations and molecular dynamics (MD) simulations. The interaction energies of the ions with graphene were computed using density functional theory (DFT). It was found that the ions show strong interaction with a graphene cluster and the overwhelming portion of the interaction energy is the ion-π orbital interaction. The large orbital interaction can be ascribed to the two contributions arising from the ion-induced polarization of graphene and the charge transfer between ion and graphene. Their different contribution degrees reveal that the polarization effect plays a main role on the orbital interaction for ion adsorption. Comparatively, for Na/Cl atom adsorption, the charge transfer shows large part to the orbital interaction with weak atom-induced polarization. The obtained interaction energies were applied to develop new interaction potentials between ion and graphene, and then MD simulations were used to study the interfacial adsorption behavior of Na⁺ /Cl^- aqueous solution onto the graphene surface. Due to enhanced ion-π interactions, Na⁺ /Cl^- cooperatively demonstrates a strong ion adsorption layer through direct contact with the hydrophobic graphene surface. Our simulation result presents a new understanding of ion-graphene interactions.