Identifying and Analyzing Intermolecular Bonding Interactions in van der Waals Molecules

Noncovalent and Covalent O-H···O Interactions in PPh3O Cocrystals: A Correlation Study Involving QTAIM, SAPT, NBO, and IBSI Methods

PPh3O.hemihydrate polymorphs and 11 assorted PPh3O cocrystals collectively constitute a reliable stock to pursue a systematic analysis aiming to investigate the impacts of some vital issues on the TPPO.H-bond donor aggregates. The issues highlighted herein are (i) effect of varying acidity of H-bond donors on the degeneracy of lone pairs of the H-bond acceptor (PPh3O), (ii) effectiveness of the |V(r)|/G(r) and H(r)/ρ(r) parameters as a covalency metric, (iii) 3c-4e bonding in the covalent PPh3O.nitric acid cocrystal, (iv) salient features of H-bond interaction energy and an interplay of its components, (v) an intrinsic bond strength scale for the PPh3O cocrystals, and (vi) reliable empirical relations between several bond descriptors for a quick estimation of interaction energy. To be specific about point (vi), we have propounded two promising avenues for a fast semiquantitative calculation of interaction energy from an endearing nonenergetic parameter, viz., bond length: dO-H···O → ρBCP (MAPE = 2.36%) → ESAPT0 (MAPE = 9.26%), and dO-H···O → IBSI (MAPE = 1.87%) → ESAPT0 (MAPE = 9.66%). All the aforesaid issues have been explored in detail through the QTAIM, NBO, and IBSI analyses (M06-2X-D3/def2-TZVP level), as well as by the SAPT study at the SAPT0/aug-cc-pVDZ platform. The statistically valid correlation studies can be particularly conducive for practical purposes as a transformative extension of the established facts into postulates for the unknown cocrystals.

Bader's Topological Bond Path Does Not Necessarily Indicate Stabilizing Interaction-Proof Studies Based on the Ng@[3n]cyclophane Endohedral Complexes

According to Bader's quantum theory of atoms in molecules (QTAIM), the simultaneous presence of a bond path and the corresponding bond critical point between any two atoms is both a necessary and sufficient condition for the atoms to be bonded to one another. In principle, this means that this pair of atoms should make a stabilizing contribution to the molecular system. However, the multitude of so-called counterintuitive bond paths strongly suggests that this statement is not necessarily true. Particularly 'troublesome' are endohedral complexes, in which encapsulation-enforced proximity between the trapped guest (e.g., an atom) and the host's cage system usually 'produces' many counterintuitive bond paths. In the author's opinion, the best evidence to demonstrate the repulsive nature of the intra-cage guest⋯host interaction is the use of some trapping systems containing small escape channels and then showing that the initially trapped entity spontaneously escapes outside the host's cage during geometry optimization of the initially built guest@host endohedral complex. For this purpose, a group of 24 Ng@[3n]cyclophane (3≤n≤6) endohedral complexes is used. As a result, arguments are presented showing that Bader's topological bond path does not necessarily indicate a stabilizing interaction.

Tools for overcoming reliance on energy-based measures in chemistry: a tutorial review

The vast majority of literature in the chemical sciences describes fundamental chemical and physical phenomena using scalar measures, such as the energy, even though many phenomena are beyond the scope of scalar-based considerations. This problem exists no matter how accurately the associated energies are calculated. The solution that is explained in this work is to remove the reliance on scalar quantum chemical measures and instead utilize the vector-based and full symmetry-breaking nature of next generation quantum theory of atoms in molecules (NG-QTAIM). The connection with experiment on neutral chiral molecules is explained. A selection of non-energy-based explanations are provided: the functioning of molecular devices, why the cis-effect is the exception rather than the rule, stereochemical phenomena including chiral discrimination, quantifying chiral character of formally achiral molecules, mixed S and R stereoisomer character and the effect of an applied electric field. Current and future developments along with suggestions for future avenues of investigation are discussed. This tutorial review provides the practical details required to implement NG-QTAIM for a range of phenomena that are not accessible with energy-based measures. Step-by-step worked examples are included with data sets and instructions for use of commercial and open-source software along with examples of how to interpret the results.

Methylammonium Tetrel Halide Perovskite Ion Pairs and Their Dimers: The Interplay between the Hydrogen-, Pnictogen- and Tetrel-Bonding Interactions

The structural stability of the extensively studied organic-inorganic hybrid methylammonium tetrel halide perovskite semiconductors, MATtX₃ (MA = CH₃NH₃⁺; Tt = Ge, Sn, Pb; X = Cl, Br, I), arises as a result of non-covalent interactions between an organic cation (CH₃NH₃⁺) and an inorganic anion (TtX₃^-). However, the basic understanding of the underlying chemical bonding interactions in these systems that link the ionic moieties together in complex configurations is still limited. In this study, ion pair models constituting the organic and inorganic ions were regarded as the repeating units of periodic crystal systems and density functional theory simulations were performed to elucidate the nature of the non-covalent interactions between them. It is demonstrated that not only the charge-assisted N-H···X and C-H···X hydrogen bonds but also the C-N···X pnictogen bonds interact to stabilize the ion pairs and to define their geometries in the gas phase. Similar interactions are also responsible for the formation of crystalline MATtX₃ in the low-temperature phase, some of which have been delineated in previous studies. In contrast, the Tt···X tetrel bonding interactions, which are hidden as coordinate bonds in the crystals, play a vital role in holding the inorganic anionic moieties (TtX₃^-) together. We have demonstrated that each Tt in each [CH₃NH₃⁺•TtX₃^-] ion pair has the capacity to donate three tetrel (σ-hole) bonds to the halides of three nearest neighbor TtX₃^- units, thus causing the emergence of an infinite array of 3D TtX₆^4- octahedra in the crystalline phase. The TtX₄^4- octahedra are corner-shared to form cage-like inorganic frameworks that host the organic cation, leading to the formation of functional tetrel halide perovskite materials that have outstanding optoelectronic properties in the solid state. We harnessed the results using the quantum theory of atoms in molecules, natural bond orbital, molecular electrostatic surface potential and independent gradient models to validate these conclusions.

Controlling the Diradical Character of Thiele Like Compounds

Organic diradicals play an important role in many fields of chemistry, biochemistry, and materials science. In this work, by means of high-level theoretical calculations, we have investigated the effect of representative chemical substituents in p-quinodimethane (pQDM) and Thiele's hydrocarbons with respect to the singlet-triplet energy gap, a feature characterizing their diradical character. We show how the nature of the substituents has a very important effect in controlling the singlet-triplet energy gap so that several compounds show diradical features in their ground electronic state. Importantly, steric effects appear to play the most determinant role for pQDM analogues, with minor effects of the substituents in the central ring. For Thiele like compounds, we found that electron-withdrawing groups in the central ring favor the quinoidal form with a low or almost null diradical character, whereas electron-donating group substituents favor the aromatic-diradical form if the electron donation does not exceed 6-π electrons. In this case, if there is an excess of electron donation, the diradical character is reduced. The electronic spectrum of these compounds is also calculated, and we predict that the most intense bands occur in the visible region, although in some cases characteristic electronic transition in the near-IR region may appear.

High-resolution structural study on pyridin-3-yl ebselen and its N-methylated tosylate and iodide derivatives

The crystal structure of the pyridine-substituted benzisoselenazolinone 2-(pyridin-3-yl)-2,3-dihydro-1,2-benzoselenazol-3-one (C₁₂H₈N₂OSe, 2), related to the antioxidant ebselen [systematic name: 2-phenyl-1,2-benzoselenazol-3(2H)-one, 1], is characterized by strong intermolecular N...Se(-N) chalcogen bonding, where the N...Se distance of 2.3831 (6) Å is well within the sum of the van der Waals radii for N and Se (3.34 Å). This strong interaction results in significant lengthening of the internal N-Se distance, consistent with significant population of the Se-N σ* antibonding orbital. Much weaker intermolecular O...Se chalcogen bonding occurs between the amide-like O atom in 2 and the less polarized C-Se bond in this structure. Charge density analysis of 2 using multipole refinement of high-resolution data allowed the electrostatic surface potential for 2 to be mapped, and clearly reveals the σ-hole at the extension of the Se-N bond as an area of positive electrostatic potential. Topological analysis of the electron-density distribution in 2 was carried out within the Quantum Theory of Atoms in Molecules (QTAIM) framework and revealed bond paths and (3,-1) bond critical points (BCPs) for the N...Se-N moiety consistent with a closed-shell interaction; however, the potential energy term is suggestive of electron sharing. Analysis of the electron localization function (ELF) for the strong N...Se and the weak O...Se chalcogen-bonding interactions in the structure of 2 suggest significant electron sharing in the former interaction, and a largely electrostatic interaction in the latter. Conversion of 2 to its N-methylated derivatives by reaction with methyl iodide [1-methyl-3-(3-oxo-2,3-dihydro-1,2-benzoselenazol-2-yl)pyridin-1-ium iodide, C₁₃H₁₁N₂OSe⁺·I^-] and methyl tosylate [1-methyl-3-(3-oxo-2,3-dihydro-1,2-benzoselenazol-2-yl)pyridin-1-ium toluenesulfonate trihydrate, C₁₃H₁₁N₂OSe⁺·C₇H₇O₃S^-·3H₂O] removes the possibility of N...Se chalcogen bonding and instead structures are obtained where the iodide and tosylate counter-ions fulfill the role of chalcogen-bond acceptors, with a strong I^-...Se interaction in the iodide salt and a weaker p-Tol-SO₃^-...Se interaction in the tosylate salt.

Hydrogen bond networks of dimethylsulfoxide (DMSO) pentamer

Understanding of clusters of dimethylsulfoxide (DMSO) is important in several applications in Chemistry. Despite its importance, very few studies of DMSO clusters, (DMSO)_n, have been reported in comparison to systems such as water clusters or methanol clusters. In order to provide further understanding of DMSO clusters, we investigated the structures and non-covalent interactions of the (DMSO)_n, n=5. Therefore, the potential energy surface (PES) of the DMSO pentamer has been examined using classical molecular dynamics. The structures generated using classical molecular dynamics are further optimized at the PW6B95D3/aug-cc-pVDZ level of theory. To comprehend the non-covalent bondings in the DMSO pentamer, we carried out a quantum theory of atoms in molecule (QTAIM) analysis. In addition, the effects of temperature on the structural stability is investigated between 20 and 500K. It comes out that seven different kind of non-covalent bondings can be found in DMSO pentamers.

Direction of theoretical and experimental investigation into the mechanism of n-HA/Si-PA-SC@Ag as a bio-based heterogeneous catalyst in the reduction reactions

In the present study, a natural-based heterogeneous catalyst is synthesized. For this purpose, nano-hydroxyapatite (n-HA) is prepared, silica-modified and functionalized with phthalimide. Finally, Ag²⁺ was immobilized onto n-HA/Si-PA-SC and reduced to Ag nanoparticles by Bellis perennis flowers extract. n-HA/Si-PA-SC@Ag characterized by TGA, FTIR, SEM/EDX, XRD, TEM, BET and ICP-AES techniques. Moreover, metal-ligand interactions in n-HA/Si-PA-SC@Ag complex models were assessed to make a quantitative representation for the immobilization behavior of Ag NPs on the surface of n-HA/Si-PA-SC through quantum chemistry computations. Furthermore, the performance of n-HA/Si-PA-SC@Ag was studied in the nitroarene, methylene blue and congo red reductions. Finally, the recyclability study as well as Ag-leaching verified that, n-HA/Si-PA-SC@Ag was stable and reused-up to four times without losing its activity.

Self-Assembling Behaviour of Perylene, Perylene Diimide, and Thionated Perylene Diimide Deciphered through Non-Covalent Interactions

The π-conjugated supramolecular polymers (SMP) have gained vast popularity in materials chemistry and biomedicine due to their spectacular self-assembling behaviour. A detailed account of the electronic structure and bonding through quantum theory of atoms-in-molecules, non-covalent interactions, and energy decomposition analysis (EDA) in the oligomers of perylene, perylene diimide (PDI), and thionated-PDI (t-PDI) is presented. The oligomers of all three molecules show a slip angle of θ≈62° thus forming H-aggregates. The stacking pattern in perylene oligomers prefers a slip-stacked brick-layer order, while the bulkier PDI and t-PDI prefer a parallel step-wise pattern in their oligomers. Successive addition of monomers leads to a consequent rise in the association energy, although to a much greater extent in PDI and t-PDI than in perylene. While the major contribution to this association energy comes from the dispersion interactions in all three systems, the steric interactions in t-PDI quench the cooperativity in its SMP formation. A detailed analysis of the non-covalent interactions reveals the presence of π-π, π-hole⋅⋅⋅O=C, and π-hole⋅⋅⋅S=C electrostatic interactions playing a crucial role in the self-assembly process, which can be further implemented on developing force field-based methods for understanding the self-assembling mechanism in higher degree of oligomers.

Probing the Directionality of S···O/N Chalcogen Bond and Its Interplay with Weak C-H···O/N/S Hydrogen Bond Using Molecular Electrostatic Potential

The directionality of the chalcogen bond (Ch-bond) formed by S and its interplay with other weak interactions have important chemical and biological implications. Here, dimers made of CH₃-S-X and O/N containing nucleophiles are studied and found to be stabilized by coexisting S···O/N and C-H···O/N interactions. Based on experimentally accessible electron density and molecular electrostatic potentials (MESPs), we showed that reciprocity between S···O/N and C-H···O/N interactions in the stability of cumulative molecular interaction (ΔE) was dependent on the strength of the σ-hole on S (V_s,max). Direct correlation between ΔE of dimers with V_s,max of S supports the electrostatic nature of the Ch-bond. Such interplay of the Ch-bond is necessary for its directionality in complex nucleophiles (carbonyl groups) with multiple electron-rich centers, which is explained using MESP. A correlation between the MESP minima in the π-region and the strength of the S-π interaction explains the directional selectivity of the Ch-bond.

The Cis-Effect Explained Using Next-Generation QTAIM

We used next-generation QTAIM (NG-QTAIM) to explain the cis-effect for two families of molecules: C₂X₂ (X = H, F, Cl) and N₂X₂ (X = H, F, Cl). We explained why the cis-effect is the exception rather than the rule. This was undertaken by tracking the motion of the bond critical point (BCP) of the stress tensor trajectories T_σ(s) used to sample the U_σ-space cis- and trans-characteristics. The T_σ(s) were constructed by subjecting the C1-C2 BCP and N1-N2 BCP to torsions ± θ and summing all possible T_σ(s) from the bonding environment. During this process, care was taken to fully account for multi-reference effects. We associated bond-bending and bond-twisting components of the T_σ(s) with cis- and trans-characteristics, respectively, based on the relative ease of motion of the electronic charge density ρ(r_b). Qualitative agreement is found with existing experimental data and predictions are made where experimental data is not available.

Determining Repulsion in Cyclophane Cages

Superphane, i.e., [2.2.2.2.2.2](1,2,3,4,5,6)cyclophane, is a very convenient molecule in studying the nature of guest⋯host interactions in endohedral complexes. Nevertheless, the presence of as many as six ethylene bridges in the superphane molecule makes it practically impossible for the trapped entity to escape out of the superphane cage. Thus, in this article, I have implemented the idea of using the superphane derivatives with a reduced number of ethylene linkers, which leads to the [2n] cyclophanes where n<6. Seven such cyclophanes are then allowed to form endohedral complexes with noble gas (Ng) atoms (He, Ne, Ar, Kr). It is shown that in the vast majority of cases, the initially trapped Ng atom spontaneously escapes from the cyclophane cage, creating an exohedral complex. This is the best proof that the Ng⋯cyclophane interaction in endohedral complexes is indeed highly repulsive, i.e., destabilizing. Apart from the ‘sealed’ superphane molecule, endohedral complexes are only formed in the case of the smallest He atom. However, it has been shown that in these cases, the Ng⋯cyclophane interaction inside the cyclophane cage is nonbonding, i.e., repulsive. This highly energetically unfavorable effect causes the cyclophane molecule to ‘swell’.

Rational Computational Design of Systems Exhibiting Strong Halogen Bonding Involving Fluorine in Bicyclic Diamine Derivatives

Perhaps the most controversial and rare aspect of the halogen bonding interaction is the potential of fluorine in compounds to serve as a halogen bond donor. In this note, we provide clear and convincing examples of hypothetical molecules in which fluorine is strongly halogen bonding in a metastable state. Of particular note is a polycyclic system inspired by Selectfluor, which has been controversially proposed to engage in halogen bonding.

Synthesis, X-ray, Hirshfeld, and AIM Studies on Zn(II) and Cd(II) Complexes with Pyridine Ligands

The synthesis and crystal structures of three heteroleptic complexes of Zn(II) and Cd(II) with pyridine ligands (ethyl nicotinate (EtNic), N,N-diethylnicotinamide (DiEtNA), and 2-amino-5-picoline (2Ampic) are presented. The complex [Zn(EtNic)2Cl2] (1) showed a distorted tetrahedral coordination geometry with two EtNic ligand units and two chloride ions as monodentate ligands. Complexes [Zn(DiEtNA)(H2O)4(SO4)]·H2O (2) and [Cd(OAc)2(2Ampic)2] (3) had hexa-coordinated Zn(II) and Cd(II) centers. In the former, the Zn(II) was coordinated with three different monodentate ligands, which were DiEtNA, H2O, and SO42−. In 3, the Cd(II) ion was coordinated with two bidentate acetate ions and two monodentate 2Ampic ligand units. The supramolecular structures of the three complexes were elucidated using Hirshfeld analysis. In 1, the most important interactions that governed the molecular packing were O···H (15.5–15.6%), Cl···H (13.6–13.8%), Cl···C (6.3%), and C···H (10.3–10.6%) contacts. For complexes 2 and 3, the H···H, O···H, and C···H contacts dominated. Their percentages were 50.2%, 41.2%, and 7.1%, respectively, for 2 and 57.1%, 19.6%, and 15.2%, respectively, for 3. Only in complex 3, weak π-π stacking interactions between the stacked pyridines were found. The Zn(II) natural charges were calculated using the DFT method to be 0.8775, 1.0559, and 1.2193 for complexes 1–3, respectively. A predominant closed-shell character for the Zn–Cl, Zn–N, Zn–O, Cd–O, and Cd–N bonds was also concluded from an atoms in molecules (AIM) study.

Dimethylsulfoxide (DMSO) Clusters Dataset: DFT Relative Energies, Non-Covalent Interactions, and Cartesian Coordinates

Theoretical understanding of dimethylsulfoxide (DMSO) liquid depends on the understanding of the DMSO clusters. In this work, we provide the structures and the energetics of the DMSO clusters. The structures have been generated using ABCluster and further optimized at the MP2/aug-cc-pVDZ level of theory. The final structures have been optimized at two different levels of theory: PW6B95D3/aug-cc-pVDZ and ωB97XD/aug-cc-pVDZ. The Cartesian coordinates of the structures optimized at the MP2/aug-cc-pVDZ level of theory are also reported. The relative energies of the structures can be used to locate the most favorable structures of the DMSO clusters. The Cartesian coordinates of the structures can be used for further investigations on DMSO clusters. In addition, we report the data related to the quantum theory of atoms in molecule (QTAIM) analysis of the investigated clusters. The QTAIM data reported in this work can be used to understand and determine the nature of non-covalent interactions in DMSO clusters. For further reading and discussion on the data reported here, please report to the original manuscript Malloum and Conradie (2022) [1].

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).

Next generation quantum theory of atoms in molecules for the design of emitters exhibiting thermally activated delayed fluorescence with laser irradiation

The effect of a static electric (E)-field and an unchirped and chirped laser pulse field on the cycl[3.3.3]azine molecule was investigated using next-generation quantum theory of atoms in molecules (NG-QTAIM). Despite the magnitude of the E-field of the laser pulses being an order of magnitude lower than for the static E-field, the variation of the energy gap between the lowest lying singlet (S₁ ) and triplet (T₁ ) excited states was orders of magnitude greater for the laser pulse than for the static E-field. Insights into the response of the electronic structure were captured by NG-QTAIM, where differences in the inverted singlet-triplet gap due to the laser pulses were significant larger compared to those induced by the static E-field. The response of the S₁ and T₁ excited states, as determined by NG-QTAIM, switched discontinuously between weak and strong chemical character for the static E-field. In contrast, the response to the laser pulses, determined by NG-QTAIM, is to induce a continuous range of chemical character, indicating the unique ability of the laser pulses to induce polarization effects in the form of "mixed" bond types. Our analysis demonstrates that NG-QTAIM is a useful tool for understanding the response to laser irradiation of the lowest-lying singlet S₁ and triplet T₁ excited states of emitters exhibiting thermally activated delayed fluorescence. The chirped laser pulse led to more frequent instances of the desired outcome of an inverted singlet-triplet gap than the unchirped pulse, indicating its usefulness as a tool to design more efficient organic light-emitting diode devices.

Three types of noncovalent interactions studied between pyrazine and XF

Three types noncovalent interactions (type I, II and III) between pyrazine (C₄H₄N₂) and XF (X = F, Cl, Br, and I) have been discovered at the MP2/aug-cc-pVTZ level. TypeI is σ-hole interaction between the positive site on the halogen X of XF and the negative site on one of the pyrazine nitrogens. Type II is counterintuitive σ-hole interaction driven by polarization between the positive site on the halogen X of XF and a portion of the pyrazine ring. Type III is an interaction between the lateral regions of the halogen X of XF and the position of the pyrazine ring. Through comparing the calculated interaction energy, we can know that the type II and type III interactions are weaker than the corresponding type I interactions, and type III interactions are weaker than the corresponding type II interactions in C₄H₄N₂-XF complexes. SAPT analysis shows that the electrostatic energy are the major source of the attraction for the type I (σ-hole) interactions while the type III interactions are mainly dispersion energy. For the type II (counterintuitive σ-hole) interactions in C₄H₄N₂-XF (X = F and Cl) complexes, electrostatic energy are the major source of the attraction, while in C₄H₄N₂-XF (X = Br and I) complexes, the electrostatic term, induction and dispersion play equally important role in the total attractive interaction. NBO analysis, AIM theory, and conceptual DFT are also being utilized.

Structures, binding energies and non-covalent interactions of furan clusters

Understanding of the furan solvent is subjected to the knowledge of the structures of the furan clusters and interactions taking place therein. Although, furan clusters can be very important to determine the dynamics and the properties of the furan solvent, there has been only a few investigations reported on furan dimer. In this work, we have explored the potential energy surfaces (PESs) of the furan clusters using two incremental levels of theory. Structures have been initially generated using classical molecular dynamics followed by full optimization at the MP2/aug-cc-pVDZ level of theory. The results show that the most stable structure of the furan dimer has a stacking configuration while that of the trimer has a cyclic configuration. We have noted that the structures of the furan tetramer have no definite configurations. In addition, we have performed a quantum theory of atoms in molecule (QTAIM) analysis to identify all possible non-covalent interactions of the furan clusters. The results show that six different types of non-covalent interactions can be identified in furan clusters. We have noted that the CH⋯C and CH⋯O hydrogen bondings are the strongest non-covalent interactions while the H⋯H bonding interaction is found to be the weakest. Furthermore, we have assessed the performance of ten DFT functionals in calculating the binding energies of the furan clusters. The ten DFT functionals (M05, M05-2X, M06, M06-2X, M08HX, PBE0, ωB97XD, PW6B95D3, APFD and MN15) have been benchmarked to DLPNO-CCSD(T)/CBS. The functionals M05-2X and M06 are recommended for further affordable investigations of the furan clusters.

New insights of QTAIM and stress tensor to finding non-competitive/competitive torquoselectivity of cyclobutene

This study aims to investigate the phenomenon of torquoselectivity through three thermal cyclobutene ring-opening reactions (N1-N3). This research focuses on the nature of the chemical bond, electronic reorganization, predicting non-competitive or competitive reactions, and torquoselectivity preference within Quantum Theory of Atoms in Molecules (QTAIM) and stress tensor frameworks. Various theoretical analyses for these reactions, such as metallicity ξ(r_b), ellipticity ε, total local energy density H(r_b), stress tensor polarizability ℙ_σ, stress tensor eigenvalue λ_3σ, and bond-path length, display differently for non-competitive and competitive reactions as well as for the conrotatory preferences either it is the transition state outward conrotatory (TS_out) or transition state inward conrotatory (TS_in) directions by presenting degeneracy or non-degeneracy in their results. The ellipticity profile provides the motion of the bond critical point locations due to the different substituents of cyclobutene. In agreement with experimental results, examinations demonstrated that N1 is a competitive reaction and N2-N3 are non-competitive reactions with TS_out and TS_in preference directions, respectively. The concordant results of QTAIM and stress tensor scalar and vectors with experimental results provide a better understanding of reaction mechanisms.

Does the Presence of a Bond Path Really Mean Interatomic Stabilization? The Case of the Ng@Superphane (Ng = He, Ne, Ar, and Kr) Endohedral Complexes

Using a fairly structurally flexible and, therefore, very suitable for this type of research, superphane molecule, we demonstrate that the inclusion of a noble gas atom (Ng = He, Ne, Ar, and Kr) inside it and, thus, the formation of the Ng@superphane endohedral complex, leads to its ‘swelling’. Positive values of both the binding and strain energies prove that encapsulation and in turn ‘swelling’ of the superphane molecule is energetically unfavorable and that the Ng⋯C interactions in the interior of the cage are destabilizing, i.e., repulsive. Additionally, negative Mayer Bond Orders indicate the antibonding nature of Ng⋯C contacts. This result in combination with the observed Ng⋯C bond paths shows that the presence of a bond path in the molecular graph does not necessarily prove interatomic stabilization. It is shown that the obtained conclusions do not depend on the computational methodology, i.e., the method and the basis set used. However, on the contrary, the number of bond paths may depend on the methodology. This is yet another disadvantageous finding that does not favor the treatment of bond paths on molecular graphs as indicators of chemical bonds. The Kr@superphane endohedral complex features one of the longest C–C bonds ever reported (1.753 Å).

Comparative study of molecular docking, structural, electronic, vibrational and hydrogen bonding interactions on 4-hydroxy benzo hydrazide (4HBH) and its newly designed derivative [(E)–N′-((1H-Pyrrol-2-YL)methylene) –4-hydroxy benzo hydrazide and its isomers (I, II and III)] (potential inhibitors) for COVID-19 protease

Studies have shown that hydrazides and thier derivatives are used for pharmaceutical and medicinal purposes. At present, the whole world is suffering for COVID-19 virus. There are some vaccines or medicines available to treat this disease all over the world. Today the one fourth of the world’s population is under lockdown condition. In this scenario, scientists from the whole world are doing different types of research on this disease. Being a molecular modeller, this inspires us to design new types of species (may be drugs) which may be capable for COVID-19 Protease. In the present effort, we have performed docking studies of title compounds with COVID-19 protein (6LU7) for anti-COVID-19 activity. A comparative quantum chemical calculations of molecular geometries (bond lengths and bond angles) of 4-Hydroxy Benzo Hydrazide (4HBH) and its newly designed derivatve [(E)-N′-((1H-Pyrrol-2-YL)Methylene) –4-Hydroxy Benzo Hydrazide and its isomers (I, II and III)] in the ground state have also been carried out due to its biological importance and compared with the similer type of compound found in literature i.e. benzohydrazide. The optimized geometry and wavenumber of the vibrational bands of the molecules have been calculated by density functional theory (DFT) using Becke’s three-parameters hybrid functional (B3LYP/CAM-B3LYP) with 6–311G (d, p) as the basis set. Vibrational wavenumbers are compared with the observed FT-IR and FTRaman spectra of 4-Hydroxy Benzo Hydrazide. TDDFT calculations are also done on the same level of theory and a theoretical UV-vis spectrum of title molecules are also drawn. HOMO-LUMO analysis has been done to describe the way the molecule interacts with other species. Natural bond orbitals (NBO) analysis has been carried out to inspect the intra- and inter- molecular hydrogen-bonding, conjugative and hyper conjugative interactions and their second order stabilization energy. Nonlinear optical (NLO) analysis has been performed to study the non-linear optical properties of the molecule by computing the first hyperpolarizability. The variation of thermodynamic properties with temperature has been studied. QATIM analysis shows that hydrogen bonding occurs in 4HBH, isomer II and III respectively.

Chirality without Stereoisomers: Insight from the Helical Response of Bond Electrons

The association between molecular chirality and helical characteristics known as the chirality-helicity equivalence is determined for the first time by quantifying a chirality-helicity measure consistent with photoexcitation circular dichroism experiments. This is demonstrated using a formally achiral S_N 2 reaction and a chiral S_N 2 reaction. Both the achiral and chiral S_N 2 reactions possess significant values of the chirality-helicity measure and display a Walden inversion, i. e. an inversion of the chirality between the reactant and product. We also track the chirality-helicity measure along the reaction path and discover the presence of chirality at the transition state and two intermediate structures for both reactions. We demonstrate the need for the chirality-helicity measure to differentiate between steric effects due to eclipsed conformations and chiral behaviors in formally achiral species. We explain the significance of this work for asymmetric synthetic reactions including the intermediate structures where the Cahn-Ingold-Prelog (CIP) rules cannot be used.

Theoretical study on the noncovalent interactions involving triplet diphenylcarbene

The properties of some types of noncovalent interactions formed by triplet diphenylcarbene (DPC³) have been investigated by means of density functional theory (DFT) calculations and quantum theory of atoms in molecule (QTAIM) studies. The DPC³···LA (LA = AlF₃, SiF₄, PF₅, SF₂, ClF) complexes have been analyzed from their equilibrium geometries, binding energies, and properties of electron density. The triel bond in the DPC³···AlF₃ complex exhibits a partially covalent nature, with the binding energy - 65.7 kJ/mol. The tetrel bond, pnicogen bond, chalcogen bond, and halogen bond in the DPC³···LA (LA = SiF₄, PF₅, SF₂, ClF) complexes show the character of a weak closed-shell noncovalent interaction. Polarization plays an important role in the formation of the studied complexes. The strength of intermolecular interaction decreases in the order LA = AlF₃ > ClF > SF₂ > SiF₄ > PF₅. The electron spin density transfers from the radical DPC³ to ClF and SF₂ in the formation of halogen bond and chalcogen bond, but for the DPC³···AlF₃/SiF₄/PF₅ complexes, the transfer of electron spin density is minimal.

Halogen-Halogen Nonbonded Interactions

Halogen-halogen nonbonded interactions were studied for methyl halides and phenyl halides using both B3LYP and MP2 along with 6-311+G* and aug-cc-pVTZ. With the methyl halides, the linear approach was found to lead to little stabilization, whereas the "90°" approach gave 1-2 kcal/mol. This modest stabilization was due to long-range electron correlation effects. The lowest-energy arrangement had the molecules side-by-side, with the major stabilization being derived from halogen-hydrogen interactions. The results for methyl bromide were quite similar. Chlorobenzene dimer with the 90° orientation gave a small stabilization energy, but the best arrangement had the two benzene rings oriented over each other. The meta orientation of the chlorines had a lower energy than ortho or para. The dimerization energy was larger than that for two benzene rings sitting directly above each other, suggesting that whereas Cl···Cl interaction is not very important, the effect of the halogen on the electron distribution does have an effect. This suggests that much of the crystallographic results for these compounds may not be due to halogen-halogen interactions but rather the interaction between the substituted benzene rings along with crystal forces.

The Crystal Structure and Intermolecular Interactions in Fenamic Acids-Acridine Complexes

In order to improve pharmaceutical properties of drugs, complexes are synthesized as combinations with other chemical substances. The complexes of fenamic acid and its derivatives, such as mefenamic-, tolfenamic- and flufenamic acid, with acridine were obtained and the X-ray structures were discussed. Formation of the crystals is determined by the presence of the intermolecular O–H^…N hydrogen bond that occur between fenamic acids and acridine. Intermolecular interactions stabilizing the crystals such as π^…π stacking, C–H^…X (X = O, Cl) intermolecular hydrogen bonds as well as C–H^…π and other dispersive interactions were analyzed by theoretical methods: the quantum theory of atoms in molecules (QTAIM) and noncovalent interaction (NCI) approaches.

Orbital-Free Quantum Crystallographic View on Noncovalent Bonding: Insights into Hydrogen Bonds, π⋅⋅⋅π and Reverse Electron Lone Pair⋅⋅⋅π Interactions

A detailed analysis of a complete set of the local potentials that appear in the Euler equation for electron density is carried out for noncovalent interactions in the crystal of a uracil derivative using experimental X-ray charge density. The interplay between the quantum theory of atoms in molecules and crystals and the local potentials and corresponding inner-crystal electronic forces of electrostatic and kinetic origin is explored. Partitioning of crystal space into atomic basins and atomic-like potential basins led us to the definite description of interatomic interaction and charge transfer. Novel physically grounded bonding descriptors derived within the orbital-free quantum crystallography provided the detailed examination of π-stacking and intricate C=O⋅⋅⋅π interactions and nonclassical hydrogen bonds present in the crystal. The donor-acceptor character of these interactions is revealed by analysis of Pauli and von Weizsäcker potentials together with well-known functions, e. g., deformation electron density and electron localization function. In this way, our analysis throws light on aspects of these closed-shell interactions hitherto hidden from the description.

The Structures, Molecular Orbital Properties and Vibrational Spectra of the Homo- and Heterodimers of Sulphur Dioxide and Ozone. An Ab Initio Study

The structures of a number of dimers of sulphur dioxide and ozone were optimized by means of a series of ab initio calculations. The dimer species were classified as either genuine energy minima or transition states of first or higher order, and the most probable structures consistent with the experimental data were confirmed. The molecular orbitals engaged in the interactions resulting in adduct formation were identified and relations between the orbitals of the dimers of the valence isoelectronic monomer species were examined. The vibrational spectra of the most probable structures were computed and compared with those reported in the literature, particularly with spectra observed in cryogenic matrices. The calculations were extended to predict the properties of a number of possible heterodimers formed between sulphur dioxide and ozone.

Predicted structure and selectivity of 3d transition metal complexes with glutamic N,N-bis(carboxymethyl) acid

Structure and selectivity of 3d transition metal complexes with glutamic N,N-bis(carboxymethyl) acid are analyzed and predicted from DFT calculations.

What Is the Main Feature Distinguishing the Through-Space Interactions in Cyclophanes from Their Aliphatic Analogues?

Classical cyclophanes with two benzene rings have been compared with cyclophanes with one benzene ring replaced with an aliphatic part and aliphatic compounds, which are cyclophane analogues. Analysis of geometry, atomic charges, and aromatic and steric energy and investigation of intramolecular noncovalent interactions and charge mobility show that there is no special feature that distinguishes the classical cyclophanes from aliphatic analogues, so the definition of cyclophanes can be extended to other compounds.

Interatomic exchange-correlation interaction energy from a measure of quantum theory of atoms in molecules topological bonding: A diatomic case

The potential relations between the measure of topological interatomic bonding-integrals of electron density with respect to internuclear axis over the corresponding quantum theory of atoms in molecules (QTAIM)-defined interatomic surface (IAS)-and interatomic exchange-correlation contributions from the interacting quantum atoms approach are discussed. The quantum chemical computations of 38 equilibrium diatomic systems at different levels of theory (HF, MP2, MP4SDQ, and CCSD) are invoked to support abstract considerations. Parameters of excellent correlations between IAS integrals and interatomic exchange-correlation energy are found by the optimization. The performance of these trends depends on the accuracy of the electronic correlation treatment. The resulting trends are a unique feature of equilibrium states, whereas more complicated dependencies are explored for several systems at non-equilibrium conditions. The relations of established trends with other IAS-based estimations of strength of bonding interactions between topological atoms and issues explored for multiatomic systems are briefly discussed.

Diversity and uniformity in anion-π complexes of thiocyanate with aromatic, olefinic and quinoidal π-acceptors

Despite the progress in the study of anion-π interactions, there are still inconsistencies in the use of this term and the experimental data about factors affecting the strength of such bonding are limited. To shed light on these issues, we explored supramolecular associations between NCS^- anions and a series of aromatic, olefinic or quinoidal π-acceptors. Combined experimental and computational studies revealed that all these complexes were formed by an attraction of the anion to the face of the π-system, and the arrangements of thiocyanate followed the areas of the most positive potentials on the surfaces of the π-acceptors. The stabilities of the complexes increased with the π-acceptor strength (reflected by their reduction potentials), and were essentially independent of the magnitudes of the maximum electrostatic potentials on their surfaces. The complexes showed intense absorption bands in the UV-Vis range, and the energies of these bands were correlated with the difference of the redox potentials of the anions and π-acceptors. Such features, as well as results of atoms-in-molecules and non-covalent index analyses suggested that besides electrostatics, molecular orbital interactions play a substantial role in the formation of these complexes. The unified trends in variations of the characteristics of the complexes between thiocyanate and a variety of π-acceptors point to their common nature. To embrace diversity and uniformity of the anion-π associates, we suggest (following the halogen bond's definition) that anion-π bonding occurs when there is evidence of a net attraction between the anions and the face of the electrophilic π-system.