Scrutinizing ion-π and ion-σ interactions using the noncovalent index and energy decomposition analysis

Towards a predictive model for polymer solubility using the noncovalent interaction index: polyethylene as a case study

In this work we present the development of a novel, quantitative solubility descriptor based on the non-covalent interaction index. It is presented as a more insightful alternative to Hansen's solubility parameters and the COSMO model to assess and predict polymer solubility in different solvents. To this end, we studied the solvation behaviour as a function of the chain length of a single chain of arguably the most simple polymer, polyethylene, in anisole (solvent) and methanol (poor solvent) via molecular dynamics simulations. It was found that in anisole the solute maximized its interface with the solvent, whereas in methanol the macromolecule formed rod-like structures by folding on itself once the chain length surpassed a certain barrier. We assessed this behaviour - which can be related to solubility - quantitatively and qualitatively via well-known descriptors, namely the solvation free energy, and the solvent accessible surface area. In addition, we propose the non-covalent interaction (NCI) index as a versatile descriptor, providing information on the strength, as well as the nature, of the solute-solvent interactions, the solute's intramolecular interactions and on the solute's conformation, both qualitatively and quantitatively. Finally, as a quantitative measure for solubility, defined in this context as the solute's tendency to maximize its interactions with the solvent, we propose two new NCI-based descriptors: the relative integrated NCI density and the integrated NCI difference. The former represents the quantitative difference in solute-solvent interactions between a fully extended coil and the actual conformation during simulation and the latter the quantitative difference between the intermolecular (solute-solvent) and the intramolecular (in the solute) non-covalent interactions. The easy interpretation and calculation of these novel quantities open up the possibility of fast, reliable and insightful high-throughput screening of different (anti)solvent and solute combinations.

Unveiling the Impact of Aggregation on Optical Anisotropy of Triazaacephenanthrylene Single Crystals. A Combined Quantum Crystallography and Conceptual Density Functional Theory Approach

Triazaacephenanthrylene (TAAP) triclinic single crystals show substantial optical anisotropy of absorption and fluorescence. The maximum effect can be correlated with the direction perpendicular to the plane of chromophores connected in a head-to-tail manner via weak dispersive interactions. This phenomenon is uncommon as usually the existence of postulated π···π interactions between the molecules forming dimers or stacks cause quenching of fluorescence. Herein we present a comprehensive study of inter- and intramolecular interactions in the crystal of TAAP enriched with the investigation of aromaticity. Our results show that intramolecular interactions stabilize the overall conformation of the molecule whereas dispersive forces determine the aggregation between TAAP molecules. In fact, there is no conventional π···π interaction between the molecules in the dimer. Instead, we observed a close contact between the lone pair of the bridgehead N10B atom and π-deficient pyrazine ring from an adjacent molecule. Optical anisotropy in TAAP crystals was directly correlated with the alignment of the molecular transition dipole moments caused by specific molecular self-assembly.

An assessment of the random-phase approximation functional and characteristics analysis for noncovalent cation-π interactions

The binding energy is of great importance in understanding the formation and stability of noncovalent interactions. However, the determination of the binding energy with high precision and efficiency in medium- and long-range noncovalent interactions is still challenging for quantum chemistry. Here, we assess the performance of random-phase approximation (RPA), a fully non-local fifth-rung of the Jacob ladder functional, in determining the binding energy of cation-π systems (cation = Li⁺, Na⁺, Be²⁺, Mg²⁺, Al⁺, and NH₄⁺; π = C₆H₆), which, to the best of our knowledge, has not been investigated. Using experimental results as the benchmark, we systematically compared the RPA method to the other ab initio methods (DFT/B3LYP, MP2, CCSD(T), and QCISD(T)) both in calculation accuracy and efficiency. From the perspective of accuracy, RPA is the best among these approaches, followed by the CCSD(T) and QCISD(T) methods. DFT/B3LYP and MP2 provide the worst accuracy. In addition, the computational efficiency of RPA is much faster than that of CCSD(T) and QCISD(T). We believe that RPA is a robust method for the precise description of medium- and long-range noncovalent interactions and is capable of providing benchmarking data. The interaction strength and interaction nature of cation-π systems are further analyzed by atoms in molecules (AIM) and the color-mapped reduced density gradient (RDG) isosurface, which are consistent with the characteristics of a typical cation-π interaction.

Multi-center covalency: revisiting the nature of anion–π interactions

Plot of the delocalization index, δ(X⁻,Ω), scaled by the internuclear distance, R_X,Ω, versus the exchange–correlation potential energy, V_XC(X⁻,Ω), for anion–π complexes.