Data to understand the nature of non-covalent interactions in the thiophene clusters.

We have reported herein the data to understand the nature and number of non-covalent interactions that stabilize the structures of the thiophene clusters. In addition, we have also provided the optimized Cartesian coordinates of all the structures of the investigated thiophene clusters. Initially, the geometries have been generated using the ABCluster code which performs a global optimization to locate local and global minima structures of molecular clusters. The located geometries have been optimized at the MP2/aug-cc-pVDZ level of theory using Gaussian 16 suite of programs. To understand the nature of non-covalent interactions, we have performed a quantum theory of atoms in molecules (QTAIM) analysis on all the structures of the thiophene dimer. Furthermore, the QTAIM analysis has been performed also on the most stable structure of the thiophene trimer and tetramer. We have used the AIMAll program to perform the QTAIM analysis. The data reported in this paper contains the critical points, the bonds paths and their related properties, for each investigated structures. Besides, the data contains the optimized Cartesian coordinates of all the investigated structures of the thiophene clusters. This can be use for any further investigations involving thiophene clusters. For further information and analysis, the reader is referred to the original related research article (Malloum and Conradie, 2022).

Specifications Table

Value of the Data

We have reported the critical points, bond paths and their related properties for the structures of the thiophene dimer. These data will be useful to understand the nature and number of non-covalent interactions that stabilize the structures of the thiophene clusters.

The optimized Cartesian coordinates of the thiophene clusters will save computational time for all further investigations involving thiophene clusters.

Using this data one would be able to determine all non-covalent bonds in thiophene clusters and therefore, will provide proper understanding of non-covalent bonds in aromatic and heterocyclic molecular clusters.

The optimized Cartesian coordinates of the structures of the thiophene clusters will be useful for investigations involving solvation processes and/or ions transfer in thiophene solvent.

The data provided in this paper will be of great interest for investigations using quantum cluster equilibrium theory to compute the thermodynamics properties of the liquid thiophene.

Data Description

We have reported in Fig. 1 the structures of the thiophene dimer, their relative energies, critical points, and bond paths. Each sub-caption of the structure reports the name and the relative energy in kcal/mol as calculated at the MP2/aug-cc-pVDZ level of theory. For each structure, critical points (Bond critical points, BCP, are reported in green color, Ring critical points, RCP, are in red color while cage critical points, CCPs are in blue color) and bond paths (covalent bonds in solid lines and non-covalent bonds in dash lines) are also given in the figure. The properties of the critical points of all the structures of the thiophene dimer, thiophene trimer, and those of the most stable structures of the thiophene tetramer are reported in Excel files as supplementary material. Each of the tables in the Excel file has seven columns which describe respectively, the name of the bond critical point (name), the two atoms involved in the bond critical point (Atoms), the electron density at the bond critical point (), the Laplacian of the electron density at the bond critical point (), the ellipticity of the bond (Ellipticity), the electron kinetic energy density at the bond critical point (K) and the difference between the bond path length with a corresponding geometric bond length (BPL-GBL_I). In the supplementary material, we have also reported the Cartesian coordinates of all the investigated structures of the thiophene dimer, trimer and tetramer as optimized at the MP2/aug-cc-pVDZ level of theory.

Fig. 1

Critical points and bond paths of all the investigated structures of the thiophene dimer. Bond critical points (BCP) are reported in green color, Ring critical points (RCP) are in red color while cage critical points (CCPs) are in blue color. Structures of Thio2_1 and Thio2_10 are also reported in [1].

Experimental Design, Materials and Methods

The most important task in exploring the potential energy surfaces (PESs) of molecular clusters is the location of possible stable structures. Thus, to obtain the Cartesian coordinates of the structures reported in the supplementary material, we explored thoroughly the PESs of thiophene clusters from dimer to tetramer using classical molecular dynamics as implemented in ABCluster [2], [3]. ABCluster generates possible structures and classifies them following their relative energies based on classical potential energy . The classical potential energy implemented in ABCluster is constituted of Lenard-Jones and electrostatics interactions:Where and are the indices of the molecules, and are the indices of the atoms in molecules and , respectively. is the distance between atom and . The accuracy of the exploration of a given PES using ABCluster is based on the setting of three parameters: the scout limit , the size of the population of trial solutions and the maximum cycle number . For the potential energy surfaces of thiophene clusters, we used , and . For each cluster size , we generated 500 geometries using ABCluster. ABCluster is based on the bee colony algorithm. The authors of ABCluster have stated that a scout limit between 3 and 5 is enough to achieve accurate exploration of a given PES. It is worth noting that ABCluster has been used in our previous works to explore PESs of molecular clusters [4], [5], [6], [7], [8]. For further reading on ABCluster, the reader is referred to our previous works and the original works of Zhang and Dolg [2], [3]. In practical use, ABCluster generates several geometries that are similar one to another. Thus, from the generated geometries, we selected those that are different one from another. The selected geometries are fully optimized at the MP2/aug-cc-pVDZ level of theory. Optimizations have been performed using Gaussian 16 suite of program. We used the tight option of the optimization for high accuracy. To ensure accurate location of stable minima, frequencies calculations have been also performed at the same level of theory. All the located structures have positive frequencies highlighting their stability. After optimizations, some structures have been optimized to the same isomer. The redundant structures have been eliminated from the final structures reported in the supplementary material.

After locating the structures, we used their Gaussian formatted checkpoint files to perform (to generate the data related to) the quantum theory of atoms in molecule (QTAIM) analysis. The QTAIM analysis is performed to provide data necessary to understand the nature of non-covalent bondings in thiophene clusters. QTAIM analyses the topology of the electron density, , of a given molecule to locate its critical points. These critical points are obtained by vanishing the first order derivatives of the electron density . Depending on the sign of the second order derivatives, there are usually four critical points used in QTAIM analysis: (3, -3) atom critical point (ACP); (3, -1) bond critical point (BCP); (3, 1) ring critical point (RCP); and (3, 3) cage critical point (CCP). A bond path connects two atom critical points along which lies a bond critical points. Bader [9] stated that the analysis of the topology of the electron density and the properties of bond critical points can provide universal description of chemical bondings. When the second order derivatives of the electron density (the Laplacian of the electron density ) is positive (respectively negative), it indicates that the electron is locally depleted (respectively concentrated) [10]. The values of and at the BCP can be used to distinguish between different types of bonding (covalent, ionic bonding, hydrogen bonding, and van der Waals interactions) [11], [12]. The range proposed for and at a bond critical point for a hydrogen bond to exist, is 0.002-0.035 for and 0.024-0.139 for [13]. A positive value of at a bond critical point indicates a non-covalent bonding, while a negative value of at a bond critical point indicates a covalent bonding [14]. Generally, the larger the value of the electron density at a bond critical point, the stronger the corresponding bonding [15], [16]. Using the description above, one can now analyse the data provided as excel file in the supplementary material to understand the non-covalent interactions in thiophene clusters.

CRediT authorship contribution statement

Alhadji Malloum: Conceptualization, Methodology, Validation, Formal analysis, Investigation, Data curation, Writing – original draft, Visualization. Jeanet Conradie: Resources, Visualization, Writing – review & editing, Supervision, Funding acquisition, Project administration.

Declaration of Competing Interest

The authors declare that they have no known competing financial interests or personal relationships which have, or could be perceived to have, influenced the work reported in this article.

Subject	Chemistry
Specific subject area	Physical and Theoretical Chemistry
Type of data	TableFigure
How data were acquired	Data were obtain using Gaussian 16 computational chemistry program. Data for quantum theory of atoms in molecule analysis (QTAIM) are calculated using AIMAll program.
Data format	RawAnalyzed
Parameters for datacollection	Raw data (optimized Cartesian coordinates of the structures) are extracted directly from the Gaussian output files. Analyzed data (description of critical points and their related figures) were obtained from the AIMAll program.
Description of datacollection	Optimization of the geometries have been performed using the resources of the Center of High Performance Computing (CHPC), South Africa. The QTAIM analysis has been performed in our laboratory (Physical Chemistry Laboratory of the Department of Chemistry).
Data source location	Institution: Department of Chemistry, University of the Free StateCity/Town/Region: BloemfonteinCountry: South Africa
Data accessibility	With the article
Related researcharticle	Alhadji Malloum and Jeanet Conradie, Non-Covalent Interactions in Small Thiophene Clusters, J. Mol. Liq. 347 (2022) 118301 [1]. https://doi.org/10.1016/j.molliq.2021.118301