Characterization of a Closed-Shell Fluorine−Fluorine Bonding Interaction in Aromatic Compounds on the Basis of the Electron Density

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.

Dispersion, Rehybridization, and Pentacoordination: Keys to Understand Clustering of Boron and Aluminum Hydrides and Halides

The structure, stability, and bonding characteristics of dimers and trimers involving BX₃ and AlX₃ (X = H, F, Cl) in the gas phase, many of them explored for the first time, were investigated using different DFT (B3LYP, B3LYP/D3BJ, and M06-2X) and ab initio (MP2 and G4) methods together with different energy decomposition formalisms, namely, many-body interaction-energy and localized molecular orbital energy decomposition analysis. The electron density of the clusters investigated was analyzed with QTAIM, electron localization function, NCIPLOT, and adaptive natural density partitioning approaches. Our results for triel hydride dimers and Al₂X₆ (X = F, Cl) clusters are in good agreement with previous studies in the literature, but in contrast with the general accepted idea that B₂F₆ and B₂Cl₆ do not exist, we have found that they are predicted to be weakly bound systems if dispersion interactions are conveniently accounted for in the theoretical schemes used. Dispersion interactions are also dominant in both homo- and heterotrimers involving boron halide monomers. Surprisingly, B₃F₉ and B₃Cl₉ C_3v cyclic trimers, in spite of exhibiting rather strong B-X (X = F, Cl) interactions, were found to be unstable with respect to the isolated monomers due to the high energetic cost of the rehybridization of the B atom, which is larger than the two- and three-body stabilization contributions when the cyclic is formed. Another important feature is the enhanced stability of both homo- and heterotrimers in which Al is the central atom because Al is systematically pentacoordinated, whereas this is not the case when the central atom is B, which is only tri- or tetra-coordinated.

Synthesis and Structure Elucidation of Novel Spirooxindole Linked to Ferrocene and Triazole Systems via [3 + 2] Cycloaddition Reaction

In the present work, a novel heterocyclic hybrid of a spirooxindole system was synthesized via the attachment of ferrocene and triazole motifs into an azomethine ylide by [3 + 2] cycloaddition reaction protocol. The X-ray structure of the heterocyclic hybrid (1″R,2″S,3R)-2″-(1-(3-chloro-4-fluorophenyl)-5-methyl-1H-1,2,3-triazole-4-carbonyl)-5-methyl-1″-(ferrocin-2-yl)-1″,2″,5″,6″,7″,7a″-hexahydrospiro[indoline-3,3″-pyrrolizin]-2-one revealed very well the expected structure, by using different analytical tools (FTIR and NMR spectroscopy). It crystallized in the triclinic-crystal system and the P-1-space group. The unit cell parameters are a = 9.1442(2) Å, b = 12.0872(3) Å, c = 14.1223(4) Å, α = 102.1700(10)°, β = 97.4190(10)°, γ = 99.1600(10)°, and V = 1484.81(7) ų. There are two molecules per unit cell and one formula unit per asymmetric unit. Hirshfeld analysis was used to study the molecular packing of the heterocyclic hybrid. H···H (50.8%), H···C (14.2%), Cl···H (8.9%), O···H (7.3%), and N···H (5.1%) are the most dominant intermolecular contacts in the crystal structure. O···H, N···H, H···C, F···H, F···C, and O···O are the only contacts that have the characteristic features of short and significant interactions. AIM study indicated predominant covalent characters for the Fe–C interactions. Also, the electron density (ρ(r)) at the bond critical point correlated inversely with the Fe–C distances.

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

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.

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

Molecular, Supramolecular Structures Combined with Hirshfeld and DFT Studies of Centrosymmetric M(II)-azido {M=Ni(II), Fe(II) or Zn(II)} Complexes of 4-Benzoylpyridine

The supramolecular structures of the three metal (II) azido complexes [Fe(4bzpy)4(N3)2]; 1, [Ni(4bzpy)4(N3)2]; 2 and [Zn(4bzpy)2(N3)2]n; 3 with 4-benzoylpyridine (4bzpy) were presented. All complexes contain hexa-coordinated divalent metal ions with a slightly distorted octahedral MN6 coordination sphere. Complexes 1 and 2 are monomeric with terminal azido groups while 3 is one-dimensional coordination polymer containing azido groups with μ(1,1) and μ(1,3) bridging modes of bonding. Hirshfeld analysis was used to quantitatively determine the different contacts affecting the molecular packing in the studied complexes. The most common interactions are the polar O…H and N…H interactions and the hydrophobic C…H contacts. The charges at the M(II) sites are calculated to be 1.004, 0.847, and 1.147 e for complexes 1–3, respectively. The degree of asymmetry is the highest in the case of the terminal azide in complexes 1 and 2 while was found the lowest in the μ(1,1) and μ(1,3) azide bonding modes in the Zn(II) complex 3. These facts were further explained in terms of atoms in molecules (AIM) topological parameters.

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.

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.

Improving phase-transfer catalysis by enhancing non-covalent interactions

A theoretical study of the interactions established between an alkaloid quinine-derived PTC and different anions of interest was performed. Ion pairing competes with an intermolecular hydrogen bond between the PT counteranion and potential HB donors.

Existence of untypical halogen-involving interactions in crystal packings: a statistical and first-principles study

There is a common perception by the scientific community that a halogen-involving interaction forms when the distance between the donor atom and the acceptor atom is less than the sum of their van der Waals (vdW) radii.

Interplay of π-stacking and inter-stacking interactions in two-component crystals of neutral closed-shell aromatic compounds: periodic DFT study

This paper bridges the gap between high-level ab initio computations of gas-phase models of 1 : 1 arene–arene complexes and calculations of the two-component (binary) organic crystals using atom–atom potentials. The studied crystals consist of electron-rich and electron-deficient compounds, which form infinite stacks (columns) of heterodimers. The sublimation enthalpy of crystals has been evaluated by DFT periodic calculations, while intermolecular interactions have been characterized by Bader analysis of the periodic electronic density. The consideration of aromatic compounds without a dipole moment makes it possible to reveal the contribution of quadrupole–quadrupole interactions to the π-stacking energy. These interactions are significant for heterodimers formed by arenes with more than 2 rings, with absolute values of the traceless quadrupole moment (Q_zz) larger than 10 D Å. The further aggregation of neighboring stacks is due to the C–H⋯F interactions in arene/perfluoroarene crystals. In crystals consisting of arene and an electron-deficient compound such as pyromellitic dianhydride, aggregation occurs due to the C–H⋯O interactions. The C–H⋯F and C–H⋯O inter-stacking interactions make the main contribution to the sublimation enthalpy, which exceeds 150 kJ mol⁻¹ for the two-component crystals formed by arenes with more than 2 rings.

The interplay of π-stacking and inter-stacking interactions in two-component organic crystals without conventional hydrogen bonds.

Synthesis, X-ray Crystal Structure and Antimicrobial Activity of Unexpected Trinuclear Cu(II) Complex from s-Triazine-Based Di-Compartmental Ligand via Self-Assembly

The synthesis and X-ray crystal structure of the trinuclear [Cu3(HL)(Cl)2(NO3)(H2O)5](NO3)2 complex of the s-triazine-based di-compartmental ligand, 2-methoxy-4,6-bis(2-(pyridin-2-ylmsethylene)hydrazinyl)-1,3,5-triazine (H2L), are presented. The Cu1 and Cu2 are penta-coordinated with CuN3ClO coordination environment, distorted square pyramidal coordination geometry while Cu3 is hexa-coordinated with CuN2O4 coordination sphere, and distorted octahedral geometry. The complex crystallized in the primitive P-1 triclinic crystal system with two molecular units per unit cell. Its packing is dominated by the O–H (35.5%) and Cl–H (8.8%) hydrogen bonding interactions as well as the π–π stacking (2.3%) and anion–π-stacking interactions (3.7%). The different coordination interactions were analyzed using atoms in molecules (AIM) theory, and the number of charge transferences from the ligand group to Cu(II) were determined using natural bond orbital calculations. The effect of the free ligand and its Cu(II) complex on the tested pathogenic microbes (Staphylococcus aureus, S. epidermidis, Enterococcus faecalis, Escherichia coli, Salmonella typhi and Pseudomonas aeruginosa) and one fungal isolate (Candida albicans) is presented. Both have wide spectrum antimicrobial activity against the selected microorganism. It is observed that the free ligand at 180 µg/mL was more effective than its Cu(II) complex and showed close results compared to the positive control gentamicin. At higher concentrations (1 mg/mL), the Cu(II) complex was found to be more active against S. epidermidis, E. coli and C. albicans than the lower concentration. The minimum inhibitory concentration (MIC) and minimum bactericidal concentration (MBC) values are also lower for the Cu(II) complex than the free ligand.

Intramolecular H-Bond Is Formed in 2-Fluorophenol and 2-Fluorothiophenol, but It May Not Be the Main Pathway of the JFH Coupling Constant Transmission

The intramolecular CF···HX (X = O or S) H-bond and J_FH spin-spin coupling constants (SSCCs) in 2-fluorophenol and 2-fluorothiophenol were investigated experimentally by ¹H and ¹⁹F NMR and theoretically in the framework of the natural bond orbital analysis. In contrast with recent findings from the literature, the results obtained in this work showed that an intramolecular H-bond is formed in the cis conformers and that it has an electrostatic origin. Such an intramolecular electrostatic H-bond is the interaction that rules the conformational preferences on these compounds. Moreover, the natural J-coupling analysis indicated that the J_HF SSCC in 2-fluorophenol has its origin on LP(F)/σ_OH steric interactions and should be labeled as ^4TSJ_HF. On the other hand, the analogue SSCC for 2-fluorothiophenol has its origin on the intramolecular H-bond LP(F) → σ_SH^* hyperconjugative interaction and should be described as a ^1hJ_FH SSCC.

Synthesis, Molecular and Supramolecular Structures of New Cd(II) Pincer-Type Complexes with s-TriazineCore Ligand

The manuscript described the synthesis and characterization of the new [Cd(BDMPT)2](ClO4)2; 1 and [Cd2(MBPT)2(H2O)2Cl](ClO4)3.4H2O ; 2s-triazine pincer-type complexes, where BDMPT and MBPT are 2,4-bis(3,5-dimethyl-1H-pyrazol-1-yl)-6-methoxy-1,3,5-triazine and 2-methoxy-4,6-bis(2-(pyridin-2-ylmsethylene)hydrazinyl)-1,3,5-triazinerespectively.The synthesized complexes were characterized using Fourier-transform infrared spectroscopy (FTIR), 1H and 13C NMR spectroscopy, and the single-crystal X-ray diffraction technique.The homoleptic mononuclear complex (1)contains a hexa-coordinated Cd(II) center with two tridentate N-pincer ligand (BDMPT) with a highly distorted octahedral coordination environment located as an intermediate case between the octahedron and trigonal prism. The heteroleptic dinuclear complex (2) contains two hepta-coordinated Cd(II) coordination spheres where each Cd(II) is coordinated with one pentadentate pincer N-chelate (MBPT), one water, and one bridged chloride ligand connecting the two metal ions. The different intermolecular interactions in the studied complexes were quantified using Hirshfeld analysis. Their thermal stabilities and FTIR spectra were compared with the corresponding free ligands. The strength and nature of Cd–N, Cd–O, and Cd–Cl coordination interactions were discussed in light of atoms in molecules calculations (AIM). The M(II)–BDMPT and M(II)–MBPT interaction energies revealed that such sterically hindered ligands have higher affinity toward large-size metal ions (M =Cd) compared to smaller ones (M= Ni or Mn).

Nature of halogen-centered intermolecular interactions in crystal growth and design: Fluorine-centered interactions in dimers in crystalline hexafluoropropylene as a prototype

The wide occurrence of halogen-centered noncovalent interactions in crystal growth and design prompted this study, which includes a mini review of recent advances in the field. Particular emphasis is placed on providing compelling theoretical evidence of the formation of these interactions between sites of positive electrostatic potential, as well as between sites of negative electrostatic potential, localized on the electrostatic surfaces of the bound fluorine atoms in a prototypical system, hexafluoropropylene (C₃ F₆ ), upon its interaction with another same molecule to form (C₃ F₆ )₂ dimers. The existence of σ- and π-hole interactions is shown for the stable dimers. Even so, weakly bound interactions locally responsible in holding the molecular fragments together cannot and should not be overlooked since they are partly responsible for determining the overall geometry of the crystal. The results of combined quantum theory of atoms in molecules, molecular electrostatic surface potential, and reduced density gradient noncovalent interaction analyses showed that these latter interactions do indeed play a role in the stability and growth of crystalline C₃ F₆ itself and the (C₃ F₆ )₂ dimers. A symmetry adapted perturbation theory energy decomposition analysis leads to the conclusion that a great majority of the (C₃ F₆ )₂ dimers examined are the consequence of dispersion (and electrostatics), with nonnegligible contribution from polarization, which together competes with an exchange repulsion component to determine the equilibrium geometries. In a few structures of the (C₃ F₆ )₂ dimer, the fluorine is found to serve as a six-center five-bond donor/acceptor, as found for carbon in other systems (Malischewski and Seppelt, Angew. Chem. Int. Ed. 2017, 56, 368). © 2019 Wiley Periodicals, Inc.

On the Uselessness of Bond Paths Linking Distant Atoms and on the Violation of the Concept of Privileged Exchange Channels

We refer to frequently used determinants suggesting dominant interactions between distant atoms in various dimers. First of all, we show, against the still-prevailling opinion, that, in general, bond paths have nothing in common with dominant intermolecular interactions and therefore they are useless in such cases. Quite the contrary, reliable information about dominant intermolecular interactions can be obtained by means of electrostatic potential maps, which very convincingly explain mutual orientation of molecules in a dimer. For the first time, numerous examples of interactions that violate both the concept of privileged exchange channels proposed by Pendás and his collaborators as well as inequalities obtained by Tognetti and Joubert for the β parameter related to secondary interactions are presented. The possible cause of this violation is suggested. We also show that the so-called counterintuitive bond paths result from quite natural behavior of the electron density gradient vector, i. e. searching for those areas of space that are characterized by large values of electron density or the most expanded its distributions.

Is the Fluorine in Molecules Dispersive? Is Molecular Electrostatic Potential a Valid Property to Explore Fluorine-Centered Non-Covalent Interactions?

Can two sites of positive electrostatic potential localized on the outer surfaces of two halogen atoms (and especially fluorine) in different molecular domains attract each other to form a non-covalent engagement? The answer, perhaps counterintuitive, is yes as shown here using the electronic structures and binding energies of the interactions for a series of 22 binary complexes formed between identical or different atomic domains in similar or related halogen-substituted molecules containing fluorine. These were obtained using various computational approaches, including density functional and ab initio first-principles theories with M06-2X, RHF, MP2 and CCSD(T). The physical chemistry of non-covalent bonding interactions in these complexes was explored using both Quantum Theory of Atoms in Molecules and Symmetry Adapted Perturbation Theories. The surface reactivity of the 17 monomers was examined using the Molecular Electrostatic Surface Potential approach. We have demonstrated inter alia that the dispersion term, the significance of which is not always appreciated, which emerges either from an energy decomposition analysis, or from a correlated calculation, plays a structure-determining role, although other contributions arising from electrostatic, exchange-repulsion and polarization effects are also important. The 0.0010 a.u. isodensity envelope, often used for mapping the electrostatic potential is found to provide incorrect information about the complete nature of the surface reactive sites on some of the isolated monomers, and can lead to a misinterpretation of the results obtained.

Modulating intramolecular chalcogen bonds in aromatic (thio)(seleno)phene-based derivatives

Intramolecular chalcogen interactions have been studied for four different derivatives of compounds within two different families, S or Se, to evaluate the effect of these IMChBs in the stability of the interacting and non-interacting systems.

Solving the enigma of weak fluorine contacts in the solid state: a periodic DFT study of fluorinated organic crystals

The nature and strength of weak interactions with organic fluorine in the solid state are revealed by periodic density functional theory (periodic DFT) calculations coupled with experimental data on the structure and sublimation thermodynamics of crystalline organofluorine compounds. To minimize other intermolecular interactions, several sets of crystals of perfluorinated and partially fluorinated organic molecules are considered. This allows us to establish the theoretical levels providing an adequate description of the metric and electron-density parameters of the C–F⋯F–C interactions and the sublimation enthalpy of crystalline perfluorinated compounds. A detailed comparison of the C–F⋯F–C and C–H⋯F–C interactions is performed using the relaxed molecular geometry in the studied crystals. The change in the crystalline packing of aromatic compounds during their partial fluorination points to the structure-directing role of C–H⋯F–C interactions due to the dominant electrostatic contribution to these contacts. C–H⋯F–C and C–H⋯O interactions are found to be identical in nature and comparable in energy. The factors that determine the contribution of these interactions to the crystal packing are revealed. The reliability of the results is confirmed by considering the superposition of the electrostatic potential and electron density gradient fields in the area of the investigated intermolecular interactions.

The nature and strength of weak C–H⋯F–C and C–F⋯F–C interactions and their role in organofluorine molecular crystals were studied using periodic DFT coupled with CSD data mining and experimental sublimation enthalpies.

Combined QTAIM and ETS-NOCV investigation of the interactions in ClnM[PhB(NtBu)2] complexes with M = Si & Ge (n = 0), As & Sb (n = 1), Te & Po (n = 2)

In this work, the nature of the chemical interactions between the metalloid atom (M = Si, Ge, As, Sb, Te, Po) and the nitrogen atoms in the bora-amidinate (bam) complexes (Cl_nM[PhB(N^tBu)₂]) are investigated, mainly via density-based indices. The descriptors used are derived using the quantum theory of atoms in molecules and natural orbitals for chemical valence approaches. It is shown that the strongest interaction is achieved with silicon. Indeed, it is generally the lightest metalloid in a particular group of the periodic table (i.e., Si, As, and Te for groups 14-16, respectively) that exhibits the strongest bond in the bam complex. This suggests that the atomic radius of the metalloid is a useful parameter for predicting the bonding strength. Extended transition state (ETS) decomposition results indicate that the interactions are more electrostatic than due to orbital interactions.

Bond paths between distant atoms do not necessarily indicate dominant interactions

The goal of the article is to revive discussion on the interpretation of bond paths linking distant atoms, particularly tracing weak interactions in dimers. According to the Pendás' concept of privileged exchange channel, a bond path is formed between this pair of competing atoms, which is associated with larger value of the exchange energy. We point out that, due to the short-range nature of the exchange energy, bond paths linking distant atoms clearly become doubtful indicators of dominant intermolecular interactions, particularly if some other characteristics (geometric, spectroscopic, based on electrostatic parameters, etc.) indicate other intermolecular interactions as dominant. Several such cases are thoroughly investigated. We show that electrostatic parameters are much more reliable indicators of dominant intermolecular interactions than bond paths. Then, we pay attention that the presence of ("unexpected", i.e., not necessarily indicating dominant intermolecular interactions) bond paths between pairs of atoms featuring highly expanded charge distributions can be easily explained by visual exploration of isodensity contour plots. As always pointing in the direction of the steepest increase, the gradient vector of the electron density favors areas of its high values gaining higher exchange energy, yet being blind to highly electron deficient areas nearby, which, however, can quite often be involved in dominant intermolecular interactions as strongly suggested by many other bonding analysis. We also suggest that an interatomic component of Hellmann-Feynman force would most likely be the most reliable indicator of attractive or repulsive character of individual interatomic interaction. © 2018 Wiley Periodicals, Inc.

Can Combined Electrostatic and Polarization Effects Alone Explain the F···F Negative-Negative Bonding in Simple Fluoro-Substituted Benzene Derivatives? A First-Principles Perspective

The divergence of fluorine-based systems and significance of their nascent non-covalent chemistry in molecular assemblies are presented in a brief review of the field. Emphasis has been placed to show that type-I and -II halogen-centered F···F long-ranged intermolecular distances viable between the entirely negative fluorine atoms in some fluoro-substituted dimers of C6H6 can be regarded as the consequence of significant non-covalent attractive interactions. Such attractive interactions observed in the solid-state structures of C6F6 and other similar fluorine-substituted aromatic compounds have frequently been underappreciated. While these are often ascribed to crystal packing effects, we show using first-principles level calculations that these are much more fundamental in nature. The stability and reliability of these interactions are supported by their negative binding energies that emerge from a supermolecular procedure using MP2 (second-order Møller-Plesset perturbation theory), and from the Symmetry Adapted Perturbation Theory, in which the latter does not determine the interaction energy by computing the total energy of the monomers or dimer. Quantum Theory of Atoms in Molecules and Reduced Density Gradient Non-Covalent Index charge-density-based approaches confirm the F···F contacts are a consequence of attraction by their unified bond path (and bond critical point) and isosurface charge density topologies, respectively. These interactions can be explained neither by the so-called molecular electrostatic surface potential (MESP) model approach that often demonstrates attraction between sites of opposite electrostatic surface potential by means of Coulomb’s law of electrostatics, nor purely by the effect of electrostatic polarization. We provide evidence against the standalone use of this approach and the overlooking of other approaches, as the former does not allow for the calculation of the electrostatic potential on the surfaces of the overlapping atoms on the monomers as in the equilibrium geometry of a complex. This study thus provides unequivocal evidence of the limitation of the MESP approach for its use in gaining insight into the nature of reactivity of overlapped interacting atoms and the intermolecular interactions involved.

Revealing Factors Influencing the Fluorine-Centered Non-Covalent Interactions in Some Fluorine-Substituted Molecular Complexes: Insights from First-Principles Studies

We examine the equilibrium structure and properties of six fully or partially fluorinated hydrocarbons and several of their binary complexes using computational methods. In the monomers, the electrostatic surface of the fluorine is predicted to be either entirely negative or weakly positive. However, its lateral sites are always negative. This enables the fluorine to display an anisotropic distribution of charge density on its electrostatic surface. While this is the electrostatic surface scenario of the fluorine atom, its negative sites in some of these monomers are shown to have the potential to engage in attractive engagements with the negative site(s) on the same atom in another molecule of the same type, or a molecule of a different type, to form bimolecular complexes. This is revealed by analyzing the results of current state-of-the-art computational approaches such as DFT, together with those obtained from the quantum theory of atoms in molecules, molecular electrostatic surface potential and symmetry adapted perturbation theories. We demonstrate that the intermolecular interaction energy arising in part from the universal London dispersion, which has been underappreciated for decades, is an essential factor in explaining the attraction between the negative sites, although energy arising from polarization strengthens the extent of the intermolecular interactions in these complexes.

Hydride-Triel Bonds

In this article, we present the results of our comprehensive studies of 72 dimers of the R3XXH⋯YR3Y type (X = Si, Ge; Y = B, Al, Ga; R^X  = H, Cl, Me; R^Y  = H, F, Cl, Me) and featuring hydride-triel bonds (i.e., charge-inverted hydrogen bonds). Influence of X and Y atoms as well as R^X and R^Y substituents on various properties of these dimers is investigated in detail. In particular the strength of the H⋯Y hydride-triel bonds is paid a close attention and it is shown that hydride-triel bonds can be strong enough to considerably determine structure and properties of molecular systems. In addition, properties of the investigated dimers are largely governed by the charge transfer from the Lewis base to the Lewis acid, which is particularly important if more bulky and polarizable R^Y and Y atoms are present in the YR3Y molecule. Several excellent linear (R² close to 1) and exponential correlations between pairs of diverse parameters are presented. Few instances are discussed where somewhat unexpected bond paths exist between two atoms featuring partial negative charges (e.g., between hydride hydrogen and halogen and between lateral sides of two halogens) showing that in some cases a bond path prefers to link two closely spaced electron-rich atoms instead of two atoms that are expected to form a bond. © 2018 Wiley Periodicals, Inc.

Structural and computational understanding of weak interactions in “bridge-flipped” isomeric tetrafluoro-bis-benzylideneanilines

Structural and computational insights into inter-molecular interactions in isomeric bridge-flipped tetrafluoro-bis-benzylideneanilines.

Energetic Analysis of Conjugated Hydrocarbons Using the Interacting Quantum Atoms Method

A number of aromatic, antiaromatic, and nonaromatic organic molecules was analyzed in terms of the contributions to the electronic energy defined in the quantum theory of atoms in molecules and the interacting quantum atoms method. Regularities were found in the exchange and electrostatic interatomic energies showing trends that are closely related to those of the delocalization indices defined in the theory. In particular, the CC interaction energies between bonded atoms allow to rationalize the energetic stabilization associated with the bond length alternation in conjugated polyenes. This approach also provides support to Clar's sextet rules devised for aromatic systems. In addition, the H⋯H bonding found in some of the aromatic molecules studied was of an attractive nature, according to the stabilizing exchange interaction between the bonded H atoms. © 2017 Wiley Periodicals, Inc.

Quantitative analysis of solid-state diversity in trifluoromethylated phenylhydrazones

The cooperative roles of various structural motifs associated with the presence of different intermolecular interactions in the formation of molecular crystals are investigated in a series of trifluoromethylated phenylhydrazones. Out of the six compounds analysed, two exhibit three-dimensional structural similarities with geometrically equivalent building blocks, while a third exists as two polymorphic forms crystallized from ethanol solutions at low temperature (277 K) and room temperature (298 K), respectively. The compounds were characterized via single-crystal and powder X-ray diffraction techniques and differential scanning calorimetry. In the absence of any strong hydrogen bonding, the supramolecular constructs are primarily stabilized via molecular pairs with a high dispersion-energy contribution, due to the presence of molecular stacking along the molecular backbone along with C-H...π interactions in the solid state, in preference to an electrostatic contribution. The interaction energies for the most stabilizing molecular building blocks are in the range -29 to -43 kJ mol^-1. In addition, weak N-H...F, C-H...F and N-H...C interactions and F...F, F...C, F...N and C...N contacts act as secondary motifs, providing additional stability to the crystal packing. The overall molecular arrangements are carefully analysed in terms of their nature and energetics, and the roles of different molecular pairs towards the crystal structure are delineated. A topological study using the quantum theory of atoms in molecules was used to characterize all the atomic interactions in the solid state. It established the presence of (3, -1) bond critical points and the closed-shell nature of all the interactions.