Related Experiment Video
Updated: Mar 25, 2026

Computation of Atmospheric Concentrations of Molecular Clusters from ab initio Thermochemistry
Published on: April 8, 2020
Specific energy contributions from competing hydrogen-bonded structures in six polymorphs of phenobarbital
Thomas Gelbrich1, Doris E Braun1, Ulrich J Griesser1
1Institute of Pharmacy, University of Innsbruck, Innrain 52c, 6020 Innsbruck, Austria.
Background:
In solid state structures of organic molecules, identical sets of H-bond donor and acceptor functions can result in a range of distinct H-bond connectivity modes. Specifically, competing H-bond structures (HBSs) may differ in the quantitative proportion between one-point and multiple-point H-bond connections. For an assessment of such HBSs, the effects of their internal as well as external (packing) interactions need to be taken into consideration. The semi-classical density sums (SCDS-PIXEL) method, which enables the calculation of interaction energies for molecule-molecule pairs, was used to investigate six polymorphs of phenobarbital (Pbtl) with different quantitative proportions of one-point and two-point H-bond connections.
Results:
The structures of polymorphs V and VI of Pbtl were determined from single crystal data. Two-point H-bond connections are inherently inflexible in their geometry and lie within a small PIXEL energy range (-45.7 to -49.7 kJ mol(-1)). One-point H-bond connections are geometrically less restricted and subsequently show large variations in their dispersion terms and total energies (-23.1 to -40.5 kJ mol(-1)). The comparison of sums of interaction energies in small clusters containing only the strongest intermolecular interactions showed an advantage for compact HBSs with multiple-point connections, whereas alternative HBSs based on one-point connections may enable more favourable overall packing interactions (i.e. V vs. III). Energy penalties associated with experimental intramolecular geometries relative to the global conformational energy minimum were calculated and used to correct total PIXEL energies. The estimated order of stabilities (based on PIXEL energies) is III > I > II > VI > X > V, with a difference of just 1.7 kJ mol(-1) between the three most stable forms.
Conclusions:
For an analysis of competing HBSs, one has to consider the contributions from internal H-bond and non-H-bond interactions, from the packing of multiple HBS instances and intramolecular energy penalties. A compact HBS based on multiple-point H-bond connections should typically lead to more packing alternatives and ultimately to a larger number of viable low-energy structures than a competing one-point HBS (i.e. dimer vs. catemer). Coulombic interaction energies associated with typical short intermolecular C-H···O contact geometries are small in comparison with dispersion effects associated with the packing complementary molecular shapes.Graphical abstractCompeting H-bond motifs can differ markedly in their energy contributions.
More Related Videos
05:51Isotopic Effect in Double Proton Transfer Process of Porphycene Investigated by Enhanced QM/MM Method
Published on: July 19, 2019
08:49Incorporating Target Protein Structure Flexibility and Dynamics in Computational Drug Discovery Using Ensemble-Based Docking Analysis
Published on: June 20, 2025
Related Concept Videos
Acidity and Basicity of Alcohols and Phenols
Physical Properties of Alcohols and Phenols
Alcohols possess a higher boiling point than aliphatic hydrocarbons of similar...
Resonance and Hybrid Structures
Resonance Structures and Resonance Hybrids
The Lewis structure of a nitrite anion (NO2−) may actually be drawn in two different ways, distinguished by the locations of the N–O and N=O bonds.
Conformations of Butane
Crystal Field Theory - Tetrahedral and Square Planar Complexes
Crystal field theory (CFT) is applicable to molecules in geometries other than octahedral. In octahedral complexes, the lobes of the dx2−y2 and dz2 orbitals point directly at the ligands. For tetrahedral complexes, the d orbitals remain in place, but with only four ligands located between the axes. None of the orbitals points directly at the tetrahedral ligands. However, the dx2−y2 and dz2 orbitals (along the Cartesian axes) overlap with the ligands less than the dxy,...
Structure of Benzene: Molecular Orbital Model