Enhancing hydrogen positions in X-ray structures of transition metal hydride complexes with dynamic quantum crystallography.

Hirshfeld atom refinement (HAR) is a method which enables the user to obtain more accurate positions of hydrogen atoms bonded to light chemical elements using X-ray data. When data quality permits, this method can be extended to hydrogen-bonded transition metals (TMs), as in hydride complexes. However, addressing hydrogen thermal motions with HAR, particularly in TM hydrides, presents a challenge. At the same time, proper description of thermal vibrations can be vital for determining hydrogen positions correctly. In this study, we employ tools such as SHADE3 and Normal Mode Refinement (NoMoRe) to estimate anisotropic displacement parameters (ADPs) for hydrogen atoms during HAR and IAM refinements performed for seven structures of TM (Fe, Ni, Cr, Nb, Rh and Os) and metalloid (Sb) hydride complexes for which both the neutron and the X-ray structures have been determined. A direct comparison between neutron and HAR/SHADE3/NoMoRe ADPs reveals that the similarity between neutron hydrogen ADPs and those estimated with NoMoRe or SHADE3 is significantly higher than when hydrogen ADPs are refined with HAR. Regarding TM-H bond lengths, traditional HAR exhibits a slight advantage over the other methods. However, combining NoMoRe/SHADE3 with HAR results in a minor decrease in agreement with neutron TM-H bond lengths. For the Cr complex, for which high-resolution X-ray data were collected, an investigation of resolution-related effects was possible.

Theoretically derived thermodynamic properties can be improved by the refinement of low-frequency modes against X-ray diffraction data.

Herein, a framework for the estimation of the thermodynamic properties of molecular crystals via the refinement of frequencies from density functional theory calculations against X-ray diffraction data is presented. The framework provides an efficient approach to including the contribution of acoustic modes in the thermodynamic properties. The obtained heat capacities for urea, the α- and ß-glycine polymorphs, benzoic acid, and 4'-hydroxyacetophenone are in good agreement with those from adiabatic calorimetry.

Further Validation of Quantum Crystallography Approaches.

Quantum crystallography is a fast-developing multidisciplinary area of crystallography. In this work, we analyse the influence of different charge density models (i.e., the multipole model (MM), Hirshfeld atom refinement (HAR), and the transferable aspherical atom model (TAAM)), modelling of the thermal motion of hydrogen atoms (anisotropic, isotropic, and with the aid of SHADE or NoMoRe), and the type of radiation used (Mo Kα and Cu Kα) on the final results. To achieve this aim, we performed a series of refinements against X-ray diffraction data for three model compounds and compared their final structures, geometries, shapes of ADPs, and charge density distributions. Our results were also supported by theoretical calculations that enabled comparisons of the lattice energies of these structures. It appears that geometrical parameters are better described (closer to the neutron values) when HAR is used; however, bonds to H atoms more closely match neutron values after MM or TAAM refinement. Our analysis shows the superiority of the NoMoRe method in the description of H-atom ADPs. Moreover, the shapes of the ADPs of H atoms, as well as their electron density distributions, were better described with low-resolution Cu Kα data in comparison to low-resolution Mo Kα data.

A combined model of electron density and lattice dynamics refined against elastic diffraction data. Thermodynamic properties of crystalline L-alanine.

Thermodynamic stability is an essential property of crystalline materials, and its accurate calculation requires a reliable description of the thermal motion - phonons - in the crystal. Such information can be obtained from periodic density functional theory (DFT) calculations, but these are costly and in some cases insufficiently accurate for molecular crystals. This deficiency is addressed here by refining a lattice-dynamics model, derived from DFT calculations, against accurate high-resolution X-ray diffraction data. For the first time, a normal-mode refinement is combined with the refinement of aspherical atomic form factors, allowing a comprehensive description and physically meaningful deconvolution of thermal motion and static charge density in the crystal. The small and well diffracting L-alanine system was used. Different lattice-dynamics models, with or without phonon dispersion, and derived from different levels of theory, were tested, and models using spherical and aspherical form factors were compared. The refinements indicate that the vibrational information content in the 23 K data is too small to study lattice dynamics, whereas the 123 K data appear to hold information on the acoustic and lowest-frequency optical phonons. These normal-mode models show slightly larger refinement residuals than their counterparts using atomic displacement parameters, and these features are not removed by considering phonon dispersion in the model. The models refined against the 123 K data, regardless of their sophistication, give calculated heat capacities for L-alanine within less than 1 cal mol-1 K-1 of the calorimetric measurements, in the temperature range 10-300 K. The findings show that the normal-mode refinement method can be combined with an elaborate description of the electron density. It appears to be a promising technique for free-energy determination for crystalline materials at the expense of performing a single-crystal elastic X-ray diffraction determination combined with periodic DFT calculations.

X-ray diffraction data as a source of the vibrational free-energy contribution in polymorphic systems.

In this contribution we attempt to answer a general question: can X-ray diffraction data combined with theoretical computations be a source of information about the thermodynamic properties of a given system? Newly collected sets of high-quality multi-temperature single-crystal X-ray diffraction data and complementary periodic DFT calculations of vibrational frequencies and normal mode vectors at the Γ point on the yellow and white polymorphs of di-methyl 3,6-di-chloro-2,5-di-hydroxy-terephthalate are combined using two different approaches, aiming to obtain thermodynamic properties for the two compounds. The first approach uses low-frequency normal modes extracted from multi-temperature X-ray diffraction data (normal coordinate analysis), while the other uses DFT-calculated low-frequency normal mode in the refinement of the same data (normal mode refinement). Thermodynamic data from the literature [Yang et al. (1989), Acta Cryst. B45, 312-323] and new periodic ab initio DFT supercell calculations are used as a reference point. Both approaches tested in this work capture the most essential features of the systems: the polymorphs are enantiotropically related, with the yellow form being the thermodynamically stable system at low temperature, and the white form at higher temperatures. However, the inferred phase transition temperature varies between different approaches. Thanks to the application of unconventional methods of X-ray data refinement and analysis, it was additionally found that, in the case of the yellow polymorph, anharmonicity is an important issue. By discussing contributions from low- and high-frequency modes to the vibrational entropy and enthalpy, the importance of high-frequency modes is highlighted. The analysis shows that larger anisotropic displacement parameters are not always related to the polymorph with the higher vibrational entropy contribution.

Dynamic quantum crystallography: lattice-dynamical models refined against diffraction data. II. Applications to L-alanine, naphthalene and xylitol.

In the first paper of this series [Hoser & Madsen (2016). Acta Cryst. A72, 206-214], a new approach was introduced which enables the refinement of frequencies of normal modes obtained from ab initio periodic computations against single-crystal diffraction data. In this contribution, the performance of this approach is tested by refinement against data in the temperature range from 23 to 205 K on the molecular crystals of L-alanine, naphthalene and xylitol. The models, which are lattice-dynamical models derived at the Γ point of the Brillouin zone, are able to describe the atomic vibrations of L-alanine and naphthalene to a level where the residual densities are similar to those obtained from the independent atom model. For the more flexible molecule xylitol, larger deviations are found. Hydrogen ADPs (anisotropic displacement parameters) derived from the models are in similar or better agreement with neutron diffraction results than ADPs obtained by other procedures. The heat capacity calculated after normal mode refinement for naphthalene is in reasonable agreement with the heat capacity obtained from calorimetric measurements (to less than 1 cal mol-1 K-1 below 300 K), with deviations at higher temperatures indicating anharmonicity. Standard uncertainties and correlation of the refined parameters have been derived based on a Monte Carlo procedure. The uncertainties are quite small and probably underestimated.

Dynamic quantum crystallography: lattice-dynamical models refined against diffraction data. I. Theory.

This study demonstrates and tests the refinement of a lattice-dynamical model derived from periodic ab initio calculations at the Γ point against elastic diffraction data (X-ray or neutron). Refinement of only a handful of parameters is sufficient to obtain a similar agreement with the data as the conventional crystallographic model using anisotropic displacement parameters. By refinement against X-ray data, H displacement parameters are obtained which compare favourably with those from neutron diffraction experiments. The approach opens the door for evaluating thermodynamic properties, and for refinement against multi-temperature data, against inelastic diffraction data, spectroscopic information and thermal diffuse scattering data.

Yes, one can obtain better quality structures from routine X-ray data collection.

Single-crystal X-ray diffraction structural results for benzidine dihydrochloride, hydrated and protonated N,N,N,N-peri(dimethylamino)naphthalene chloride, triptycene, dichlorodimethyltriptycene and decamethylferrocene have been analysed. A critical discussion of the dependence of structural and thermal parameters on resolution for these compounds is presented. Results of refinements against X-ray data, cut off to different resolutions from the high-resolution data files, are compared to structural models derived from neutron diffraction experiments. The Independent Atom Model (IAM) and the Transferable Aspherical Atom Model (TAAM) are tested. The average differences between the X-ray and neutron structural parameters (with the exception of valence angles defined by H atoms) decrease with the increasing 2θmax angle. The scale of differences between X-ray and neutron geometrical parameters can be significantly reduced when data are collected to the higher, than commonly used, 2θmax diffraction angles (for Mo Kα 2θmax > 65°). The final structural and thermal parameters obtained for the studied compounds using TAAM refinement are in better agreement with the neutron values than the IAM results for all resolutions and all compounds. By using TAAM, it is still possible to obtain accurate results even from low-resolution X-ray data. This is particularly important as TAAM is easy to apply and can routinely be used to improve the quality of structural investigations [Dominiak (2015 ▸). LSDB from UBDB. University of Buffalo, USA]. We can recommend that, in order to obtain more adequate (more accurate and precise) structural and displacement parameters during the IAM model refinement, data should be collected up to the larger diffraction angles, at least, for Mo Kα radiation to 2θmax = 65° (sin θmax/λ < 0.75 Å(-1)). The TAAM approach is a very good option to obtain more adequate results even using data collected to the lower 2θmax angles. Also the results of translation-libration-screw (TLS) analysis and vibrational entropy values are more reliable for 2θmax > 65°.

Single crystal X-ray diffraction, spectroscopic and mass spectrometric studies of furanocoumarin peucedanin.

The structure of peucedanin, isolated from Peucedanum tauricum Bieb. (Apiaceae), has been established using single crystal X-ray diffraction. This furanocoumarin isolated from the light petroleum extract of P. tauricum fruits was characterized by high resolution EI-MS, sATR-FTIR and 2D NMR spectroscopic techniques. The EI-MS showed the typical fragmentation pattern of methoxyfuranocoumarins. Extensive 1D (1H and 13C) as well as 2D NMR data enabled complete assignment of the carbon atoms in the peucedanin molecule. The FTIR data confirms intermolecular hydrogen bonding between peucedanin molecules in polar solvents. Peucedanin crystallises in the R-3 space group from the trigonal system with one molecule in the asymmetric part of the unit cell. The crystal lattice of peucedanin consists of the molecules arranged in separate columns. They are related by two fold screw axes and centres of symmetry. Interestingly, peucedanin columns form two channels per unit cell with a diameter of 7.5angstrom going through the crystal lattice in the Z-direction. These channels are filled with disordered water molecules, which are surrounded by hydrophobic methyl groups and are located exactly at the centres of the channels. The peucedanin molecules are stacked in a single column with the opposite orientation of the neighbouring molecules. These results could be interesting in further application of this molecule, for example in biological tests of its activity.

Statistical analysis of multipole-model-derived structural parameters and charge-density properties from high-resolution X-ray diffraction experiments.

A comprehensive analysis of various properties derived from multiple high-resolution X-ray diffraction experiments is reported. A total of 13 charge-density-quality data sets of α-oxalic acid dihydrate (C2H2O4·2H2O) were subject to Hansen-Coppens-based modelling of electron density. The obtained parameters and properties were then statistically analysed yielding a clear picture of their variability across the different measurements. Additionally, a computational approach (CRYSTAL and PIXEL programs) was utilized to support and examine the experimental findings. The aim of the study was to show the real accuracy and interpretation limits of the charge-density-derived data. An investigation of raw intensities showed that most of the reflections (60-70%) fulfil the normality test and the lowest ratio is observed for weak reflections. It appeared that unit-cell parameters are determined to the order of 10(-3) Š(for cell edges) and 10(-2) ° (for angles), and compare well with the older studies of the same compound and with the new 100 K neutron diffraction data set. Fit discrepancy factors are determined within a 0.5% range, while the residual density extrema are about ±0.16 (3) e Å(-3). The geometry is very well reproducible between different data sets. Regarding the multipole model, the largest errors are present on the valence shell charge-transfer parameters. In addition, symmetry restrictions of multipolar parameters, with respect to local coordinate systems, are well preserved. Standard deviations for electron density are lowest at bond critical points, being especially small for the hydrogen-bonded contacts. The same is true for kinetic and potential energy densities. This is also the case for the electrostatic potential distribution, which is statistically most significant in the hydrogen-bonded regions. Standard deviations for the integrated atomic charges are equal to about 0.1 e. Dipole moments for the water molecule are comparable with the ones presented in various earlier studies. The electrostatic energies should be treated rather qualitatively. However, they are quite well correlated with the PIXEL results.

Towards the best model for H atoms in experimental charge-density refinement.

The consequences of different treatments of H atoms in experimental charge-density studies are discussed. Geometric and topological parameters obtained after applying four different H-atom models in multipolar refinement on high-resolution X-ray data only were compared with the results obtained for a reference joint high-resolution X-ray/neutron refinement. The geometry and the topological critical point and integrated parameters closest to the reference values were obtained after a mixed refinement (high-order refinement of heavy atoms, low-angle refinement of H atoms and elongation of the X-H distance to the average neutron bond lengths) supplemented by an estimation of the anisotropic thermal motions of H atoms using the SHADE program. Such a procedure works very well even for strong hydrogen bonds. The worst fit to the reference results for both critical point and integrated parameters was obtained when only the standardization to the average neutron X-H distances was applied. The non-H-atom parameters are also systematically influenced by the H-atom modeling. In order to compare topological and integrated properties calculated for H and non-H atoms in multipolar refinement when there are no neutron data, the same treatment of H atoms (ideally the mixed refinement + estimated anisotropic atomic displacement parameters for H atoms) should be applied.