Towards accurate and precise positions of hydrogen atoms bonded to heavy metal atoms

Hirshfeld atom refinement and dynamical refinement of hexagonal ice structure from electron diffraction data

Reaching beyond the commonly used spherical atomic electron density model allows one to greatly improve the accuracy of hydrogen atom structural parameters derived from X-ray data. However, the effects of atomic asphericity are less explored for electron diffraction data. In this work, Hirshfeld atom refinement (HAR), a method that uses an accurate description of electron density by quantum mechanical calculation for a system of interest, was applied for the first time to the kinematical refinement of electron diffraction data. This approach was applied here to derive the structure of ordinary hexagonal ice (Ih). The effect of introducing HAR is much less noticeable than in the case of X-ray refinement and it is largely overshadowed by dynamical scattering effects. It led to only a slight change in the O-H bond lengths (shortening by 0.01 Å) compared with the independent atom model (IAM). The average absolute differences in O-H bond lengths between the kinematical refinements and the reference neutron structure were much larger: 0.044 for IAM and 0.046 Å for HAR. The refinement results changed considerably when dynamical scattering effects were modelled - with extinction correction or with dynamical refinement. The latter led to an improvement of the O-H bond length accuracy to 0.021 Å on average (with IAM refinement). Though there is a potential for deriving more accurate structures using HAR for electron diffraction, modelling of dynamical scattering effects seems to be a necessary step to achieve this. However, at present there is no software to support both HAR and dynamical refinement.

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.

Experimental charge density and topological analysis of tetraaquabis(hydrogenmaleato)nickel(II): a comparison with Hirshfeld atom refinement

Experimental charge density analysis is conducted on the coordination compound tetraaquabis(hydrogenmaleato)nickel(II), which exhibits a short intramolecular hydrogen bond. Through topological analysis, the nature of Ni-O bonds is concluded to be intermediate between ionic and covalent, but mainly presenting an ionic character, while the short hydrogen bond is classified as covalent in nature. The compound was also analysed after Hirshfeld atom refinement performed using NoSpherA2. A topological analysis was conducted on the molecular wavefunction and the results are compared with those obtained from experiment. In general, there is good agreement between the refinements, and the chemical bonds involving H atoms are in better agreement with what is expected from neutron data after HAR than they are after multipole refinement.

Hirshfeld Atom Refinement of Metal-Organic Complexes: Treatment of Hydrogen Atoms Bonded to Transition Metals

Hydrogen positions in hydrides play a key role in hydrogen storage materials and high-temperature superconductors. Our recently published study of five crystal structures of transition-metal-bound hydride complexes showed that using aspherical atomic scattering factors for Hirshfeld atom refinement (HAR) resulted in a systematic elongation of metal-hydrogen bonds compared to using spherical scattering factors with the Independent Atom Model (IAM). Even though only standard-resolution X-ray data was used, for the highest-quality data, we obtained excellent agreement between the X-ray and the neutron-derived bond lengths. We present an extended version of this study including 10 crystal structures of metal-organic complexes containing hydrogen atoms bonded to transition-metal atoms for which both X-ray and neutron data are available. The neutron structures were used as a benchmark, and the X-ray structures were refined by applying Hirshfeld atom refinement using various basis sets and DFT functionals in order to investigate the influence of the technical aspects on the length of metal-hydrogen bonds. The result of including relativistic effects in the Hamiltonian and using a cluster of multipoles simulating interactions with a crystal environment during wave function calculations was examined. The effect of the data quality on the final result was also evaluated. The study confirms that a high quality of experimental data is the key factor allowing us to obtain significant improvement in transition metal (TM)-hydrogen bond lengths from HAR in comparison with the IAM. Individual adjustments and better choices of the basis set can improve hydrogen positions. Average differences between TM-H bond lengths obtained with various DFT functionals upon including relativistic effects or between double-ζ and triple-ζ basis sets were not statistically significant. However, if all bonds formed by H atoms were considered, significant differences caused by different refinement strategies were observed. Finally, we examined the refinement of atomic thermal motions. Anisotropic refinement of hydrogen thermal motions with HAR was feasible only in some cases, and isotropically refined hydrogen thermal motions were in similar agreement with neutron values whether obtained with HAR or with the IAM.

How do density functionals affect the Hirshfeld atom refinement?

In this work, the effect of mixing different amounts of Hartree-Fock (HF) exchange with hybrid density functionals applied to the Hirshfeld atom refinement (HAR) of urea and oxalic acid dihydrate is explored. Together, the influence of using different basis sets, methods (including MP2 and HF) and cluster sizes (to model bulk effects) is studied. The results show that changing the amount of HF exchange, no matter the level of theory, has an impact almost exclusively on the H atom refinement parameters. Contrary to pure quantum mechanical calculations where good geometries are obtained with intermediate HF exchange mixtures, in the HAR the best match with neutron diffraction reference values is not necessarily found for these admixtures. While the non-hydrogen covalent bond lengths are insensitive to the combination of method or basis set employed, the X-H bond lengths always increase proportionally to the HF exchange for the analysed systems. This outcome is opposite to what is normally observed from geometry optimisations, i.e., shorter bonds are obtained with greater HF exchange. Additionally, the thermal ellipsoids tend to shrink with larger HF exchange, especially for the H atoms involved in strong hydrogen bonding. Thus, it may be the case that the development of density functionals or basis sets suitable for quantum crystallography should take a different path than those fitted for quantum chemistry calculations.

Accurate crystal structure of ice VI from X-ray diffraction with Hirshfeld atom refinement

Water is an essential chemical compound for living organisms, and twenty of its different crystal solid forms (ices) are known. Still, there are many fundamental problems with these structures such as establishing the correct positions and thermal motions of hydrogen atoms. The list of ice structures is not yet complete as DFT calculations have suggested the existence of additional and - to date - unknown phases. In many ice structures, neither neutron diffraction nor DFT calculations nor X-ray diffraction methods can easily solve the problem of hydrogen atom disorder or accurately determine their anisotropic displacement parameters (ADPs). Here, accurate crystal structures of H2O, D2O and mixed (50%H2O/50%D2O) ice VI obtained by Hirshfeld atom refinement (HAR) of high-pressure single-crystal synchrotron and laboratory X-ray diffraction data are presented. It was possible to obtain O-H/D bond lengths and ADPs for disordered hydrogen atoms which are in good agreement with the corresponding single-crystal neutron diffraction data. These results show that HAR combined with X-ray diffraction can compete with neutron diffraction in detailed studies of polymorphic forms of ice and crystals of other hydrogen-rich compounds. As neutron diffraction is relatively expensive, requires larger crystals which can be difficult to obtain and access to neutron facilities is restricted, cheaper and more accessible X-ray measurements combined with HAR can facilitate the verification of the existing ice polymorphs and the quest for new ones.

Influence of modelling disorder on Hirshfeld atom refinement results of an organo-gold(I) compound

Details of the validation of disorder modelling with Hirshfeld atom refinement (HAR) for a previously investigated organo-gold(I) compound are presented here. The impact of refining disorder on HAR results is discussed using an analysis of the differences of dynamic structure factors. These dynamic structure factor differences are calculated from thermally smeared quantum mechanical electron densities based on wavefunctions that include or exclude electron correlation and relativistic effects. When disorder is modelled, the electron densities stem from a weighted superposition of two (or more) different conformers. Here this is shown to impact the relative importance of electron correlation and relativistic effect estimates expressed by the structure factor magnitudes. The role of disorder modelling is also compared with the effect of the treatment of hydrogen anisotropic displacement parameter (ADP) values and atomic anharmonicity of the gold atom. The analysis of ADP values of gold and disordered carbon atoms showed that the effect of disorder significantly altered carbon ADP values and did not influence those of the gold atom.

Fragmentation and transferability in Hirshfeld atom refinement

Hirshfeld atom refinement (HAR) is one of the most effective methods for obtaining accurate structural parameters for hydrogen atoms from X-ray diffraction data. Unfortunately, it is also relatively computationally expensive, especially for larger molecules due to wavefunction calculations. Here, a fragmentation approach has been tested as a remedy for this problem. It gives an order of magnitude improvement in computation time for larger organic systems and is a few times faster for metal-organic systems at the cost of only minor differences in the calculated structural parameters when compared with the original HAR calculations. Fragmentation was also applied to polymeric and disordered systems where it provides a natural solution to problems that arise when HAR is applied. The concept of fragmentation is closely related to the transferable aspherical atom model (TAAM) and allows insight into possible ways to improve TAAM. Hybrid approaches combining fragmentation with the transfer of atomic densities between chemically similar atoms have been tested. An efficient handling of intermolecular interactions was also introduced for calculations involving fragmentation. When applied in fragHAR (a fragmentation approach for polypeptides) as a replacement for the original approach, it allowed for more efficient calculations. All of the calculations were performed with a locally modified version of Olex2 combined with a development version of discamb2tsc and ORCA. Care was taken to efficiently use the power of multicore processors by simple implementation of load-balancing, which was found to be very important for lowering computational time.

Hirshfeld atom refinement based on projector augmented wave densities with periodic boundary conditions

Hirshfeld atom refinement (HAR) is an X-ray diffraction refinement method that, in numerous publications, has been shown to give H-atom bond lengths in close agreement with neutron diffraction derived values. Presented here is a first evaluation of an approach using densities derived from projector augmented wave (PAW) densities with three-dimensional periodic boundary conditions for HAR. The results show an improvement over refinements that neglect the crystal environment or treat it classically, while being on a par with non-periodic approximations for treating the solid-state environment quantum mechanically. A suite of functionals were evaluated for this purpose, showing that the SCAN and revSCAN functionals are most suited to these types of calculation.

New venues in electron density analysis

We provide a comprehensive overview of the chemical information within the electron density: how to extract information, but also how to obtain and how to assess the quality of the...

Refinement of K[HgI3]·H2O using non-spherical atomic form factors

The crystal structure model for potassium tri-iodido-mercurate(II) monohydrate, K[HgI3]·H2O, based on single-crystal data, was reported 50 years ago [Nyqvist & Johansson (1971 ▸). Acta Chem. Scand. 25, 1615-1629]. We have now redetermined this structure with X-ray diffraction data at 0.70 Å resolution collected at 153 K using Ag Kα radiation. Combined quantum mechanical methods (ORCA) and computation of non-spherical scattering form factors (NoSpherA2) allowed the refinement of the shape of the water mol-ecule with anisotropic H atoms, despite the presence of heavy elements in the crystal. The refined shape of the water mol-ecule via this Hirshfeld refinement is close to that determined for liquid water by neutron diffraction experiments. Moreover, the Laplacian of the electron density clearly shows how electron density accumulates along the O-H σ-valence bonds in the water mol-ecule.