research letters
Hirshfeld atom
and dynamical of hexagonal ice structure from electron diffraction dataaBiological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury 101, Warszawa, Warszawa 02-089, Poland, and bCentre of New Technologies, University of Warsaw, S. Banacha 2c, Warsaw 02-097, Poland
*Correspondence e-mail: ml.chodkiewicz@uw.edu.pl, kwozniak@chem.uw.edu.pl
This article is part of a collection of articles on Quantum Crystallography, and commemorates the 100th anniversary of the development of Quantum Mechanics.
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 h). The effect of introducing HAR is much less noticeable than in the case of X-ray 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 results changed considerably when dynamical scattering effects were modelled – with extinction correction or with dynamical 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.
(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 of electron diffraction data. This approach was applied here to derive the structure of ordinary hexagonal ice (I1. Introduction
Ice is one of the most abundant materials on Earth, common also in many other celestial bodies of the solar system. It is also one of the most important and the most studied materials. The structure of ordinary hexagonal ice (Ih) has long been established by X-ray diffraction (Barnes, 1929) which led to the determination of oxygen positions. Neutron diffraction measurements of heavy ice with powder (Wollan et al., 1949) and single-crystal (Peterson & Levy, 1957) diffraction later allowed the determination of hydrogen atom positions. Neutron measurements revealed hydrogen atom disorder consistent with a model proposed by Pauling (1935), with two half-occupied hydrogen sites lying at lines joining neighbouring oxygen atoms. Further single-crystal neutron diffraction studies (Kuhs & Lehmann, 1981, 1983) have confirmed the unusually long O—H distances: about 1.004 Å when the harmonic model of atomic displacement was used [similar results from powder neutron diffraction were also obtained by Fortes et al. (2015)]. The results from the study by Kuhs & Lehmann (1981, 1983) are used as a reference in this work since they are of very high precision (O—H bond length standard deviation not exceeding 0.6 mÅ). An electron density distribution study with X-ray diffraction (Goto et al., 1990) included anharmonic treatment of atomic displacement for oxygen and hydrogen atom positions. Recently, the structure of ice has been studied with 3D electron diffraction (3D ED) (Martynowycz & Gonen, 2019). 3D ED is an umbrella term (see e.g. Gemmi et al., 2019) for several techniques for collecting ED data in 3D The technique allows for from even very small samples, but data interpretation is hampered by strong dynamical scattering effects. A method dedicated to the treatment of the dynamical scattering effects in 3D ED was developed by Palatinus et al. (2015a) which significantly improves the accuracy of structure parameters (Palatinus et. al, 2015b) compared with the kinematical refinement.
Reference neutron measurement values of the two non-equivalent O—H bond lengths in hexagonal ice are very similar: 1.0040 (6) and 1.0036 (3) Å. Both X-ray and ED-derived structures exhibit O—H bond lengths that are considerably different from the neutron measurement values. Short (below 0.9 Å), in the case of X-rays, and long in the case of 3D ED (1.05 and 1.13 Å), bond lengths may be at least partially attributed to the applied spherically symmetric models of atomic electron density (X-ray) and electrostatic potential (3D ED). Such a model [independent atom model (IAM)] is commonly employed in X-ray and electron crystallographic
It neglects the effects of chemical bonding and interatomic interactions in the solid state. This results in poor structural parameters for the hydrogen atoms – covalent bonds to hydrogen are too short in the case of X-rays (by about 0.1 Å) and it is usually not possible to determine the anisotropic ADPs. In practice, there are two commonly used aspherical models which can be used for deriving (freely refining) the structural parameters of hydrogen: transferable aspherical atom model (TAAM) and Hirshfeld atom (HAR).TAAM is based on the assumption of atomic density transferability between atoms in similar chemical environments and uses the Hansen–Coppens multipole model (Hansen & Coppens, 1978) for the parameterization of atomic densities (Pichon-Pesme et al., 1995). TAAM leads to much more accurate structural parameters than IAM, as shown previously for several databanks of atomic electron density parameters (Zarychta et al., 2007; Domagała et al., 2012; Nassour et al., 2017; Dittrich et al., 2004, 2013; Volkov et al., 2007; Dominiak et al., 2007; Bąk et al., 2011; Jarzembska & Dominiak, 2012; Jha et al., 2020).
HAR (Jayatilaka & Dittrich, 2008; Capelli et al., 2014) uses electron densities calculated for the system of interest which makes it more computationally demanding than TAAM, but avoids assumptions related to transferability and takes into account intermolecular interactions which, in principle, make it more accurate than TAAM. Atomic densities in HAR are obtained by Hirshfeld partition (Hirshfeld, 1977) of the total electron density. Other partitions were also tested (Chodkiewicz et al., 2020) but none were clearly best. Though most of the HAR refinements were performed using wavefunction calculations for isolated systems, it was also possible to use a periodic wavefunction model (Wall, 2016; Ruth et al., 2022). HAR on X-ray diffraction data produced hydrogen bond lengths very similar to those obtained from neutron diffraction experiments (e.g. Capelli et al., 2014; Woińska et al., 2014, 2016; Fugel et al., 2018; Sanjuan-Szklarz et al., 2020; Ruth et al., 2022); however, the difference could be quite significant in the case of hydrogen bonded to heavy elements (Woińska et al., 2021).
TAAM had not yet been used for an ice structure study. HAR had been used only once to obtain an ice VI structure from X-ray data (Chodkiewicz et al., 2022) which lead to a much more accurate structure than that from IAM refinement.
Although the advantages of these models are pretty well established for X-ray diffraction, the effect of atomic asphericity in ED is much less explored. The introduction of 3D electron crystallography (Kolb et al., 2007) started a period of rapid development of ED techniques applicable to determination.
Electrostatic potentials determine electron scattering, which leads to rather different scattering factors for neutral and charged atoms (see e.g. Saha et al., 2022). Application of a partial charge model (Yonekura & Maki-Yonekura, 2016) allowed for improved fitting statistics.
While both TAAM and HAR form factors are based on atomic electron densities with partial atomic charges, they also model the asphericity of the atomic charge distribution. This feature was already modelled 20 years ago by Zhong et al. (2002), who used the superposition of potentials calculated for small fragments to model electrostatic potentials in a larger system. They observed significant differences in scattering factors when compared with free atom superposition, especially for low-resolution reflections. Aspherical scattering factors for electrons were parameterized by Zheng et al. (2009) by parameterizing scattering from p and d orbitals; however, its practical application requires knowledge of the orbital orientation and occupancy.
A et al., 2020). Data for a small molecule (carbamazepine) have been used. Refinements within kinematical approximation were performed against both experimental and theoretically modelled data (to avoid effects originating from dynamical scattering). The conclusions from this work were drawn mostly on the basis of model data since the experimental data gave much higher discrepancies in bond lengths and ADPs for both TAAM and IAM, and some trends visible for the model data were probably obscured by experimental errors or dynamical scattering effects. Overall, in the case of the against model data, a better agreement between the structural parameters and better have been observed for TAAM (compared with IAM), with significantly better ADPs for non-hydrogen atoms and slightly better hydrogen atom positions (which were already quite good with IAM, the average root mean square deviation (RMSD) was 0.02 versus 0.01 Å for TAAM). RMSDs for X—H bond lengths from against experimental data were similar for TAAM and IAM refinements: 0.075 Å for TAAM and 0.078 Å for IAM. TAAM was also used for against experimental ED data for β-glycine (Jha et al., 2021). It was shown to give better and more reasonable X—H bond lengths than IAM.
using aspherical atom form factors against electron scattering data has been performed with TAAM (GruzaIn this work, we test HAR for 3D ED. The techniques are combined for the first time to assess the effect of modelling the asphericity of the atomic electron densities and the corresponding electrostatic potentials by the HAR approach on the kinematical h). A combination of HAR and dynamical would be better suited for such an analysis, but there is no software yet that enables both techniques. Therefore, dynamical scattering effects in this work are modelled only in the case of IAM.
of 3D ED data. The method is applied to hexagonal ice (I2. Measurements
Hexagonal ice crystals eventually appear whenever we prepare grids for data collection at the temperature of liquid nitrogen by absorbing moisture from the surroundings. A glow discharged lacey carbon 200 mesh Cu grid was used for data collection. MicroED samples for ice VI crystals were prepared by directly applying the grid to liquid nitrogen. A Thermo Fisher Scientific Glacios cryo transmission electron microscope (TEM) equipped with a field emission gun and operated at 200 kV and −192°C was used for data collection. The microscope was equipped with a Thermo Fisher Scientific CETA-D detector, an autoloader with a twelve-grid holder and EPU-D software for automated data collection. A 50 µm condenser aperture, spot size 11 and gun lens 8 were set and the diffraction datasets were collected under the parallel illumination condition with a very low dose. The crystal was continuously rotated from −55 to +55° with a tilt speed of 1° s−1 and frames were recorded after each rotation of 0.5°. The single frame exposure time was 0.5 s. A total of 219 image frames were collected. Data reduction was carried out in CrysAlis (Rigaku Corporation). In total, 54 virtual frames were made. Data collection details are given in Table S1 of the supporting information. Raw diffraction images and associated files documenting data reduction various steps are available online (Chodkiewicz et al., 2024) at https://doi.org/10.18150/IVF5BA (Repository for Open Data, Interdisciplinary Centre for Mathematical and Computational Modelling, University of Warsaw, Warsaw, Poland).
3. Refinements
3.1. Kinematical refinements
3.1.1. In Olex2
For both HAR and IAM, three refinements were performed: one with all the reflections, one with some of the reflections omitted (presumably those which deviate most from experimental values, discussed in Section 4), and one with the extinction parameter refined and all reflections included.
IAM Olex2 (version 1.5; Dolomanov et al., 2009) was performed without I/σ(I) cut-off and using the parameterization by Peng (1999).
inIn the case of HAR, the electron atomic form factors fe(s) were calculated from the X-ray form factors fx(s) using the Mott−Bethe formula:
where m0 is the ε0 is the permitivity under vacuum and Z is the atomic number.
Since the ice Ih structure is disordered, calculation of atomic form factors for HAR required incorporation of electron densities from multiple local conformations which could occur in real crystals. Atomic densities obtained with Hirshfeld partition of the total electron density were then averaged over the conformations. We have used clusters of five water molecules (a `central' water molecule and its four nearest neighbours) to generate such local conformations. There are two different conformations of the central molecule: (1) containing H1 and H2; and (2) containing two H1 atoms, in addition each of the four neighbours of the central molecule has three possible conformations. Conformations of the neighbours were not modelled, instead one of them was assumed to provide good representation of the interaction between the central molecule and its neighbours (see Fig. S1). Other possible representations would differ in the position of the hydrogen atom which is not directly involved in the hydrogen bond to the central molecule and, as a consequence, probably only had a minor influence on modelling the effects of the intermolecular interactions. Two possible configurations of water molecules were used in the calculations.
A second-order Møller–Plesset perturbation theory (MP2) calculation with a cc-pVTZ basis set was used for the calculation of the electron densities. Quantum mechanical calculations were performed with Gaussian (Frisch et al., 2016). Details of the HAR procedure are given in the supporting information.
HAR was performed using a locally modified version of Olex2 (Dolomanov et al., 2009). A program based on a development version of the DiSCaMB library (Chodkiewicz et al., 2018) was used to generate files with atomic form factors in tsc format (Kleemiss et al., 2021; Midgley et al., 2019). Such files are then imported into Olex2 and used in the conducted with olex2.refine.
3.1.2. In Jana2020
There were two kinematical IAM refinements performed in Jana2020 (version 1.3.36; Petříček et al., 2023): one with optimization of the extinction parameter and one without. The refinements were performed without using an I/σ(I) cut-off and with electron form factor parameterization from International Tables of Crystallography (Vol. C, ch. 4.3.1.1; Colliex et al., 2006).
3.2. Dynamical IAM refinement
Similar to the kinematical Jana2020, no I/σ(I) cut-off was applied and the same source of form factors was used. The structure obtained from IAM in Olex2 was used as a starting point. All reflections with the maximal relative excitation error (RSg) equal to 0.4 were included in the The number of merged frames was set to six and the step between frames was set to four to obtain virtual frames. A sigma-weighting scheme was applied. Including the dynamical scattering description allowed the inclusion of all strong reflections in the R1 was calculated for the reflections with I > 3σ(I).
in4. Results and discussion
The results of the refinements are reported in Table 1 and additional information is provided in Table S2. R1 values for refinements with Jana2020 are reported for reflections with I > 3σ(I) and for Olex2 refinements for I > 2σ(I).
|
4.1. Kinematical refinements
While HAR gives a slightly better R1 agreement factor than IAM (19.34 versus 21.06%) and slightly cleaner but very similar residual density maps (Fig. S2), both kinematical HAR and IAM lead to quite similar structures (see Fig. 1). The average absolute difference in O—H bond lengths [compared with the neutron structure (Kuhs & Lehmann, 1981) at 60 K] is 0.044 and 0.046 Å for IAM and HAR, respectively. Bonds in structures from HAR are slightly shorter (by about 0.01 Å) than those from IAM, which is an expected effect; similar differences were obtained from IAM and TAAM refinements against theoretical data (Gruza et al., 2020). One of the O—H bonds is shorter than the other by as much as 0.09 Å, whereas they are almost the same [1.004 (6) and 1.0036 (13) Å] in the case of neutron data. A similar feature can be observed in the case of a 3D ED derived structure (Martynowycz & Gonen, 2019) (1.13 and 1.05 Å). In this case, the bond lengths were much larger than from the IAM reported here [1.03 (3) and 0.943 (11) Å]. Certainly the effect of HAR on bond length (0.01 Å) was very small in comparison with the observed inaccuracies.
Equivalent isotropic atomic displacement parameters (Ueq) for hydrogen atoms from IAM and HAR are much smaller than the corresponding neutron diffraction values (on average 0.018 versus 0.027 Å2) and also smaller but to a lesser extent for oxygen atoms (0.012 versus 0.014 Å2). A comparison of calculated and experimental structure factors revealed (Fig. 2, red diamonds) that some intense reflections had significantly lower intensities than the model examples, which might be related to dynamical scattering effects. On average, dynamical scattering increases the intensity of weaker reflections, whereas the strong reflections become less intense. Similarly, decreasing the ADPs would make high-angle intensities larger. As a result, the dynamical scattering effect can be partially modelled by artificial decreasing of ADPs and this might be the reason why the ADPs are too small.
Refinements with some of the strongest reflections omitted (see Fig. 2, the labelled reflections and the supporting information) and refinements with extinction optimization (as a way to model effects of dynamical scattering) were also performed (see Table 1 for a summary of all the refinements). The reflections with the highest difference between observed and calculated intensities were chosen manually based on visual inspection. They lead to an increase of the ADPs close to the neutron values (on average 0.025 Å2 for omitted reflections, 0.0275 Å2 for with extinction included and 0.027 Å2 for neutron measurements). Yet the problem with bond lengths remained at a similar scale.
In addition to kinematical refinements with Olex2, kinematical refinements with Jana2020 were performed as they were more appropriate for comparison with dynamical refinements also performed with Jana2020. Contrary to Olex2, Jana2020 does not use a SHELX-type weighting scheme. Such a scheme was used in the refinements described so far. In addition, Olex2 refinements were based on |F|2 while Jana2020 on |F|. O—H bond lengths from kinematical with Jana2020 were longer than those from Olex2, especially O—H2 (1.088 vs 1.030 Å). Hydrogen atom ADPs were refined isotropically, otherwise non-positive definite ADPs would appear for one of the hydrogen atoms.
4.2. Dynamical IAM refinement
Dynamical IAM R factors, R1 decreased by more than 4.75 percentage points relative to kinematical IAM However, during the procedure, some reflections were automatically omitted (treated as outliers), namely 100, 101, 102 and 103. They were included in the kinematical and were far from being well reproduced (see Fig. 2). The large discrepancy between the model and observed intensities for other strong reflections disappeared (Fig. 2). The bond lengths also significantly improved, i.e. the average absolute difference from reference neutron values dropped from 0.064 to 0.021 Å. However, it was still not possible to refine all hydrogen atoms with anisotropic ADPs.
lowered4.3. Discussion
4.3.1. Observations from previous refinements
An effect of atomic asphericity on the hydrogen atom parameters was quite large in the case of X-ray structure R factors). For discussion of the current results, it was convenient to focus on the accuracy of the lengths of covalent bonds to hydrogen (X—H), since the other parameters – ADPs and R factors – can be strongly influenced by dynamical scattering (Palatinus et al., 2015a; Klar et al., 2023). Inclusion of the asphericity of atomic electron densities in the models usually led to substantial elongation of the X—H bonds (usually by more than 0.1 Å) in the case of X-ray structure but it is predicted to be modest in the case of 3D ED. IAM against quantum mechanically simulated 3D ED intensities for carbamazepine (Gruza et al., 2020) led to slightly long X—H bonds [0.024 Å on average for 0.83 Å resolution and only by 0.006 Å for 0.6 Å resolution (RMSDs of 0.024 and 0.012 Å, respectively)]. Experimental verification of this prediction, however, was hindered by the presence of dynamical scattering effects in 3D ED which is typically not accounted for. The differences from against simulated data were far smaller than those from refinements against experimental data for the same system. Using the neutron structure as a reference gave an RMSD of 0.078 Å for IAM and nearly the same (0.075 Å) for TAAM X—H bond lengths. Even higher differences occurred in the case of β-glycine (Jha et al., 2021). While there was no reference structure to use for comparison in this case, some of the X—H bond lengths were clearly too large, e.g. for C—H it reached 1.32 (8) Å for IAM and 1.17 (5) Å for TAAM whereas the average from the neutron measurements was 1.097 Å (Allen & Bruno, 2010).
It may have also lead to a significant lowering of the agreement factors (4.3.2. Discussion of the current results
In the current study, X—H bonds from IAM are on average only 0.01 Å longer than those from HAR. The differences between IAM and aspherical model values were larger in the previous studies – the bonds obtained with TAAM were on average 0.036 Å longer for carbamazepine (Gruza et al., 2020; 3D ED experimental data) and 0.055 Å longer for β-glycine (Jha et al., 2021) in comparison with the IAM values. In the case of ice, these differences are of the order of standard uncertainties (0.01 and 0.03 Å in the case of HAR) for bond lengths, but the trend seemed to be quite consistent; all X—H bonds from IAM were longer in the case of ice and glycine.
Although we were able to observe the asphericity-related trends in bond lengths, which were predicted using simulated data, the related effects were relatively small. They were much smaller than the effect of modelling dynamical scattering effects with the F| in Jana2020 and on |F|2 in Olex2). Dynamical largely improved the and accuracy of the structure. At present, it is not possible to perform HAR with dynamical and it is also difficult to fully assess its influence on structure with ED data. While there is no additional theoretical development needed for implementing dynamical HAR, no software (that we are aware of) currently supports such refinement.
of extinction parameter or changes related to using different least squares software (probably as a result of using a different weighting scheme and/or different objective function – on |5. Conclusions
Kinematical HAR against 3D ED data has been performed for the first time. HAR brings clear improvement in the case of X-ray diffraction, yet in the case studied here – ice Ih – its advantages were still overshadowed by errors introduced by a lack of modelling of the dynamical scattering not being accounted for.
In the case of the investigated structure, it was possible to freely refine hydrogen atom positions and anisotropic ADPs for both HAR and IAM in kinematical refinement.
The average absolute differences in O—H bond lengths between the reference neutron structure and both for HAR and for IAM were between 0.035 and 0.053 Å, on average 0.044 Å for IAM and 0.046 Å for HAR. In the case of X-ray refinements, such differences for IAM are usually one order of magnitude larger than for HAR, here they are very similar. The two O—H bond lengths from the ED experiment differed quite considerably (by about 0.06–0.07 Å), yet they were almost identical in the case of neutron data. It is expected that application of aspherical atomic electron density models would lead to shorter X—H bonds and indeed bonds from HAR were shorter than those from IAM, but only by 0.01 Å on average. The effects of using HAR instead of IAM in kinematical refinements were quite small in comparison with the inaccuracies of the structures obtained (the 0.01 Å difference in bond lengths versus at least 0.035 Å error in bond lengths). Dynamical was only possible in the case of IAM and considerably improved the resulting structure, the bond lengths became closer to neutron values (0.021 Å on average) and the large discrepancies between observed and calculated intensities were radically reduced.
HAR may potentially produce better structural models from 3D ED than IAM; however, it did not show any advantage over IAM in the structure studied and the differences between HAR- and IAM-derived structures were relatively small. Modelling of the dynamical scattering effects seems to be necessary to take advantage of a more advanced HAR model.
Supporting information
https://doi.org/10.1107/S2052252524006808/vq5006sup1.cif
contains datablocks Jana2020_IAM, JANA2020_IAM_dynamic, Jana2020_IAM_exti, HAR, har_exti, har_omit, Olex_IAM, Olex_iam_exti, Olex_iam_omit. DOI:Supporting tables and figures. DOI: https://doi.org/10.1107/S2052252524006808/vq5006sup2.pdf
H2O | F(000) = 40 |
Mr = 18 | Dx = 1.002 Mg m−3 |
Hexagonal, P63/mmc | Electron radiation, λ = 0.0251 Å |
Hall symbol: -P 6c;-2c | Cell parameters from 1533 reflections |
a = 4.3815 (13) Å | θ = 0.2–1.7° |
c = 7.183 (4) Å | µ = 0 mm−1 |
V = 119.42 (8) Å3 | T = 81 K |
Z = 4 | Plate |
Esperanto-CrysAlisPro-abstract goniometer imported esperanto images diffractometer | Rint = 0.072 |
Radiation source: transmission electron microscope | θmax = 1.6°, θmin = 0.2° |
Transmission electron microscope monochromator | h = −9→9 |
2005 measured reflections | k = −9→9 |
265 independent reflections | l = −14→14 |
199 reflections with I > 3σ(I) |
Refinement on F | 1 constraint |
R[F2 > 2σ(F2)] = 0.210 | All H-atom parameters refined |
wR(F2) = 0.239 | Weighting scheme based on measured s.u.'s w = 1/(σ2(F) + 0.0001F2) |
S = 9.24 | (Δ/σ)max = 0.018 |
265 reflections | Δρmax = 0.18 e Å−3 |
9 parameters | Δρmin = −0.39 e Å−3 |
0 restraints |
x | y | z | Uiso*/Ueq | Occ. (<1) | |
O1 | 0.333333 | 0.666667 | 0.5616 (6) | 0.008 (2) | |
H1 | 0.450 (4) | 0.550 (4) | 0.510 (7) | 0.026 (7)* | 0.5 |
H2 | 0.333333 | 0.666667 | 0.713 (12) | 0.025 (14)* | 0.5 |
O1—H1 | 0.96 (3) | H1—H1iii | 0.78 (3) |
O1—H1i | 0.96 (3) | H2—H2iv | 0.53 (12) |
O1—H1ii | 0.96 (3) | ||
H1—O1—H2 | 112.72 | H1i—O1—H1 | 106.04 |
Symmetry codes: (i) −y+1, x−y+1, z; (ii) −x+y, −x+1, z; (iii) y, x, −z+1; (iv) −y+1, −x+1, −z+3/2. |
D—H···A | D—H | H···A | D···A | D—H···A |
O1—H1···O1iii | 0.96 (3) | 1.73 (3) | 2.680 (2) | 175 (5) |
O1—H2···O1iv | 1.09 (8) | 1.62 (8) | 2.706 (7) | ** |
Symmetry codes: (iii) y, x, −z+1; (iv) −y+1, −x+1, −z+3/2. |
H2O | F(000) = 40 |
Mr = 18 | Dx = 1.002 Mg m−3 |
Hexagonal, P63/mmc | Electron radiation, λ = 0.0251 Å |
Hall symbol: -P 6c;-2c | Cell parameters from 1533 reflections |
a = 4.3815 (13) Å | θ = 0.2–1.7° |
c = 7.183 (4) Å | µ = 0 mm−1 |
V = 119.42 (8) Å3 | T = 81 K |
Z = 4 | Plate |
Esperanto-CrysAlisPro-abstract goniometer imported esperanto images diffractometer | 623 reflections with I > 3σ(I) |
Radiation source: transmission electron microscope | θmax = 1.6°, θmin = 0.2° |
Transmission electron microscope monochromator | h = −9→9 |
5180 measured reflections | k = −9→9 |
1106 independent reflections | l = −14→14 |
Refinement on F | 1 constraint |
R[F2 > 2σ(F2)] = 0.103 | All H-atom parameters refined |
wR(F2) = 0.129 | Weighting scheme based on measured s.u.'s w = 1/(σ2(F) + 0.0001F2) |
S = 2.54 | (Δ/σ)max < 0.001 |
1106 reflections | Δρmax = 0.10 e Å−3 |
64 parameters | Δρmin = −0.25 e Å−3 |
0 restraints |
x | y | z | Uiso*/Ueq | Occ. (<1) | |
O1 | 0.333333 | 0.666667 | 0.5622 (2) | 0.0107 (3) | |
H1 | 0.4554 (12) | 0.5446 (12) | 0.5107 (18) | 0.0169 (16)* | 0.5 |
H2 | 0.333333 | 0.666667 | 0.707 (2) | 0.015 (3)* | 0.5 |
O1—H1 | 0.998 (8) | O1—H2 | 1.037 (18) |
O1—H1i | 0.998 (8) | H1—H1iii | 0.694 (10) |
O1—H1ii | 0.998 (8) | H2—H2iv | 0.62 (2) |
H1—O1—H2 | 111.77 | H1i—O1—H1 | 107.08 |
Symmetry codes: (i) −y+1, x−y+1, z; (ii) −x+y, −x+1, z; (iii) y, x, −z+1; (iv) −y+1, −x+1, −z+3/2. |
D—H···A | D—H | H···A | D···A | D—H···A |
O1—H1···O1iii | 0.998 (8) | 1.687 (7) | 2.6829 (10) | 176.3 (12) |
O1—H2···O1iv | 1.037 (18) | 1.661 (18) | 2.698 (4) | ** |
Symmetry codes: (iii) y, x, −z+1; (iv) −y+1, −x+1, −z+3/2. |
H1.999998O1.000002 | F(000) = 40 |
Mr = 18 | Dx = 1.002 Mg m−3 |
Hexagonal, P63/mmc | Electron radiation, λ = 0.0251 Å |
Hall symbol: -P 6c;-2c | Cell parameters from 1533 reflections |
a = 4.3815 (13) Å | θ = 0.2–1.7° |
c = 7.183 (4) Å | µ = 0 mm−1 |
V = 119.42 (8) Å3 | T = 81 K |
Z = 4 | Plate |
Esperanto-CrysAlisPro-abstract goniometer imported esperanto images diffractometer | Rint = 0.072 |
Radiation source: transmission electron microscope | θmax = 1.6°, θmin = 0.2° |
Transmission electron microscope monochromator | h = −9→9 |
2005 measured reflections | k = −9→9 |
265 independent reflections | l = −14→14 |
199 reflections with I > 3σ(I) |
Refinement on F | All H-atom parameters refined |
R[F2 > 2σ(F2)] = 0.143 | Weighting scheme based on measured s.u.'s w = 1/(σ2(F) + 0.0001F2) |
wR(F2) = 0.166 | (Δ/σ)max = 0.003 |
S = 6.45 | Δρmax = 0.19 e Å−3 |
265 reflections | Δρmin = −0.38 e Å−3 |
10 parameters | Extinction correction: SHELXL-2017/1 (Sheldrick, 2015), Fc*=kFc[1+0.001xFc2λ3/sin(2θ)]-1/4 |
0 restraints | Extinction coefficient: 430 (70) |
1 constraint |
x | y | z | Uiso*/Ueq | Occ. (<1) | |
O1 | 0.333333 | 0.666667 | 0.5618 (4) | 0.0206 (19) | |
H1 | 0.443 (3) | 0.557 (3) | 0.521 (4) | 0.023 (4)* | 0.5 |
H2 | 0.333333 | 0.666667 | 0.706 (7) | 0.027 (9)* | 0.5 |
O1—H1 | 0.885 (16) | O1—H2 | 1.04 (5) |
O1—H1i | 0.885 (16) | H1—H1iii | 0.91 (2) |
O1—H1ii | 0.885 (16) | H2—H2iv | 0.63 (7) |
H1—O1—H2 | 109.4 | H1i—O1—H1 | 109.54 |
Symmetry codes: (i) −y+1, x−y+1, z; (ii) −x+y, −x+1, z; (iii) y, x, −z+1; (iv) −y+1, −x+1, −z+3/2. |
D—H···A | D—H | H···A | D···A | D—H···A |
O1—H1···O1iii | 0.885 (16) | 1.796 (16) | 2.6808 (15) | 180 (3) |
O1—H2···O1iv | 1.04 (5) | 1.67 (5) | 2.704 (5) | ** |
Symmetry codes: (iii) y, x, −z+1; (iv) −y+1, −x+1, −z+3/2. |
H2O | Dx = 1.002 Mg m−3 |
Mr = 18 | Electron radiation, λ = 0.0251 Å |
Hexagonal, P63/mmc | Cell parameters from 1533 reflections |
a = 4.3815 (13) Å | θ = 0.2–1.7° |
c = 7.183 (4) Å | µ = 0.000 mm−1 |
V = 119.42 (8) Å3 | T = 81 K |
Z = 4 | Plate |
F(000) = 12.160 |
Esperanto-CrysAlisPro-abstract goniometer imported esperanto images diffractometer | θmax = 1.6°, θmin = 0.2° |
2005 measured reflections | h = −9→9 |
265 independent reflections | k = −9→9 |
215 reflections with I ≥ 2u(I) | l = −14→14 |
Rint = 0.072 |
Refinement on F2 | 0 restraints |
Least-squares matrix: full | 0 constraints |
R[F2 > 2σ(F2)] = 0.193 | All H-atom parameters refined |
wR(F2) = 0.387 | w = 1/[σ2(Fo2) + (0.171P)2] where P = (Fo2 + 2Fc2)/3 |
S = 1.77 | (Δ/σ)max = 0.0002 |
265 reflections | Δρmax = 0.19 e Å−3 |
13 parameters | Δρmin = −0.41 e Å−3 |
Refinement. Aspherical atom refinement in Olex2 with DiSCaMB generated form factors Hirshfeld Atom Refinement using Gaussian 16 generated electron density calculated at MP2/cc-pVTZ level of theory, two possible configurations of a water molecule were taken into account, both modelled with a cluster of 5 water molecules i.e. including 1 neighbors of the central molecule of water to mimic crystal field atomic electron scattering factors calculated from X-ray form factors with Mott-Bethe formula |
x | y | z | Uiso*/Ueq | Occ. (<1) | |
O1 | 0.333333 | 0.666667 | 0.4378 (3) | 0.0118 (6) | |
H1 | 0.103 (3) | 0.5513 (16) | 0.481 (2) | 0.018 (2) | 0.500000 |
H2 | 0.333333 | 0.666667 | 0.296 (4) | 0.019 (5) | 0.499980 |
U11 | U22 | U33 | U12 | U13 | U23 | |
O1 | 0.0116 (6) | 0.0116 (6) | 0.0122 (10) | 0.0058 (3) | −0.000000 | 0.000000 |
H1 | 0.008 (4) | 0.023 (5) | 0.019 (7) | 0.004 (2) | 0.007 (4) | 0.003 (2) |
H2 | 0.026 (9) | 0.026 (9) | 0.003 (9) | 0.013 (4) | −0.000000 | 0.000000 |
O1—H1i | 0.930 (11) | O1—H2 | 1.02 (3) |
O1—H1 | 0.930 (11) | H1—H1iii | 0.82 (2) |
O1—H1ii | 0.930 (11) | H2—H2iv | 0.66 (6) |
H1—O1—H1i | 109.3 (11) | H2—O1—H1ii | 109.6 (11) |
H1ii—O1—H1i | 109.3 (10) | H2—O1—H1i | 109.6 (11) |
H1ii—O1—H1 | 109.3 (11) | H1iii—H1—O1ii | 179 (4) |
H2—O1—H1 | 109.6 (11) | H2iv—H2—O1 | 180.0 |
Symmetry codes: (i) −y+1, x−y+1, z; (ii) −x+y, −x+1, z; (iii) −x, −x+y, −z+1; (iv) −y+1, −x+1, −z+1/2. |
H2O | Dx = 1.002 Mg m−3 |
Mr = 18 | Electron radiation, λ = 0.0251 Å |
Hexagonal, P63/mmc | Cell parameters from 1533 reflections |
a = 4.3815 (13) Å | θ = 0.2–1.7° |
c = 7.183 (4) Å | µ = 0.000 mm−1 |
V = 119.42 (8) Å3 | T = 81 K |
Z = 4 | Plate |
F(000) = 12.160 |
Esperanto-CrysAlisPro-abstract goniometer imported esperanto images diffractometer | θmax = 1.6°, θmin = 0.2° |
2005 measured reflections | h = −9→9 |
265 independent reflections | k = −9→9 |
215 reflections with I ≥ 2u(I) | l = −14→14 |
Rint = 0.072 |
Refinement on F2 | 0 restraints |
Least-squares matrix: full | 0 constraints |
R[F2 > 2σ(F2)] = 0.114 | All H-atom parameters refined |
wR(F2) = 0.349 | w = 1/[σ2(Fo2) + (0.2P)2] where P = (Fo2 + 2Fc2)/3 |
S = 1.38 | (Δ/σ)max = 0.001 |
265 reflections | Δρmax = 0.08 e Å−3 |
14 parameters | Δρmin = −0.21 e Å−3 |
Refinement. Aspherical atom refinement in Olex2 with DiSCaMB generated form factors Hirshfeld Atom Refinement using Gaussian 16 generated electron density calculated at MP2/cc-pVTZ level of theory, two possible configurations of a water molecule were taken into account, both modelled with a cluster of 5 water molecules i.e. including 1 neighbors of the central molecule of water to mimic crystal field atomic electron scattering factors calculated from X-ray form factors with Mott-Bethe formula |
x | y | z | Uiso*/Ueq | Occ. (<1) | |
O1 | 0.333333 | 0.666667 | 0.4377 (2) | 0.0176 (8) | |
H1a | 0.100 (2) | 0.5501 (12) | 0.4793 (19) | 0.027 (2) | 0.500000 |
H1b | 0.333333 | 0.666667 | 0.296 (3) | 0.028 (5) | 0.499980 |
U11 | U22 | U33 | U12 | U13 | U23 | |
O1 | 0.0172 (8) | 0.0172 (8) | 0.0184 (11) | 0.0086 (4) | −0.000000 | 0.000000 |
H1a | 0.015 (4) | 0.031 (4) | 0.031 (7) | 0.0077 (19) | 0.010 (4) | 0.0048 (19) |
H1b | 0.035 (8) | 0.035 (8) | 0.015 (10) | 0.018 (4) | −0.000000 | 0.000000 |
O1—H1ai | 0.933 (9) | O1—H1b | 1.01 (2) |
O1—H1a | 0.933 (9) | H1a—H1aiii | 0.817 (17) |
O1—H1aii | 0.933 (9) | H1b—H1biv | 0.67 (5) |
H1a—O1—H1ai | 110.3 (8) | H1b—O1—H1aii | 108.6 (9) |
H1aii—O1—H1ai | 110.3 (7) | H1b—O1—H1ai | 108.6 (9) |
H1aii—O1—H1a | 110.3 (8) | H1aiii—H1a—O1ii | 177 (3) |
H1b—O1—H1a | 108.6 (9) | H1biv—H1b—O1 | 180.0 |
Symmetry codes: (i) −y+1, x−y+1, z; (ii) −x+y, −x+1, z; (iii) −x, −x+y, −z+1; (iv) −y+1, −x+1, −z+1/2. |
H2O | Dx = 1.002 Mg m−3 |
Mr = 18 | Electron radiation, λ = 0.0251 Å |
Hexagonal, P63/mmc | Cell parameters from 1533 reflections |
a = 4.3815 (13) Å | θ = 0.2–1.7° |
c = 7.183 (4) Å | µ = 0.000 mm−1 |
V = 119.42 (8) Å3 | T = 81 K |
Z = 4 | Plate |
F(000) = 12.160 |
Esperanto-CrysAlisPro-abstract goniometer imported esperanto images diffractometer | θmax = 1.6°, θmin = 0.4° |
2005 measured reflections | h = −9→9 |
255 independent reflections | k = −9→9 |
205 reflections with I ≥ 2u(I) | l = −14→14 |
Rint = 0.072 |
Refinement on F2 | 0 restraints |
Least-squares matrix: full | 0 constraints |
R[F2 > 2σ(F2)] = 0.106 | All H-atom parameters refined |
wR(F2) = 0.296 | w = 1/[σ2(Fo2) + (0.1566P)2 + 0.0022P] where P = (Fo2 + 2Fc2)/3 |
S = 1.05 | (Δ/σ)max = 0.0004 |
255 reflections | Δρmax = 0.08 e Å−3 |
13 parameters | Δρmin = −0.15 e Å−3 |
Refinement. Aspherical atom refinement in Olex2 with DiSCaMB generated form factors Hirshfeld Atom Refinement using Gaussian 16 generated electron density calculated at MP2/cc-pVTZ level of theory, two possible configurations of a water molecule were taken into account, both modelled with a cluster of 5 water molecules i.e. including 1 neighbors of the central molecule of water to mimic crystal field atomic electron scattering factors calculated from X-ray form factors with Mott-Bethe formula |
x | y | z | Uiso*/Ueq | Occ. (<1) | |
O1 | 0.333333 | 0.666667 | 0.4380 (2) | 0.0144 (5) | |
H1 | 0.106 (4) | 0.5529 (19) | 0.479 (2) | 0.025 (3) | 0.500000 |
H2 | 0.333333 | 0.666667 | 0.300 (3) | 0.025 (5) | 0.499980 |
U11 | U22 | U33 | U12 | U13 | U23 | |
O1 | 0.0143 (6) | 0.0143 (6) | 0.0147 (8) | 0.0071 (3) | −0.000000 | 0.000000 |
H1 | 0.018 (5) | 0.031 (6) | 0.022 (6) | 0.009 (3) | 0.006 (5) | 0.003 (2) |
H2 | 0.033 (9) | 0.033 (9) | 0.010 (9) | 0.017 (4) | −0.000000 | 0.000000 |
O1—H1i | 0.911 (13) | O1—H2 | 0.99 (2) |
O1—H1 | 0.911 (13) | H1—H1iii | 0.86 (3) |
O1—H1ii | 0.911 (13) | H2—H2iv | 0.72 (5) |
H1—O1—H1i | 110.2 (9) | H2—O1—H1ii | 108.7 (11) |
H1ii—O1—H1i | 110.2 (8) | H2—O1—H1i | 108.7 (11) |
H1ii—O1—H1 | 110.2 (9) | H1iii—H1—O1ii | 178 (3) |
H2—O1—H1 | 108.7 (11) | H2iv—H2—O1 | 180.0 |
Symmetry codes: (i) −y+1, x−y+1, z; (ii) −x+y, −x+1, z; (iii) −x, −x+y, −z+1; (iv) −y+1, −x+1, −z+1/2. |
H2O | Dx = 1.002 Mg m−3 |
Mr = 18 | Electron radiation, λ = 0.02510 Å |
Hexagonal, P63/mmc | Cell parameters from 1533 reflections |
a = 4.3815 (13) Å | θ = 0.2–1.7° |
c = 7.183 (4) Å | µ = 0.000 mm−1 |
V = 119.42 (8) Å3 | T = 81 K |
Z = 4 | Plate |
F(000) = 12.160 |
Esperanto-CrysAlisPro-abstract goniometer imported esperanto images diffractometer | θmax = 1.6°, θmin = 0.2° |
2005 measured reflections | h = −9→9 |
265 independent reflections | k = −9→9 |
215 reflections with I ≥ 2u(I) | l = −14→14 |
Rint = 0.072 |
Refinement on F2 | 0 restraints |
Least-squares matrix: full | 0 constraints |
R[F2 > 2σ(F2)] = 0.211 | All H-atom parameters refined |
wR(F2) = 0.400 | w = 1/[σ2(Fo2) + (0.2P)2] where P = (Fo2 + 2Fc2)/3 |
S = 1.62 | (Δ/σ)max = 0.001 |
265 reflections | Δρmax = 0.74 e Å−3 |
13 parameters | Δρmin = −1.45 e Å−3 |
x | y | z | Uiso*/Ueq | Occ. (<1) | |
O1 | 0.333333 | 0.666667 | 0.4378 (3) | 0.0116 (6) | |
H1a | 0.099 (3) | 0.5496 (15) | 0.482 (2) | 0.018 (2) | 0.500000 |
H1b | 0.333333 | 0.666667 | 0.294 (4) | 0.018 (5) | 0.500000 |
U11 | U22 | U33 | U12 | U13 | U23 | |
O1 | 0.0114 (6) | 0.0114 (6) | 0.0119 (10) | 0.0057 (3) | −0.000000 | 0.000000 |
H1a | 0.008 (4) | 0.022 (5) | 0.020 (7) | 0.004 (2) | 0.008 (4) | 0.004 (2) |
H1b | 0.025 (9) | 0.025 (9) | 0.003 (9) | 0.013 (4) | −0.000000 | 0.000000 |
O1—H1ai | 0.943 (11) | O1—H1b | 1.03 (3) |
O1—H1a | 0.943 (11) | H1a—H1aiii | 0.80 (2) |
O1—H1aii | 0.943 (11) | H1b—H1biv | 0.64 (6) |
H1a—O1—H1ai | 109.5 (11) | H1b—O1—H1aii | 109.5 (11) |
H1aii—O1—H1ai | 109.5 (10) | H1b—O1—H1ai | 109.5 (11) |
H1aii—O1—H1a | 109.5 (11) | H1aiii—H1a—O1ii | 180 (4) |
H1b—O1—H1a | 109.5 (11) | H1biv—H1b—O1 | 180.0 |
Symmetry codes: (i) −y+1, x−y+1, z; (ii) −x+y, −x+1, z; (iii) −x, −x+y, −z+1; (iv) −y+1, −x+1, −z+1/2. |
H2O | Dx = 1.002 Mg m−3 |
Mr = 18 | Electron radiation, λ = 0.02510 Å |
Hexagonal, P63/mmc | Cell parameters from 1533 reflections |
a = 4.3815 (13) Å | θ = 0.2–1.7° |
c = 7.183 (4) Å | µ = 0.000 mm−1 |
V = 119.42 (8) Å3 | T = 81 K |
Z = 4 | Plate |
F(000) = 12.160 |
Esperanto-CrysAlisPro-abstract goniometer imported esperanto images diffractometer | θmax = 1.6°, θmin = 0.2° |
2005 measured reflections | h = −9→9 |
265 independent reflections | k = −9→9 |
215 reflections with I ≥ 2u(I) | l = −14→14 |
Rint = 0.072 |
Refinement on F2 | 0 restraints |
Least-squares matrix: full | 0 constraints |
R[F2 > 2σ(F2)] = 0.109 | All H-atom parameters refined |
wR(F2) = 0.340 | w = 1/[σ2(Fo2) + (0.2P)2] where P = (Fo2 + 2Fc2)/3 |
S = 1.35 | (Δ/σ)max = −0.0003 |
265 reflections | Δρmax = 0.25 e Å−3 |
14 parameters | Δρmin = −0.65 e Å−3 |
x | y | z | Uiso*/Ueq | Occ. (<1) | |
O1 | 0.333333 | 0.666667 | 0.4378 (2) | 0.0174 (8) | |
H1a | 0.096 (2) | 0.5479 (12) | 0.4802 (18) | 0.029 (2) | 0.500000 |
H1b | 0.333333 | 0.666667 | 0.294 (3) | 0.028 (5) | 0.500000 |
U11 | U22 | U33 | U12 | U13 | U23 | |
O1 | 0.0171 (8) | 0.0171 (8) | 0.0181 (11) | 0.0085 (4) | −0.000000 | 0.000000 |
H1a | 0.015 (4) | 0.033 (4) | 0.033 (7) | 0.0076 (18) | 0.013 (4) | 0.0063 (19) |
H1b | 0.036 (8) | 0.036 (8) | 0.012 (9) | 0.018 (4) | −0.000000 | 0.000000 |
O1—H1ai | 0.951 (8) | O1—H1b | 1.03 (2) |
O1—H1a | 0.951 (8) | H1a—H1aiii | 0.781 (16) |
O1—H1aii | 0.951 (8) | H1b—H1biv | 0.64 (4) |
H1a—O1—H1ai | 110.2 (8) | H1b—O1—H1aii | 108.7 (8) |
H1aii—O1—H1ai | 110.2 (7) | H1b—O1—H1ai | 108.7 (8) |
H1aii—O1—H1a | 110.2 (8) | H1aiii—H1a—O1ii | 177 (3) |
H1b—O1—H1a | 108.7 (8) | H1biv—H1b—O1 | 180.0 |
Symmetry codes: (i) −y+1, x−y+1, z; (ii) −x+y, −x+1, z; (iii) −x, −x+y, −z+1; (iv) −y+1, −x+1, −z+1/2. |
H2O | Dx = 1.002 Mg m−3 |
Mr = 18 | Electron radiation, λ = 0.0251 Å |
Hexagonal, P63/mmc | Cell parameters from 1533 reflections |
a = 4.3815 (13) Å | θ = 0.2–1.7° |
c = 7.183 (4) Å | µ = 0.000 mm−1 |
V = 119.42 (8) Å3 | T = 81 K |
Z = 4 | Plate |
F(000) = 12.160 |
Esperanto-CrysAlisPro-abstract goniometer imported esperanto images diffractometer | θmax = 1.6°, θmin = 0.4° |
2005 measured reflections | h = −9→9 |
255 independent reflections | k = −9→9 |
205 reflections with I ≥ 2u(I) | l = −14→14 |
Rint = 0.072 |
Refinement on F2 | 0 constraints |
Least-squares matrix: full | Primary atom site location: iterative |
R[F2 > 2σ(F2)] = 0.103 | All H-atom parameters refined |
wR(F2) = 0.288 | w = 1/[σ2(Fo2) + (0.1584P)2 + 0.0023P] where P = (Fo2 + 2Fc2)/3 |
S = 1.01 | (Δ/σ)max = 0.001 |
255 reflections | Δρmax = 0.08 e Å−3 |
13 parameters | Δρmin = −0.15 e Å−3 |
0 restraints |
x | y | z | Uiso*/Ueq | Occ. (<1) | |
O1 | 0.333333 | 0.666667 | 0.4380 (2) | 0.0142 (5) | |
H1 | 0.100 (3) | 0.5499 (17) | 0.479 (2) | 0.026 (3) | 0.500000 |
H2 | 0.333333 | 0.666667 | 0.298 (3) | 0.025 (5) | 0.499980 |
U11 | U22 | U33 | U12 | U13 | U23 | |
O1 | 0.0141 (5) | 0.0141 (5) | 0.0144 (8) | 0.0070 (3) | −0.000000 | 0.000000 |
H1 | 0.016 (5) | 0.032 (6) | 0.024 (6) | 0.008 (2) | 0.009 (5) | 0.005 (2) |
H2 | 0.034 (9) | 0.034 (9) | 0.008 (8) | 0.017 (4) | −0.000000 | 0.000000 |
O1—H1i | 0.935 (12) | O1—H2 | 1.00 (2) |
O1—H1 | 0.935 (12) | H1—H1iii | 0.81 (2) |
O1—H1ii | 0.935 (12) | H2—H2iv | 0.69 (5) |
H1—O1—H1i | 110.4 (9) | H2—O1—H1ii | 108.5 (10) |
H1ii—O1—H1i | 110.4 (8) | H2—O1—H1i | 108.5 (10) |
H1ii—O1—H1 | 110.4 (9) | H1iii—H1—O1ii | 177 (3) |
H2—O1—H1 | 108.5 (10) | H2iv—H2—O1 | 180.0 |
Symmetry codes: (i) −y+1, x−y+1, z; (ii) −x+y, −x+1, z; (iii) −x, −x+y, −z+1; (iv) −y+1, −x+1, −z+1/2. |
Acknowledgements
We are grateful to the Rigaku Corporation for sharing their version of the CrysAlis software with us which allowed us to process the ED data.
Conflict of interest
The authors declare no competing interests.
Data availability
The authors confirm that the data supporting the findings of this study are available within the article and its supporting information.
Funding information
Support of this work by the National Science Centre, Poland, through OPUS 21 (grant No. DEC-2021/41/B/ST4/03010 awarded to KW) is gratefully acknowledged. This work was partly carried out at the Biological and Chemical Research Centre, University of Warsaw, established within the project co-financed by European Union from the European Regional Development Fund under the Operational Programme Innovative Economy, 2007–2. The work was accomplished at the TEAM TECH Core Facility for Crystallographic and Biophysical Research to support the development of medicinal products sponsored by the Foundation for Polish Science. The Authors acknowledge the Ministry of Science and Higher Education programme `Excellence Initiative – Research University (2020–2026) (IDUB)' awarded to the University of Warsaw, in particular, the 1.4.1. Activity `Strengthening the core-facility potential on the Ochota Campus' and `Infrastructure for Cryomicroscopy and Electron Diffraction Core Facility' and I.4.2 (grant No. BOB-IDUB-622–20/2021). Funding for the Renovation and Development of Research Infrastructure – development and maintenance of infrastructure – is acknowledged.
References
Allen, F. H. & Bruno, I. J. (2010). Acta Cryst. B66, 380–386. Web of Science CrossRef CAS IUCr Journals Google Scholar
Bąk, J. M., Domagała, S., Hübschle, C., Jelsch, C., Dittrich, B. & Dominiak, P. M. (2011). Acta Cryst. A67, 141–153. Web of Science CrossRef IUCr Journals Google Scholar
Barnes, W. H. (1929). Proc. R. Soc. Lond. A Math. Phys. Sci. 125, 670–693. CAS Google Scholar
Capelli, S. C., Bürgi, H.-B., Dittrich, B., Grabowsky, S. & Jayatilaka, D. (2014). IUCrJ, 1, 361–379. Web of Science CSD CrossRef CAS PubMed IUCr Journals Google Scholar
Chodkiewicz, M. L., Gajda, R., Lavina, B., Tkachev, S., Prakapenka, V. B., Dera, P. & Wozniak, K. (2022). IUCrJ, 9, 573–579. Web of Science CSD CrossRef ICSD CAS PubMed IUCr Journals Google Scholar
Chodkiewicz, M. L., Migacz, S., Rudnicki, W., Makal, A., Kalinowski, J. A., Moriarty, N. W., Grosse-Kunstleve, R. W., Afonine, P. V., Adams, P. D. & Dominiak, P. M. (2018). J. Appl. Cryst. 51, 193–199. Web of Science CrossRef CAS IUCr Journals Google Scholar
Chodkiewicz, M. L., Olech, B., Jha, K. K., Dominiak, P. & Woźniak, K. (2024). Repository for Open Data (RepOD), https://doi.org/10.18150/IVF5BA. Google Scholar
Chodkiewicz, M. L., Woińska, M. & Woźniak, K. (2020). IUCrJ, 7, 1199–1215. Web of Science CSD CrossRef CAS PubMed IUCr Journals Google Scholar
Colliex, C., Cowley, J. M., Dudarev, S. L., Fink, M., Gjønnes, J., Hilderbrandt, R., Howie, A., Lynch, D. F., Peng, L. M., Ren, G., Ross, A. W., Smith, V. H. Jr, Spence, J. C. H., Steeds, J. W., Wang, J., Whelan, M. J. & Zvyagin, B. B. (2006). International Tables for Crystallography. Vol. C, edited by E. Prince, pp. 259–429. Chester: International Union of Crystallography. Google Scholar
Dittrich, B., Hübschle, C. B., Pröpper, K., Dietrich, F., Stolper, T. & Holstein, J. J. (2013). Acta Cryst. B69, 91–104. CrossRef CAS IUCr Journals Google Scholar
Dittrich, B., Koritsánszky, T. & Luger, P. (2004). Angew. Chem. Int. Ed. 43, 2718–2721. Web of Science CSD CrossRef CAS Google Scholar
Dolomanov, O. V., Bourhis, L. J., Gildea, R. J., Howard, J. A. K. & Puschmann, H. (2009). J. Appl. Cryst. 42, 339–341. Web of Science CrossRef CAS IUCr Journals Google Scholar
Domagała, S., Fournier, B., Liebschner, D., Guillot, B. & Jelsch, C. (2012). Acta Cryst. A68, 337–351. Web of Science CrossRef IUCr Journals Google Scholar
Dominiak, P. M., Volkov, A., Li, X., Messerschmidt, M. & Coppens, P. (2007). J. Chem. Theory Comput. 3, 232–247. Web of Science CrossRef CAS PubMed Google Scholar
Fortes, A. D. (2015). Powder Diffr. 30, 149–157. Web of Science CrossRef ICSD CAS Google Scholar
Frisch, M. J., Trucks, G. W., Schlegel, H. B., Scuseria, G. E., Robb, M. A., Cheeseman, J. R., Scalmani, G., Barone, V., Petersson, G. A., Nakatsuji, H., Li, X., Caricato, M., Marenich, A. V., Bloino, J., Janesko, B. G., Gomperts, R., Mennucci, B., Hratchian, H. P., Ortiz, J. V., Izmaylov, A. F., Sonnenberg, J. L., Williams-Young, D., Ding, F., Lipparini, F., Egidi, F., Goings, J., Peng, B., Petrone, A., Henderson, T., Ranasinghe, D., Zakrzewski, V. G., Gao, J., Rega, N., Zheng, G., Liang, W., Hada, M., Ehara, M., Toyota, K., Fukuda, R., Hasegawa, J., Ishida, M., Nakajima, T., Honda, Y., Kitao, O., Nakai, H., Vreven, T., Throssell, K., Montgomery, J. A. Jr, Peralta, J. E., Ogliaro, F., Bearpark, M. J., Heyd, J. J., Brothers, E. N., Kudin, K. N., Staroverov, V. N., Keith, T. A., Kobayashi, R., Normand, J., Raghavachari, K., Rendell, A. P., Burant, J. C., Iyengar, S. S., Tomasi, J., Cossi, M., Millam, J. M., Klene, M., Adamo, C., Cammi, R., Ochterski, J. W., Martin, R. L., Morokuma, K., Farkas, O., Foresman, J. B. & Fox, D. J. (2016). Gaussian 16, Revision C. 0, Gaussian, Inc., Wallingford CT, 2016. Google Scholar
Fugel, M., Jayatilaka, D., Hupf, E., Overgaard, J., Hathwar, V. R., Macchi, P., Turner, M. J., Howard, J. A. K., Dolomanov, O. V., Puschmann, H., Iversen, B. B., Bürgi, H.-B. & Grabowsky, S. (2018). IUCrJ, 5, 32–44. Web of Science CSD CrossRef CAS PubMed IUCr Journals Google Scholar
Gemmi, M., Mugnaioli, E., Gorelik, T. E., Kolb, U., Palatinus, L., Boullay, P., Hovmöller, S. & Abrahams, J. P. (2019). ACS Cent. Sci. 5, 1315–1329. Web of Science CrossRef CAS PubMed Google Scholar
Goto, A., Hondoh, T. & Mae, S. (1990). J. Chem. Phys. 93, 1412–1417. CrossRef ICSD CAS Web of Science Google Scholar
Gruza, B., Chodkiewicz, M. L., Krzeszczakowska, J. & Dominiak, P. M. (2020). Acta Cryst. A76, 92–109. Web of Science CrossRef IUCr Journals Google Scholar
Hansen, N. K. & Coppens, P. (1978). Acta Cryst. A34, 909–921. CrossRef CAS IUCr Journals Web of Science Google Scholar
Hirshfeld, F. L. (1977). Theor. Chim. Acta, 44, 129–138. CrossRef CAS Web of Science Google Scholar
Jarzembska, K. N. & Dominiak, P. M. (2012). Acta Cryst. A68, 139–147. Web of Science CrossRef CAS IUCr Journals Google Scholar
Jayatilaka, D. & Dittrich, B. (2008). Acta Cryst. A64, 383–393. Web of Science CrossRef CAS IUCr Journals Google Scholar
Jha, K. K., Gruza, B., Chodkiewicz, M. L., Jelsch, C. & Dominiak, P. M. (2021). J. Appl. Cryst. 54, 1234–1243. Web of Science CSD CrossRef CAS IUCr Journals Google Scholar
Jha, K. K., Gruza, B., Kumar, P., Chodkiewicz, M. L. & Dominiak, P. M. (2020). Acta Cryst. B76, 296–306. Web of Science CSD CrossRef IUCr Journals Google Scholar
Klar, P. B., Krysiak, Y., Xu, H., Steciuk, G., Cho, J., Zou, X. & Palatinus, L. (2023). Nat. Chem. 15, 848–855. Web of Science CSD CrossRef CAS PubMed Google Scholar
Kleemiss, F., Dolomanov, O. V., Bodensteiner, M., Peyerimhoff, N., Midgley, L., Bourhis, L. J., Genoni, A., Malaspina, L. A., Jayatilaka, D., Spencer, J. L., White, F., Grundkötter-Stock, B., Steinhauer, S., Lentz, D., Puschmann, H. & Grabowsky, S. (2021). Chem. Sci. 12, 1675–1692. Web of Science CSD CrossRef CAS Google Scholar
Kolb, U., Gorelik, T., Kübel, C., Otten, M. T. & Hubert, D. (2007). Ultramicroscopy, 107, 507–513. Web of Science CrossRef PubMed CAS Google Scholar
Kuhs, W. F. & Lehmann, M. S. (1981). Nature, 294, 432–434. CrossRef CAS Web of Science Google Scholar
Kuhs, W. F. & Lehmann, M. S. (1983). J. Phys. Chem. 87, 4312–4313. CrossRef CAS Web of Science Google Scholar
Martynowycz, M. W. & Gonen, T. (2019). ChemRxiv, https://doi.org/10.26434/chemrxiv.8298641.v1. Google Scholar
Midgley, L., Bourhis, L. J., Dolomanov, O., Peyerimhoff, N. & Puschmann, H. (2019). arXiv:1911.08847. Google Scholar
Nassour, A., Domagala, S., Guillot, B., Leduc, T., Lecomte, C. & Jelsch, C. (2017). Acta Cryst. B73, 610–625. Web of Science CrossRef IUCr Journals Google Scholar
Palatinus, L., Corrêa, C. A., Steciuk, G., Jacob, D., Roussel, P., Boullay, P., Klementová, M., Gemmi, M., Kopeček, J., Domeneghetti, M. C., Cámara, F. & Petříček, V. (2015b). Acta Cryst. B71, 740–751. Web of Science CrossRef IUCr Journals Google Scholar
Palatinus, L., Petříček, V. & Corrêa, C. A. (2015a). Acta Cryst. A71, 235–244. Web of Science CrossRef IUCr Journals Google Scholar
Pauling, L. (1935). J. Am. Chem. Soc. 57, 2680–2684. CrossRef CAS Google Scholar
Peng, L. M. (1999). Micron, 30, 625–648. Web of Science CrossRef CAS Google Scholar
Peterson, S. W. & Levy, H. A. (1957). Acta Cryst. 10, 70–76. CrossRef ICSD CAS IUCr Journals Web of Science Google Scholar
Petříček, V., Palatinus, L., Plášil, J. & Dušek, M. (2023). Z. Kristallogr. Cryst. Mater. 238, 271–282. Google Scholar
Pichon-Pesme, V., Lecomte, C. & Lachekar, H. (1995). J. Phys. Chem. 99, 6242–6250. CrossRef CAS Web of Science Google Scholar
Ruth, P. N., Herbst-Irmer, R. & Stalke, D. (2022). IUCrJ, 9, 286–297. Web of Science CSD CrossRef CAS PubMed IUCr Journals Google Scholar
Saha, A., Nia, S. S. & Rodríguez, J. A. (2022). Chem. Rev. 122, 13883–13914. Web of Science CrossRef CAS PubMed Google Scholar
Sanjuan-Szklarz, W. F., Woińska, M., Domagała, S., Dominiak, P. M., Grabowsky, S., Jayatilaka, D., Gutmann, M. & Woźniak, K. (2020). IUCrJ, 7, 920–933. Web of Science CSD CrossRef CAS PubMed IUCr Journals Google Scholar
Volkov, A., Messerschmidt, M. & Coppens, P. (2007). Acta Cryst. D63, 160–170. Web of Science CrossRef CAS IUCr Journals Google Scholar
Wall, M. E. (2016). IUCrJ, 3, 237–246. Web of Science CrossRef CAS PubMed IUCr Journals Google Scholar
Woińska, M., Chodkiewicz, M. L. & Woźniak, K. (2021). Chem. Commun. 57, 3652–3655. Google Scholar
Woińska, M., Grabowsky, S., Dominiak, P. M., Woźniak, K. & Jayatilaka, D. (2016). Sci. Adv. 2, e1600192. Web of Science PubMed Google Scholar
Woińska, M., Jayatilaka, D., Spackman, M. A., Edwards, A. J., Dominiak, P. M., Woźniak, K., Nishibori, E., Sugimoto, K. & Grabowsky, S. (2014). Acta Cryst. A70, 483–498. Web of Science CSD CrossRef IUCr Journals Google Scholar
Wollan, E. O., Davidson, W. L. & Shull, C. G. (1949). Phys. Rev. 75, 1348. CrossRef Web of Science Google Scholar
Yonekura, K. & Maki-Yonekura, S. (2016). J. Appl. Cryst. 49, 1517–1523. Web of Science CrossRef CAS IUCr Journals Google Scholar
Zarychta, B., Pichon-Pesme, V., Guillot, B., Lecomte, C. & Jelsch, C. (2007). Acta Cryst. A63, 108–125. Web of Science CrossRef CAS IUCr Journals Google Scholar
Zheng, J.-C., Wu, L. & Zhu, Y. (2009). J. Appl. Cryst. 42, 1043–1053. Web of Science CrossRef CAS IUCr Journals Google Scholar
Zhong, S., Dadarlat, V. M., Glaeser, R. M., Head-Gordon, T. & Downing, K. H. (2002). Acta Cryst. A58, 162–170. Web of Science CrossRef CAS IUCr Journals Google Scholar
This is an open-access article distributed under the terms of the Creative Commons Attribution (CC-BY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.