Is it possible to derive quantitative information on polarization of electron density from the multipolar model?

Protonated nucleobases are not fully ionized in their chloride salt crystals and form metastable base pairs further stabilized by the surrounding anions

This paper presents experimental charge-density studies of cytosinium chloride, adeninium chloride hemihydrate and guaninium dichloride crystals based on ultra-high-resolution X-ray diffraction data and extensive theoretical calculations. The results confirm that the cohesive energies of the studied systems are dominated by contributions from intermolecular electrostatic interactions, as expected for ionic crystals. Electrostatic interaction energies (Ees) usually constitute 95% of the total interaction energy. The Ees energies in this study were several times larger in absolute value when compared, for example, with dimers of neutral nucleobases. However, they were not as large as some theoretical calculations have predicted. This was because the molecules appeared not to be fully ionized in the studied crystals. Apart from charge transfer from chlorine to the protonated nucleobases, small but visible charge redistribution within the nucleobase cations was observed. Some dimers of singly protonated bases in the studied crystals, namely a cytosinium-cytosinium trans sugar/sugar edge pair and an adeninium-adeninium trans Hoogsteen/Hoogsteen edge pair, exhibited attractive interactions (negative values of Ees) or unusually low repulsion despite identical molecular charges. The pairs are metastable as a result of strong hydrogen bonding between bases which overcompensates the overall cation-cation repulsion, the latter being weakened due to charge transfer and molecular charge-density polarization.

Charge density of the biologically active molecule (2-oxo-1,3-benzoxazol-3(2H)-yl)acetic acid

(2-Oxo-1,3-benzoxazol-3(2H)-yl)acetic acid is a member of a biologically active class of compounds. Its molecular structure in the crystal has been determined by X-ray diffraction, and its gas phase structure was obtained by quantum chemical calculations at the B3LYP/6-311++G(d,p) level of theory. In order to understand the dynamics of the molecule, two presumably soft degrees of freedom associated with the relative orientation of the planar benzoxazolone system and its substituent at the N atom were varied systematically. Five conformers have been identified as local minima on the resulting two-dimensional potential energy surface within an energy window of 27 kJ mol(-1). The energetically most favourable minimum closely matches the conformation observed in the crystal. Based on high-resolution diffraction data collected at low temperature, the experimental electron density of the compound was determined. Comparison with the electron density established by theory for the isolated molecule allowed the effect of intermolecular interactions to be addressed, in particular a moderately strong O-H...O hydrogen bond with a donor...acceptor distance of 2.6177 (9) Å: the oxygen acceptor is clearly polarized in the extended solid. The hydrogen bond connects consecutive molecules to chains, and the pronounced charge separation leads to stacking between neighburs with antiparallel dipole moments perpendicular to the chain direction.

Host perturbation in a β-hydroquinone clathrate studied by combined X-ray/neutron charge-density analysis: implications for molecular inclusion in supramolecular entities

X-ray/neutron (X/N) diffraction data measured at very low temperature (15 K) in conjunction with ab initio theoretical calculations were used to model the crystal charge density (CD) of the host-guest complex of hydroquinone (HQ) and acetonitrile. Due to pseudosymmetry, information about the ordering of the acetonitrile molecules within the HQ cavities is present only in almost extinct, very weak diffraction data, which cannot be measured with sufficient accuracy even by using the brightest X-ray and neutron sources available, and the CD model of the guest molecule was ultimately based on theoretical calculations. On the other hand, the CD of the HQ host structure is well determined by the experimental data. The neutron diffraction data provide hydrogen anisotropic thermal parameters and positions, which are important to obtain a reliable CD for this light-atom-only crystal. Atomic displacement parameters obtained independently from the X-ray and neutron diffraction data show excellent agreement with a |ΔU| value of 0.00058 Å(2) indicating outstanding data quality. The CD and especially the derived electrostatic properties clearly reveal increased polarization of the HQ molecules in the host-guest complex compared with the HQ molecules in the empty HQ apohost crystal structure. It was found that the origin of the increased polarization is inclusion of the acetonitrile molecule, whereas the change in geometry of the HQ host structure following inclusion of the guest has very little effect on the electrostatic potential. The fact that guest inclusion has a profound effect on the electrostatic potential suggests that nonpolarizable force fields may be unsuitable for molecular dynamics simulations of host-guest interaction (e.g., in protein-drug complexes), at least for polar molecules.