Abstract
In this paper, we report the synthesis, characterization and theoretical studies of two new cocrystals of 1,5 bis(3-pyridyl)-3,4-diaza-2,4-hexadiene (3-bpdh) with two simple organic acids; oxalic and succinic with 1:1 stoichiometric ratio. These cocrystals, [(3-bpdh)(oxalic acid)] (
These structures have been studied theoretically using RB3LYP hybrid density functional method with 6–31G(d) basis set. DFT calculations have been used to predict the geometry of the cocrystal structures. Furthermore, molecular electrostatic potential (MESP) maps of the compounds have been generated. Theoretical results show good agreement with the experimental data and these data show hydrogen bonding has played an important role in supramolecular structure formation. The cocrystal design was done based on the perfect selection of two distinct hydrogen donor and acceptor molecules.
Introduction
A cocrystal is composed of two or more different components that are generally solid at room temperature in a crystalline lattice. The cocrystal structures are formed via non-ionic and non-covalent intermolecular interactions such as van der Waals forces and hydrogen bonding. Therefore, it may be possible to fine-tune physical properties by exercising precise control over the supramolecular assembly by crystal engineering [1]. Crystal engineering strategies would be improved physicochemical properties of chemical compounds or generate new properties for utilitarian applications without altering the inherent individual properties such as pharmaceutical and electronical [2–4], nano scale architectures [5], host–guest systems [6], topochemical reactions [7], solid-state reactivity [8] and interpenetrated network [9]. So, the cocrystal structures represent a significant synthetic challenge and a hot topic with considerable attentions in recently twenty years [10–13] and several theoretical studies have been performed on these structural molecules [14–19]. In this research, synthesis and characterization and theoretical studies of two novel cocrystals; [(3-bpdh)(oxalic acid)] (
Experimental
Reagent
All reagents and solvents used for the synthesis and analysis were commercially available, purchased from Merck and used as received without further purifications.
Methods
Infrared spectra were measured on a Perkin-Elmer 597 and Nicolet 510P spectrophotometers in the range 400–4000 cm–1 in KBr discs technique. The microanalyses were carried out using a Heraeus CHN Rapid analyzer. The melting points were measured on an electro-thermal 9100 apparatus. 1H-NMR spectra were measured with a BRUKER DRX-500 AVANCE spectrometer at 500 MHz. Thermal behavior was measured with PL-STA 1500 apparatus.
Preparation of [(3-bpdh)(oxalic acid)] (1)
1,5 bis(3-pyridyl)-3,4-diaza-2,4-hexadiene (3-bpdh) (0.12 g, 0.5 mmol) and oxalic acid (0.063 g, 0.5 mmol) were mixed and ground manually using a pestle and mortar then 5 mL of ethanol was added dropwise to obtain a pasty solution and stirred for 1 h at room temperature. The prepared solution was allowed to evaporate slowly for 2 days at room temperature (ca. 25°C). The yellowish crystals of the desired product precipitated, which were washed by acetone and dried in air.
Yield: 0.15 g, (89%). Mp: 195°C; FT-IR (KBr; cm–1): 1051 (m), 1120 (w), 1183 (m), 1240 (m), 1294 (s), 1326 (w), 1367 (m), 1408 (w), 1479 (w), 1610 (s), 1709 (s), 2396 (b), 2922 (w), 3014 (w), 3069 (w), 3426 (m), 3500 (w), 3615 (w), 3660 (w). Elemental Analysis calculated for C16H16N4O4(
Preparation of [(3-bpdh)(succinic acid)] (2)
This compound was prepared by mixing 1,5 bis(3-pyridyl)-3,4-diaza-2,4-hexadiene (3-bpdh) (0.12 g, 0.5 mmol) and succinic acid (0.059 g, 0.5 mmol). The mixture was ground with a mortar and pestle and kneaded, during which ethanol (5 mL) was added dropwise to obtain a pasty solution for 1 h at room temperature. The yellowish crystals of cocrystal
Yield: 0.15 g, (87%). Mp: 157°C; 1H NMR (500 MHz, DMSO-d6, ppm): 2.3 (s, 6 H), 2.42 (s, 4 H), 7.49–7.52 (d.q, 2 H), 8.26–8.29. (d.t, 2 H), 8.65–8.67 (d.d, 2 H), 9.1 (d, 2 H), 12.17 (s, 2 H); FT-IR (KBr; cm–1): 1043 (m), 1078 (w), 1122 (m), 1182 (m), 1329 (s), 1368 (w), 1411 (m), 1490 (w), 1608 (s), 1704 (m), 2362 (b), 2936 (w), 3068 (w), 3183 (m), 3500 (w), 3615 (w), 3660 (w). Elemental Analysis Calculated for C18H20N4O4(
Crystal structure determination
Crystallographic measurements were made at 293(2) K for [(3-bpdh)(oxalic acid)] (
Computational
Geometrical characterization of the synthesized compounds were performed in the gas phase at the level of density functional theory (DFT) using Becke’s three parameter exchange-functional (B3) combined with gradient-corrected correlation functional of Lee, Yang and Parr (LYP) (B3LYP) at 6–31G(d) level of theory. Starting geometries for calculation were taken from the X-ray crystal structure of the related cocrystals. To validate the optimization of the structures, frequency calculations were performed and the results showed no negative (imaginary) frequencies.
Results
Spectroscopic characterization
The infrared spectra of these two structures (solid in KBr pellets) show some harmonic vibrational frequencies in the range of 400–4000 cm–1. The vibrational analysis enables the characterization of a given molecular structure. The peak of the C = N and C = O groups are observed as strong bands centered at 1610 and 1709 cm–1 for [(3-bpdh)(oxalic acid)](
The relatively medium absorption band at 2922 cm–1 and 2936 cm–1 is due to the C–H modes of compounds
The 1H-NMR spectrum in DMSO-d6 show characteristic signals at δ= 7.49, 8.26, 8.65 and 9.11 ppm for proton in the aromatic ring. Single broad peak at δ= 12.17 ppm related to O–H acid and the other signals at 2.31 and 2.42 ppm correspond to CH3 and CH2 groups of 3-bpdh and succinic acid, respectively.
Crystal structure
Cocrystal

a) ORTEP plot b) a perspective view of hydrogen bonded 1D assembly c) 2D view assembly of [(3-bpdh)(oxalic acid)] (
Crystal data and structure refinement for compounds [(3-bpdh)(oxalic acid)]
Weak hydrogen bonding between hydrogen atoms of aryl rings and oxygen atoms of oxalic acid (C (Aryl) –H ... O) form a two dimensional (2D) structure (Fig. 1c). As (Fig. 2a) shows; two other weak hydrogen bonding between hydrogen atoms of methyl groups and oxygen atoms of oxalic acid from adjacent layers (C (Me) –H(12) ... O(4)) and (C (Me) –H(6) ... O(3)) form a 3D structure of cocrystal

3D views of the cocrystal assemblies of cocrystal a) [(3-bpdh)(oxalic acid)] (
Cocrystal (

a) ORTEP plot b) a perspective view of hydrogen bonded 1D assembly c) 2D assembly view of [(3-bpdh)(succinic acid)] (
As (Fig. 2b) shows, some weak hydrogen bonding is observed between methyl group and oxygen atom of two adjacent layers (C (Me) –H ... O). These interactions lead to form three dimensional (3D) supramolecular structure of cocrystal
Hydrogen bond length (Å)/angel (°) in the structure of [(3-bpdh)(oxalic acid)]
The comparison between X-ray powder diffraction (XRD) patterns of the simulated X-ray crystallography and experimental XRD pattern was done and good matches were observed (Fig. 4a, b).

XRD Pattern of a1) simulated pattern based on single crystal data, a2) synthesized powder of [(3-bpdh)(oxalic acid)] (
The theoretical studies enable the characterization of these molecular cocrystal structures

Optimized structure of cocrystal structures a)[(3-bpdh)(oxalic acid)] (
Experimental (X-ray) and calculated [RB3LYP/6–31 ... G(d)] selected bond lengths (Å), angles (°) and torsion angels (°) for [(DABTZ)(bpo)]
MESP displays three-dimensional visual regions for the electrophilic attack of charged point-like reagents on molecules and chemical reactivity of them [29].
The MESP diagrams (Fig. 6a, b) for these cocrystals were obtained based on the DFT optimized results. As these diagrams show the electrophilic sites (negative electrostatic potential in red color) such as carboxylic groups make hydrogen bond acceptor while electron deficient groups (positive electrostatic potential in blue color) such as nitrogen atoms in 3-bpdh makes hydrogen bond donor regions.

Molecular electrostatic potential map of cocrystals a)[(3-bpdh)(oxalic acid)] (
Effectually, MEP contour map provides a simple way to predict how these cocrystals could interact, geometrically. The binding of an attacking molecule such as drug to a target depends on the variation in electrostatic potential of them. As depicted in (Fig. 6) the negative regions are partial on the C = O double bonds of acids and nitrogen atoms of pyridine rings in 3-bpdh ligand, and positive ones on O–H bond of carboxylic groups.
The thermal decomposition behavior of these two new cocrystals [(3-bpdh)(oxalic acid)] (

TGA and DTA diagrams of two cocrystals a) [(3-bpdh)(oxalic acid)] (
Thermal decomposition of cocrystal
All of the DTA curves were merged into one picture and illustrated in (Fig. 8). DTA diagram of oxalic acid shows two endothermic process at 102 and 200C (Fig. 8a1). These peaks and another endothermic peak of 3-bpdh ligand at 97°C (Fig. 8a2) were removed after co-crystallization. Results revealed that the peaks of cocrystal

DTA diagrams of a1) oxalic acid; b1) succinic acid; a2 and b2) 3-bpdh; a3) [(3-bpdh)(oxalic acid)] (
The succinic acid (Fig. 8b1) and 3-bpdh ligand (Fig. 8b2) exhibited endothermic peaks at 190, 250 and 97°C, respectively.
The new peaks in the DTA curves of cocrystal
Herein, we report the synthesize of two new 1:1 stoichiometry ratio cocrystals of 3-bpdh as hydrogen acceptor with two distinct simple acids. Hydrogen bonding controls the packing of these two cocrystals. They are characterized and formulated as [(3-bpdh)(oxalic acid)] (1) and [(3-bpdh)(succinic acid)] (2) according to the obtained results from theoretical and experimental studies (DFT, elemental, thermal, 1H-NMR, FTIR spectroscopies, and single X-ray crystallography). It is gratifying to see that there is good agreement between experimental and theoretical results.
Footnotes
Acknowledgments
Support of this investigation by the Iran National Science Foundation is gratefully acknowledged.
CCDC 1437628 and 1437645 contain the supplementary crystallographic data for [(3-bpdh)(oxalic acid)] (
.
