Structure and Stability of the Ionic Liquid Clusters [EMIM]n[BF4]n+1- (n = 1-9): Implications for Electrochemical Separations

A many-body energy decomposition analysis (MB-EDA) scheme based on a target state optimization self-consistent field (TSO-SCF) method

In this paper, we combine an energy decomposition analysis (EDA) scheme with many-body expansion (MBE) to develop a MB-EDA method to study the cooperative and anti-cooperative effects in molecular cluster systems. Based on the target state optimization self-consistent field (TSO-SCF) method, the intermolecular interaction energy can be decomposed into five chemically meaningful terms, i.e., electrostatic, exchange, polarization, charge transfer and dispersion interaction energies. MB-EDA can decompose each of these terms in MBE. This MB-EDA has been applied to 3 different cluster systems: water clusters, ionic liquid clusters, and acetonitrile-methane clusters. This reveals that electrostatic, exchange, and dispersion interactions are highly pairwise additive in all systems. In water and ionic liquid clusters, the many-body effects are significant in both polarization and charge transfer interactions, but are cooperative and anti-cooperative, respectively. For acetonitrile-methane clusters, which do not involve hydrogen bonds or charge-charge Coulombic interactions, the many-body effects are quite small. The chemical origins of different many-body effects are deeply analyzed. The MB-EDA method has been implemented in Qbics (https://qbics.info) and can be a useful tool for understanding the many-body behavior in molecular aggregates at the quantum chemical level of theory.

Electrochemical Imaging of Precisely-Defined Redox and Reactive Interfaces

Understanding the diverse electrochemical reactions occurring at electrode-electrolyte interfaces (EEIs) is a critical challenge to developing more efficient energy conversion and storage technologies. Establishing a predictive molecular-level understanding of solid electrolyte interphases (SEIs) is challenging due to the presence of multiple intertwined chemical and electrochemical processes occurring at battery electrodes. Similarly, chemical conversions in reactive electrochemical systems are often influenced by the heterogeneous distribution of active sites, surface defects, and catalyst particle sizes. In this mini review, we highlight an emerging field of interfacial science that isolates the impact of specific chemical species by preparing precisely-defined EEIs and visualizing the reactivity of their individual components using single-entity characterization techniques. We highlight the broad applicability and versatility of these methods, along with current state-of-the-art instrumentation and future opportunities for these approaches to address key scientific challenges related to batteries, chemical separations, and fuel cells. We establish that controlled preparation of well-defined electrodes combined with single entity characterization will be crucial to filling key knowledge gaps and advancing the theories used to describe and predict chemical and physical processes occurring at EEIs and accelerating new materials discovery for energy applications.

Tracing Vortex Clustering in a Superconductor by the Magnetic Flux Distribution

By investigating spatial configurations of the intermediate mixed state in an intertype superconductor, it is shown that vortex clustering can be characterized by the sample averaged distribution of the penetrating magnetic field. The clustering is manifested in the two-peak structure of the distribution. The second peak indicates a spot a material occupies in the phase diagram of superconductivity types. The conclusions are general and do not depend on details of the model.

Target State Optimized Density Functional Theory for Electronic Excited and Diabatic States

A flexible self-consistent field method, called target state optimization (TSO), is presented for exploring electronic excited configurations and localized diabatic states. The key idea is to partition molecular orbitals into different subspaces according to the excitation or localization pattern for a target state. Because of the orbital-subspace constraint, orbitals belonging to different subspaces do not mix. Furthermore, the determinant wave function for such excited or diabatic configurations can be variationally optimized as a ground state procedure, unlike conventional ΔSCF methods, without the possibility of collapsing back to the ground state or other lower-energy configurations. The TSO method can be applied both in Hartree-Fock theory and in Kohn-Sham density functional theory (DFT). The density projection procedure and the working equations for implementing the TSO method are described along with several illustrative applications. For valence excited states of organic compounds, it was found that the computed excitation energies from TSO-DFT and time-dependent density functional theory (TD-DFT) are of similar quality with average errors of 0.5 and 0.4 eV, respectively. For core excitation, doubly excited states and charge-transfer states, the performance of TSO-DFT is clearly superior to that from conventional TD-DFT calculations. It is shown that variationally optimized charge-localized diabatic states can be defined using TSO-DFT in energy decomposition analysis to gain both qualitative and quantitative insights on intermolecular interactions. Alternatively, the variational diabatic states may be used in molecular dynamics simulation of charge transfer processes. The TSO method can also be used to define basis states in multistate density functional theory for excited states through nonorthogonal state interaction calculations. The software implementing TSO-DFT can be accessed from the authors.

Conformational energies of microsolvated Na+ clusters with protic and aprotic solvents from GFNn-xTB methods

Performance of contemporary tight-binding semiempirical GFNn-xTB methods for the conformational energies of singly charged sodium clusters Na⁺ (S)_n (n = 4-8) with 3 protic and 8 aprotic solvents is examined against the reference RI-MP2/CBS method. The median Pearson correlation coefficients of ρ = 0.84 (GFN2-xTB) and ρ = 0.82 (GFN1-xTB) do not give the clear preference to any tested approach. GFN1-xTB method demonstrates more stable performance than its GFN2-xTB successor with the average mean absolute errors (MAEs)/mean signed errors (MSEs) of 1.2/0.2 and 2.3/1.6 kcal mol^-1 , respectively. Conformational energies produced by the computationally efficient DFT functional PBE and double-ζ basis set complemented with -D3(BJ) dispersion correction are suitable for the preliminary sampling (median ρ = 0.93), but should be used with a caution for the calculations of the average ensemble properties (MAE/MSE = 1.7/1.1 kcal mol^-1 ). Higher-ranking PBE0-D3(BJ) and ωB97M-V with triple-ζ basis sets yield significantly lower MAEs/MSEs of 0.55/0.20 and 0.51/0.23 kcal mol^-1 , respectively.

Surface Properties of Ionic Liquids: A Mass Spectrometric View Based on Soft Cluster-Induced Desorption

Desorption/ionization induced by neutral clusters (DINeC) in combination with mass spectrometry (MS) was used for the investigation of the molecular composition of the surface of ionic liquids (IL). Based on the surface sensitivity of DINeC-MS, accumulation of either cations or anions was discriminated on the surface of bulk IL depending on the molecular structure of the IL components. In particular, cations with long alkyl chains aggregate on the surface, but this tendency is more reduced the larger the respective anion is; in the case of larger anions and smaller cations, it can be even reversed. For thin layers of IL, the ratio between cations and anions as detected in the mass spectra was found to be further influenced by the substrate surface.

Efficient evaluation of electrostatic potential with computerized optimized code

The evaluation of molecular electrostatic potential (ESP) is a performance bottleneck for many computational chemical tasks like restrained ESP charge fitting or quantum mechanics/molecular mechanics simulations. In this paper, an efficient algorithm for the evaluation of ESP is proposed. It regroups the expression in terms of primitive Gaussian type orbitals (GTOs) with identical angular momentum types and nuclei centers. Each term is calculated using a computerized optimized code. This algorithm was integrated into the wavefunction analysis program Multiwfn and was tested on several large systems. In the cases of dopamine and remdesivir, the performance of this algorithm was comparable to or better than some popular state-of-the-art codes. For meta1-organic framework-5, where the number of GTOs and ESP points is 4840 and 259 262, respectively, our code could finish the evaluation in 1874 seconds on ordinary hardware. It also exhibits good parallelization scaling. The source code of this algorithm is freely available and can become a useful tool for computational chemists.

Vortex Interactions and Clustering in Thin Superconductors

Cluster formation is a focus of interdisciplinary research in both chemistry and physics. Here we discuss the exotic example of this phenomenon in the vortex matter of a thin superconductor. In superconducting films, the clustering takes place because of particular properties of the vortex interactions in the crossover or intertype regime between superconductivity types I and II. These interactions are controlled by the two parameters that are responsible for the crossover, Ginzburg-Landau parameter κ, which specifies the superconducting material of the film, and film thickness d, which controls effects due to stray magnetic fields outside the sample. We demonstrate that their competition gives rise to a complex spatial dependence of the interaction potential between vortices, favoring the formation of chainlike vortex clusters.

Gas-Phase Binding Energies and Dissociation Dynamics of 1-Alkyl-3-Methylimidazolium Tetrafluoroborate Ionic Liquid Clusters

Ionic liquids (ILs) have become increasingly popular due to their useful and unique properties, yet there are still many unanswered questions regarding their fundamental interactions. In particular, details regarding the nature and strength of the intrinsic cation-anion interactions and how they influence the macroscopic properties of ILs are still largely unknown. Elucidating the molecular-level details of these interactions is essential to the development of better models for describing ILs and enabling the purposeful design of ILs with properties tailored for specific applications. Current uses of ILs are widespread and diverse and include applications for energy storage, electrochemistry, designer/green solvents, separations, and space propulsion. To advance the understanding of the energetics, conformations, and dynamics of gas-phase IL clustering relevant to space propulsion, threshold collision-induced dissociation approaches are used to measure the bond dissociation energies (BDEs) of the 2:1 clusters of 1-alkyl-3-methylimidazolium cations and tetrafluoroborate, [2C_nmim:BF₄]⁺. The cation, [C_nmim]⁺, is varied across the series, 1-ethyl-3-methylimidazolium [C₂mim]⁺, 1-butyl-3-methylimidazolium [C₄mim]⁺, 1-hexyl-3-methylimidazolium [C₆mim]⁺, and 1-octyl-3-methylimidazolium [C₈mim]⁺, to examine the structural and energetic effects of the size of the 1-alkyl substituent on binding. Complementary electronic structure calculations are performed to determine the structures and energetics of the [C_nmim]⁺ and [BF₄]^- ions and their binding preferences in the (C_nmim:BF₄) ion pairs and [2C_nmim:BF₄]⁺ clusters. Several levels of theory, B3LYP, B3LYP-GD3BJ, and M06-2X, using the 6-311+G(d,p) basis set for geometry optimizations and frequency analyses and the 6-311+G(2d,2p) basis set for energetics, are benchmarked to examine their abilities to properly describe the nature of the binding interactions and to reproduce the measured BDEs. The modest structural variation among these [C_nmim]⁺ cations produces only minor structural changes and variation in the measured BDEs of the [2C_nmim:BF₄]⁺ clusters. Present findings indicate that the dominant cation-anion interactions involve the 3-methylimidazolium moieties and that these clusters are sufficiently small that differences in packing effects associated with the variable length of the 1-alkyl substituents are not yet significant.