<title>Introduction</title> <p>Bonding in organic compounds and nuances inherent in crystal packing engender boundless diversity in the arrangements of hydrogen bonded organic solid states. Due to its relative strength and its directional nature, the hydrogen bond has drawn much attention as a structural element in the design of crystalline phases [##UREF##0##1##–##UREF##3##5##]. The paradox that hydrogen bonding is important as both an element of structural diversity and design vanishes given that diversity depends on maximizing the number of hydrogen bonding options available to molecules whereas design focuses on controlling the direction and minimizing the number of hydrogen bonds. The current approach to probing relationships between molecular structure and packing involves substitutional modifications to an interesting parent molecule with limitations on hydrogen bond structural diversity [##UREF##4##6##–##UREF##6##8##].</p> <p>Hydrogen bonds optimally positioned, <bold>1</bold> (##FIG##0##Figure 1##), can lead to infinite polymers (tape) [##UREF##7##9##–##UREF##8##10##], cyclic n-mers [##UREF##1##2##,##REF##12224954##11##–##REF##12431093##12##], or dimers. For example, imidazolecarboxamidines <bold>2</bold> and imidazoles <bold>3</bold> have similar hydrogen bonding options offering sp<sup>2</sup>-NH hydrogen bond donors and sp<sup>2</sup>-N atom hydrogen bond acceptors. However <bold>3</bold> directs the hydrogen bond donor and acceptor approximately linearly [##UREF##9##13##] versus the ~90° intramolecular dihedral angle in <bold>2</bold>. Large angles between hydrogen bond donor and acceptor, as in <bold>3</bold>, should predispose solid state tapes [##UREF##9##13##], whereas small angles should prefer dimers. Planar motifs allow the favorable linearity in the hydrogen bond angles [##UREF##10##14##–##UREF##11##16##]. Intermediate dihedral angles ought to straddle the two crystalline phase motifs in the production of rings larger than dimers. The hydrogen bond motifs of amides encapsulate this concept in that tape and dimer dominate the solid state [##UREF##2##3##,##UREF##6##8##,##UREF##12##17##]. Crystalline phase hydrogen bonding near the parametric tipping point [##UREF##13##18##] between 0° and 180° might likely afford a variety of hydrogen bonded motifs.</p> <p>To simplify the interplay between directionality and the motif of the crystalline phase n-mer we studied a molecule with one hydrogen bond donor and one acceptor. Motif diversity increases sharply with more donors or acceptors [##UREF##1##2##]. The intermolecular hydrogen bond between the amidine sp<sup>2</sup>-NH and the imidazole sp<sup>2</sup>-N was maintained in all crystalline phases examined. Apparently the sp<sup>2</sup>-N atom in the imidazole accepts hydrogen bonds better than the sp<sup>2</sup>-N atom in the amidine group, and the non-involvement of amidine, as the weaker electron donor, agrees with previous studies of competitive solid state hydrogen bonding [##UREF##14##19##]. Rotamer/tautomer <bold>4</bold> was not observed in the crystalline phases, providing further control and predictability.</p> <p>Crystallization requires non-equilibrium conditions to progress [##REF##15701033##20##–##UREF##15##21##]; however, predictions are usually modeled based on thermodynamic considerations. With the notion that packing enthalpy mandates the crystalline phase [##UREF##16##22##], the question asked by this study was: how do the stabilities of the hydrogen bonded n-mers compare energetically in the absence of packing?</p>
<title>Results and Discussion</title> <p>The molecules in ##TAB##0##Table 1## were synthesized by combining imidazoles with commercially available carbodiimides as in ##FIG##1##Figure 2##. Even though the synthesis is easy, these molecules are very rare in the chemical literature. The products were crystallized under various conditions. In one of 21 syntheses, (R<sup>1</sup> = NH<sub>2</sub>) a more complex molecule than the carboxamidine was isolated due to the inclusion of two carbodiimide moieties in the product; see ##SUPPL##0##Supporting Information File 1##. Two other crystalline phases, with R<sup>1</sup> = NH<sub>2</sub>, favored hydrogen bond tapes with the involvement of R<sup>1 </sup>hydrogen bonding. We spent little time on these structures; they are not included in the current study. However, the preparations of these are included in the ##SUPPL##0##Supporting Information File 1##. Hydrogen bonding in the remaining substances was categorized into four groups: cyclic <italic>C</italic><sub>i</sub> dimers, cyclic <italic>C</italic><sub>1</sub> dimers, cyclic n-mers and infinite tape. Without the addition of a hydrogen bond donor for R<sup>1</sup>, the tape motif appears to be unlikely. Some effort was made to find polymorphs. Vials of material were arrayed in a variety of solvents and the unit cells were indexed. The same solvent-free crystal structures or crystals not suitable for diffraction (disordered or too small) were obtained.</p> <p>Hydrogen bonded dimers possessing inversion centers, <italic>C</italic><sub>i</sub> dimers, comprised the most popular bonding motif found in the solid state imidazolecarboxamidines. To investigate the nature of this energetic preference we compared by computation the gas-phase stabilities of these dimers versus the other observed hydrogen bonded motifs. Counter to intuition, in the absence of packing interactions the <italic>C</italic><sub>i</sub> dimers in their crystalline phase atomic coordinates were calculated to be considerably less stable than the less popular structural motifs.</p> <p>Literature on organic solid states contains much dialog regarding the minimization of Z' (molecules in the asymmetric unit) by associating structural or conformational enantiomers. Wallach's hypothesis foreshadowed this dialog: the racemic crystalline phase is more dense and more stable than the analogous optically pure crystalline phase [##UREF##17##23##]; although, exceptions have been noted [##UREF##18##24##]. Musing about this issue, Brock and Dunitz state, “Inversion centers are especially favorable for crystal packing because they diminish like-like interactions and are uniquely compatible with translation.” [##UREF##19##25##] Symmetry is a powerful component in packing: 83% of the entries in the Cambridge Structural Database that do not symmetry-relate molecules possess pseudosymmetry within 0.5 Å [##UREF##20##26##]. Molecules capable of either chiral or achiral space groups prefer the latter with concomitant minimization of Z'. The infrequency of chiral space groups in the CSD (~1:9) [##UREF##21##27##] may manifest a bias in the data toward small Z' possessing inversion symmetry.</p> <p>Packing could select a particular motif because dispersion forces factor in the construction of the organic solid state [##REF##26626389##28##–##UREF##22##29##]. If modern incantations of Wallach's hypothesis apply, the preferred solid state motif of <bold>2</bold> is likely the <italic>C</italic><sub>i</sub> dimer and structures that successfully compete with the <italic>C</italic><sub>i</sub> dimer should have increased stabilities from identifiable atomic parameters.</p> <p>The dihedral angles, θ, between the imidazole and the amidine moieties characterize the solid-state conformation and are reported in ##TAB##0##Table 1##. This parameter is defined in ##FIG##2##Figure 3## by the amidine N, C atoms and imidazole N, C atoms. With all else equal, molecules with θ of equal value but of opposite sign are conformational enantiomers.</p> <p>Even though only five structures did <italic>not</italic> crystallize as hydrogen bonded <italic>C</italic><sub>i</sub> dimers, comparing the atomic parameters of these to the atomic parameters of the <italic>C</italic><sub>i</sub> dimers is instructive. The directionality of the hydrogen bonds in this family of molecules approximates the hypothetical tipping point between atomic parameters that favor the infinite hydrogen bonded tape motif and the dimer. Despite the fact that the <italic>C</italic><sub>i</sub> dimer was the most common, the <italic>C</italic><sub>1</sub> dimer, trimer, tetramer and tape are calculated below to have more stable hydrogen bonding. The <italic>C</italic><sub>i</sub> dimer also tended to have the calculated least stable θ dihedral angles. Compensative packing must render the <italic>C</italic><sub>i</sub> dimer competitive.</p> <p>Cursory examination of molecular models shows that the imidazole moieties could stack with the R<sup>1</sup> substituents pointing either in the same or, as in structure <bold>5c</bold> in ##FIG##3##Figure 4##, in opposite directions. R<sup>1</sup> substituents pointing in opposite directions were the most popular, occurring in 12 of the 14 dimers. In light of Wallach's rule, an obvious advantage of this arrangement is the possibility that the dimeric units possess an inversion center and afford the molecules the assumed advantage of pairing two conformational enantiomers. The molecules that crystallized as dimers of conformational enantiomers with Z' = 1 had +/− pairs of single valued θ that varied between absolute values of 54° and 114°.</p> <p>The only polymorph found in this study, dimer <bold>7c2</bold>, nearly missed the inversion center (root mean square difference from a perfect inversion center of the C and N atoms = 0.03 Å, θ = 91.0°, θ' = −91.8°, Z' = 2). This is a common situation [##UREF##20##26##]. For perspective on this result, <italic>C</italic><sub>1</sub> dimers <bold>8a</bold> and <bold>8b</bold> missed the inversion center by RMS differences ~8.0 Å, a distance similar to the dimensions of the unit cell. The space groups encountered in this study were without exception achiral so the dimers with θ values: 91.8 and −91.0 were also present in the crystalline phase of <bold>7c2</bold>. Due to its structural proximity to <italic>C</italic><sub>i</sub>, for the purpose of taxonomy, <bold>7c2</bold> was classed as a <italic>C</italic><sub>i</sub> dimer. Differences between <bold>7c2</bold> and polymorph <bold>7c1</bold>, a true <italic>C</italic><sub>i</sub> dimer, are discussed below.</p> <p>When the dimer crystallized with the two R<sup>1</sup> substituents pointing in the same direction, θ and θ' within the dimer had different values of the same sign. Nature did not use a <italic>C</italic><sub>2</sub> operation to symmetrize these values. Only two molecules, <bold>8a</bold> and <bold>8b</bold> crystallized as <italic>C</italic><sub>1</sub> dimers, thus limiting any generalizations about the range of θ in these cases. The <italic>C</italic><sub>1</sub> dimers paired their aromatic substituents at R<sup>2</sup> in <bold>8a</bold> and <bold>8b</bold> and at R<sup>3</sup> in <bold>8b</bold>. Optimizing π-stacking, hydrogen bonding and θ likely allowed these two <italic>C</italic><sub>1</sub> dimer solid states in lieu of the otherwise ubiquitous <italic>C</italic><sub>i</sub> dimer.</p> <p>There were three other imidazolecarboxamidines in this study that did not crystallize as dimers: trimer <bold>6b</bold>, tetramer <bold>9a</bold> and tape <bold>11a</bold> (see ##TAB##0##Table 1## and ##FIG##4##Figure 5##–##FIG##6##Figure 7##). Like the <italic>C</italic><sub>1</sub> dimers <bold>8a</bold> and <bold>8b</bold>, structures <bold>6b</bold>, <bold>9a</bold> and <bold>11a</bold> tended to possess more stable calculated θ angles and hydrogen bonds than those found in the <italic>C</italic><sub>i</sub> dimers. Structures <bold>6b</bold> and <bold>9a</bold> are interesting in their putative ontological relationship to the dimers. One molecule in the trimer asymmetric unit was disordered. The two ordered molecules were analogous to an open <italic>C</italic><sub>i</sub> dimer with θ angles numerically close but of opposite sign. The best solution of the disorder modeled two molecules with large θ of opposite signs with unequal levels of occupancy. Tetramer <bold>9a</bold> does not suffer from this ambiguity; it is approximately an open dimer of <italic>C</italic><sub>i</sub> dimers with large θ of alternating sign. The molecules in the trimer and tetramer are unrelated by symmetry; Z' = 3 and 4 respectively.</p> <p>A comparison of the calculated stabilities of the sets of hydrogen-bound n-mers to determine if any hydrogen bonded motif enjoyed an energetic advantage follows in the text below. Two approaches based on the principle of practical model chemistries [##UREF##23##30##] queried the stabilities of the crystalline phase n-mers in the absence of packing.</p> <p>The potential energy of the imidazolecarboxamidine as a function of θ (##FIG##7##Figure 8##) was investigated by performing relaxed scans with <italic>Gaussian (G03)</italic> [##UREF##24##31##] at rhf/6-311+g(d,p) on hypothetical model monomer <bold>2a</bold> (R<sup>1</sup>, R<sup>2</sup> = H and R<sup>3</sup> = CH<sub>3</sub>). ##FIG##7##Figure 8## plots the potential energy of <bold>2a</bold> as a function of θ. The experimental crystal structure θ parameters of the molecules in ##TAB##0##Table 1##, categorized by hydrogen bond motif are included on the graph. At θ angles near 0° or 180° steric factors should increase the energies associated with θ and bring into the question whether using <bold>2</bold> as a model for the θ energy in all structures is useful. However, from θ = 50–130° a variety of steric environments are present in <bold>5</bold>–<bold>11</bold>. For example the steric nature of ##TAB##0##Table 1## entries <bold>11</bold> and <bold>7</bold> could possibly constrain θ to ~90°, but they do not; θ for <bold>7</bold> is near 90°, but θ for <bold>11</bold> is relatively small. Independent of sterics, the <italic>C</italic><sub>i</sub> dimers on average clustered in the high-energy area around θ ~90–100° in ##FIG##7##Figure 8## whereas the other crystalline phases tended to have θ parameters associated with lower energies corresponding to the gains in π bonding as θ approached planarity. The dashed arrow shows how the energetic content associated with θ of <bold>7c2</bold> changed when switched to the true <italic>C</italic><sub>i</sub> polymorph, <bold>7c1</bold>. The red icons represent <italic>C</italic><sub>i</sub> dimers with θ of exceptional stability and a parameter in a <italic>C</italic><sub>1</sub> dimer that is relatively unstable; these are discussed below.</p> <p>A more holistic calculation that allowed gross comparisons of the stabilities of all solid-state n-mers in the absence of packing interactions yielded paradoxical conclusions similar to the preceding calculations presented in ##FIG##7##Figure 8##. These calculations are more holistic in that more than one parameter is the focus of the calculation and the results are paradoxical because the more popular motif, the <italic>C</italic><sub>i</sub> dimer, is again calculated to be less stable.</p> <p>The steps of this calculation are a bit complex; a flow chart is presented in ##FIG##8##Figure 9##. Step 1: crystallographic information files (cif) were written as <italic>Gaussian (G03)</italic> input files, thus removing the material from the crystalline phase and bringing it into the gas-phase. Step 2: The atoms corresponding to R<sup>1</sup>, R<sup>2</sup> and R<sup>3</sup> were replaced with hydrogen atoms while preserving the relative positions of the remaining heavy atoms, this gives a set of structures corresponding to <bold>2b</bold>, R<sup>1</sup>, R<sup>2</sup>, R<sup>3</sup>, = H, that differ only in hydrogen bonded motif and atomic position. For the size of the molecules under study, accurately calculating dispersion forces in the clusters would have entailed an unreasonable high level of theory [##UREF##25##32##]. Step 3: The NHN hydrogen bond lengths, dihedral angle θ, four inter-imidazole-ring bond angles and one Cartesian coordinate per molecule were frozen. The remaining atomic parameters were optimized at the rhf/6-311+g(d,p) level of theory. It is important to optimize the C-H and N-H bond lengths to remove crystallographic errors generated by the algorithmic assignment of H atom positions. Step 4: Assurance that the solid state coordinates were not severely perturbed by optimization was gained from C- and N-atom RMS differences between the X-ray structure coordinates and those of the corresponding optimized structures; the RMS differences were calculated using gOpenMol. An RMS difference value of 0.04 Å was tolerated (0.02 Å average RMS difference). The two high values near 0.035 Å were not consequential. Step 5: Basis set superposition error (BSSE) biases the calculation of the hydrogen bond energies [##UREF##26##33##–##UREF##27##34##]; therefore, counterpoise correction was applied. The per-molecule, mostly-strong-local-energetic contributions to the stabilities of the n-mers were accessible by simply dividing the energies from these calculations by n. This calculation should include the effects of θ and hydrogen bonding.</p> <p>Calculating the analogous stability of the linear hydrogen bonded tape, <bold>11a</bold>, was more complicated. The corresponding linear dimer, trimer, tetramer and pentamer of <bold>11a</bold> were subjected to the above method. From the slope of the energy/n vs. n relationship, the per-molecule energy of the linear hydrogen bonded tape was calculated. The effect of the non-hydrogen bonded termini was further diminished by extrapolating the curve to n = 100.</p> <p>Since all values of n-mer/n are associated with the same molecular formula <bold>2b</bold>, the energies per molecule allow fair comparison of the energies due to the pi-energy effect of θ and hydrogen bonding in the n-mers in the absence of packing. ##FIG##9##Figure 10## summarizes the results of 20 calculations of this type by graphing the calculated, gas-phase, n-mer/n energies against the packing densities (molecular mass ⋉ Z/cell volume).</p> <p>With all else equal, packing density should correlate only grossly with solid state intermolecular interactions [##UREF##16##22##,##UREF##18##24##,##UREF##20##26##,##UREF##28##35##]. Stark differences in hydrogen bonding should enhance this correlation because the strength and directionality of hydrogen bonding can perturb packing in organic solids. Kitaigorodskii posits that organic molecules in crystalline phases fill space nearly as efficiently as close-packed spheres ~0.74 [##UREF##28##35##]. This occurs when the dimples and bumps of one molecule spatially correspond with the bumps and dimples of a lattice mate. Structures reliant on hydrogen bonds could violate this general rule by decreasing the packing coefficient due to the directionality of hydrogen bonds [##UREF##29##36##]. ##FIG##9##Figure 10## shows that the structures with aromatic substituents (X's in the graph) tended to be denser. In general, aromatic organics are denser than aliphatic organics [##UREF##30##37##]. This is likely due to the fact that bonds are shorter for sp<sup>2</sup>/sp<sup>2</sup> atoms than for sp<sup>3</sup>/sp<sup>3</sup> atoms. The effect of aromaticity on density is likely enhanced because these double aromatic substituents in these small structures accounted for much of the molecular mass.</p> <p>In ##FIG##9##Figure 10##, calculations again find that the <italic>C</italic><sub>i</sub> dimers (frozen crystalline phase coordinates) are least stable; the circles are all high on the Y axis in ##FIG##9##Figure 10##. There are four points in ##FIG##9##Figure 10## that contradict the trend, one high-energy <italic>C</italic><sub>1</sub> structure and three low-energy <italic>C</italic><sub>i</sub> structures.</p> <p>One point in ##FIG##9##Figure 10## that contradicts the hypothesis that lattice-free <italic>C</italic><sub>i</sub> dimers are least stable is the red diamond corresponding to <bold>8b</bold>, an unstable <italic>C</italic><sub>1</sub> dimer. However, this molecule has the most aromatic groups and is the densest non-<italic>C</italic><sub>i</sub> dimer. The four aromatic groups in <bold>8b</bold> interact extensively which is readily apparent upon examination of the packing. The method of the calculations summarized by ##FIG##9##Figure 10## replaced the aromatic substituents with hydrogen atoms. The difference between <bold>8a</bold> and <bold>8b</bold> is <italic>i</italic>Pr versus Tol at the amidine N atoms (R<sup>3</sup>). Perhaps surface area-dependent dispersive interactions in the nucleation process of C<sub>1</sub> dimer, <bold>8b</bold> perturbed hydrogen bonding away from optimum. Aromatic stacking is quoted anywhere between 2 and 0.5 kcal/mol so a scenario in which the eight aromatic-interactions in dimer <bold>8b</bold> perturbed the energies of the hydrogen bonds is very reasonable.</p> <p>The red circles in ##FIG##9##Figure 10## represent <italic>C</italic><sub>i</sub> dimer <bold>9b1</bold>–<bold>9b3</bold> cocrystallized 1:1 with EtOAc, Toluene, and Et<sub>2</sub>O respectively. The three red out-of-place circles to the left of the graph in ##FIG##7##Figure 8## also belong to <bold>9b1</bold>–<bold>9b3</bold>. Molecule <bold>9b</bold> was the only one in this study to crystallize with solvent. Further attempts to obtain <bold>9b</bold> solvent-free resulted in 1:1 inclusions of CH<sub>3</sub>CN, isopropyl ether, and chlorobenzene which were not analyzed completely. Solvent appears in ~15% of neutral organics in the CSD and has been attributed to interrupted crystallization processes [##UREF##31##38##]. The presence of solvent molecules in <bold>9b</bold> neatly allows rationalization of the high densities of <bold>9b1</bold>–<bold>9b3</bold> and the low, per molecule, gas-phase energies, and the low-energy θ parameters. From these observations we make the following hypothesis. <italic>The solution states of</italic>
<bold><italic>9b</italic></bold>
<italic>did not completely desolvate during the nucleation process thereby maintaining near solution-state conformation. This allowed the low calculated energies, and allowed the low-energy θ parameters. Disordered solvent molecules in the crystalline phase of</italic>
<bold><italic>9</italic></bold>
<italic>occupied voids at the faces of the two hydrogen-bound imidazole rings thereby increasing the densities of</italic>
<bold><italic>9</italic></bold>
<italic>and preserving near solution-state conformation for the C</italic><italic><sub>i</sub></italic><italic> dimer of </italic><bold><italic>9</italic></bold>.</p> <p>More significance of calculations of θ and of the hydrogen bond motifs was attained by plotting the energies in ##FIG##9##Figure 10## against the energies in ##FIG##7##Figure 8##. If all the points in ##FIG##10##Figure 11## had fallen on a diagonal line, the energies in ##FIG##9##Figure 10## would have been ascribable to the θ parameter of ##FIG##7##Figure 8## with no motif-dependent, energetic differences from intermolecular hydrogen bonding. However, ##FIG##10##Figure 11## shows that tetramer <bold>9a</bold> (triangles) optimized hydrogen bonds better (smaller X axis values) than trimer <bold>6b</bold> (open squares) even though the tetramer had less optimum θ parameters (larger average Y axis values) than the trimer. The sign of θ alternates around the tetrameric ring, thus point symmetry possibilities for the reduction of Z' in this structure were <italic>C</italic><sub>2</sub> and <italic>S</italic><sub>4</sub>, but not <italic>C</italic><sub>i</sub>. Axiomatically, wide XH-Y angles stabilize hydrogen bonds [##UREF##32##39##]; crystal structures prefer ~linear hydrogen bonds like those in ##FIG##5##Figure 6## for tetramer <bold>9a</bold>, over the less linear hydrogen bonds in <bold>5c</bold> and <bold>11a</bold> in ##FIG##4##Figure 5## and ##FIG##6##Figure 7## [##UREF##10##14##–##UREF##11##16##].</p> <p>##FIG##10##Figure 11## also supports the hypothesis that inter-dimer dispersive forces edited the optimum solution state conformation and n-mer molecularity. The low-energy <italic>C</italic><sub>i</sub> dimers <bold>9b1</bold>–<bold>3</bold> (red circles in ##FIG##10##Figure 11##) corresponded to the <italic>C</italic><sub>i</sub> dimers with the lowest-energy θ angles in ##FIG##7##Figure 8##. Little perturbation of the solution state occurred upon nucleation because EtOAc, toluene and ether solvent molecules filled the voids and preserved the optimum solution-state conformation. In each case, the lacuna in the lattices occupied by solvent allows for much disorder on the part of the solvent. These three points really do not at all contradict the hypothesis that bringing <italic>C</italic><sub>i</sub> symmetric dimers of <bold>2</bold> from solution into the solid state results in low-energy packing and high-energy local interactions.</p> <p>The density of the crystalline phase of <bold>7c</bold> decreased in going from truly <italic>C</italic><sub>i</sub> symmetric <bold>7c1</bold> to near <italic>C</italic><sub>i</sub> symmetric <bold>7c2</bold> in ##FIG##9##Figure 10##. ##FIG##10##Figure 11## analyzes this change in terms of hydrogen bonding and θ-derived energies. Attaining <italic>C</italic><sub>i</sub> symmetry is accompanied by stabilization of θ-derived energies (~0.74 kcal/mol) and destabilization of hydrogen bonding (0.80 kcal/mol). Within error these effects all but cancel and this result is consonant with the view that polymorphism is the result of a subtle balance of orthogonal forces associated with certain atomic parameters [##UREF##33##40##]. A relatively large increase in density accompanied the attainment of true <italic>C</italic><sub>i</sub> symmetry. The caveat here is of course that only one polymorph was found and thus generalizations will have to wait for a study of another system.</p> <p>The <italic>C</italic><sub>i</sub> dimers in general had less stable hydrogen bonds and non-optimum θ angles which put them in the upper right corner of the graph in ##FIG##10##Figure 11##. The <italic>C</italic><sub>1</sub> dimer also suffers due to high-energy hydrogen bonds (right side of ##FIG##10##Figure 11##) but this dimer can optimize one of the two structures in the asymmetric unit as in <bold>8b</bold>, or both fairly well as in <bold>8a</bold>, thereby minimizing the average local interaction energy. The high-energy molecular component in <italic>C</italic><sub>1</sub> dimer <bold>8b</bold> (high-energy red diamond in ##FIG##10##Figure 11##) corresponded to the highest-energy θ angle in ##FIG##7##Figure 8## that lost π resonance between the amidine and the imidazole moiety. As discussed earlier, this structure is likely a result of extensive π-stacking.</p>
<title>Results and Discussion</title> <p>The molecules in ##TAB##0##Table 1## were synthesized by combining imidazoles with commercially available carbodiimides as in ##FIG##1##Figure 2##. Even though the synthesis is easy, these molecules are very rare in the chemical literature. The products were crystallized under various conditions. In one of 21 syntheses, (R<sup>1</sup> = NH<sub>2</sub>) a more complex molecule than the carboxamidine was isolated due to the inclusion of two carbodiimide moieties in the product; see ##SUPPL##0##Supporting Information File 1##. Two other crystalline phases, with R<sup>1</sup> = NH<sub>2</sub>, favored hydrogen bond tapes with the involvement of R<sup>1 </sup>hydrogen bonding. We spent little time on these structures; they are not included in the current study. However, the preparations of these are included in the ##SUPPL##0##Supporting Information File 1##. Hydrogen bonding in the remaining substances was categorized into four groups: cyclic <italic>C</italic><sub>i</sub> dimers, cyclic <italic>C</italic><sub>1</sub> dimers, cyclic n-mers and infinite tape. Without the addition of a hydrogen bond donor for R<sup>1</sup>, the tape motif appears to be unlikely. Some effort was made to find polymorphs. Vials of material were arrayed in a variety of solvents and the unit cells were indexed. The same solvent-free crystal structures or crystals not suitable for diffraction (disordered or too small) were obtained.</p> <p>Hydrogen bonded dimers possessing inversion centers, <italic>C</italic><sub>i</sub> dimers, comprised the most popular bonding motif found in the solid state imidazolecarboxamidines. To investigate the nature of this energetic preference we compared by computation the gas-phase stabilities of these dimers versus the other observed hydrogen bonded motifs. Counter to intuition, in the absence of packing interactions the <italic>C</italic><sub>i</sub> dimers in their crystalline phase atomic coordinates were calculated to be considerably less stable than the less popular structural motifs.</p> <p>Literature on organic solid states contains much dialog regarding the minimization of Z' (molecules in the asymmetric unit) by associating structural or conformational enantiomers. Wallach's hypothesis foreshadowed this dialog: the racemic crystalline phase is more dense and more stable than the analogous optically pure crystalline phase [##UREF##17##23##]; although, exceptions have been noted [##UREF##18##24##]. Musing about this issue, Brock and Dunitz state, “Inversion centers are especially favorable for crystal packing because they diminish like-like interactions and are uniquely compatible with translation.” [##UREF##19##25##] Symmetry is a powerful component in packing: 83% of the entries in the Cambridge Structural Database that do not symmetry-relate molecules possess pseudosymmetry within 0.5 Å [##UREF##20##26##]. Molecules capable of either chiral or achiral space groups prefer the latter with concomitant minimization of Z'. The infrequency of chiral space groups in the CSD (~1:9) [##UREF##21##27##] may manifest a bias in the data toward small Z' possessing inversion symmetry.</p> <p>Packing could select a particular motif because dispersion forces factor in the construction of the organic solid state [##REF##26626389##28##–##UREF##22##29##]. If modern incantations of Wallach's hypothesis apply, the preferred solid state motif of <bold>2</bold> is likely the <italic>C</italic><sub>i</sub> dimer and structures that successfully compete with the <italic>C</italic><sub>i</sub> dimer should have increased stabilities from identifiable atomic parameters.</p> <p>The dihedral angles, θ, between the imidazole and the amidine moieties characterize the solid-state conformation and are reported in ##TAB##0##Table 1##. This parameter is defined in ##FIG##2##Figure 3## by the amidine N, C atoms and imidazole N, C atoms. With all else equal, molecules with θ of equal value but of opposite sign are conformational enantiomers.</p> <p>Even though only five structures did <italic>not</italic> crystallize as hydrogen bonded <italic>C</italic><sub>i</sub> dimers, comparing the atomic parameters of these to the atomic parameters of the <italic>C</italic><sub>i</sub> dimers is instructive. The directionality of the hydrogen bonds in this family of molecules approximates the hypothetical tipping point between atomic parameters that favor the infinite hydrogen bonded tape motif and the dimer. Despite the fact that the <italic>C</italic><sub>i</sub> dimer was the most common, the <italic>C</italic><sub>1</sub> dimer, trimer, tetramer and tape are calculated below to have more stable hydrogen bonding. The <italic>C</italic><sub>i</sub> dimer also tended to have the calculated least stable θ dihedral angles. Compensative packing must render the <italic>C</italic><sub>i</sub> dimer competitive.</p> <p>Cursory examination of molecular models shows that the imidazole moieties could stack with the R<sup>1</sup> substituents pointing either in the same or, as in structure <bold>5c</bold> in ##FIG##3##Figure 4##, in opposite directions. R<sup>1</sup> substituents pointing in opposite directions were the most popular, occurring in 12 of the 14 dimers. In light of Wallach's rule, an obvious advantage of this arrangement is the possibility that the dimeric units possess an inversion center and afford the molecules the assumed advantage of pairing two conformational enantiomers. The molecules that crystallized as dimers of conformational enantiomers with Z' = 1 had +/− pairs of single valued θ that varied between absolute values of 54° and 114°.</p> <p>The only polymorph found in this study, dimer <bold>7c2</bold>, nearly missed the inversion center (root mean square difference from a perfect inversion center of the C and N atoms = 0.03 Å, θ = 91.0°, θ' = −91.8°, Z' = 2). This is a common situation [##UREF##20##26##]. For perspective on this result, <italic>C</italic><sub>1</sub> dimers <bold>8a</bold> and <bold>8b</bold> missed the inversion center by RMS differences ~8.0 Å, a distance similar to the dimensions of the unit cell. The space groups encountered in this study were without exception achiral so the dimers with θ values: 91.8 and −91.0 were also present in the crystalline phase of <bold>7c2</bold>. Due to its structural proximity to <italic>C</italic><sub>i</sub>, for the purpose of taxonomy, <bold>7c2</bold> was classed as a <italic>C</italic><sub>i</sub> dimer. Differences between <bold>7c2</bold> and polymorph <bold>7c1</bold>, a true <italic>C</italic><sub>i</sub> dimer, are discussed below.</p> <p>When the dimer crystallized with the two R<sup>1</sup> substituents pointing in the same direction, θ and θ' within the dimer had different values of the same sign. Nature did not use a <italic>C</italic><sub>2</sub> operation to symmetrize these values. Only two molecules, <bold>8a</bold> and <bold>8b</bold> crystallized as <italic>C</italic><sub>1</sub> dimers, thus limiting any generalizations about the range of θ in these cases. The <italic>C</italic><sub>1</sub> dimers paired their aromatic substituents at R<sup>2</sup> in <bold>8a</bold> and <bold>8b</bold> and at R<sup>3</sup> in <bold>8b</bold>. Optimizing π-stacking, hydrogen bonding and θ likely allowed these two <italic>C</italic><sub>1</sub> dimer solid states in lieu of the otherwise ubiquitous <italic>C</italic><sub>i</sub> dimer.</p> <p>There were three other imidazolecarboxamidines in this study that did not crystallize as dimers: trimer <bold>6b</bold>, tetramer <bold>9a</bold> and tape <bold>11a</bold> (see ##TAB##0##Table 1## and ##FIG##4##Figure 5##–##FIG##6##Figure 7##). Like the <italic>C</italic><sub>1</sub> dimers <bold>8a</bold> and <bold>8b</bold>, structures <bold>6b</bold>, <bold>9a</bold> and <bold>11a</bold> tended to possess more stable calculated θ angles and hydrogen bonds than those found in the <italic>C</italic><sub>i</sub> dimers. Structures <bold>6b</bold> and <bold>9a</bold> are interesting in their putative ontological relationship to the dimers. One molecule in the trimer asymmetric unit was disordered. The two ordered molecules were analogous to an open <italic>C</italic><sub>i</sub> dimer with θ angles numerically close but of opposite sign. The best solution of the disorder modeled two molecules with large θ of opposite signs with unequal levels of occupancy. Tetramer <bold>9a</bold> does not suffer from this ambiguity; it is approximately an open dimer of <italic>C</italic><sub>i</sub> dimers with large θ of alternating sign. The molecules in the trimer and tetramer are unrelated by symmetry; Z' = 3 and 4 respectively.</p> <p>A comparison of the calculated stabilities of the sets of hydrogen-bound n-mers to determine if any hydrogen bonded motif enjoyed an energetic advantage follows in the text below. Two approaches based on the principle of practical model chemistries [##UREF##23##30##] queried the stabilities of the crystalline phase n-mers in the absence of packing.</p> <p>The potential energy of the imidazolecarboxamidine as a function of θ (##FIG##7##Figure 8##) was investigated by performing relaxed scans with <italic>Gaussian (G03)</italic> [##UREF##24##31##] at rhf/6-311+g(d,p) on hypothetical model monomer <bold>2a</bold> (R<sup>1</sup>, R<sup>2</sup> = H and R<sup>3</sup> = CH<sub>3</sub>). ##FIG##7##Figure 8## plots the potential energy of <bold>2a</bold> as a function of θ. The experimental crystal structure θ parameters of the molecules in ##TAB##0##Table 1##, categorized by hydrogen bond motif are included on the graph. At θ angles near 0° or 180° steric factors should increase the energies associated with θ and bring into the question whether using <bold>2</bold> as a model for the θ energy in all structures is useful. However, from θ = 50–130° a variety of steric environments are present in <bold>5</bold>–<bold>11</bold>. For example the steric nature of ##TAB##0##Table 1## entries <bold>11</bold> and <bold>7</bold> could possibly constrain θ to ~90°, but they do not; θ for <bold>7</bold> is near 90°, but θ for <bold>11</bold> is relatively small. Independent of sterics, the <italic>C</italic><sub>i</sub> dimers on average clustered in the high-energy area around θ ~90–100° in ##FIG##7##Figure 8## whereas the other crystalline phases tended to have θ parameters associated with lower energies corresponding to the gains in π bonding as θ approached planarity. The dashed arrow shows how the energetic content associated with θ of <bold>7c2</bold> changed when switched to the true <italic>C</italic><sub>i</sub> polymorph, <bold>7c1</bold>. The red icons represent <italic>C</italic><sub>i</sub> dimers with θ of exceptional stability and a parameter in a <italic>C</italic><sub>1</sub> dimer that is relatively unstable; these are discussed below.</p> <p>A more holistic calculation that allowed gross comparisons of the stabilities of all solid-state n-mers in the absence of packing interactions yielded paradoxical conclusions similar to the preceding calculations presented in ##FIG##7##Figure 8##. These calculations are more holistic in that more than one parameter is the focus of the calculation and the results are paradoxical because the more popular motif, the <italic>C</italic><sub>i</sub> dimer, is again calculated to be less stable.</p> <p>The steps of this calculation are a bit complex; a flow chart is presented in ##FIG##8##Figure 9##. Step 1: crystallographic information files (cif) were written as <italic>Gaussian (G03)</italic> input files, thus removing the material from the crystalline phase and bringing it into the gas-phase. Step 2: The atoms corresponding to R<sup>1</sup>, R<sup>2</sup> and R<sup>3</sup> were replaced with hydrogen atoms while preserving the relative positions of the remaining heavy atoms, this gives a set of structures corresponding to <bold>2b</bold>, R<sup>1</sup>, R<sup>2</sup>, R<sup>3</sup>, = H, that differ only in hydrogen bonded motif and atomic position. For the size of the molecules under study, accurately calculating dispersion forces in the clusters would have entailed an unreasonable high level of theory [##UREF##25##32##]. Step 3: The NHN hydrogen bond lengths, dihedral angle θ, four inter-imidazole-ring bond angles and one Cartesian coordinate per molecule were frozen. The remaining atomic parameters were optimized at the rhf/6-311+g(d,p) level of theory. It is important to optimize the C-H and N-H bond lengths to remove crystallographic errors generated by the algorithmic assignment of H atom positions. Step 4: Assurance that the solid state coordinates were not severely perturbed by optimization was gained from C- and N-atom RMS differences between the X-ray structure coordinates and those of the corresponding optimized structures; the RMS differences were calculated using gOpenMol. An RMS difference value of 0.04 Å was tolerated (0.02 Å average RMS difference). The two high values near 0.035 Å were not consequential. Step 5: Basis set superposition error (BSSE) biases the calculation of the hydrogen bond energies [##UREF##26##33##–##UREF##27##34##]; therefore, counterpoise correction was applied. The per-molecule, mostly-strong-local-energetic contributions to the stabilities of the n-mers were accessible by simply dividing the energies from these calculations by n. This calculation should include the effects of θ and hydrogen bonding.</p> <p>Calculating the analogous stability of the linear hydrogen bonded tape, <bold>11a</bold>, was more complicated. The corresponding linear dimer, trimer, tetramer and pentamer of <bold>11a</bold> were subjected to the above method. From the slope of the energy/n vs. n relationship, the per-molecule energy of the linear hydrogen bonded tape was calculated. The effect of the non-hydrogen bonded termini was further diminished by extrapolating the curve to n = 100.</p> <p>Since all values of n-mer/n are associated with the same molecular formula <bold>2b</bold>, the energies per molecule allow fair comparison of the energies due to the pi-energy effect of θ and hydrogen bonding in the n-mers in the absence of packing. ##FIG##9##Figure 10## summarizes the results of 20 calculations of this type by graphing the calculated, gas-phase, n-mer/n energies against the packing densities (molecular mass ⋉ Z/cell volume).</p> <p>With all else equal, packing density should correlate only grossly with solid state intermolecular interactions [##UREF##16##22##,##UREF##18##24##,##UREF##20##26##,##UREF##28##35##]. Stark differences in hydrogen bonding should enhance this correlation because the strength and directionality of hydrogen bonding can perturb packing in organic solids. Kitaigorodskii posits that organic molecules in crystalline phases fill space nearly as efficiently as close-packed spheres ~0.74 [##UREF##28##35##]. This occurs when the dimples and bumps of one molecule spatially correspond with the bumps and dimples of a lattice mate. Structures reliant on hydrogen bonds could violate this general rule by decreasing the packing coefficient due to the directionality of hydrogen bonds [##UREF##29##36##]. ##FIG##9##Figure 10## shows that the structures with aromatic substituents (X's in the graph) tended to be denser. In general, aromatic organics are denser than aliphatic organics [##UREF##30##37##]. This is likely due to the fact that bonds are shorter for sp<sup>2</sup>/sp<sup>2</sup> atoms than for sp<sup>3</sup>/sp<sup>3</sup> atoms. The effect of aromaticity on density is likely enhanced because these double aromatic substituents in these small structures accounted for much of the molecular mass.</p> <p>In ##FIG##9##Figure 10##, calculations again find that the <italic>C</italic><sub>i</sub> dimers (frozen crystalline phase coordinates) are least stable; the circles are all high on the Y axis in ##FIG##9##Figure 10##. There are four points in ##FIG##9##Figure 10## that contradict the trend, one high-energy <italic>C</italic><sub>1</sub> structure and three low-energy <italic>C</italic><sub>i</sub> structures.</p> <p>One point in ##FIG##9##Figure 10## that contradicts the hypothesis that lattice-free <italic>C</italic><sub>i</sub> dimers are least stable is the red diamond corresponding to <bold>8b</bold>, an unstable <italic>C</italic><sub>1</sub> dimer. However, this molecule has the most aromatic groups and is the densest non-<italic>C</italic><sub>i</sub> dimer. The four aromatic groups in <bold>8b</bold> interact extensively which is readily apparent upon examination of the packing. The method of the calculations summarized by ##FIG##9##Figure 10## replaced the aromatic substituents with hydrogen atoms. The difference between <bold>8a</bold> and <bold>8b</bold> is <italic>i</italic>Pr versus Tol at the amidine N atoms (R<sup>3</sup>). Perhaps surface area-dependent dispersive interactions in the nucleation process of C<sub>1</sub> dimer, <bold>8b</bold> perturbed hydrogen bonding away from optimum. Aromatic stacking is quoted anywhere between 2 and 0.5 kcal/mol so a scenario in which the eight aromatic-interactions in dimer <bold>8b</bold> perturbed the energies of the hydrogen bonds is very reasonable.</p> <p>The red circles in ##FIG##9##Figure 10## represent <italic>C</italic><sub>i</sub> dimer <bold>9b1</bold>–<bold>9b3</bold> cocrystallized 1:1 with EtOAc, Toluene, and Et<sub>2</sub>O respectively. The three red out-of-place circles to the left of the graph in ##FIG##7##Figure 8## also belong to <bold>9b1</bold>–<bold>9b3</bold>. Molecule <bold>9b</bold> was the only one in this study to crystallize with solvent. Further attempts to obtain <bold>9b</bold> solvent-free resulted in 1:1 inclusions of CH<sub>3</sub>CN, isopropyl ether, and chlorobenzene which were not analyzed completely. Solvent appears in ~15% of neutral organics in the CSD and has been attributed to interrupted crystallization processes [##UREF##31##38##]. The presence of solvent molecules in <bold>9b</bold> neatly allows rationalization of the high densities of <bold>9b1</bold>–<bold>9b3</bold> and the low, per molecule, gas-phase energies, and the low-energy θ parameters. From these observations we make the following hypothesis. <italic>The solution states of</italic>
<bold><italic>9b</italic></bold>
<italic>did not completely desolvate during the nucleation process thereby maintaining near solution-state conformation. This allowed the low calculated energies, and allowed the low-energy θ parameters. Disordered solvent molecules in the crystalline phase of</italic>
<bold><italic>9</italic></bold>
<italic>occupied voids at the faces of the two hydrogen-bound imidazole rings thereby increasing the densities of</italic>
<bold><italic>9</italic></bold>
<italic>and preserving near solution-state conformation for the C</italic><italic><sub>i</sub></italic><italic> dimer of </italic><bold><italic>9</italic></bold>.</p> <p>More significance of calculations of θ and of the hydrogen bond motifs was attained by plotting the energies in ##FIG##9##Figure 10## against the energies in ##FIG##7##Figure 8##. If all the points in ##FIG##10##Figure 11## had fallen on a diagonal line, the energies in ##FIG##9##Figure 10## would have been ascribable to the θ parameter of ##FIG##7##Figure 8## with no motif-dependent, energetic differences from intermolecular hydrogen bonding. However, ##FIG##10##Figure 11## shows that tetramer <bold>9a</bold> (triangles) optimized hydrogen bonds better (smaller X axis values) than trimer <bold>6b</bold> (open squares) even though the tetramer had less optimum θ parameters (larger average Y axis values) than the trimer. The sign of θ alternates around the tetrameric ring, thus point symmetry possibilities for the reduction of Z' in this structure were <italic>C</italic><sub>2</sub> and <italic>S</italic><sub>4</sub>, but not <italic>C</italic><sub>i</sub>. Axiomatically, wide XH-Y angles stabilize hydrogen bonds [##UREF##32##39##]; crystal structures prefer ~linear hydrogen bonds like those in ##FIG##5##Figure 6## for tetramer <bold>9a</bold>, over the less linear hydrogen bonds in <bold>5c</bold> and <bold>11a</bold> in ##FIG##4##Figure 5## and ##FIG##6##Figure 7## [##UREF##10##14##–##UREF##11##16##].</p> <p>##FIG##10##Figure 11## also supports the hypothesis that inter-dimer dispersive forces edited the optimum solution state conformation and n-mer molecularity. The low-energy <italic>C</italic><sub>i</sub> dimers <bold>9b1</bold>–<bold>3</bold> (red circles in ##FIG##10##Figure 11##) corresponded to the <italic>C</italic><sub>i</sub> dimers with the lowest-energy θ angles in ##FIG##7##Figure 8##. Little perturbation of the solution state occurred upon nucleation because EtOAc, toluene and ether solvent molecules filled the voids and preserved the optimum solution-state conformation. In each case, the lacuna in the lattices occupied by solvent allows for much disorder on the part of the solvent. These three points really do not at all contradict the hypothesis that bringing <italic>C</italic><sub>i</sub> symmetric dimers of <bold>2</bold> from solution into the solid state results in low-energy packing and high-energy local interactions.</p> <p>The density of the crystalline phase of <bold>7c</bold> decreased in going from truly <italic>C</italic><sub>i</sub> symmetric <bold>7c1</bold> to near <italic>C</italic><sub>i</sub> symmetric <bold>7c2</bold> in ##FIG##9##Figure 10##. ##FIG##10##Figure 11## analyzes this change in terms of hydrogen bonding and θ-derived energies. Attaining <italic>C</italic><sub>i</sub> symmetry is accompanied by stabilization of θ-derived energies (~0.74 kcal/mol) and destabilization of hydrogen bonding (0.80 kcal/mol). Within error these effects all but cancel and this result is consonant with the view that polymorphism is the result of a subtle balance of orthogonal forces associated with certain atomic parameters [##UREF##33##40##]. A relatively large increase in density accompanied the attainment of true <italic>C</italic><sub>i</sub> symmetry. The caveat here is of course that only one polymorph was found and thus generalizations will have to wait for a study of another system.</p> <p>The <italic>C</italic><sub>i</sub> dimers in general had less stable hydrogen bonds and non-optimum θ angles which put them in the upper right corner of the graph in ##FIG##10##Figure 11##. The <italic>C</italic><sub>1</sub> dimer also suffers due to high-energy hydrogen bonds (right side of ##FIG##10##Figure 11##) but this dimer can optimize one of the two structures in the asymmetric unit as in <bold>8b</bold>, or both fairly well as in <bold>8a</bold>, thereby minimizing the average local interaction energy. The high-energy molecular component in <italic>C</italic><sub>1</sub> dimer <bold>8b</bold> (high-energy red diamond in ##FIG##10##Figure 11##) corresponded to the highest-energy θ angle in ##FIG##7##Figure 8## that lost π resonance between the amidine and the imidazole moiety. As discussed earlier, this structure is likely a result of extensive π-stacking.</p>
<title>Conclusion</title> <p>In a family of molecules this study found 15 crystalline phases in which the hydrogen bonded motif was <italic>C</italic><sub>i</sub> symmetric, Z' = 1 and five other structures in which the components of the hydrogen bonded motif were not related by improper symmetry. Other than the tape motif in which the hydrogen bonded components followed a screw axis, the other four structures had Z' ≥ 2. The packing of <italic>C</italic><sub>i</sub> symmetric pairs stabilized this most popular motif even though the <italic>C</italic><sub>i</sub> motif was destabilized relative to other motifs by conformational and hydrogen bonding effects. Our analysis of this small data set separated local versus dispersive contributions to stability. In Gavezzotti's statistical search of the CSD for energy-edited symmetry preferences in Z' = 2 vs. Z' = 1 structures, a relationship was found between solid state symmetry and the stabilities of pair-wise interactions [##UREF##20##26##]. In related work, Steed <italic>et al.</italic> found that stereogenic atoms included in normally <italic>C</italic><sub>i</sub> symmetric hydrogen bonded dimers, increases Z' from 1 to 2, indicating that rotational symmetry is less propitious in the minimization of Z' than inversion [##UREF##34##41##]. Very related to the current report is Wheeler's discovery that heterochiral isosteric molecules conserve the solid state inversion-symmetric motifs of their racemic analogues [##UREF##35##42##–##UREF##36##43##]. These results relate to Wallach's hypothesis: either packing interactions are optimal when <italic>C</italic><sub>i</sub>-symmetric units nucleate or <italic>C</italic><sub>i</sub>-symmetric interactions are more stable in solution and hence get included in the crystalline phase. Regarding less condensed states, optically pure gas phase methyl lactate favors the tetramer over dimers more than the racemic mixture [##UREF##37##44##]. The lack of the energetically competitive heterochiral dimer in the optically pure mixture could have produced that result.</p> <p>With all things equal, apparently packing prefers centro-symmetric pairs. Here, structures other than the <italic>C</italic><sub>i</sub> dimer required more stability from hydrogen bonding and conformation to compensate for non-optimal packing forces. Here, less dense crystalline phases resulted when packing forces yielded control of conformation and hydrogen bond motif to local, more directional forces.</p> <p>Should all molecules capable of hydrogen bonding build lattices of <italic>C</italic><sub>i</sub> symmetric aggregates with Z' = 1 at the expense of local interactions? No, local and dispersive forces can also synergize to construct the solid state. Crystal structure databases are likely mosaics containing molecular families with structural aspects that compete and synergize to varying degrees. More work has to be put toward a holistic understanding of interplay between packing and solid state synthons that are usually the function of strong local interactions [##UREF##38##45##]. However when there is competition between local interactions and the dispersive interactions, this work suggests that small <italic>C</italic><sub>i</sub> symmetric units have a slight thermodynamic packing advantage which could be the basis for Wallach's rule.</p> <p>The results bring into question predictive methods based on energy minimization and their level of accuracy, especially in the prediction of hydrogen bonding options that are proximal in energy [##REF##26626389##28##,##UREF##39##46##–##UREF##40##47##]. Prediction of the organic solid state is challenged by the fact that many crystalline phases likely result from marginal differences in large opposing effects. While the problem of calculating marginal differences in large energetic contributions to the organic solid state has been the subject of previous conjecture, this study is one example of the nature of the problem, unveiled and dissected.</p> <p>The results provide a caveat for mining crystal structure databases and translating structural popularity to energy. Such searches should be as structurally broad as possible. For example, tendencies of a particular dihedral angle to adopt a certain average value could be gleaned by looking at many crystal structures. However the applicability of the parameter outside the solid state may be severely and systematically diminished by symmetry-edited packing effects. The most common dihedral angle is the highest-energy dihedral angle from ##FIG##7##Figure 8## with a fairly large data set. This study detailed relevant interactions in how such an observation could occur. Similar statements can be made about the most popular hydrogen bond motif followed by caveats regarding the use of crystal-structure derived atomic parameters to broadly characterize hydrogen bonding energies.</p>
<p>Substitutional changes to imidazolecarboxamidine that preserved intermolecular hydrogen bonding in the solid state were used to study the relationship between packing and the hydrogen bond motif. Various motifs competed, but the most common imidazolecarboxamidine crystalline phase was a <italic>C</italic><sub>i</sub> symmetric dimer that established inversion centers by associating enantiomeric tautomers. Counter to intuition, the calculated gas-phase energies per molecule of the solid state atomic coordinates of the <italic>C</italic><sub>i</sub> dimer motifs were higher than those of the <italic>C</italic><sub>1</sub> dimer, trimer, tetramer and tape motifs, while the packing densities of <italic>C</italic><sub>i</sub> dimers were found to be higher. This result was interpreted as an enhanced ability of the <italic>C</italic><sub>i</sub> dimers to pack. If other motifs competed, the hydrogen bonds and conformations should be lower in energy than the <italic>C</italic><sub>i</sub> dimer. The results detail the effect of packing on the conformation in these molecules. The results are interpreted as a rough measure of the energetic compromise between packing and the energies related to the coordinates involving one dihedral angle and hydrogen bonding. The results establish a connection between solution and solid phase conformation.</p>
<title>Supporting Information</title>
<p>We thank the National Science Foundation, DGE9870691 for a predoctoral fellowship to PGW and CHE-9702287 for financial support. The beginning stages of this research were funded by NASA-EPSCoR of Kentucky. We thank the University of Kentucky Information Technology Computing Center for CPU time.</p>
<fig id="F1" position="float"><label>Figure 1</label><caption><p><bold>1:</bold> An intuitive prediction regarding the relationship between crude hydrogen bond donor/acceptor directionality and hydrogen bond molecularity. <bold>2:</bold> The amidine and imidazole moieties are not coplanar due to steric interactions; the CN–CN dihedral angle (50–90°) directs hydrogen bond vectors (dashed arrows). <bold>3:</bold> Imidazole prefers the tape. <bold>4:</bold> This tautomer/rotamer is not observed.</p></caption></fig> <fig id="F2" position="float"><label>Figure 2</label><caption><p>Facile syntheses of imidazole carboxamidines from commercial imidazoles and carbodiimides furnished a series of crystalline phases with related hydrogen bonding.</p></caption></fig> <fig id="F3" position="float"><label>Figure 3</label><caption><p>The NCNC dihedral angle, θ, between the hydrogen bond donors and acceptors, was assigned values between +180° and −180°. Structures with opposite signs are conformational enantiomers.</p></caption></fig> <fig id="F4" position="float"><label>Figure 4</label><caption><p>Stereoview of Dimer <bold>5c</bold>. This dimer stacked imidazole rings with R<sup>1</sup> pointing in opposite directions.</p></caption></fig> <fig id="F5" position="float"><label>Figure 5</label><caption><p>Stereoview of trimer <bold>6b</bold>.</p></caption></fig> <fig id="F6" position="float"><label>Figure 6</label><caption><p>Stereoview of tetramer <bold>9a</bold>.</p></caption></fig> <fig id="F7" position="float"><label>Figure 7</label><caption><p>Stereoview of linear hydrogen bond tape <bold>11a</bold>.</p></caption></fig> <fig id="F8" position="float"><label>Figure 8</label><caption><p>The calculated (rhf/6-311+g(d,p)) potential energy (kcal/mol) of <italic>N,N</italic>'-dimethyl-1H-imidazole-1-carboxamidine, <bold>2</bold>, R<sup>1</sup>,R<sup>2</sup>=H, R<sup>3</sup>=Me versus dihedral angle θ (degrees). Icons represent crystal structure θ values: (∘) = <italic>C</italic><sub>i</sub> dimer, (◊) = <italic>C</italic><sub>1</sub> dimer, (□) = trimer, (▵) = tetramer, (■) = tape. The red icons are in apparent contradiction of the trend: the <italic>C</italic><sub>2</sub> dimer has high-energy θ; these are discussed further.</p></caption></fig> <fig id="F9" position="float"><label>Figure 9</label><caption><p>A Flow chart for the calculation of the energies of the n-mers minus the effects of packing and substituent interactions. See text.</p></caption></fig> <fig id="F10" position="float"><label>Figure 10</label><caption><p>Icons correspond to those in ##FIG##7##Figure 8##. Crosses indicate structures with aromatic groups. The calculated (rhf/6-311+g(d,p)) energy (kcal/mol) with counterpoise BSSE correction of <bold>2b</bold> (##FIG##0##Figure 1##: R<sup>1</sup>−R<sup>3</sup> = H) in the particular crystal structure atomic coordinates versus empirical packing density. Broad conclusions: The crystalline phases with aromatic groups were denser. Gas phase hydrogen bonds in the <italic>C</italic><sub>i</sub> dimers were less stable but <italic>C</italic><sub>i</sub> solid states were denser; the red icons are obvious exceptions. See text.</p></caption></fig> <fig id="F11" position="float"><label>Figure 11</label><caption><p>The icon legend is identical to ##FIG##7##Figure 8## and ##FIG##9##Figure 10##. The Y-axis from ##FIG##7##Figure 8## energies (θ only) and the X-axis from ##FIG##9##Figure 10## energies are compared. The high and low energies of θ are mostly responsible for the two crystalline phases that lie outside the prediction that the gas-phase <italic>C</italic><sub>i</sub> dimer should be high-energy. Dashed line has slope = 1.</p></caption></fig>
<table-wrap id="T1" position="float"><label>Table 1</label><caption><p>The imidazolecarboxamidines synthesized and crystallized for this study. Given are the hydrogen bond molecularity, the imidazole/amidine dihedral angle (θ), and the space group. Further information regarding data collection can be found in the ##SUPPL##1##Supporting Information File 2##.</p></caption><table frame="hsides" rules="groups"><tr><td align="left" colspan="2" valign="middle" rowspan="1">                 Amidine<break/>────────────<break/>Imidazole</td><td align="left" valign="middle" rowspan="1" colspan="1"><break/><bold>a</bold></td><td align="left" valign="middle" rowspan="1" colspan="1"><break/><bold>b</bold></td><td align="left" valign="middle" rowspan="1" colspan="1"><break/><bold>c</bold></td></tr><tr><td align="left" colspan="5" valign="middle" rowspan="1"><hr/></td></tr><tr><td align="left" valign="middle" rowspan="1" colspan="1"><bold>5</bold></td><td align="center" valign="middle" rowspan="1" colspan="1"></td><td align="left" valign="middle" rowspan="1" colspan="1">no crystal</td><td align="left" valign="middle" rowspan="1" colspan="1">dimer<break/>±81.5°, <italic>P</italic>-1</td><td align="left" valign="middle" rowspan="1" colspan="1">dimer<break/>±87.0°, <italic>P</italic>2<sub>1</sub>/c<break/>hexane or EtOAc</td></tr><tr><td align="left" valign="middle" rowspan="1" colspan="1"><bold>6</bold></td><td align="center" valign="middle" rowspan="1" colspan="1"></td><td align="left" valign="middle" rowspan="1" colspan="1">no crystal</td><td align="left" valign="middle" rowspan="1" colspan="1">trimer<break/>−102.6°, 108.3°, (104.3°, −98.6°), <italic>P</italic>2<sub>1</sub>/c</td><td align="left" valign="middle" rowspan="1" colspan="1">dimer<break/>±66.7°, <italic>C</italic>2/c<break/>hexane or EtOAc</td></tr><tr><td align="left" valign="middle" rowspan="1" colspan="1"><bold>7</bold></td><td align="center" valign="middle" rowspan="1" colspan="1"></td><td align="left" valign="middle" rowspan="1" colspan="1">dimer<break/>±91.8°, <italic>P</italic>-1<break/>hexane</td><td align="left" valign="middle" rowspan="1" colspan="1">dimer<break/>±100.7°, <italic>P</italic>2<sub>1</sub>/n<break/>MeOH/EtOAc</td><td align="left" valign="middle" rowspan="1" colspan="1">(1) dimer<break/>±108.5°, <italic>P</italic>2<sub>1</sub>/n, hexane<break/>------------<break/>(2) dimer polymorph<break/>−91.8°, 91.0°, <italic>P</italic>2<sub>1</sub>/c</td></tr><tr><td align="left" valign="middle" rowspan="1" colspan="1"><bold>8</bold></td><td align="center" valign="middle" rowspan="1" colspan="1"></td><td align="left" valign="middle" rowspan="1" colspan="1"><italic>C</italic><sub>1</sub> dimer<break/>56.8°, 67.8°, <italic>P</italic>-1<break/>hexane or EtOAc</td><td align="left" valign="middle" rowspan="1" colspan="1"><italic>C</italic><sub>1</sub> dimer<break/>95.5°, 60.5°, <italic>P</italic>2<sub>1</sub>/n</td><td align="left" valign="middle" rowspan="1" colspan="1">no crystal</td></tr><tr><td align="left" valign="middle" rowspan="1" colspan="1"><bold>9</bold></td><td align="center" valign="middle" rowspan="1" colspan="1"></td><td align="left" valign="middle" rowspan="1" colspan="1">tetramer<break/>101.2°, −100.1°, 103.0°, −106.0°, <italic>P</italic>2<sub>1</sub>/n<break/>hexane</td><td align="left" valign="middle" rowspan="1" colspan="1">(1)−(3) dimer<break/>±54.3°, <italic>P</italic>2<sub>1</sub>/n<break/>±53.0°, <italic>P</italic>2<sub>1</sub>/n<break/>±54.1°, <italic>P</italic>2<sub>1</sub>/n<break/>(1) EtOAc, (2) C<sub>6</sub>H<sub>5</sub>CH<sub>3</sub>, (3) Et<sub>2</sub>O all 1:1</td><td align="left" valign="middle" rowspan="1" colspan="1">dimer<break/>±88.0°, <italic>P</italic>2<sub>1</sub>/n</td></tr><tr><td align="left" valign="middle" rowspan="1" colspan="1"><bold>10</bold></td><td align="center" valign="middle" rowspan="1" colspan="1"></td><td align="left" valign="middle" rowspan="1" colspan="1">no crystal</td><td align="left" valign="middle" rowspan="1" colspan="1">dimer<break/>±91.1°, <italic>P</italic>-1</td><td align="left" valign="middle" rowspan="1" colspan="1">dimer<break/>±73.3°, <italic>P</italic>-1<break/>hexane or EtOAc</td></tr><tr><td align="left" valign="middle" rowspan="1" colspan="1"><bold>11</bold></td><td align="center" valign="middle" rowspan="1" colspan="1"></td><td align="left" valign="middle" rowspan="1" colspan="1">tape<break/>63.2°, <italic>P</italic>2<sub>1</sub>/c</td><td align="left" valign="middle" rowspan="1" colspan="1">dimer<break/>±99.7°, <italic>P</italic>2<sub>1</sub>/n</td><td align="left" valign="middle" rowspan="1" colspan="1">dimer<break/>±98.7°, <italic>P</italic>-1</td></tr></table></table-wrap>
<supplementary-material content-type="local-data" id="SD1"><label>File 1</label><caption><p>Experimental Section. The Experimental Section describes the synthesis and purification of all substances described, and general experimental procedures.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="SD2"><label>File 2</label><caption><p>Crystal Data and Structure Refinement Information</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="SD3"><label>File 3</label><caption><p>Sample <sup>1</sup>H NMR Spectra for <bold>10b</bold>, and <bold>14b</bold></p></caption></supplementary-material> <supplementary-material content-type="local-data" id="SD4"><label>File 4</label><caption><p>Crystallographic Information Files (CIF) for all the compounds reported in this work. File Names: Xy_descriptor.cif. Xy is the compound number and letter designation found in Table I and/or described in the Experimental Section. The descriptor is a short tag indicating hydrogen bond motif or symmetry.</p></caption></supplementary-material>
<table-wrap-foot><fn id="TFN1"><p>If not noted, crystals were from EtOAc.</p></fn></table-wrap-foot>
<graphic xlink:href="Beilstein_J_Org_Chem-04-23-g002"/> <graphic xlink:href="Beilstein_J_Org_Chem-04-23-g003"/> <inline-graphic xlink:href="Beilstein_J_Org_Chem-04-23-i001.jpg"/> <inline-graphic xlink:href="Beilstein_J_Org_Chem-04-23-i002.jpg"/> <inline-graphic xlink:href="Beilstein_J_Org_Chem-04-23-i003.jpg"/> <inline-graphic xlink:href="Beilstein_J_Org_Chem-04-23-i004.jpg"/> <inline-graphic xlink:href="Beilstein_J_Org_Chem-04-23-i005.jpg"/> <inline-graphic xlink:href="Beilstein_J_Org_Chem-04-23-i006.jpg"/> <inline-graphic xlink:href="Beilstein_J_Org_Chem-04-23-i007.jpg"/> <inline-graphic xlink:href="Beilstein_J_Org_Chem-04-23-i008.jpg"/> <inline-graphic xlink:href="Beilstein_J_Org_Chem-04-23-i009.jpg"/> <inline-graphic xlink:href="Beilstein_J_Org_Chem-04-23-i010.jpg"/> <graphic xlink:href="Beilstein_J_Org_Chem-04-23-g004"/> <graphic xlink:href="Beilstein_J_Org_Chem-04-23-g005"/> <graphic xlink:href="Beilstein_J_Org_Chem-04-23-g006"/> <graphic xlink:href="Beilstein_J_Org_Chem-04-23-g007"/> <graphic xlink:href="Beilstein_J_Org_Chem-04-23-g008"/> <graphic xlink:href="Beilstein_J_Org_Chem-04-23-g009"/> <graphic xlink:href="Beilstein_J_Org_Chem-04-23-g010"/> <graphic xlink:href="Beilstein_J_Org_Chem-04-23-g011"/> <graphic xlink:href="Beilstein_J_Org_Chem-04-23-g012"/>
<media mime-subtype="pdf" mimetype="application" xlink:href="Beilstein_J_Org_Chem-04-23-s001.pdf" xlink:type="simple" id="d39e1771" position="anchor"/> <media mime-subtype="pdf" mimetype="application" xlink:href="Beilstein_J_Org_Chem-04-23-s002.pdf" xlink:type="simple" id="d39e1778" position="anchor"/> <media mime-subtype="pdf" mimetype="application" xlink:href="Beilstein_J_Org_Chem-04-23-s003.pdf" xlink:type="simple" id="d39e1793" position="anchor"/> <media mime-subtype="zip" mimetype="application" xlink:href="Beilstein_J_Org_Chem-04-23-s004.zip" xlink:type="simple" id="d39e1800" position="anchor"/>
[{"label": ["1"], "surname": ["Steiner"], "given-names": ["T"], "source": ["Angew Chem, Int Ed"], "year": ["2002"], "volume": ["41"], "fpage": ["48"], "lpage": ["76"], "pub-id": ["10.1002/1521-3773(20020104)41:1<48::AID-ANIE48>3.0.CO;2-U"]}, {"label": ["2"], "surname": ["Baures", "Rush", "Wiznycia", "Desper", "Helfrich", "Beatty"], "given-names": ["P W", "J R", "A V", "J", "B A", "A M"], "source": ["Cryst Growth Des"], "year": ["2002"], "volume": ["2"], "fpage": ["653"], "lpage": ["664"], "pub-id": ["10.1021/cg025549j"]}, {"label": ["3"], "surname": ["Etter"], "given-names": ["M C"], "source": ["Acc Chem Res"], "year": ["1990"], "volume": ["23"], "fpage": ["120"], "lpage": ["126"], "pub-id": ["10.1021/ar00172a005"]}, {"label": ["5"], "surname": ["Krische", "Lehn"], "given-names": ["M J", "J M"], "article-title": ["The utilization of persistent H-bonding motifs in the self-assembly of supramolecular architectures"], "source": ["Structure and Bonding"], "year": ["2000"], "volume": ["96"], "publisher-loc": ["Berlin/Heidelberg, Germany"], "publisher-name": ["Springer"], "fpage": ["3"], "lpage": ["29"], "pub-id": ["10.1007/3-540-46591-X_1"]}, {"label": ["6"], "surname": ["Zerkowski", "MacDonald", "Seto", "Wierda", "Whitesides"], "given-names": ["J A", "J C", "C T", "D A", "G M"], "source": ["J Am Chem Soc"], "year": ["1994"], "volume": ["116"], "fpage": ["2382"], "lpage": ["2391"], "pub-id": ["10.1021/ja00085a018"]}, {"label": ["7"], "surname": ["Zerkowski", "Whitesides"], "given-names": ["J A", "G A"], "source": ["J Am Chem Soc"], "year": ["1994"], "volume": ["116"], "fpage": ["4298"], "lpage": ["4304"], "pub-id": ["10.1021/ja00089a019"]}, {"label": ["8"], "surname": ["Gallant", "Phan Viet", "Wuest"], "given-names": ["M", "M T", "J D"], "source": ["J Org Chem"], "year": ["1991"], "volume": ["56"], "fpage": ["2284"], "lpage": ["2286"], "pub-id": ["10.1021/jo00007a007"]}, {"label": ["9"], "surname": ["Mel\u00e9ndez", "Hamilton", "Murto", "Heino"], "given-names": ["R E", "A D", "J", "E L"], "source": ["Top Curr Chem"], "year": ["1998"], "volume": ["198"], "fpage": ["97"], "lpage": ["129"]}, {"label": ["10"], "surname": ["Schwiebert", "Chin", "MacDonald", "Whitesides"], "given-names": ["K E", "D N", "J C", "G M"], "source": ["J Am Chem Soc"], "year": ["1996"], "volume": ["118"], "fpage": ["4018"], "lpage": ["4029"], "pub-id": ["10.1021/ja952836l"]}, {"label": ["13"], "surname": ["Cammers", "Parkin"], "given-names": ["A", "S"], "source": ["CrystEngComm"], "year": ["2004"], "volume": ["6"], "fpage": ["168"], "lpage": ["172"], "pub-id": ["10.1039/b405096g"]}, {"label": ["14"], "surname": ["Taylor", "Kennard"], "given-names": ["R", "O"], "source": ["Acc Chem Res"], "year": ["1984"], "volume": ["17"], "fpage": ["320"], "lpage": ["326"], "pub-id": ["10.1021/ar00105a004"]}, {"label": ["16"], "surname": ["Steiner", "Desiraju"], "given-names": ["T", "G R"], "source": ["Chem Commun"], "year": ["1998"], "fpage": ["891"], "lpage": ["892"], "pub-id": ["10.1039/a708099i"]}, {"label": ["17"], "surname": ["Leiserowitz", "Schmidt"], "given-names": ["L", "G M J"], "source": ["J Chem Soc A"], "year": ["1969"], "fpage": ["2372"], "lpage": ["2382"], "pub-id": ["10.1039/j19690002372"]}, {"label": ["18"], "surname": ["Gladwell"], "given-names": ["M"], "source": ["The Tipping Point (How Little Things Can Make a Big Difference)"], "year": ["2000"], "publisher-loc": ["New York, NY"], "publisher-name": ["Little, Brown and Company"]}, {"label": ["19"], "surname": ["Bertolasi", "Gilli", "Ferretti", "Gilli"], "given-names": ["V", "P", "V", "G"], "source": ["Acta Crystallogr, Sect B: Struct Sci"], "year": ["2001"], "volume": ["57"], "fpage": ["591"], "lpage": ["598"], "pub-id": ["10.1107/S010876810100814X"]}, {"label": ["21"], "surname": ["Das", "Banerjee", "Mondal", "Howard", "Boese", "Desiraju"], "given-names": ["D", "R", "R", "J A K", "R", "G R"], "source": ["Chem Commun"], "year": ["2006"], "fpage": ["555"], "lpage": ["557"], "pub-id": ["10.1039/b514076e"]}, {"label": ["22"], "surname": ["Dunitz", "Filippini", "Gavezzotti"], "given-names": ["J D", "G", "A"], "source": ["Helv Chim Acta"], "year": ["2000"], "volume": ["83"], "fpage": ["2317"], "lpage": ["2335"], "pub-id": ["10.1002/1522-2675(20000906)83:9<2317::AID-HLCA2317>3.0.CO;2-R"]}, {"label": ["23"], "surname": ["Wallach"], "given-names": ["O"], "source": ["Liebigs Ann Chem"], "year": ["1895"], "volume": ["286"], "fpage": ["90"], "lpage": ["143"], "pub-id": ["10.1002/jlac.18952860105"]}, {"label": ["24"], "surname": ["Reutzel-Edens", "Russell", "Yu"], "given-names": ["S M", "V A", "L"], "source": ["J Chem Soc, Perkin Trans 2"], "year": ["2000"], "volume": ["5"], "fpage": ["913"], "lpage": ["924"], "pub-id": ["10.1039/a909259e"]}, {"label": ["25"], "surname": ["Brock", "Dunitz"], "given-names": ["C P", "J D"], "source": ["Chem Mater"], "year": ["1994"], "volume": ["6"], "fpage": ["1118"], "lpage": ["1127"], "pub-id": ["10.1021/cm00044a010"]}, {"label": ["26"], "surname": ["Gavezzotti"], "given-names": ["A"], "source": ["CrystEngComm"], "year": ["2008"], "volume": ["10"], "fpage": ["389"], "lpage": ["398"], "pub-id": ["10.1039/b714349d"]}, {"label": ["27"], "surname": ["Brock", "Schweizer", "Dunitz"], "given-names": ["C P", "W B", "J D"], "source": ["J Am Chem Soc"], "year": ["1991"], "volume": ["113"], "fpage": ["9811"], "lpage": ["9820"], "pub-id": ["10.1021/ja00026a015"]}, {"label": ["29"], "surname": ["Dunitz", "Gavezzotti"], "given-names": ["J D", "A"], "source": ["Cryst Growth Des"], "year": ["2005"], "volume": ["5"], "fpage": ["2180"], "lpage": ["2189"], "pub-id": ["10.1021/cg050098z"]}, {"label": ["30"], "surname": ["Foresman", "Frisch"], "given-names": ["J B", "A"], "source": ["Exploring Chemistry with Electronic Structure Methods"], "year": ["1996"], "publisher-loc": ["Pittsburgh, PA"], "publisher-name": ["Gaussian, Inc."], "fpage": ["7"]}, {"label": ["31"], "collab": ["Gaussian 03"], "source": ["Revision C.02"], "year": ["2004"], "publisher-loc": ["Wallingford, CT"], "publisher-name": ["Gaussian, Inc."]}, {"label": ["32"], "surname": ["Philp", "Robinson"], "given-names": ["D", "J M A"], "source": ["J Chem Soc, Perkin Trans 2"], "year": ["1998"], "fpage": ["1643"], "lpage": ["1650"], "pub-id": ["10.1039/a800931g"]}, {"label": ["33"], "surname": ["Muguet", "Robinson"], "given-names": ["F F", "G W"], "source": ["J Chem Phys"], "year": ["1995"], "volume": ["102"], "fpage": ["3648"], "lpage": ["3654"], "pub-id": ["10.1063/1.468593"]}, {"label": ["34"], "surname": ["Simon", "Duran", "Dannenberg"], "given-names": ["S", "M", "J J"], "source": ["J Chem Phys"], "year": ["1996"], "volume": ["105"], "fpage": ["11024"], "lpage": ["11031"], "pub-id": ["10.1063/1.472902"]}, {"label": ["35"], "surname": ["Kitaigorodskii"], "given-names": ["A I"], "source": ["Molecular Crystals and Molecules"], "year": ["1973"], "publisher-loc": ["New York"], "publisher-name": ["Academic Press"]}, {"label": ["36"], "surname": ["Dunitz", "Gavezzotti"], "given-names": ["J D", "A"], "source": ["Angew Chem, Int Ed"], "year": ["2005"], "volume": ["44"], "fpage": ["1766"], "lpage": ["1787"], "pub-id": ["10.1002/anie.200460157"]}, {"label": ["37"], "surname": ["Immirzi", "Perini"], "given-names": ["A", "B"], "source": ["Acta Crystallogr, Sect A"], "year": ["1977"], "volume": ["A33"], "fpage": ["216"], "lpage": ["218"], "pub-id": ["10.1107/S0567739477000448"]}, {"label": ["38"], "surname": ["Mondal", "Howard"], "given-names": ["R", "J A K"], "source": ["CrystEngComm"], "year": ["2005"], "volume": ["7"], "fpage": ["462"], "lpage": ["464"], "pub-id": ["10.1039/b509170p"]}, {"label": ["39"], "surname": ["Adalsteinsson", "Maulitz", "Bruice"], "given-names": ["H", "A H", "T C"], "source": ["J Am Chem Soc"], "year": ["1996"], "volume": ["118"], "fpage": ["7689"], "lpage": ["7693"], "pub-id": ["10.1021/ja954267n"]}, {"label": ["40"], "surname": ["Brehmer", "Weber", "Cano"], "given-names": ["T H", "E", "F H"], "source": ["J Phys Org Chem"], "year": ["2000"], "volume": ["13"], "fpage": ["63"], "lpage": ["74"], "pub-id": ["10.1002/(SICI)1099-1395(200001)13:1<63::AID-POC194>3.0.CO;2-E"]}, {"label": ["41"], "surname": ["Anderson", "Afarinkia", "Yu", "Goeta", "Steed"], "given-names": ["K M", "K", "H-W", "A E", "J W"], "source": ["Cryst Growth Des"], "year": ["2006"], "volume": ["6"], "fpage": ["2109"], "lpage": ["2113"], "pub-id": ["10.1021/cg0603265"]}, {"label": ["42"], "surname": ["Wheeler", "Grove", "Davis", "Kassel"], "given-names": ["K A", "R C", "R E", "W S"], "source": ["Angew Chem, Int Ed"], "year": ["2008"], "volume": ["47"], "fpage": ["78"], "lpage": ["81"], "pub-id": ["10.1002/anie.200704007"]}, {"label": ["43"], "surname": ["Lineberry", "Benjamin", "Davis", "Kassel", "Wheeler"], "given-names": ["A M", "E T", "R E", "W S", "K A"], "source": ["Cryst Growth Des"], "year": ["2008"], "volume": ["8"], "fpage": ["612"], "lpage": ["619"], "pub-id": ["10.1021/cg700789y"]}, {"label": ["44"], "surname": ["Adler", "Borho", "Reiher", "Suhm"], "given-names": ["T B", "N", "M", "M A"], "source": ["Angew Chem, Int Ed"], "year": ["2006"], "volume": ["45"], "fpage": ["3440"], "lpage": ["3445"], "pub-id": ["10.1002/anie.200600380"]}, {"label": ["45"], "surname": ["Bowers", "Hopkins", "Yap", "Wheeler"], "given-names": ["J R", "G W", "G P A", "K A"], "source": ["Cryst Growth Des"], "year": ["2005"], "volume": ["5"], "fpage": ["727"], "lpage": ["736"], "pub-id": ["10.1021/cg0497391"]}, {"label": ["46"], "surname": ["Karamertzanis", "Price"], "given-names": ["P G", "S L"], "source": ["J Chem Theory Comput"], "year": ["2006"], "volume": ["2"], "fpage": ["1184"], "lpage": ["1199"], "pub-id": ["10.1021/ct600111s"]}, {"label": ["47"], "surname": ["Chaka", "Zaniewski", "Youngs", "Tessier", "Klopman"], "given-names": ["A M", "R", "W", "C", "G"], "source": ["Acta Crystallogr, Sect B: Struct Sci"], "year": ["1996"], "volume": ["52"], "fpage": ["165"], "lpage": ["183"], "pub-id": ["10.1107/S0108768195006987"]}]
Beilstein J Org Chem. 2008 Jul 7; 4:23