Next Article in Journal
Methanolic Extract of Piper sarmentosum Attenuates Obesity and Hyperlipidemia in Fructose-Induced Metabolic Syndrome Rats
Next Article in Special Issue
Measuring Shared Electrons in Extended Molecular Systems: Covalent Bonds from Plane-Wave Representation of Wave Function
Previous Article in Journal
Encapsulation of Bioactive Phytochemicals in Plant-Based Matrices and Application as Additives in Meat and Meat Products
Previous Article in Special Issue
Combining Molecular Dynamic Information and an Aspherical-Atom Data Bank in the Evaluation of the Electrostatic Interaction Energy in Multimeric Protein-Ligand Complex: A Case Study for HIV-1 Protease
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hydrogen-Bonded and Halogen-Bonded: Orthogonal Interactions for the Chloride Anion of a Pyrazolium Salt

1
Institute of Inorganic Chemistry, RWTH Aachen University, Landoltweg 1, 52074 Aachen, Germany
2
Institute of Molecular Science, Shanxi University, Taiyuan 030006, China
*
Author to whom correspondence should be addressed.
Molecules 2021, 26(13), 3982; https://doi.org/10.3390/molecules26133982
Submission received: 31 May 2021 / Revised: 25 June 2021 / Accepted: 25 June 2021 / Published: 29 June 2021

Abstract

:
In the hydrochloride of a pyrazolyl-substituted acetylacetone, the chloride anion is hydrogen-bonded to the protonated pyrazolyl moiety. Equimolar co-crystallization with tetrafluorodiiodobenzene (TFDIB) leads to a supramolecular aggregate in which TFDIB is situated on a crystallographic center of inversion. The iodine atom in the asymmetric unit acts as halogen bond donor, and the chloride acceptor approaches the σ -hole of this TFDIB iodine subtending an almost linear halogen bond, with Cl···I = 3.1653 ( 11 ) Å and Cl···I–C = 179.32(6)°. This contact is roughly orthogonal to the N–H···Cl hydrogen bond. An analysis of the electron density according to Bader’s Quantum Theory of Atoms in Molecules confirms bond critical points (bcps) for both short contacts, with ρ bcp = 0.129 for the halogen and 0.321 e Å 3 for the hydrogen bond. Our halogen-bonded adduct represents the prototype for a future class of co-crystals with tunable electron density distribution about the σ -hole contact.

1. Introduction

Starting with his first publications in the 1920s, Linus Pauling took a never-ceasing interest in both theoretical [1,2] and experimental [3,4] insights in chemical bonds. This fascination is reflected by his famous textbook The Nature of the Chemical Bond [5] and peaked in him being awarded the Nobel prize in chemistry in 1954 “for his research into the nature of the chemical bond and its application to the elucidation of the structure of complex substances” [6]. Pauling achieved the aforementioned experimental elucidation of chemical structures predominantly by X-ray methods, e.g., single-crystal X-ray diffraction (SCXRD). In The Nature of the Chemical Bond, Pauling also recognized the strongly electrostatic nature of hydrogen bonds.
A halogen’s ability to form favourable interactions with nucleophiles is called a halogen bond [7]. Similar to hydrogen bonds, they form between an electron donating and an electron accepting site; in the case of halogen bonds the halogen is the electron accepting partner (Figure 1) [7,8].
By definition, the halogen is the donor of the halogen bond. Because of the localization of the electron deficiency on the opposing site of the σ -bond to R it is called the σ -hole [9]. First crystallographic evidence for a halogen bond was given by the Hassel group in 1954: they described the 1:1 adduct of 1,4-dioxane with Br2 [10]. Although the term halogen bond was not used until introduced by Dumas et al. in 1978, this short O···Br contact of approx. 2.71  Å paved the way for forthcoming research regarding the character of these interactions. Since then, halogen bonds have been studied and utilized extensively in various fields including macromolecular and supramolecular chemistry, theory concerning chemical bonds and crystal engineering.
More recently, halogen bonds and related interactions have gained increasing attention from theory [11,12,13,14], and the σ -hole model has also been supported by experimental charge density studies [15,16,17,18,19]. In view of their occurrence in various fields, including macromolecular and supramolecular chemistry, halogen bonds have also been applied in the design of extended solids, i.e., in crystal engineering [20]. It has also been shown that molecules with enforced intramolecular hydrogen bonds to a halogen bond donor form stronger halogen bonds to adjacent halides [21].
For the design of extended structures, we often employ heterobifunctional molecules. Their different sites can selectively interact to form well-ordered assemblies, mostly based on coordinative bonds [22]. As our ditopic molecules usually exhibit functional groups with nitrogen lone pairs, e.g., nitriles [23,24], pyridines [25] or pyrazoles [26,27], they may alternatively engage as nucleophiles in halogen bonds (Figure 2) [28].
We have recently shown that pyrazolyl substituted acetylacetones are heterobifunctional molecules suitable for the construction of well-ordered bimetallic coordination polymers [27]. In this study we do not report metal coordination but rather hydrogen and halogen bonds with a hydrochloride. We discuss the synthesis and characterization of 3-(3,5-dimethyl-1-phenyl-1H-pyrazol-4-yl)acetylacetone (HacacPhPz, 1) and the cocrystal formed by its hydrochloride with tetrafluorodiiodobenzene (TFDIB) 2, a popular building block in crystal engineering [20].

2. Results and Discussion

2.1. Experimental Crystal Structures

We will first address the heterobifunctional molecule HacacPhPz itself. Like our previous reported compounds, it is obtained by the reaction of the appropriate substituted hydrazine derivative, e.g., phenylhydrazine, with tetraacetylethane [26,27,29]. HacacPhPz (1) crystallizes in the orthorhombic space group P b c a with Z = 16 and Z = 2 ; one of the independent molecules is shown in Figure 3.
Like most acetylacetones substituted in 3-position, 1 predominantly exists as the enol tautomer in solid state as well as in solution (Figure S3).
In the crystalline solid at 100 K , the enol hydrogen is clearly localized and is detected as local electron density maximum in a difference Fourier map (Figure 4). This assignment is corroborated by the C–O and C–C bond lengths in the acetylacetone moiety: in the case of the first molecule in the asymmetric residue depicted in Figure 4, O1–C2 ( 1.311 ( 2 )  Å) is significantly longer than O2–C4 ( 1.278 ( 2 )  Å), and in the carbon backbone C2–C3 ( 1.395 ( 3 )  Å) is shorter than C3–C4 ( 1.424 ( 3 )  Å).
In the following comparison, references are accompanied by their CCDC refcodes. Unsubstituted acetylacetone (refcode LIWPIQ [31]) crystallizes as the enol in the orthorhombic space group P n m a with the crystallographic mirror plane perpendicular to the least squares plane through the molecule, running through C3 and the hydrogen attached to it. For reasons of symmetry, disorder of the enol H with half occupancy is enforced. Boese et al. noted that a stable refinement of the structure could be conducted in the non-centrosymmetric subgroup P n a 2 1 , with two local electron density maxima for the disordered enol H; occupancies resulted as equal within error. The interatomic C–O ( 1.2909 ( 13 )  Å) and C–C ( 1.4028 ( 16 )  Å) distances in unsubstituted acetylacetone represent average values of the bond lengths found in 1. The C=O bond distance for the keto group in 1 ( 1.278 ( 2 )  Å) is slightly longer than corresponding distances in small ketones such as methyl ethyl ketone (MEK, LASLAU [32], 1.2128 ( 19 )  Å) or pentanone (ZEJFEZ01 [33], 1.216 ( 4 )  Å). In contrast, the C–C bond adjacent to the keto side in 1 ( 1.395 ( 3 )  Å) is shorter than in the non-conjugated compounds ( 1.485 ( 2 )  Å for LASLAU, 1.493 ( 5 )  Å for ZEJFEZ01). We terminate our comparison with the literature by referring to a very recent article by Martinez et al. [34] in which the authors could correlate the enol H position, i.e., acetylacetone tautomerism with the presence of halogen bonds.
An overlay of both molecules in the asymmetric unit (Figure 5) reveals that the independent residues do not adopt the same conformation. Their planar subunits acetylacetone, pyrazolyl and phenyl subtend different dihedral angles.
In the packing in 1, no particularly close directional interactions occur. Acetylacetone moieties and phenyl substituents aggregate in separate regions; C–H···O interactions of approx. 2.5  Å represent the shortest intermolecular contacts.
In order to probe the influence of the N-substitution on the ability to form halogen bonds, co-crystallization of 1 with TFDIB was attempted. We have not been able to isolate any adduct, presumably due to steric congestion about the nucleophilic N1, but this site was still expected to be accessible by a proton. Therefore, hydrochloric acid was added and crystals of the solid HacacPhPz·HCl·0·5TFDIB (2) were obtained. This product crystallizes in the monoclinic space group P 2 1 / c with Z = 4 , i.e., with a TFDIB moiety on a crystallographic centre of inversion (Figure 6).
The chloride anion Cl1 bridges two protonated H2acacPhPz+ cations to a central TFDIB moiety via two roughly orthogonal short contacts. The first interaction involves a classical N–H···Cl hydrogen bond, with a donor···acceptor distance of 2.97  Å; additional geometric details have been compiled in the caption of Figure 6. We recall that Pauling has emphasized the strongly electrostatic nature of hydrogen bonds as early as 1939. In the second short contact, the chloride ion acts as halogen bond acceptor: in agreement with σ -hole theory [7,8], it approaches the iodine atom in the direction opposite to the C–I bond; the arrangement Cl1···I1–C17 is almost linear. Only a few examples of this kind have been described to date [35,36,37,38,39,40], and Table 1 provides an overview of the most relevant geometry parameters.
Both linear and roughly orthogonal interactions are known for interhalogene ions. According to the popular VSEPR (Valence Shell Electron Pair Repulsion, aka Gillespie-Nyholm) interpretation [41,42] a linear disposition for two bonding partners and three lone pairs can be expected in formal 10 valence electron species such as the central iodine in I3 or the recently described ClF2 ions [43]. In contrast, for a formal 8 electron moiety such as ClF2+ [44] or all the examples given in Table 1 in which a Cl anion formally employs two out of four lone pairs, one for accepting the hydrogen bond and one as a nucleophile towards iodine, a bent arrangement can be predicted.

2.2. Theoretical Evaluation of the Halogen Bond

The Hirshfeld surface [45] about the chlorine anion is shown in Figure 7 and clearly reflects the short contacts mentioned above as red contact regions.
Additional insight concerning the electron density distribution in 2 was obtained via a single point calculation, followed by an analysis of the resulting electron density by Bader’s Quantum Theory of Atoms in Molecules [47]. Figure 8 shows a trajectory plot of the gradient of the electron density from a similar view direction as the Hirshfeld surface in Figure 7; both the classical N–H···Cl hydrogen bond and the Cl···I halogen bond are associated with essentially linear bond paths and ( 3 , 1 ) critical points (bcps).
Table 2 summarizes relevant characteristics of the electron density in the bcps of the short contacts and its derived properties.
We are not aware of experimental charge density studies on I···Cl halogen bonds and therefore cannot compare ρ bcp for I1···Cl1 to other structures of the same class but we have encountered similar values for the shortest Cl···Cl contacts investigated by high resolution X-ray diffraction [49,50]. Geometric and electronic parameters of the classical N–H···Cl hydrogen bond may be compared to the situation in bis(2-chloropyridinium) tetrachloridozincate; for the latter compound, the electron density was established experimentally [51]. Table 2 shows that the calculated properties for the hydrogen bond in 2 and our previous experimental observations for a contact of similar geometry match well, albeit the experimental values are obviously associated with appreciable standard uncertainties. In addition to distance criteria and ρ bcp , energy densities in the bond critical point have proven useful to categorize short contacts. Only the hydrogen bond is associated with a negative total energy density E, the sum of the (positive) kinetic energy density G and the (negative) potential energy density V. The halogen bond, in contrast, is characterized by a small positive E bcp , typically encountered for weak closed shell interactions [52].
The Laplacian of the electron density in the region of the short contacts is depicted in Figure 9. A slight polarization of the Cl anion towards the hydrogen can be perceived. The iodine atom is associated with two well visible valence charge concentrations perpendicular to the σ and the halogen bond.
The short contact between the chloride anion and the iodine atom of the fluorinated aromatic ring does not only reflect the geometry of a σ -hole interaction but also matches the electrostatics of such a contact: Figure 10 shows the electrostatic potential for 2. The positive region on the iodine atom opposite to its covalent bond to the ring carbon atom is clearly visible.

3. Conclusions and Outlook

Hydrogen bonds and halogen bonds contribute to the structure of crystalline 2. This balance does not come by accident and may become a useful motif in crystal engineering. Hydrochlorides and -bromides of a wide range of basic compounds are readily available. The interaction of their halide anion X with a halogen bond donor represents a rarely exploited way to tune and investigate both the geometry and the charge density of halogen bonds. The formation of a classical hydrogen bond between the protonated site and X will often occur, and the halide engaged in this interaction can be expected to act as a good halogen bond acceptor, both with respect to its electronic and steric properties. The possibility to vary the parent base, the hydrogen halide, and the halogen bond donor molecule makes this approach attractive for the systematic study of halogen bonds. In favorable cases, the crystalline products will not only allow conventional structural studies and theoretical calculations but can also become the target of experimental charge density studies.

4. Experimental Section

4.1. Computational Details

Prior to the single point calculation, C–H and N–H distances were idealized to values consistent with results from neutron diffraction [53]. The single point calculation was performed with Gaussian [54]; the MIDIX basis set [55] was used. In the calculation complete residues were taken into account and therefore a slightly expanded asymmetric unit was used. It comprised the residues depicted in Figure S5, i.e., a full rather than a half TFDIB molecule. The electron density associated with the single point calculation was analyzed with AIMAll [48] and Multiwfn [56] and interpreted according to Bader’s QTAIM [47]. The kinetic energy density G and the ratio between kinetic energy density and electron density, G / ρ in the bcp, were derived as suggested by Abramov [57], and the potential energy density V was obtained according to the local virial theorem [58,59].

4.2. Materials and Methods

All chemicals were used without further purification. Magnetic resonance spectra were measured using a Bruker Avance II UltrashieldT11 plus 400 instrument (400 M Hz , referenced to SiMe4). Infrared spectra were recorded using a Nicolet Avatar 360 E.S.P. spectrometer in potassium bromide windows. Elemental analyses were performed using a Heraeus CHNO-Rapid VarioEL. 3,4-Diacetylhexane-2,5-dione (tetraacetylethane, TAE) was prepared according to known procedures [26,29,60]. Intensity data was collected with a Bruker D8 goniometer equipped with an APEX CCD area detector and an Incoatec microsource (Mo-K α radiation, λ = 0.71073 Å, multilayer optics). Temperature was maintained by using an Oxford Cryostream 700 instrument, Oxfordshire, UK. Data was integrated with SAINT [61] and corrected for absorption by multi-scan methods [62]. The structures were solved by intrinsic phasing [63] and refined by full matrix least squares procedures against F 2 , as implemented in SHELXL-18. Crystal data, data collection parameters and refinement results have been compiled in Table S1. CIFs have been deposited under CCDC No. 2086575 (1) and 2086574 (2). Powder diffraction experiments were performed on flat samples at room temperature using a STOE STADI P diffractometer with Guinier geometry (Cu-K α 1 , λ = 1.54059 Å, Johann germanium monochromator, STOE image plate detector IP-PSD, 0.005° step width in 2 θ ).

4.3. Synthesis of 3-(3,5-Dimethyl-1-phenyl-1 H-pyrazol-4-yl)acetylacetone, HacacPhPz, 1

TAE ( 0.991   g , 5.0   mmol , 1.0   eq . ) was dispersed in absolute ethanol (50 m L ) under N2 atmosphere and heated to reflux. When all TAE was dissolved, phenylhydrazine ( 0.49   m L , 5.0   mmol , 1.0   eq . ) was added dropwise and the solution was kept at reflux for an additional 5 h . The solvent was removed under reduced pressure yielding a yellow oil. Column chromatography (silica, n-pentane/EtOAc, 3:1, v:v; R f = 0.71 ) yields a pale yellow solid. Yield: 523 m g (39%). Colorless plate-shaped crystals suitable for SCXRD were obtained after recrystallization from n-hexane. 1 H NMR (CDCl3, 400 M Hz ): δ 16.89 (s, 1H), 7.51–7.42 (m, 4H), 7.41–7.33 (m, 1H), 2.18 (s, 6H), 1.94 (s, 6H). 13 C{ 1 H} NMR (CDCl3, 100 M Hz ): δ 192.83, 148.73, 139.94, 138.16, 129.27, 127.63, 124.68, 115.34, 104.20, 23.87, 12.33, 11.52. CHN: anal. calcd. for C16H18N2O2: C: 71.1%, H: 6.7%, N: 10.4%; found: C: 70.6%, H: 6.6%, N: 10.2%. HRMS-ESI ( m / z ): [M + H] + calcd. for C16H19N2O2+: 271.14410; found: 271.14409. mp: 63.5   C . While some similarities can be derived from the simulated and experimental powder pattern, phase purity could not be confirmed (Figure S1).

4.4. Synthesis of HacacPhPz·HCl·0·5TFDIB, 2

HacacPhPz ( 13.6   m g , 0.05   mmol , 2.0   eq . ) and TFDIB ( 10.1   m g , 0.025   mmol , 1.0   eq . ) were dissolved in CHCl3 ( 0.25   m L ) each. The two solutions were combined and placed in a bigger vial containing aq. HCl ( 0.91   m L , 5.4   M ) for slow vapor diffusion. Colorless crystals were obtained after slow evaporation of the solvent at room temperature. Yield: 13.4   m g (76%). The powder pattern of ground 2 does not match the simulation of the experimentally established single crystal structure (Figure S1, bottom). In order to investigate whether 2 decomposes upon grinding, a moist pH-indicator paper was placed in the head space of a bulk sample of 2. The color change of this indicator within several hours proved the release of a gaseous acid, most likely HCl. We do, however, not observe complete loss of HCl and re-generation of the constituents 1 and TFDIB but formation of an unknown crystalline product (Figure S2). On the scale of single crystals rather than ground powder, the loss of HCl is slower, in agreement with the matching microanalytical data.

Supplementary Materials

The following are available online, Figure S1: Simulated and experimental powder patterns of 1 (top) and 2 (bottom), Figure S2: Experimental powder patterns of 2 and simulated patterns of 1 and TFDIB [64], Figure S3: 1 H (top) and 13 C{ 1 H} (bottom) NMR spectra of 1 measured in CDCl3 (*) at room temperature, Figure S4: HSQC (red) and HMBC (green) NMR spectra of 1 measured in CDCl3 (*) at room temperature, Figure S5: Structure fragment used for the single point calculation discussed in the main text, Section 2.2, Table S1: Crystal data and refinement results for SCXRD data for 1 and 2 measured at T = 100   K , Table S2: Topological properties of interactions at their bond critical point ( 3 , 1 ) of 2.

Author Contributions

Conceptualization, S.v.T. and U.E.; methodology, U.E.; software, R.W.; validation, S.v.T and D.B.; formal analysis, S.v.T., R.W. and U.E.; investigation, S.v.T. and D.B.; resources, U.E.; data curation, S.v.T. and U.E.; writing—original draft preparation, S.v.T., D.B., R.W. and U.E.; writing—review and editing, S.v.T., D.B., R.W. and U.E.; visualization, S.v.T., D.B., R.W. and U.E.; supervision, S.v.T. and U.E.; project administration, U.E.; funding acquisition, U.E. All authors have read and agreed to the published version of the manuscript.

Funding

An RWTH fellowship for doctoral students to S.v.T. is gratefully acknowledged.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article or Supplementary Materials.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
bcpbond critical point
HacacPhPz3-(3,5-dimethyl-1-phenyl-1H-pyrazol-4-yl)acetylacetone
QTAIMQuantum Theory of Atoms in Molecules
SCXRDsingle-crystal X-ray diffraction
TAEtetraacetylethane
TFDIBtetrafluorodiiodobenzene

References

  1. Pauling, L. The Dynamic Model of the Chemical Bond and its Application to the Structure of Benzene. J. Am. Chem. Soc. 1926, 48, 1132–1143. [Google Scholar] [CrossRef]
  2. Pauling, L. The theoretical prediction of the physical properties of many electron atoms and ions. Mole refraction, diamagnetic susceptibility, and extension in space. Proc. R. Soc. Lond. A 1927, 114, 181–211. [Google Scholar] [CrossRef] [Green Version]
  3. Pauling, L. The Crystal Structure of Magnesium Stannide. J. Am. Chem. Soc. 1923, 45, 2777–2780. [Google Scholar] [CrossRef]
  4. Pauling, L.; Dickinson, R.G. The Crystal Structure of Uranyl Nitrate Hexahydrate. J. Am. Chem. Soc. 1924, 46, 1615–1622. [Google Scholar] [CrossRef]
  5. Pauling, L. The Nature of the Chemical Bond, 3rd ed.; Cornell University Press: Ithaca, NY, USA, 1960. [Google Scholar]
  6. Nobel Prize Committee. The Nobel Prize in Chemistry 1954; Nobel Prize Committee: Stockholm, Sweden, 1954. [Google Scholar]
  7. Cavallo, G.; Metrangolo, P.; Milani, R.; Pilati, T.; Priimagi, A.; Resnati, G.; Terraneo, G. The Halogen Bond. Chem. Rev. 2016, 116, 2478–2601. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Bent, H.A. Structural chemistry of donor-acceptor interactions. Chem. Rev. 1968, 68, 587–648. [Google Scholar] [CrossRef]
  9. Clark, T.; Hennemann, M.; Murray, J.S.; Politzer, P. Halogen bonding: The sigma-hole. Proceedings of “Modeling interactions in biomolecules II”, Prague, September 5th–9th, 2005. J. Mol. Model. 2007, 13, 291–296. [Google Scholar] [CrossRef]
  10. Hassel, O.; Hvoslef, J. The Structure of Bromine 1,4-Dioxonate. Acta Chim. Scand. 1954, 8, 873. [Google Scholar] [CrossRef] [Green Version]
  11. Politzer, P.; Murray, J.S. Halogen bonding: An interim discussion. ChemPhysChem 2013, 14, 278–294. [Google Scholar] [CrossRef] [PubMed]
  12. Wolters, L.P.; Schyman, P.; Pavan, M.J.; Jorgensen, W.L.; Bickelhaupt, F.M.; Kozuch, S. The many faces of halogen bonding: A review of theoretical models and methods. Comput. Mol. Sci. 2014, 4, 523–540. [Google Scholar] [CrossRef]
  13. Kirshenboim, O.; Kozuch, S. How to Twist, Split and Warp a sigma-Hole with Hypervalent Halogens. J. Phys. Chem. A 2016, 120, 9431–9445. [Google Scholar] [CrossRef]
  14. Murray, J.S.; Politzer, P. Interaction and Polarization Energy Relationships in sigma-Hole and pi-Hole Bonding. Crystals 2020, 10, 76. [Google Scholar] [CrossRef] [Green Version]
  15. Forni, A. Experimental and theoretical study of the Br...N halogen bond in complexes of 1,4-dibromotetrafluorobenzene with dipyridyl derivatives. J. Phys. Chem. A 2009, 113, 3403–3412. [Google Scholar] [CrossRef]
  16. Wang, R.; Kalf, I.; Englert, U. Insight into trifluoromethylation—Experimental electron density for Togni reagent I. RSC Adv. 2018, 8, 34287–34290. [Google Scholar] [CrossRef] [Green Version]
  17. Wang, R.; Hartnick, D.; Englert, U. Short is strong: Experimental electron density in a very short N···I halogen bond. Z. Kristallogr. Cryst. Mater. 2018, 233, 733–744. [Google Scholar] [CrossRef]
  18. Wang, R.; George, J.; Potts, S.K.; Kremer, M.; Dronskowski, R.; Englert, U. The many flavours of halogen bonds—Message from experimental electron density and Raman spectroscopy. Acta Crystallogr. 2019, C75, 1190–1201. [Google Scholar] [CrossRef] [Green Version]
  19. Otte, F.; Kleinheider, J.; Hiller, W.; Wang, R.; Englert, U.; Strohmann, C. Weak yet Decisive: Molecular Halogen Bond and Competing Weak Interactions of Iodobenzene and Quinuclidine. J. Am. Chem. Soc. 2021, 143, 4133–4137. [Google Scholar] [CrossRef]
  20. Nemec, V.; Lisac, K.; Bedeković, N.; Fotović, L.; Stilinović, V.; Cinčić, D. Crystal engineering strategies towards halogen-bonded metal—Organic multi-component solids: Salts, cocrystals and salt cocrystals. CrystEngComm 2021, 23, 3063–3083. [Google Scholar] [CrossRef]
  21. Riel, A.M.S.; Decato, D.A.; Sun, J.; Massena, C.J.; Jessop, M.J.; Berryman, O.B. The intramolecular hydrogen bonded-halogen bond: A new strategy for preorganization and enhanced binding. Chem. Sci. 2018, 9, 5828–5836. [Google Scholar] [CrossRef] [Green Version]
  22. Kremer, M.; Englert, U. N Donor substituted acetylacetones—Versatile ditopic ligands. Z. Kristallogr. Cryst. Mater. 2018, 233, 437–452. [Google Scholar] [CrossRef]
  23. Burrows, A.D.; Cassar, K.; Mahon, M.F.; Warren, J.E. The stepwise formation of mixed-metal coordination networks using complexes of 3-cyanoacetylacetonate. Dalton Trans. 2007, 24, 2499–2509. [Google Scholar] [CrossRef]
  24. Kondracka, M.; Englert, U. Bimetallic coordination polymers via combination of substitution-inert building blocks and labile connectors. Inorg. Chem. 2008, 47, 10246–10257. [Google Scholar] [CrossRef] [PubMed]
  25. Chen, B.; Fronczek, F.R.; Maverick, A.W. Porous Cu-Cd mixed-metal-organic frameworks constructed from Cu(Pyac)2 Bis3-(4-pyridyl)pentane-2,4-dionatocopper(II). Inorg. Chem. 2004, 43, 8209–8211. [Google Scholar] [CrossRef]
  26. Guo, Q.; Englert, U. An Acetylacetonate or a Pyrazole? Both! 3-(3,5-Dimethyl-pyrazol-4-yl)pentane-2,4-dione as a Ditopic Ligand. Cryst. Growth Des. 2016, 16, 5127–5135. [Google Scholar] [CrossRef]
  27. Van Terwingen, S.; Nachtigall, N.; Ebel, B.; Englert, U. N-Donor-Functionalized Acetylacetones for Heterobimetallic Coordination Polymers, the Next Episode: Trimethylpyrazoles. Cryst. Growth Des. 2021, 21, 2962–2969. [Google Scholar] [CrossRef]
  28. Merkens, C.; Pan, F.; Englert, U. 3-(4-Pyridyl)-2,4-pentanedione—A bridge between coordinative, halogen, and hydrogen bonds. CrystEngComm 2013, 15, 8153. [Google Scholar] [CrossRef]
  29. Mosby, W.L. The Reaction of Some 1:4-Dicarbonyl Systems with Hydrazine. J. Chem. Soc. 1957, 3997–4003. [Google Scholar] [CrossRef]
  30. Spek, A.L. Structure validation in chemical crystallography. Acta Crystallogr. 2009, D65, 148–155. [Google Scholar] [CrossRef]
  31. Boese, R.; Antipin, M.Y.; Bläser, D.; Lyssenko, K. Molecular Crystal Structure of Acetylacetone at 210 and 110 K: Is the Crystal Disorder Static or Dynamic? J. Phys. Chem. B 1998, 102, 8654–8660. [Google Scholar] [CrossRef]
  32. Yufit, D.S.; Howard, J.A.K. Low melting molecular complexes. The structures of molecular complexes of tri- and di-chloromethanes with small ketones and 1,4-dioxane. CrystEngComm 2012, 14, 2003. [Google Scholar] [CrossRef]
  33. Kalf, I.; Englert, U. The molecular conformation of pentan-3-one studied in cholic acid pentan-3-one solvate. Acta Crystallogr. 2011, C67, o206–o208. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Martinez, V.; Bedeković, N.; Stilinović, V.; Cinčić, D. Tautomeric Equilibrium of an Asymmetric beta-Diketone in Halogen-Bonded Cocrystals with Perfluorinated Iodobenzenes. Crystals 2021, 11, 699. [Google Scholar] [CrossRef]
  35. Abdalrahman, M.A.; Awwadi, F.F.; Jameson, G.B.; Landee, C.P.; Saunders, C.G.; Turnbull, M.M.; Wikaira, J.L. Effects of halogen and hydrogen bonding on defect disorder: The ladder that wasn’t there. CrystEngComm 2013, 15, 4309. [Google Scholar] [CrossRef]
  36. Zordan, F.; Purver, S.L.; Adams, H.; Brammer, L. Halometallate and halide ions: Nucleophiles in competition for hydrogen bond and halogen bond formation in halopyridinium salts of mixed halide–halometallate anions. CrystEngComm 2005, 7, 350. [Google Scholar] [CrossRef]
  37. Raatikainen, K.; Rissanen, K. Modulation of N...I and + N-H...Cl - ...I Halogen Bonding: Folding, Inclusion, and Self-Assembly of Tri- and Tetraamino Piperazine Cyclophanes. Cryst. Growth Des. 2010, 10, 3638–3646. [Google Scholar] [CrossRef]
  38. Ducharme, G.; Nesterov, E.E.; Fronczek, F.R. Private Communication (Refcode VIDHEY); CCDC: Cambridge, UK, 2018. [Google Scholar] [CrossRef]
  39. Mondal, B.; Zhang, T.; Prabhakar, R.; Captain, B.; Ramamurthy, V. Role of hydrogen bonds in molecular packing of photoreactive crystals: Templating photodimerization of protonated stilbazoles in crystalline state with a combination of water molecules and chloride ions. Photochem. Photobiol. Sci. 2014, 13, 1509–1520. [Google Scholar] [CrossRef]
  40. Fotović, L.; Stilinović, V. Halogenide anions as halogen and hydrogen bond acceptors in iodopyridinium halogenides. CrystEngComm 2020, 22, 4039–4046. [Google Scholar] [CrossRef]
  41. Gillespie, R.J.; Nyholm, R.S. Inorganic Stereochemistry. Q. Rev. Chem. Soc. 1957, 11, 339–380. [Google Scholar] [CrossRef]
  42. Gillespie, R.J.; Robinson, E.A. Electron Domains and the VSEPR Model of Molecular Geometry. Angew. Chem. Int. Ed. 1996, 35, 495–514. [Google Scholar] [CrossRef]
  43. Pröhm, P.; Schmid, J.R.; Sonnenberg, K.; Voßnacker, P.; Steinhauer, S.; Schattenberg, C.J.; Müller, R.; Kaupp, M.; Riedel, S. Improved Access to Organo-Soluble Di- and Tetrafluoridochlorate(I)/(III) Salts. Angew. Chem. Int. Ed. 2020, 59, 16002–16006. [Google Scholar] [CrossRef]
  44. Ellern, A.M.; Antipin, M.Y.; Sharabarin, A.V.; Struchkov, Y.T. Crystal structures of difluorochloronium hexafluoroniobate and hexafluorotantalate, ClF2NbF6 and ClF2TaF6. Russ. J. Inorg. Chem. 1991, 36, 2266–2270. [Google Scholar]
  45. Spackman, M.A.; Jayatilaka, D. Hirshfeld surface analysis. CrystEngComm 2009, 11, 19–32. [Google Scholar] [CrossRef]
  46. Spackman, P.R.; Turner, M.J.; McKinnon, J.J.; Wolff, S.K.; Grimwood, D.J.; Jayatilaka, D.; Spackman, M.A. CrystalExplorer: A program for Hirshfeld surface analysis, visualization and quantitative analysis of molecular crystals. J. Appl. Cryst. 2021, 54, 1006–1011. [Google Scholar] [CrossRef]
  47. Bader, R.F.W. Atoms in Molecules—A Quantum Theory; Clarendon Press: Oxford, UK, 1990. [Google Scholar]
  48. Keith, T.A. AIMAll: Version 17.01.25; TK Gristmill Software: Overland Park, KS, USA, 2017. [Google Scholar]
  49. Wang, R.; Dols, T.S.; Lehmann, C.W.; Englert, U. The halogen bond made visible: Experimental charge density of a very short intermolecular Cl...Cl donor-acceptor contact. Chem. Commun. 2012, 48, 6830–6832. [Google Scholar] [CrossRef] [PubMed]
  50. Wang, R.; Dols, T.S.; Lehmann, C.W.; Englert, U. Charge Density of Intra- and Intermolecular Halogen Contacts. Z. Anorg. Allg. Chem. 2013, 639, 1933–1939. [Google Scholar] [CrossRef]
  51. Wang, A.; Wang, R.; Kalf, I.; Dreier, A.; Lehmann, C.W.; Englert, U. Charge-Assisted Halogen Bonds in Halogen-Substituted Pyridinium Salts: Experimental Electron Density. Cryst. Growth Des. 2017, 17, 2357–2364. [Google Scholar] [CrossRef]
  52. Espinosa, E.; Alkorta, I.; Elguero, J.; Molins, E. From weak to strong interactions: A comprehensive analysis of the topological and energetic properties of the electron density distribution involving X-H...F-Y systems. J. Chem. Phys. 2002, 117, 5529–5542. [Google Scholar] [CrossRef]
  53. Allen, F.H.; Bruno, I.J. Bond lengths in organic and metal-organic compounds revisited: X-H bond lengths from neutron diffraction data. Acta Crystallogr. 2010, B66, 380–386. [Google Scholar] [CrossRef]
  54. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. GAUSSIAN 16, Revision C.01; Gaussian Inc.: Wallingford, CT, USA, 2016. [Google Scholar]
  55. Easton, R.E.; Giesen, D.J.; Welch, A.; Cramer, C.J.; Truhlar, D.G. The MIDI! basis set for quantum mechanical calculations of molecular geometries and partial charges. Theor. Chim. Acta 1996, 93, 281–301. [Google Scholar] [CrossRef]
  56. Lu, T.; Chen, F. Multiwfn: A multifunctional wavefunction analyzer. J. Comput. Chem. 2012, 33, 580–592. [Google Scholar] [CrossRef]
  57. Abramov, Y.A. On the Possibility of Kinetic Energy Density Evaluation from the Experimental Electron-Density Distribution. Acta Crystallogr. 1997, A53, 264–272. [Google Scholar] [CrossRef]
  58. Espinosa, E.; Molins, E.; Lecomte, C. Hydrogen bond strengths revealed by topological analyses of experimentally observed electron densities. Chem. Phys. Lett. 1998, 285, 170–173. [Google Scholar] [CrossRef]
  59. Espinosa, E.; Lecomte, C.; Molins, E. Experimental electron density overlapping in hydrogen bonds: Topology vs. energetics. Chem. Phys. Lett. 1999, 300, 745–748. [Google Scholar] [CrossRef]
  60. Ponomarova, V.V.; Komarchuk, V.V.; Boldog, I.; Krautscheid, H.; Domasevitch, K.V. Modular construction of 3D coordination frameworks incorporating SiF62- links: Accessing the significance of [M(pyrazole)4SiF6] synthon. CrystEngComm 2013, 15, 8280. [Google Scholar] [CrossRef]
  61. Bruker. SAINT+: Program for Reduction of Data Collected on Bruker CCD Area Detector Diffractometer; Bruker: Madison, WI, USA, 2009. [Google Scholar]
  62. Bruker. SADABS; Bruker: Madison, WI, USA, 2008. [Google Scholar]
  63. Sheldrick, G.M. SHELXT—Integrated space-group and crystal-structure determination. Acta Crystallogr. 2015, A71, 3–8. [Google Scholar] [CrossRef] [Green Version]
  64. Chaplot, S.L.; McIntyre, G.J.; Mierzejewski, A.; Pawley, G.S. The High-Temperature Phase of 1,2,4,5-Tetrafluoro-3,6-diiodobenzene and the Phase Transition. Acta Crystallogr. 1981, B37, 2210–2214. [Google Scholar] [CrossRef]

Short Biography of Authors

Molecules 26 03982 i001Steven van Terwingen studied chemistry at RWTH Aachen University from 2013 on. He finished his Master degree under supervision of Ulli Englert in 2018 and is currently working on his Ph. D. thesis in the same group. He has a keen interest in crystal chemistry including crystal engineering, catalysis and structural aspects elucidated by diffraction methods.
Molecules 26 03982 i002Daniel Brüx was born in 1999 and studies chemistry at RWTH since 2018. He is currently working on his bachelor thesis under the supervision of Steven van Terwingen in the Englert group, in which he investigates the crystal chemistry of N-aryl pyrazolyl substituted acetylacetonates. His main interests are coordination chemistry and developing new crystallization procedures.
Molecules 26 03982 i003Ruimin Wang joined RWTH with a B. Sc. degree form Nankai University, China and obtained her Ph. D. with Ulli Englert in 2000. Her main area of interest as a senior scientist is experimental electron density and the application of high resolution diffraction data for the understanding of intra- and inter-molecular interactions.
Molecules 26 03982 i004Ulli Englert has fallen in love with diffraction methods decades ago as a student in Tübingen and Pisa. His chemical crystallography group in RWTH is involved in the synthesis of building blocks for extended structures and in the characterization of the products, both in the home lab and at large research facilities.
Figure 1. The halogen X exhibits an electron deficient site (red, δ + ) in direction of the σ -bond to R. The nucleophile Y can interact with this positively charged region via its lone pair, thus forming a halogen bond.
Figure 1. The halogen X exhibits an electron deficient site (red, δ + ) in direction of the σ -bond to R. The nucleophile Y can interact with this positively charged region via its lone pair, thus forming a halogen bond.
Molecules 26 03982 g001
Figure 2. Chemical structures of heterobifunctional molecules utilized by our group for crystal engineering exhibiting nitrogen lone pairs.
Figure 2. Chemical structures of heterobifunctional molecules utilized by our group for crystal engineering exhibiting nitrogen lone pairs.
Molecules 26 03982 g002
Figure 3. Displacement ellipsoid plot [30] of one molecule in the asymmetric residue of 1 (80% probability, C bonded hydrogens omitted). Selected intramolecular distances and angles (Å, °): O1–C2 1.311(2), O2–C4 1.278(2), C2–C3 1.395(3), C3–C4 1.424(3), C2–C3–C8–C9 85.7(2), N1–N2–C11–C16 28.6(2).
Figure 3. Displacement ellipsoid plot [30] of one molecule in the asymmetric residue of 1 (80% probability, C bonded hydrogens omitted). Selected intramolecular distances and angles (Å, °): O1–C2 1.311(2), O2–C4 1.278(2), C2–C3 1.395(3), C3–C4 1.424(3), C2–C3–C8–C9 85.7(2), N1–N2–C11–C16 28.6(2).
Molecules 26 03982 g003
Figure 4. Difference Fourier contour map [30] of the acetylacetone moiety in both molecules contained in the asymmetric residue of 1. Contour lines are drawn at 0.2 e Å 3 (red: positive difference, green: negative difference, blue: zero lines).
Figure 4. Difference Fourier contour map [30] of the acetylacetone moiety in both molecules contained in the asymmetric residue of 1. Contour lines are drawn at 0.2 e Å 3 (red: positive difference, green: negative difference, blue: zero lines).
Molecules 26 03982 g004
Figure 5. Overlay plot [30] of both molecules in the asymmetric residue in 1 (black: molecule 1 as shown in Figure 3, red: molecule 2 under symmetry operator x , y , z ; C bonded hydrogens omitted).
Figure 5. Overlay plot [30] of both molecules in the asymmetric residue in 1 (black: molecule 1 as shown in Figure 3, red: molecule 2 under symmetry operator x , y , z ; C bonded hydrogens omitted).
Molecules 26 03982 g005
Figure 6. Displacement ellipsoid plot [30] of 2 (80% probability, C bonded hydrogens omitted). Selected intramolecular distances and angles (Å, °): I1···Cl1 3.1653(11), Cl1···N1 2.970(2), I1···Cl1···N1 73.99(4), C17–I1···Cl1 179.32(6), Cl1···H1N–N1 172(3). Symmetry operation: a = 1 x , y , 1 z .
Figure 6. Displacement ellipsoid plot [30] of 2 (80% probability, C bonded hydrogens omitted). Selected intramolecular distances and angles (Å, °): I1···Cl1 3.1653(11), Cl1···N1 2.970(2), I1···Cl1···N1 73.99(4), C17–I1···Cl1 179.32(6), Cl1···H1N–N1 172(3). Symmetry operation: a = 1 x , y , 1 z .
Molecules 26 03982 g006
Figure 7. Hirshfeld surface [46] around Cl1 mapped with d norm ; regions marked in red represent directions of short, those in blue of long contact distances.
Figure 7. Hirshfeld surface [46] around Cl1 mapped with d norm ; regions marked in red represent directions of short, those in blue of long contact distances.
Molecules 26 03982 g007
Figure 8. Atomic basins [47,48] in 2; intramolecular bond paths and the conventional hydrogen bond are shown as solid black lines, the halogen bond and non-classical hydrogen bonds as dashed black lines.
Figure 8. Atomic basins [47,48] in 2; intramolecular bond paths and the conventional hydrogen bond are shown as solid black lines, the halogen bond and non-classical hydrogen bonds as dashed black lines.
Molecules 26 03982 g008
Figure 9. Laplacian of the electron density in 2; positive values in blue, negative values in red, contours at ± 2 n · 10 3 a. u. ( 0 n 20 ).
Figure 9. Laplacian of the electron density in 2; positive values in blue, negative values in red, contours at ± 2 n · 10 3 a. u. ( 0 n 20 ).
Molecules 26 03982 g009
Figure 10. Electrostatic potential for 2, mapped on an isosurface of electron density ρ = 0.05 a. u. [48]; red areas are associated with a positive value (0.480 a. u.), cyan areas with negative values (−0.0675 a. u.) and green areas with an ESP (0.115 a. u.).
Figure 10. Electrostatic potential for 2, mapped on an isosurface of electron density ρ = 0.05 a. u. [48]; red areas are associated with a positive value (0.480 a. u.), cyan areas with negative values (−0.0675 a. u.) and green areas with an ESP (0.115 a. u.).
Molecules 26 03982 g010
Table 1. Comparison between important distances and angles in 2 and literature known compounds with a similar structure motif.
Table 1. Comparison between important distances and angles in 2 and literature known compounds with a similar structure motif.
Molecules 26 03982 i005
Compoundd(I···Cl)/Åd(Cl···N)/Å(C–I···Cl)/°(I···Cl···N)/°
23.1653(11)2.970(2)179.32(6)73.99(4)
BEXPOL [35]3.2113.085177.81111.72
JAQNAR [36]3.3313.137169.4992.01
RUWVUB [37]3.2233.133172.72106.15
RUWWIQ [37]3.1023.103179.00111.33
VIDHEY [38]3.4893.124169.3271.89
WOQRIF [39]3.4222.989160.3877.60
JULRIU [40]3.1223.003176.4190.99
JULSAN [40]3.2403.033171.83101.30
For multiple hits in the same structure their average value was calculated and is presented here.
Table 2. Short contacts with properties of their bond critical point (bcp) ( 3 , 1 ) in 2 and experimental data from LAVNUU [49].
Table 2. Short contacts with properties of their bond critical point (bcp) ( 3 , 1 ) in 2 and experimental data from LAVNUU [49].
BondI1···Cl1 in 2Cl1···H1N in 2Cl···H in LAVNUU [49]
ρ / e Å 3 0.1290.3210.28
2 ρ / e Å 5 1.1841.7850.6
bond path length/Å3.16542.06802.11
G/a. u.0.01100.03000.018
G / ρ /a. u.0.580.630.44
V/a. u.−0.0097−0.0415−0.030
E/a. u.0.0123−0.0115−0.012
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

van Terwingen, S.; Brüx, D.; Wang, R.; Englert, U. Hydrogen-Bonded and Halogen-Bonded: Orthogonal Interactions for the Chloride Anion of a Pyrazolium Salt. Molecules 2021, 26, 3982. https://doi.org/10.3390/molecules26133982

AMA Style

van Terwingen S, Brüx D, Wang R, Englert U. Hydrogen-Bonded and Halogen-Bonded: Orthogonal Interactions for the Chloride Anion of a Pyrazolium Salt. Molecules. 2021; 26(13):3982. https://doi.org/10.3390/molecules26133982

Chicago/Turabian Style

van Terwingen, Steven, Daniel Brüx, Ruimin Wang, and Ulli Englert. 2021. "Hydrogen-Bonded and Halogen-Bonded: Orthogonal Interactions for the Chloride Anion of a Pyrazolium Salt" Molecules 26, no. 13: 3982. https://doi.org/10.3390/molecules26133982

APA Style

van Terwingen, S., Brüx, D., Wang, R., & Englert, U. (2021). Hydrogen-Bonded and Halogen-Bonded: Orthogonal Interactions for the Chloride Anion of a Pyrazolium Salt. Molecules, 26(13), 3982. https://doi.org/10.3390/molecules26133982

Article Metrics

Back to TopTop