Direct air capture of CO­₂ – topological analysis of the experimental electron density (QTAIM) of the highly insoluble carbonate salt of a 2,6-pyridine-bis­(iminoguanidine), (PyBIGH₂)(CO₃)(H₂O)₄

Chemical bonding and all intermolecular interactions in the highly insoluble carbonate salt of a 2,6-pyridine-bis(iminoguanidine), (PyBIGH₂)(CO₃)(H₂O)₄, recently employed in the direct air capture of CO₂ via crystallization, have been analyzed within the framework of the quantum theory of atoms in molecules (QTAIM) based on the experimental electron density derived from X-ray diffraction data obtained at 20 K. Accurate hydrogen positions were included based on an analogous neutron diffraction study at 100 K. Topological features of the covalent bonds demonstrate the presence of multiple bonds of various orders within the PyBIGH₂²⁺ cation. Strong hydrogen bonds define ribbons comprising carbonate anions and water molecules. These ribbons are linked to stacks of essentially planar dications via hydrogen bonds from the guanidinium moieties and an additional one to the pyridine nitro­gen. The linking hydrogen bonds are approximately perpendicular to the anion–water ribbons. The observation of these putative interactions provided motivation to characterize them by topological analysis of the total electron density. Thus, all hydrogen bonds have been characterized by the properties of their (3,−1) bond critical points. Weaker interactions between the PyBIGH₂²⁺ cations have similarly been characterized. Integrated atomic charges are also reported. A small amount of cocrystallized hydroxide ion (∼2%) was also detected in both the X-ray and neutron data, and included in the multipole model for the electron-density refinement. The small amount of additional H⁺ required for charge balance was not detected in either the X-ray or the neutron data. The results are discussed in the context of the unusually low aqueous solubility of (PyBIGH₂)(CO₃)(H₂O)₄ and its ability to sequester atmospheric CO₂.

1. Introduction

Given the strong connection between climate change and the greenhouse gases in the atmosphere, developing new methods to reduce their concentration in the air and alleviate global warming is of major importance. Besides curbing the use of fossil fuels, one strategy for limiting the increase in the atmospheric CO­₂ concentration is based on carbon capture and storage (CCS) (Lackner, 2003; Reiner, 2016). Significant progress has been made in employing CCS at point sources of CO₂ emission, such as coal- or gas-fired power plants. However, this approach does not address the problem of diffuse sources of CO­₂, such as households and transportation, which are responsible for approximately 50% of total greenhouse gas emissions. Most climate change mitigation scenarios aiming to limit global warming to 2°C or less now include implementation of negative emissions technologies (NETs) that target net reductions of the atmospheric CO₂ concentration, currently ∼408 p.p.m. One promising approach among various NETs under consideration is direct air capture (DAC), a process that removes CO₂ from the air by engineered chemical reactions (Keith, 2009; Lackner et al., 2012; Sanz-Pérez et al., 2016; Keith et al., 2018).

In a recent study, an aqueous solution of a 2,6-pyridine-bis(iminoguanidine), PyBIG [the displayed resonance form is based on the reported crystal structure (Seipp et al., 2017; Brethomé et al., 2018)], was found to efficiently absorb CO₂ from the atmosphere and convert it into the crystalline carbonate salt (PyBIGH₂)(CO₃)(H₂O)₄ (Seipp et al., 2017; Brethomé et al., 2018). The main driver for this reaction is the extremely low aqueous solubility of the carbonate salt (K_sp = 1.0 × 10⁻⁹), comparable to CaCO₃ (3.7–8.7 × 10⁻⁹), which pushes the overall equilibrium towards the carbonate formation despite the very low concentration of CO₂ in the air. After filtration, crystalline (PyBIGH₂)(CO­₃)(H₂O)₄ is mildly heated at 120°C to release the CO₂ (which can be sent to storage) and regenerate the PyBIG sorbent for reuse in another DAC cycle. Thus, this crystallization-based approach offers the prospect for energy-efficient DAC technology, provided the synthesis of PyBIG and the overall CO₂ capture process can be optimized and scaled up cost effectively.

The previously reported X-ray crystal structure of (PyBIGH₂)(CO­₃)(H₂O)₄ showed that the hydrated carbonate salt comprises an elaborate hydrogen-bonded network involving the carbonate anion, guanidinium cations, water molecules and the pyridine N atom (Seipp et al., 2017; Brethomé et al., 2018). However, the precise geometrical parameters and energetics of the hydrogen bonds and other intermolecular interactions present in this structure, which presumably play important roles in the unusually low aqueous solubility of (PyBIGH₂)(CO­₃)(H₂O)₄, have yet to be determined. To this end, as reported in this article, we have determined precisely all hydrogen positions from neutron diffraction data, which provides an accurate geometrical description of all the hydrogen bonds present. The interaction energies of these hydrogen bonds have been estimated from a topological analysis of the electron density, as determined from extremely accurate high-resolution X-ray diffraction data. At the same time, we have characterized all of the covalent bonds and the integrated atomic charges within the framework of the quantum theory of atoms in molecules (QTAIM) (Bader, 1994).

2. Experimental

2.1. Data collection and reduction

2.1.1. X-ray experiment

Colorless crystals of (PyBIGH₂)(CO₃)(H₂O)₄ were obtained by slow reaction of an aqueous solution of PyBIG with atmospheric CO₂. A single crystal (0.31 × 0.20 × 0.15 mm) was subsequently mounted with oil on top of a thin-walled glass capillary, and cooled to 20 K with an open-flow helium cryostat (Hardie et al., 1998; Kirschbaum et al., 1999). X-ray diffraction measurements were performed with a Rigaku diffractometer equipped with a Mo rotating anode generator operating at 50 kV and 300 mA (ULTRAX-18 Mo Kα, curved graphite monochromator) and using a RAPID-II cylindrical image-plate detector. To obtain highly redundant data, runs collecting 30 × 6° ω scans were performed at χ = 0°, ϕ = 0 and 180°, and at χ = 40°, ϕ = 0, 90, 180 and 270°. These runs were augmented by collecting an analogous set with ω offset by 3°. Thus, frames were overlapped by a half-frame width to improve scaling and allow for the omission of partial and overlapping reflections. An exposure time of 180 s per image was chosen to maximize he scattering power and avoid saturation of the strongest reflections.

The collected data were indexed and reflection positions predicted using the program HKL2000 (Otwinowski & Minor, 1997). Data were integrated with the program VIIPP, applying an image-plate flood-field correction, and with background and reflection profiles averaged over the whole data set, as described previously (Zhurova et al., 1999, 2008; Zhurov & Pinkerton, 2013). Partial and overlapped reflections were rejected during the integration. The effects of absorption (μ = 0.122 mm⁻¹) and thermal diffuse scattering at 20 K were considered to be negligible. Additional outliers were identified and removed manually through equivalence comparison to minimize outlier contamination. This is particularly important for removing errors from multiple scattering, and for identifying previously unidentified partial and overlapping reflections. This resulted in 1.41% of measured data (2730 out of 193 867 reflections) being additionally rejected prior to merging and scaling of the data in the space group with the program SORTAV (Blessing, 1995, 1987, 1997). Corrections of reflection intensities for λ/2 contamination were also made (Kirschbaum et al., 1997; Gianopoulos et al., 2017). Other experimental details are listed in Table 1.

Table 1
Experimental details

CIFs for both experiments are provided in the supporting information.

	X-ray	Neutron
Empirical formula	C9H15N9^2+·CO3^2−·4H2O	C9H15N9^2+·CO3^2−·4H2O
Crystal size (mm)	0.31 × 0.20 × 0.15	2.01 × 1.50 × 0.37
Crystal shape	Plate	Plate
Wavelength (Å)	0.71073	0.40–3.39 (TOF)
Crystal system	Triclinic	Triclinic
Temperature (K)	20	100
Space group
a (Å)	8.2090(2)	8.2420 (2)
b (Å)	8.5762 (2)	8.6011 (3)
c (Å)	13.8676 (4)	13.8821 (4)
α (°)	72.591 (2)	72.792 (3)
β (°)	78.815 (2)	78.998 (3)
γ (°)	71.0422 (17)	70.789 (2)
V (Å^3), Z	875.97 (4), 2	882.92 (5), 2
Density (g cm^−3)	1.447	1.436
μ (mm^−1)	0.122	0.1506 + 0.1027λ
(sinθ/λ)max (Å^−1)	1.30	2.45
Reflections integrated	189348	44971
Rint, average data multiplicity	0.027, 9.1	0.0963, 5.6
Completeness: sinθ/λ < 0.76 Å^−1, all data (%)	99.0/80.0	91.4
Independent reflections	25805	7955
Used reflections	18335 [I > 3σ(I)]	7955
Spherical refinement
R1[F, I > 2σ(I)], wR2(F^2), GOF	0.028, 0.082, 1.055	0.034,0.065, 1.128
Δρmin/max for X-rays (e Å^−3), for neutrons (fm Å^−3)	−0.33/0.79	−1.185/1.104
Multipole refinement
No. of parameters	1232
R1[F, I > 3σ(I)], wR2(F^2), GOF	0.018, 0.020, 1.115
Δρmin/max (e Å^−3), sinθ/λ < 1.3 Å^−1	−0.186, 0.272
Weighting scheme: a, b†	0.0038, 0.0038

†

2.1.2. Neutron experiment

Data from a single-crystal plate (2.01 × 1.50 × 0.37 mm) prepared as above were obtained using the TOPAZ single-crystal neutron time-of-flight (TOF) Laue diffractometer (Jogl et al., 2011; Schultz et al., 2014) at the Spallation Neutron Source (SNS) at Oak Ridge National Laboratory. The diffractometer is equipped with 24 detectors, each with an active area of 15 × 15 cm, arranged on a near-spherical detector array tank. The initial moderator-to-sample flight path is 18 m and the sample-to-detector distances vary in the range 39–46 cm. The total path length of 18.4 m and the SNS pulse rate of 60 Hz provides a wavelength bandwidth of 3.6 Å. The crystal was mounted on a MiTeGen loop using cyano­acrylate glue and cooled to 100 K for data collection. A total of 47 crystal orientations optimized with CrystalPlan software (Zikovsky et al., 2011) were used to ensure better than 95% coverage of a hemisphere of reciprocal space. Data were displayed, auto-indexed and integrated using the suite of algorithms in Mantid (Arnold et al., 2014). The raw Bragg intensities were obtained using the three-dimensional ellipsoidal Q-space integration method (Schultz et al., 1984). Data reduction including the neutron TOF spectrum, Lorentz, and detector efficiency corrections were carried out with the ANVRED3 program (Schultz et al., 1984). A Gaussian numerical absorption correction was applied with μ = 0.1506 + 0.1027λ mm⁻¹. The reduced data were saved in SHELX HKLF2 format, in which the neutron wavelength for each reflection was recorded separately.

2.2. Refinements

The crystal structure of PyBIG carbonate was reported previously (Seipp et al., 2017; Brethomé et al., 2018) and we have preserved the setting of the unit cell and the atom numbering used in that work. Based on our experimental neutron data, the crystal structure was re-refined within the SHELXTL program suite (Sheldrick, 2015) using the previously reported structure as the starting model. All atoms were refined using anisotropic thermal motion. An initial least-squares refinement based on the X-ray data was also carried out (SHELXTL). Anisotropic thermal motion was considered for all non-hydrogen atoms, and the hydrogen atoms were refined isotropically. From this starting point, a multipole refinement based on the Hansen–Coppens pseudo-atom formalism (Hansen & Coppens, 1978) [equation (1)], as implemented in the MoPro program package (Jelsch et al., 2005), using the Volkov and co-workers relativistic data bank (Volkov et al., 2006), was performed,

where ρ_c and ρ_v are spherical core and valence densities normalized to one electron, P_c and P_v are the core and spherical valence populations, respectively, R_l represents normalized Slater-type radial functions, y_lm are real angular spherical harmonics, and P_lm refers to the multipole population of the mth term of the lth order. The κ_s and κ_l terms are expansion–contraction coefficients for the spherical and multipolar valence densities, respectively.

All `heavy' atoms were refined to the hexadecapole level, while the hydrogen atoms were refined up to dipoles plus the bond-directed quadrupole, with C—H, N—H and O—H distances constrained to the values obtained from the neutron study. In the initial stages of refinement, chemical constraints for similar atoms were applied; however, these constraints were gradually released, and the final model was refined unconstrained (24 refined multipole populations for each `heavy' atom and 4 refined multipole populations for hydrogen atoms), with the exception of κ parameters (see below). The molecular electroneutrality requirement was applied throughout for the total structure. This allowed for charge transfer among the charged species rather than constraining their formal charge. The expansion–contraction parameters κ_s and κ_l for the non-hydrogen atoms were refined in ten groups according to their chemical equivalence, while κ_s and κ_l for hydrogen atoms were set to 1.2. The final description of the anisotropic thermal motion for the hydrogen atoms was obtained from SHADE-3.1 (Madsen, 2006).

The residual map calculated after the multipole refinement still had one unidentified peak significantly above background. Examination of the neutron residual showed the same small feature along with a negative neighbor (Figs. S8 and S9). These features were identified as a small number of cocrystallized hydroxide ions. The refined occupancies were 0.017 (2) from the X-ray data (IAM model) and 0.022 (2) from neutrons. No evidence was found in the neutron data for the H⁺ required for charge balance, hence we assume that it is disordered over the available oxygen and nitro­gen sites. The final multipole refinement then included a variable occupancy for the contribution from a spherically modeled oxygen atom.

Topological analysis of the total electron density was carried out with the program packages MoPro (Jelsch et al., 2005), XDPROP (Volkov et al., 2006) and WinXPRO (Stash & Tsirelson, 2002, 2005).

2.3. Evaluation of X-ray data quality

Analysis of statistical measures of data and multipole model quality have been deposited; all suggest excellent data and an excellent model. Averaged ratios (in 0.05 Å⁻¹ bins) of observed and calculated structure factors (Fig. S1) as well as the normal probability plot (Fig. S2) indicate good model fitting for the whole sinθ/λ range. The residual electron density maps (Fig. S3) are low (ρ_min/max = −0.186/0.272 e Å⁻³, calculated for the complete data set) and featureless as confirmed by a fractal dimension plot (Fig. S4). The total electron density was non-negative everywhere.

3. Results and discussion

3.1. Structure

As reported previously (Seipp et al., 2017; Brethomé et al., 2018), the structure is made up of essentially planar PyBIGH₂²⁺ dications, carbonate anions and four water molecules. The asymmetric unit and the atom labeling are shown in Fig. 1. The anions and water molecules form a ribbon (Fig. 2), and based on distance criteria, we may already propose that this structure is strongly hydrogen bonded, as shown in the figure. The cations form stacks (Fig. 3) that are approximately perpendicular to the plane of the anion–water ribbons.

Figure 1 The asymmetric unit of PyBIG carbonate tetrahydrate, as determined from 100 K neutron diffraction data, showing the atom numbering. Displacement ellipsoids are at the 50% probability level (Macrae et al., 2008).

Figure 2 Potential hydrogen bonds in the anion–water ribbons. The figure is based on the neutron structure, with displacement ellipsoids at the 20% probability level and distances in Å (Macrae et al., 2008). Color scheme – oxygen, red; carbon, dark gray; hydrogen, light gray.

Figure 3 The stacking of PyBIGH2^2+ cations. The figure is based on the neutron structure, with displacement ellipsoids drawn at the 20% probability level (Macrae et al., 2008). Color scheme: nitrogen, blue; carbon, dark gray; hydrogen, light gray.

The anion–water ribbon is canted at an angle of ∼23.6° above and below the ac plane, and extends about 1.91 Å above and below the plane (Fig. S5). The ribbon has a maximal thickness of ∼1.49 Å (on the basis of heavy atoms), while the H52 atom is oriented nearly perpendicular (81.1°) to the mean plane of the ribbon. While each ribbon extends infinitely along the a axis, the width of each ribbon is ∼11.82 Å, about 2 Å shorter than the length of the c axis. The distance between nearest neighbors on different adjacent ribbons is approximately 4.2 Å and gives rise to a channel between neighboring ribbons. In this context, it is unsurprising that the sites of the partially occupied hydroxide ions fall in this cavity and are suggestive of a stabilizing interaction between neighboring ribbons (Fig. S5). On the basis of the neutron diffraction results, the nearest hydroxide HO⋯H distance is ∼1.90 Å, while the nearest HO—H⋯OH distance is ∼1.73 Å. The cations lie roughly above and below the ac plane containing the extended network of ribbons. When viewed along the [101] direction it becomes clear that the nearest cations are all hydrogen-bond donors to the water–anion ribbons and form linear arrays along the [101] vector, nearly in the (¼ 0 ¼) plane. Slightly further from the ribbons are cation arrays (along [101]) wherein the hydrogen-bond-accepting pyridine N5 atom is oriented towards the ribbon (Figs. S6 and S7).

The distances shown in Fig. 4 suggest strong hydrogen bonds between the guanidinium hydrogen atoms and a variety of oxygen atoms, as well as a water hydrogen bonded to the pyridine nitro­gen. Note that the cations in Fig. 4 have been truncated to emphasize the possible hydrogen-bond interactions. It is also clear from Fig. 3 that there is potential for additional interactions between the π-systems of neighboring cation sheets as they are only separated by ∼3.2 Å. The observation of these putative noncovalent interactions provided much of the motivation to characterize them by topological analysis of the total electron density.

Figure 4 Potential hydrogen bonds between the anion–water ribbon and neighboring cations. All cations have been truncated to enhance the visibility of the hydrogen bonds. The figure is based on the neutron structure, with displacement ellipsoids drawn at the 20% proabability level and distances in Å (Macrae et al., 2008). Color scheme: oxygen, red; nitrogen, blue; carbon, dark gray; hydrogen, light gray.

4. Electron density

4.1. Atomic charges

The integrated charges of the atomic basins delimited by the zero-flux surfaces and their volumes are listed in Table 2. For convenience, the charges are also reported in Fig. 5. The accuracy of the integration for each atom was confirmed by a small value of the integrated Laplacian (Lagrangian). The atomic charges sum to zero as required; however, the total charge of the cation and anion differ from the formal value of 2.0, indicating significant charge transfer. Concomitantly, all water molecules are close to neutral, and the disordered OH group contributes a small amount of negative charge. The sum of the atomic volumes is close to the unit-cell volume per asymmetric unit with an error of ∼0.2%. All oxygen atoms have significant negative charges of similar magnitudes, whether in the anion or in the water molecules. The nitro­gens are all strongly negative, and may be differentiated according to their type (NH₂ < N_pyridine < NH < N_imino). The carbon atoms have a wide range of mainly positive charges that strongly correlate with their environment (C_carbonate > C_guanidine > C_imino > C_pyridine – the latter being slightly negative). As expected, the hydrogen atoms are all strongly positive and again may be grouped by type (H₂O > NH > CH).

Table 2
Integrated atomic charges (q) and volumes (Ω)

Atom	q (e)	Ω (Å^3)	Atom	q (e)	Ω (Å^3)
PyBIG cation			Carbonate
N1	−1.14	17.06	C10	2.02	4.43
H12	0.54	2.33	O1	−1.17	17.90
H11	0.52	2.08	O2	−1.16	17.97
C1	1.33	5.32	O3	−1.14	17.08
N2	−1.19	19.93	Total	−1.45	57.38
H22	0.51	2.25
H21	0.54	2.15	Water
C2	0.38	10.28	O4	−1.12	18.46
H2	0.16	7.04	H42	0.64	1.33
N3	−0.85	12.51	H41	0.60	1.70
H3	0.57	1.83	Total	0.12	21.49
C3	0.31	7.98
N4	−0.60	14.72	Water
C4	−0.11	11.65	O5	−1.21	20.95
H4	0.16	6.12	H52	0.60	1.94
C5	−0.07	11.61	H51	0.58	2.16
H5	0.14	6.89	Total	−0.04	25.05
N5	−0.90	13.88
C6	−0.09	11.98	Water
H6	0.15	6.16	O6	−1.21	18.00
N6	−0.65	13.02	H62	0.65	1.57
C7	0.31	8.03	H61	0.60	1.75
N7	−0.79	12.29	Total	0.04	21.32
H7	0.54	1.87
C8	0.39	10.32	Water
H8	0.16	7.01	O7	−1.21	19.83
N8	−1.16	16.23	H72	0.63	1.76
H82	0.50	2.94	H71	0.61	1.74
H81	0.53	2.23	Total	0.03	23.33
C9	1.27	5.09
N9	−1.22	17.51	Hydroxyl
H92	0.54	3.26	O8	−0.01	14.82
H91	0.53	1.98
Total	1.32	275.55	Total/total	0.01	438.94
			Cell volume/2		437.99

Figure 5 Integrated atomic charges (black) and topological bond orders (red).

The deformation density in the plane of the dication is mapped in Fig. 6(a) and clearly shows a significant concentration of electron density in all of the covalent bonds, as well as the presence of lone pairs on the imino and pyridine nitro­gen atoms. The covalent bonding density is also well represented for the anion and for the water molecules. Again, the expected lone-pair regions on the oxygen atoms are also well defined.

Figure 6 Deformation density (a) in the plane of the dication and (b) in the plane of the anion. Blue contours are positive density and red ones are negative. The contour level is 0.10 e Å^−3. (c) The anion deformation density iso-surface at 0.15 e Å^−3.

More complete information on the nature of the bonding may be obtained from a topological analysis of the total electron density (Fig. 7 and Table 3). All (3,−1) critical points for the covalent bonds in the dication, the carbonate anion, and selected water molecules are indicated by yellow spheres in Figs. 7(a) and 7(b). Their characteristics are tabulated in Table 3. All covalent bonds have significant electron density at the critical point, with negative values of the Laplacian as required. Of particular interest is the extent of electron delocalization (π bonding) in the essentially planar dication. In general, all bonds in the molecular skeleton are short, with significant ellipticities at the critical points indicating important π-character. Complementary information on the nature of these bonds may be obtained from the topological bond orders as defined by n_topo = a + bλ₃ + c(λ₁ + λ₂) + dρ_CP, where ρ is the electron density at the critical point, λ_1,2,3 are obtained from the Hessian matrix, and the coefficients (a, b, c, d) were taken from the literature (Howard & Lamarche, 2003; Tsirelson et al., 2006, 2007; Bartashevich et al., 2011). The bond orders for the skeleton of the cation, which range from 1.078 to 1.382 and show close to twofold molecular symmetry in their value, further indicate the delocalized π-character of the C—C and C—N bonds (Table 3 and Fig. 5). The strongest bonds are those involving the imino atoms N4 and N6, whereas the weakest are the substituted guanidinium C—N bonds (C1—N3 and C9—N7) and the substituents of the pyridine ring (C2—C3 and C7—C8). The bond orders of the N—H bonds are significantly lower than those of C—H, corresponding to the higher positive charges on the H(N) atoms compared with H(C).

Table 3
Characteristics of covalent bond critical points in PyBIG carbonate tetrahydrate

ε = λ₁/λ₂ – 1; n_topo = a + bλ₃ + c(λ₁ + λ₂) + dρ_CP (Howard & Lamarche, 2003; Tsirelson et al., 2006, 2007; Bartashevich et al., 2011); complete tables of bond critical-point properties have been included in the supporting information.

Atom 1	Atom 2	ρ(r) (e Å^−3)	∇^2ρ(r) (e Å^−5)	Rij (Å)	ε	ntopo
PyBIG
C1	N1	2.492	−26.400	1.325	0.167	1.260
C1	N2	2.455	−30.350	1.327	0.219	1.157
C1	N3	2.372	−28.050	1.352	0.173	1.104
C2	C3	1.909	−16.050	1.468	0.102	1.084
C2	H2	1.832	−19.780	1.095	0.046	0.894
C2	N4	2.637	−31.840	1.286	0.219	1.342
C3	C4	2.139	−20.070	1.398	0.192	1.265
C3	N5	2.349	−22.860	1.347	0.113	1.169
C4	C5	2.162	−20.020	1.391	0.138	1.314
C4	H4	1.878	−20.420	1.088	0.035	0.928
C5	C6	2.168	−20.320	1.388	0.163	1.304
C5	H5	1.870	−20.250	1.086	0.022	0.927
C6	C7	2.111	−19.480	1.402	0.184	1.258
C6	H6	1.892	−20.850	1.085	0.028	0.921
C7	C8	1.899	−15.870	1.470	0.117	1.078
C7	N5	2.374	−23.750	1.345	0.112	1.183
C8	H8	1.836	−19.560	1.096	0.055	0.912
C8	N6	2.637	−31.870	1.284	0.192	1.342
C9	N7	2.352	−25.500	1.358	0.186	1.121
C9	N8	2.486	−28.670	1.322	0.209	1.215
C9	N9	2.453	−27.400	1.326	0.188	1.212
N1	H11	2.134	−31.950	1.005	0.036	0.635
N1	H12	2.129	−32.370	1.037	0.042	0.689
N2	H21	2.091	−31.280	1.029	0.037	0.659
N2	H22	2.155	−30.160	1.013	0.045	0.730
N3	H3	2.031	−31.380	1.057	0.043	0.626
N3	N4	2.375	−7.047	1.362	0.084	1.382
N6	N7	2.406	−7.880	1.357	0.074	1.367
N7	H7	2.039	−30.490	1.041	0.047	0.622
N8	H81	2.092	−30.960	1.030	0.025	0.657
N8	H82	2.203	−32.860	0.998	0.037	0.687
N9	H91	2.124	−31.530	1.015	0.043	0.663
N9	H92	2.079	−30.100	1.045	0.041	0.695
Carbonate
C10	O1	2.506	−30.380	1.289	0.111	1.309
C10	O2	2.476	−30.030	1.289	0.113	1.289
C10	O3	2.467	−29.830	1.292	0.128	1.302
Waters
O4	H41	2.274	−36.720	0.980	0.046	0.595
O4	H42	2.254	−41.790	0.971	0.024	0.550
O5	H51	2.223	−33.490	0.967	0.016	0.449
O5	H52	2.248	−37.030	0.973	0.009	0.528
O6	H61	2.205	−36.630	0.981	0.022	0.505
O6	H62	2.210	−41.440	0.980	0.005	0.526
O7	H71	2.243	−38.560	0.973	0.005	0.516
O7	H72	2.159	−37.750	0.984	0.020	0.472

Figure 7 Bond paths and critical points (yellow spheres) for (a) all covalent bonds and intermolecular hydrogen bonds in the carbonate–water ribbons, for (b) all covalent bonds in the dication and cation–cation stacking interactions, and (c) selected bond paths and critical points linking the cationic stacks and anionic ribbons. Color scheme: C, black; N, blue; O, red; H, green.

Although Fig. 6(b) implies well resolved lone pairs on all the oxygen atoms of the carbonate anion, these are actually cuts through a doughnut-like charge distribution as shown by the iso-surface plot in Fig. 6(c). As expected, the carbonate anion is strongly covalently bonded, and the charge distribution is highly polarized, the oxygen atoms carrying high negative charges and the carbon atom being strongly positive.

4.2. Closed-shell interactions

The most important noncovalent interactions in this structure are the hydrogen bonds, which range in strength from modest to strong. All bond paths have been identified from the topology of the electron density, and their (3,−1) bond critical points characterized (Table 4 and Fig. 7). Every hydrogen atom in the structure bonds to a neighboring oxygen or to the pyridine nitrogen (N5) except H22, which only bonds to the partially occupied OH group. All carbonate oxygen atoms accept three hydrogen bonds and all water oxygen atoms accept two. The only nitro­gen atom that accepts a hydrogen bond is the pyridine N5 atom.

Table 4
Characteristics of bond critical points for closed shell intermolecular interactions

Complete tables of bond critical point properties are provided in the supporting information.

Atom 1	Atom 2	Symmetry†	ρ(r) (e Å^−3)	∇^2ρ(r) (e Å^−5)	Rij (Å)‡	De (kJ mol^−1)	Acceptor group/donor group
O1	H7	46401	0.358	2.093	1.698	66.0	Carbonate/imino
O2	H11	55501	0.354	2.208	1.725	65.5	Carbonate/guan
O3	H91	46401	0.317	2.527	1.725	57.5	Carbonate/guan
O6	H3	65501	0.299	2.466	1.724	53.1	Water/imino
O3	H61	66502	0.296	2.590	1.727	52.8	Carbonate/water
O2	H81	45401	0.279	2.700	1.751	49.6	Carbonate/guan
O1	H71	56502	0.266	2.170	1.790	44.2	Carbonate/water
O7	H62	55501	0.258	2.544	1.766	44.2	Water/water
O2	H41	56502	0.251	2.074	1.798	40.6	Carbonate/water
O5	H72	65602	0.237	2.561	1.797	40.0	Water/water
O3	H21	55501	0.244	2.180	1.830	39.7	Carbonate/guan
O6	H42	55501	0.200	2.171	1.868	31.3	Water/water
O1	H51	56502	0.187	2.292	1.897	29.4	Carbonate/water
N5	H52	55501	0.201	1.608	1.963	28.8	Pyridine/water
O4	H12	65501	0.154	1.475	2.003	20.5	Water/guan
O7	H82	66602	0.108	1.718	2.115	15.5	Water/guan
O4	H92	55401	0.092	1.708	2.129	13.6	Water/guan
O7	H6	66602	0.055	0.527	2.561	4.9	Water/pyridine
O5	H5	54501	0.047	0.647	2.647	4.9	Water/pyridine
N1	N3	56502	0.051	0.562	3.256	4.8	Guan/imino
C3	C9	66602	0.057	0.463	3.280	4.7	Pyridine/guan
C7	N6	66602	0.054	0.498	3.267	4.7	Pyridine/imino
N3	C9	56602	0.049	0.458	3.398	4.1	Imino/guan
O5	C5	56602	0.042	0.454	3.216	3.6	Water/pyridine
O7	H4	64501	0.034	0.511	2.807	3.4	Water/pyridine
O6	H4	64501	0.028	0.513	2.807	3.2	Water/pyridine
O1	H8	46401	0.030	0.453	2.796	3.0	Carbonate/imino
O4	N2	64501	0.026	0.477	3.366	2.9	Water/guan
O4	N2	66502	0.021	0.277	3.550	1.8	Water/guan
O6	N8	66602	0.019	0.282	3.535	1.7	Wate/guan
O7	H8	65602	0.020	0.235	3.138	1.5	Water/imino
O1	C5	57502	0.015	0.166	3.867	1.1	Carbonate/pyridine
O1	N2	47502	0.008	0.087	4.214	0.5	Carbonate/guan

†ORTEP symmetry codes for Atom 2.
‡The shorter values for the intermolecular distances compared with those reported for the neutron study are caused by lattice contraction at lower temperature.

We have estimated the dissociation energy of each hydrogen bond based on the topological analysis and assuming the validity of the relationship by Espinosa et al. (1998, 1999). The dissociation energies of the most important hydrogen bonds range from modest (∼14 kJ mol⁻¹) to strong (∼66 kJ mol⁻¹). From the criterion (Espinosa et al., 2002), although these are still closed-shell interactions, the stronger ones are well within the so-called `transit' region, .

As shown in the figures, we may divide the hydrogen bonds into two sets, one which defines the carbonate–water ribbons and another roughly perpendicular to the first linking the anionic ribbons and cationic stacks. We recognize that this is an artificial classification as the two types have similar energies; however, we believe that it provides some insight into the design criteria for new cations to provide such insoluble materials.

As expected, the strongest hydrogen bonds involve the carbonate anion, which accepts a total of nine hydrogen bonds, five from the guanidinium groups and four from the water molecules (Fig. 8a), as well as three much weaker interactions. The correlation between the hydrogen-bond energies and the observed H⋯O contact distances (Fig. 8b) may be fitted to an exponential curve as anticipated from the derivation of the Espinosa relationship. The estimated carbonate `binding' energy from all hydrogen bonds and the three weaker interactions amount to −449.9 kJ mol<sup>−1</sup>. Notably, a large fraction of the carbonate `binding' energy (−167 kJ mol⁻¹, 37.1%) comes from hydrogen bonding to water. Clearly, the water molecules of hydration play an important role in the stability, and thereby the low aqueous solubility of (PyBIGH₂)(CO­₃)(H₂O)₄ crystals, by providing a total of −441.9 kJ mol⁻¹ in hydrogen-bonding energy. The strong hydrogen bonding of carbonate and water in these crystals is needed to partially compensate for the large free energy of dehydration of the anion (1315 kJ mol⁻¹) (Marcus, 1991) involved in the crystallization of PyBIGH₂(CO­₃)(H₂O)₄. Additionally, the lattice energy, consisting of electrostatic as well as other interactions (vide infra), must also contribute to the low aqueous solubility of these crystals.

Figure 8 Hydrogen bonding involving the carbonate anion. (a) Carbonate `binding' by five guanidinium and four water hydrogen bonds, with estimated dissociation energies in kJ mol^−1. (b) Observed correlation between hydrogen-bond energies and H⋯O contact distances; red points are from our experiment and the green line is from the work by Espinosa et al. (1998).

A complete topological analysis of the electron density reveals a number of additional bond paths, suggesting much weaker interactions. Although not rigorously justified, extrapolating the Espinosa et al. (1998, 1999) relationship suggests that dissociation energies for these additional interactions are all <5 kJ mol⁻¹ (Table 4). By definition, a bond path must begin and end at a nucleus; however, due to the close face-to-face proximity of the planar cations, many of these interactions may be better described as π–π interactions.

4.3. Lattice energy and electrostatic interactions

It is well known that high lattice energies tend to lower the solubility of crystalline compounds:

Although not the only contributor to the lattice energy [E_int, equation (2)], electrostatic interactions (E_es) tend to dominate this quantity in ionic crystals (Coppens, 1997; Volkov et al., 2006). Determination of the exchange–repulsion and dispersion to the total interaction energy are method dependent and, hence, unreliable. However, the electrostatic term may be obtained from the multipole expansion of the electron density using the methodology proposed by Volkov et al. (2004). Thus, we have determined the electrostatic crystal binding energy for (PyBIGH₂)(CO₃)(H₂O)₄ to be −583 kJ mol⁻¹. Although this may seem modest for an ionic compound, we have noted significant charge transfer between the cation and the anion, and the cationic charge is highly delocalized.

5. Conclusions

This article reports both a neutron diffraction study and a high-resolution X-ray diffraction study of a highly insoluble carbonate salt formed by crystallization of a guanidine compound with atmospheric CO₂. Both diffraction methods confirm the presence of a small amount of cocrystallized hydroxide ion. The accurately determined topological properties of the electron density characterize the delocalized nature of the bonding in the planar cation, as well as two well developed strong hydrogen-bonding schemes, one defining an anion–water ribbon, and the other essentially orthogonal to the anion–water `plane', linking the anionic ribbons to the cationic stacks. The carbonate anions are strongly hydrogen bonded in these crystals, which likely contributes to the extremely low aqueous solubility of this salt. The water molecules of hydration, hydrogen bonded to the carbonate anions and the guanidinium cations, also play an important role in the stability of these crystals. While the intermolecular interactions are dominated by strong hydrogen bonds, a number of supplementary weaker interactions have been characterized. Although their bond paths have been identified by their nuclear attractors, many may better be characterized as π–π interactions. The electrostatic contribution to the lattice energy is relatively modest due to both charge transfer and charge delocalization.

Supporting information

Computing details

Data collection: CrystalClear, (2005) for xray; SNS PyDas for neutron. Cell refinement: Otwinowski & Minor (1997) for xray; Mantidplot for neutron. Data reduction: Zhurova et al.,(1999), Zhurov et al.,(2008), Blessing, (1995,1997) for xray; Mantidplot and anvred3 for neutron. Program(s) used to refine structure: MoPro (J. Appl. Cryst. 2005, 38, 38-54) for xray; SHELXL2014/7 (Sheldrick, 2014) for neutron. Molecular graphics: 'Stash, (2002,2005), MoPro (J. Appl. Cryst. 2005, 38, 38-54)' for xray.

(xray)

Crystal data

C9H15N9^2+·CO3^2^−·4(H2O)·0.0115O	Z = 2
Mr = 381.52	F(000) = 404.18
Triclinic, P1	Dx = 1.447 Mg m^−3
Hall symbol: -P 1	Mo Kα radiation, λ = 0.71073 Å
a = 8.2090 (2) Å	Cell parameters from 189402 reflections
b = 8.5762 (2) Å	θ = 2.6–67.5°
c = 13.8676 (4) Å	µ = 0.12 mm^−1
α = 72.591 (2)°	T = 20 K
β = 78.815 (2)°	Block, colorless
γ = 71.0422 (17)°	0.31 × 0.20 × 0.15 mm
V = 875.97 (4) Å^3

Data collection

Rigaku R-axis Rapid diffractometer	18335 reflections with > 3.0σ(I)
Radiation source: rotating-anode X-ray tube	Rint = 0.023
Graphite monochromator	θmax = 67.5°, θmin = 2.6°
/w scans	h = −20→21
189348 measured reflections	k = −20→22
25805 independent reflections	l = 0→33

Refinement

Refinement on F^2	Secondary atom site location: difference Fourier map
Least-squares matrix: full	Hydrogen site location: difference Fourier map
R[F^2 > 2σ(F^2)] = 0.018	H atoms treated by a mixture of independent and constrained refinement
wR(F^2) = 0.020	w = 1/[σ^2(Fo^2) + (0.004P)^2 + 0.004P] where P = (Fo^2 + 2Fc^2)/3
S = 1.12	(Δ/σ)max = 0.002
18335 reflections	Δρmax = 0.27 e Å^−3
1232 parameters	Δρmin = −0.19 e Å^−3
197 restraints	Extinction correction: Becker-Coppens type 1 Lorentzian isotropic
Primary atom site location: structure-invariant direct methods	Extinction coefficient: 0.303 (12)

Special details

Refinement. Refinement of F² against reflections. The threshold expression of F² > 2sigma(F²) is used for calculating R-factors(gt) and is not relevant to the choice of reflections for refinement. R-factors based on F² are statistically about twice as large as those based on F, and R-factors based on ALL data will be even larger.

Fractional atomic coordinates and isotropic or equivalent isotropic displacement parameters (Ų)

	x	y	z	Uiso*/Ueq	Occ. (<1)
C10	−0.393854 (16)	0.899037 (17)	−0.191931 (13)	0.005934 (11)
O1	−0.469282 (18)	0.991765 (17)	−0.271783 (14)	0.007892 (11)
O2	−0.426140 (19)	0.757127 (18)	−0.142588 (14)	0.008795 (12)
O3	−0.285570 (17)	0.948565 (17)	−0.160417 (13)	0.007524 (10)
O4	0.75846 (2)	0.20366 (2)	0.043522 (17)	0.012323 (16)
H42	0.8342 (4)	0.1733 (9)	0.0961 (2)	0.0210 (2)
H41	0.6446 (3)	0.2041 (8)	0.0817 (4)	0.0224 (2)
O5	0.26382 (2)	0.09435 (2)	0.450433 (16)	0.010262 (13)
H52	0.2445 (9)	0.21635 (17)	0.4389 (5)	0.0193 (2)
H51	0.3375 (6)	0.0653 (8)	0.3912 (2)	0.0210 (3)
O6	0.94352 (2)	0.184288 (19)	0.199156 (16)	0.009340 (12)
H62	0.8772 (6)	0.1364 (7)	0.2614 (2)	0.0213 (2)
H61	1.0662 (2)	0.1184 (7)	0.1974 (6)	0.0225 (2)
O7	0.699718 (19)	0.098775 (19)	0.350384 (16)	0.009068 (12)
H72	0.7081 (10)	0.0295 (7)	0.42101 (17)	0.0208 (3)
H71	0.6210 (4)	0.0661 (5)	0.3213 (3)	0.0213 (2)
N1	−0.231140 (18)	0.562701 (19)	0.012676 (15)	0.008222 (12)
H12	−0.2221 (8)	0.43377 (19)	0.0365 (5)	0.0200 (2)
H11	−0.2973 (6)	0.6390 (7)	−0.0461 (2)	0.0191 (3)
C1	−0.176706 (16)	0.633633 (18)	0.068282 (13)	0.006348 (11)
N2	−0.19616 (2)	0.80022 (2)	0.041425 (16)	0.009888 (13)
H22	−0.1502 (5)	0.8524 (4)	0.0831 (3)	0.0204 (3)
H21	−0.2424 (7)	0.8700 (7)	−0.0272 (2)	0.0200 (3)
C2	0.036458 (17)	0.505275 (18)	0.284909 (13)	0.006611 (11)
H2	0.0605 (9)	0.36685 (13)	0.3040 (5)	0.0200 (2)
N3	−0.100870 (17)	0.533197 (18)	0.152569 (14)	0.006926 (11)
H3	−0.0836 (9)	0.40022 (16)	0.1731 (6)	0.0188 (2)
C3	0.091974 (16)	0.576300 (17)	0.352824 (12)	0.005675 (11)
N4	−0.042925 (17)	0.606717 (19)	0.208424 (14)	0.006717 (11)
C4	0.043718 (17)	0.752022 (18)	0.342394 (13)	0.007241 (12)
H4	−0.0362 (5)	0.8382 (6)	0.2832 (3)	0.0208 (3)
C5	0.097301 (18)	0.812333 (18)	0.409876 (14)	0.007529 (12)
H5	0.0605 (9)	0.94795 (17)	0.4049 (5)	0.0205 (2)
N5	0.187435 (16)	0.462028 (18)	0.425870 (13)	0.005677 (11)
C6	0.195735 (17)	0.696026 (17)	0.485430 (13)	0.006967 (11)
H6	0.2393 (5)	0.7373 (5)	0.5402 (2)	0.0203 (3)
N6	0.388563 (17)	0.446817 (18)	0.634204 (14)	0.006180 (11)
C7	0.238053 (16)	0.521773 (17)	0.490883 (12)	0.005520 (11)
N7	0.482825 (17)	0.330366 (17)	0.708170 (13)	0.006342 (10)
H7	0.5097 (9)	0.20005 (17)	0.7163 (5)	0.0190 (2)
C8	0.341477 (16)	0.393592 (17)	0.569772 (13)	0.006325 (11)
H8	0.3721 (9)	0.2584 (2)	0.5719 (5)	0.0204 (2)
N8	0.486444 (18)	0.560531 (18)	0.761509 (14)	0.007676 (11)
H82	0.4192 (6)	0.6373 (7)	0.7040 (2)	0.0186 (3)
H81	0.5296 (5)	0.6166 (4)	0.8038 (3)	0.0193 (3)
C9	0.532560 (16)	0.393882 (17)	0.772754 (12)	0.005579 (11)
N9	0.624964 (18)	0.282319 (18)	0.845606 (14)	0.007841 (12)
H92	0.6658 (6)	0.3229 (5)	0.8984 (2)	0.0193 (3)
H91	0.6615 (9)	0.1565 (2)	0.8475 (5)	0.0188 (2)
O8	0.9737 (11)	−0.0158 (10)	0.0804 (7)	0.0114 (8)	0.01153 (5)

Atomic displacement parameters (Ų)

	U^11	U^22	U^33	U^12	U^13	U^23
C10	0.00678 (4)	0.00536 (4)	0.00608 (6)	−0.00147 (3)	−0.00250 (4)	−0.00119 (4)
O1	0.01019 (4)	0.00694 (4)	0.00683 (7)	−0.00200 (3)	−0.00441 (4)	−0.00046 (4)
O2	0.01219 (4)	0.00653 (4)	0.00893 (7)	−0.00422 (3)	−0.00499 (4)	0.00032 (4)
O3	0.00785 (4)	0.00727 (4)	0.00861 (7)	−0.00257 (3)	−0.00360 (4)	−0.00154 (4)
O4	0.01226 (5)	0.01559 (6)	0.01024 (8)	−0.00495 (4)	−0.00291 (5)	−0.00275 (6)
H42	0.0210 (8)	0.0266 (10)	0.0168 (7)	−0.0050 (7)	−0.0104 (5)	−0.0042 (6)
H41	0.0162 (5)	0.0272 (10)	0.0208 (10)	−0.0062 (7)	0.0033 (6)	−0.0050 (8)
O5	0.01303 (5)	0.00841 (4)	0.00873 (8)	−0.00272 (4)	−0.00049 (5)	−0.00209 (5)
H52	0.0267 (10)	0.0079 (3)	0.0232 (10)	−0.0038 (8)	−0.0029 (9)	−0.0050 (8)
H51	0.0225 (10)	0.0223 (10)	0.0146 (7)	−0.0038 (8)	0.0056 (7)	−0.0072 (7)
O6	0.00844 (4)	0.00845 (4)	0.01006 (7)	−0.00148 (3)	−0.00184 (4)	−0.00129 (4)
H62	0.0209 (7)	0.0246 (10)	0.0161 (8)	−0.0111 (8)	0.0053 (6)	−0.0016 (7)
H61	0.0091 (4)	0.0255 (10)	0.0274 (10)	0.0024 (6)	−0.0022 (7)	−0.0069 (9)
O7	0.00974 (4)	0.00885 (4)	0.00883 (7)	−0.00335 (3)	−0.00255 (4)	−0.00081 (4)
H72	0.0282 (10)	0.0194 (10)	0.0114 (5)	−0.0083 (9)	−0.0043 (9)	0.0039 (7)
H71	0.0215 (7)	0.0248 (9)	0.0261 (9)	−0.0131 (6)	−0.0083 (6)	−0.0078 (7)
N1	0.01030 (4)	0.00823 (4)	0.00727 (7)	−0.00257 (3)	−0.00422 (4)	−0.00163 (4)
H12	0.0292 (10)	0.0098 (4)	0.0228 (10)	−0.0079 (8)	−0.0088 (9)	−0.0005 (9)
H11	0.0241 (9)	0.0178 (10)	0.0148 (8)	−0.0042 (8)	−0.0135 (6)	0.0020 (7)
C1	0.00705 (4)	0.00686 (4)	0.00539 (6)	−0.00169 (3)	−0.00237 (4)	−0.00118 (4)
N2	0.01420 (5)	0.00698 (4)	0.00920 (7)	−0.00321 (4)	−0.00580 (5)	−0.00003 (4)
H22	0.0280 (10)	0.0200 (9)	0.0202 (9)	−0.0097 (7)	−0.0107 (7)	−0.0073 (6)
H21	0.0279 (10)	0.0194 (10)	0.0121 (5)	−0.0054 (8)	−0.0122 (7)	0.0016 (7)
C2	0.00796 (4)	0.00648 (4)	0.00583 (6)	−0.00156 (3)	−0.00307 (4)	−0.00127 (4)
H2	0.0308 (10)	0.0074 (3)	0.0230 (10)	−0.0042 (9)	−0.0099 (8)	−0.0022 (9)
N3	0.00898 (4)	0.00650 (4)	0.00587 (7)	−0.00189 (3)	−0.00373 (4)	−0.00090 (4)
H3	0.0293 (10)	0.0089 (3)	0.0206 (10)	−0.0070 (8)	−0.0089 (9)	−0.0015 (9)
C3	0.00633 (4)	0.00570 (4)	0.00520 (6)	−0.00129 (3)	−0.00219 (4)	−0.00116 (4)
N4	0.00795 (4)	0.00704 (4)	0.00565 (7)	−0.00189 (3)	−0.00305 (4)	−0.00113 (4)
C4	0.00887 (4)	0.00559 (4)	0.00728 (7)	−0.00089 (3)	−0.00374 (4)	−0.00113 (4)
H4	0.0255 (9)	0.0160 (9)	0.0183 (9)	−0.0005 (8)	−0.0160 (6)	0.0020 (7)
C5	0.00932 (4)	0.00532 (4)	0.00818 (7)	−0.00096 (3)	−0.00367 (4)	−0.00155 (4)
H5	0.0305 (10)	0.0073 (4)	0.0262 (10)	−0.0025 (9)	−0.0104 (8)	−0.0065 (8)
N5	0.00673 (4)	0.00544 (4)	0.00516 (7)	−0.00123 (3)	−0.00240 (4)	−0.00123 (4)
C6	0.00864 (4)	0.00571 (4)	0.00720 (7)	−0.00135 (3)	−0.00336 (4)	−0.00175 (4)
H6	0.0286 (9)	0.0187 (9)	0.0209 (8)	−0.0072 (7)	−0.0137 (7)	−0.0084 (6)
N6	0.00744 (4)	0.00587 (4)	0.00565 (7)	−0.00123 (3)	−0.00309 (4)	−0.00130 (4)
C7	0.00625 (4)	0.00537 (4)	0.00521 (6)	−0.00125 (3)	−0.00220 (4)	−0.00118 (4)
N7	0.00837 (4)	0.00516 (3)	0.00603 (6)	−0.00109 (3)	−0.00369 (4)	−0.00128 (4)
H7	0.0280 (10)	0.0065 (4)	0.0230 (10)	−0.0019 (7)	−0.0089 (9)	−0.0038 (8)
C8	0.00757 (4)	0.00572 (4)	0.00613 (6)	−0.00125 (3)	−0.00316 (4)	−0.00132 (4)
H8	0.0303 (10)	0.0078 (4)	0.0254 (10)	−0.0022 (9)	−0.0097 (8)	−0.0066 (8)
N8	0.01027 (4)	0.00544 (4)	0.00814 (7)	−0.00142 (3)	−0.00385 (4)	−0.00201 (4)
H82	0.0242 (10)	0.0153 (10)	0.0160 (8)	−0.0024 (8)	−0.0140 (6)	0.0003 (8)
H81	0.0263 (10)	0.0183 (9)	0.0210 (9)	−0.0086 (7)	−0.0114 (7)	−0.0082 (6)
C9	0.00643 (4)	0.00541 (4)	0.00526 (6)	−0.00133 (3)	−0.00197 (4)	−0.00144 (4)
N9	0.00996 (4)	0.00689 (4)	0.00698 (7)	−0.00119 (3)	−0.00459 (4)	−0.00120 (4)
H92	0.0263 (10)	0.0188 (9)	0.0200 (8)	−0.0069 (8)	−0.0120 (7)	−0.0085 (6)
H91	0.0250 (10)	0.0083 (4)	0.0238 (10)	−0.0018 (9)	−0.0071 (9)	−0.0058 (7)
O8	0.015 (3)	0.009 (3)	0.011 (4)	−0.003 (2)	−0.004 (3)	−0.003 (3)

Geometric parameters (Å, º)

C10—O2	1.2886 (2)	N3—H3	1.0574 (6)
C10—O1	1.2889 (2)	C3—N5	1.3474 (2)
C10—O3	1.2918 (2)	C3—C4	1.3975 (2)
O4—H42	0.9710 (6)	C4—C5	1.3907 (2)
O4—H41	0.9801 (6)	C4—H4	1.0875 (6)
O5—H51	0.9665 (6)	C5—C6	1.3875 (2)
O5—H52	0.9734 (6)	C5—H5	1.0864 (6)
O6—H62	0.9801 (6)	N5—C7	1.3447 (2)
O6—H61	0.9811 (6)	C6—C7	1.4018 (2)
O7—H71	0.9726 (6)	C6—H6	1.0853 (6)
O7—H72	0.9839 (6)	N6—N7	1.3569 (2)
N1—C1	1.3249 (2)	N6—C8	1.2841 (2)
N1—H11	1.0051 (6)	C7—C8	1.4698 (2)
N1—H12	1.0368 (6)	N7—C9	1.3579 (2)
C1—N3	1.3517 (2)	N7—H7	1.0413 (6)
C1—N2	1.3273 (2)	C8—H8	1.0959 (6)
N2—H22	1.0130 (6)	N8—C9	1.3222 (2)
N2—H21	1.0293 (6)	N8—H82	0.9976 (6)
C2—N4	1.2861 (2)	N8—H81	1.0296 (6)
C2—C3	1.4683 (2)	C9—N9	1.3262 (2)
C2—H2	1.0947 (6)	N9—H91	1.0150 (6)
N3—N4	1.3616 (2)	N9—H92	1.0451 (6)
O2—C10—O1	120.556 (14)	C5—C4—H4	121.5 (3)
O2—C10—O3	119.441 (12)	C6—C5—C4	118.890 (13)
O1—C10—O3	120.002 (13)	C6—C5—H5	120.3 (3)
H42—O4—H41	103.7 (4)	C4—C5—H5	120.8 (3)
H51—O5—H52	105.1 (6)	C7—N5—C3	117.782 (13)
H62—O6—H61	110.9 (5)	C7—C6—C5	118.805 (13)
H71—O7—H72	107.8 (5)	C7—C6—H6	119.81 (18)
C1—N1—H11	118.5 (3)	C5—C6—H6	121.38 (18)
C1—N1—H12	120.0 (3)	N7—N6—C8	118.568 (14)
H11—N1—H12	120.9 (5)	N5—C7—C8	116.329 (12)
N3—C1—N1	118.770 (14)	N5—C7—C6	122.826 (12)
N3—C1—N2	120.749 (14)	C8—C7—C6	120.845 (12)
N1—C1—N2	120.481 (14)	N6—N7—C9	116.089 (13)
C1—N2—H22	119.76 (19)	N6—N7—H7	121.3 (3)
C1—N2—H21	118.2 (3)	C9—N7—H7	122.5 (3)
H22—N2—H21	121.4 (4)	N6—C8—C7	117.670 (13)
N4—C2—C3	119.403 (13)	N6—C8—H8	123.4 (3)
N4—C2—H2	122.8 (3)	C7—C8—H8	119.0 (3)
C3—C2—H2	117.7 (3)	C9—N8—H82	119.0 (3)
N4—N3—C1	118.390 (14)	C9—N8—H81	122.98 (18)
N4—N3—H3	121.0 (3)	H82—N8—H81	117.5 (4)
C1—N3—H3	120.6 (3)	N7—C9—N8	119.581 (12)
N5—C3—C2	115.938 (13)	N7—C9—N9	117.113 (13)
N5—C3—C4	122.971 (12)	N8—C9—N9	123.301 (13)
C2—C3—C4	121.075 (12)	C9—N9—H91	119.2 (3)
N3—N4—C2	116.315 (14)	C9—N9—H92	120.89 (19)
C3—C4—C5	118.725 (12)	H91—N9—H92	119.9 (4)
C3—C4—H4	119.8 (3)
N1—C1—N3—N4	179.015 (15)	H4—C4—C3—N5	−178.45 (19)
N1—C1—N3—H3	0.9 (7)	H4—C4—C5—C6	178.48 (14)
N1—C1—N2—H22	−177.82 (7)	H4—C4—C5—H5	−0.1 (8)
N1—C1—N2—H21	−6.4 (7)	C5—C6—C7—N5	0.05 (2)
H12—N1—C1—N3	3.8 (7)	C5—C6—C7—C8	−179.770 (12)
H12—N1—C1—N2	−176.33 (19)	C5—C4—C3—N5	0.43 (2)
H11—N1—C1—N3	174.36 (18)	H5—C5—C6—C7	178.78 (10)
H11—N1—C1—N2	−5.8 (7)	H5—C5—C6—H6	−0.8 (6)
C1—N3—N4—C2	−176.46 (2)	N5—C7—C8—N6	−178.592 (16)
N2—C1—N3—N4	−0.86 (2)	N5—C7—C8—H8	0.1 (6)
N2—C1—N3—H3	−178.97 (6)	N5—C7—C6—H6	179.60 (11)
H22—N2—C1—N3	2.1 (4)	C6—C7—C8—N6	1.24 (2)
H21—N2—C1—N3	173.5 (3)	C6—C7—C8—H8	179.93 (9)
C2—N4—N3—H3	1.6 (7)	H6—C6—C7—C8	−0.2 (4)
C2—C3—N5—C7	−178.750 (10)	N6—N7—C9—N8	0.78 (2)
C2—C3—C4—C5	178.883 (16)	N6—N7—C9—N9	−179.941 (15)
C2—C3—C4—H4	0.0 (6)	C7—C8—N6—N7	179.427 (6)
H2—C2—N4—N3	1.7 (7)	N7—N6—C8—H8	0.8 (6)
H2—C2—C3—N5	8.0 (7)	N7—C9—N8—H82	−2.2 (7)
H2—C2—C3—C4	−170.5 (5)	N7—C9—N8—H81	−174.16 (18)
N3—N4—C2—C3	−175.78 (2)	N7—C9—N9—H91	4.1 (7)
C3—N5—C7—C8	179.808 (6)	N7—C9—N9—H92	−179.15 (12)
C3—N5—C7—C6	−0.02 (2)	H7—N7—N6—C8	−4.3 (7)
C3—C4—C5—C6	−0.39 (2)	H7—N7—C9—N8	−176.4 (3)
C3—C4—C5—H5	−178.99 (11)	H7—N7—C9—N9	2.9 (7)
N4—C2—C3—N5	−174.40 (3)	C8—N6—N7—C9	178.511 (18)
N4—C2—C3—C4	7.05 (2)	N8—C9—N9—H91	−176.7 (3)
C4—C3—N5—C7	−0.22 (2)	N8—C9—N9—H92	0.1 (4)
C4—C5—C6—C7	0.16 (2)	H82—N8—C9—N9	178.6 (2)
C4—C5—C6—H6	−179.39 (8)	H81—N8—C9—N9	6.6 (4)

Hydrogen-bond geometry (Å, º)

D—H···A	D—H	H···A	D···A	D—H···A
N1—H12···O4^i	1.04 (1)	2.00 (1)	3.0067 (2)	162 (1)
N3—H3···O6^i	1.06 (1)	1.72 (1)	2.7803 (2)	177 (1)
N2—H22···O8^ii	1.01 (1)	1.74 (1)	2.637 (8)	146 (1)
O4—H42···O6	0.97 (1)	1.87 (1)	2.8084 (3)	162 (1)
O4—H42···O8	0.97 (1)	1.70 (1)	2.131 (8)	102 (1)
O5—H52···N5	0.97 (1)	1.96 (1)	2.9339 (2)	175 (1)
O6—H62···O7	0.98 (1)	1.77 (1)	2.7105 (3)	161 (1)
N1—H11···C10	1.00 (1)	2.54 (1)	3.4964 (3)	158 (1)
N1—H11···O2	1.00 (1)	1.73 (1)	2.7264 (3)	174 (1)
N2—H21···O3	1.03 (1)	1.83 (1)	2.8375 (3)	166 (1)
N8—H82···N6	1.00 (1)	2.23 (1)	2.6083 (2)	101 (1)
N9—H92···O4^iii	1.05 (1)	2.13 (1)	2.9610 (3)	135 (1)
O4—H41···C10^iv	0.98 (1)	2.55 (1)	3.4940 (3)	161 (1)
O4—H41···O2^iv	0.98 (1)	1.80 (1)	2.7693 (3)	170 (1)
O5—H51···C10^iv	0.97 (1)	2.65 (1)	3.5251 (3)	150 (1)
O5—H51···O1^iv	0.97 (1)	1.90 (1)	2.8622 (3)	176 (1)
O7—H71···O1^iv	0.97 (1)	1.79 (1)	2.7623 (2)	178 (1)
N7—H7···C10^v	1.04 (1)	2.45 (1)	3.4165 (2)	154 (1)
N7—H7···O1^v	1.04 (1)	1.70 (1)	2.7353 (2)	174 (1)
N9—H91···C10^v	1.01 (1)	2.62 (1)	3.5334 (2)	150 (1)
N9—H91···O3^v	1.01 (1)	1.73 (1)	2.7373 (2)	174 (1)
N8—H81···O2^vi	1.03 (1)	1.75 (1)	2.7613 (2)	166 (1)
O7—H72···O5^vii	0.98 (1)	1.80 (1)	2.7791 (3)	176 (1)
O6—H61···C10^viii	0.98 (1)	2.63 (1)	3.5170 (2)	150 (1)
O6—H61···O3^viii	0.98 (1)	1.73 (1)	2.6783 (2)	162 (1)
N8—H82···O7^ix	1.00 (1)	2.12 (1)	2.9302 (3)	138 (1)

Symmetry codes: (i) x−1, y, z; (ii) x−1, y+1, z; (iii) x, y, z+1; (iv) −x, −y+1, −z; (v) x+1, y−1, z+1; (vi) x+1, y, z+1; (vii) −x+1, −y, −z+1; (viii) −x+1, −y+1, −z; (ix) −x+1, −y+1, −z+1.

(neutron)

Crystal data

C9H15N9^2+·CO3^2^−·4(H2O)·0.02(OH)	Z = 2
Mr = 381.80	F(000) = 404
Triclinic, P1	Dx = 1.436 Mg m^−3
a = 8.2420 (2) Å	Neutrons radiation, λ = 0.600 Å
b = 8.6011 (3) Å	Cell parameters from 35336 reflections
c = 13.8821 (4) Å	θ = 7.2–79.0°
α = 72.792 (3)°	µ = 0.05 mm^−1
β = 78.998 (3)°	T = 100 K
γ = 70.789 (2)°	Parallel plate, colorless
V = 882.92 (5) Å^3	2.01 × 1.50 × 0.37 mm

Data collection

TOPAZ diffractometer	7146 reflections with I > 2σ(I)
Radiation source: The ORNL Spallation Neutron Scource	Rint = 0.096
time–of–flight Laue scans	θmax = 79.2°, θmin = 7.4°
Absorption correction: gaussian P. Coppens, L. Leiserowitz, D Rabinovich, Acta Cryst. (1965), 18, 1035-1038. P. R. Mallinson and K. W. Muir, J. Appl. Cryst. (1985). 18, 51-53.	h = −13→13
Tmin = 0.565, Tmax = 0.930	k = −14→14
44971 measured reflections	l = −23→23
7955 independent reflections

Refinement

Refinement on F^2	Hydrogen site location: difference Fourier map
Least-squares matrix: full	All H-atom parameters refined
R[F^2 > 2σ(F^2)] = 0.034	w = 1/[σ^2(Fo^2) + 0.0203P] where P = (Fo^2 + 2Fc^2)/3
wR(F^2) = 0.066	(Δ/σ)max = 0.001
S = 1.13	Δρmax = 1.10 e Å^−3
7955 reflections	Δρmin = −1.18 e Å^−3
452 parameters	Extinction correction: SHELXL-2014/7 (Sheldrick 2014, Fc^*=kFc[1+0.001xFc^2λ^3/sin(2θ)]^-1/4
0 restraints	Extinction coefficient: 0.0466 (13)

Special details

Geometry. All esds (except the esd in the dihedral angle between two l.s. planes) are estimated using the full covariance matrix. The cell esds are taken into account individually in the estimation of esds in distances, angles and torsion angles; correlations between esds in cell parameters are only used when they are defined by crystal symmetry. An approximate (isotropic) treatment of cell esds is used for estimating esds involving l.s. planes.

Fractional atomic coordinates and isotropic or equivalent isotropic displacement parameters (Ų)

	x	y	z	Uiso*/Ueq	Occ. (<1)
O1	−0.46937 (8)	0.99058 (8)	−0.27160 (4)	0.01276 (10)
O2	−0.42505 (8)	0.75634 (8)	−0.14350 (5)	0.01544 (11)
O3	−0.28550 (7)	0.94656 (8)	−0.16056 (5)	0.01218 (10)
C10	−0.39347 (6)	0.89770 (6)	−0.19218 (4)	0.00962 (8)
N1	−0.23176 (5)	0.56105 (5)	0.01260 (3)	0.01358 (7)
H11	−0.30091 (16)	0.63919 (16)	−0.04868 (9)	0.0246 (2)
H12	−0.22416 (18)	0.43571 (17)	0.03450 (10)	0.0283 (3)
C1	−0.17692 (6)	0.63177 (7)	0.06813 (4)	0.01067 (8)
N2	−0.19586 (5)	0.79785 (5)	0.04113 (3)	0.01747 (8)
H21	−0.24058 (18)	0.86648 (17)	−0.02803 (10)	0.0292 (3)
H22	−0.1502 (2)	0.8484 (2)	0.08238 (12)	0.0351 (3)
C2	0.03576 (6)	0.50419 (6)	0.28479 (4)	0.01021 (8)
H2	0.05986 (19)	0.36599 (17)	0.30380 (11)	0.0292 (3)
N3	−0.10168 (4)	0.53178 (5)	0.15255 (3)	0.01154 (6)
H3	−0.08508 (15)	0.40133 (15)	0.17237 (9)	0.0231 (2)
C3	0.09146 (6)	0.57500 (6)	0.35282 (4)	0.00890 (8)
N4	−0.04350 (4)	0.60514 (5)	0.20836 (3)	0.01096 (6)
C4	0.04328 (6)	0.75050 (6)	0.34224 (4)	0.01156 (8)
H4	−0.03556 (18)	0.83556 (17)	0.28260 (10)	0.0296 (3)
N5	0.18698 (4)	0.46088 (4)	0.42577 (3)	0.00884 (6)
C5	0.09722 (6)	0.81083 (7)	0.40960 (4)	0.01244 (9)
H5	0.06326 (19)	0.94620 (16)	0.40411 (11)	0.0305 (3)
N6	0.38876 (4)	0.44551 (4)	0.63384 (3)	0.00963 (6)
C6	0.19583 (6)	0.69466 (6)	0.48506 (4)	0.01081 (8)
H6	0.24078 (17)	0.73555 (17)	0.53943 (10)	0.0276 (2)
N7	0.48286 (4)	0.32902 (4)	0.70789 (3)	0.01000 (6)
H7	0.51101 (15)	0.19679 (14)	0.71631 (9)	0.0210 (2)
C7	0.23804 (6)	0.52064 (6)	0.49058 (4)	0.00840 (8)
N8	0.48736 (5)	0.55845 (5)	0.76132 (3)	0.01200 (6)
H81	0.53160 (16)	0.61538 (16)	0.80226 (10)	0.0242 (2)
H82	0.42110 (16)	0.63596 (16)	0.70202 (10)	0.0254 (2)
C8	0.34137 (6)	0.39250 (6)	0.56959 (4)	0.00939 (8)
H8	0.37312 (18)	0.25786 (16)	0.57057 (11)	0.0289 (3)
N9	0.62536 (5)	0.28047 (5)	0.84525 (3)	0.01241 (6)
C9	0.53293 (6)	0.39223 (6)	0.77244 (4)	0.00850 (8)
H91	0.66211 (15)	0.15273 (15)	0.84583 (10)	0.0239 (2)
H92	0.66677 (18)	0.31595 (18)	0.89661 (10)	0.0290 (3)
O4	0.75639 (9)	0.20278 (11)	0.04436 (6)	0.02262 (14)
H41	0.64236 (18)	0.2074 (2)	0.08215 (11)	0.0326 (3)
H42	0.83387 (18)	0.1733 (2)	0.09586 (11)	0.0330 (3)
O5	0.26444 (9)	0.09316 (8)	0.45045 (5)	0.01665 (11)
H51	0.33733 (17)	0.06369 (18)	0.39114 (10)	0.0286 (2)
H52	0.24256 (17)	0.21506 (16)	0.44004 (10)	0.0271 (2)
O6	0.94345 (8)	0.18309 (8)	0.19970 (5)	0.01548 (11)
H61	1.06664 (15)	0.12008 (17)	0.19894 (11)	0.0277 (2)
H62	0.87667 (16)	0.13713 (17)	0.26163 (10)	0.0270 (2)
O7	0.69938 (8)	0.09973 (8)	0.35046 (5)	0.01432 (11)
H71	0.62036 (16)	0.06665 (17)	0.32074 (10)	0.0259 (2)
H72	0.70552 (17)	0.03178 (17)	0.42029 (10)	0.0269 (2)
O8	0.975 (3)	−0.018 (3)	0.0776 (17)	0.007 (6)*	0.022 (3)
H88	1.078 (6)	−0.065 (6)	0.035 (4)	0.019 (9)*	0.022 (3)

Atomic displacement parameters (Ų)

	U^11	U^22	U^33	U^12	U^13	U^23
O1	0.0172 (2)	0.0113 (2)	0.0108 (2)	−0.00386 (18)	−0.00780 (18)	−0.00062 (18)
O2	0.0232 (3)	0.0112 (2)	0.0153 (3)	−0.0082 (2)	−0.0117 (2)	0.00151 (19)
O3	0.0133 (2)	0.0116 (2)	0.0134 (2)	−0.00471 (17)	−0.00655 (18)	−0.00121 (18)
C10	0.01177 (18)	0.00878 (19)	0.00928 (19)	−0.00275 (14)	−0.00515 (14)	−0.00143 (14)
N1	0.01667 (15)	0.01513 (16)	0.01062 (15)	−0.00475 (12)	−0.00651 (11)	−0.00235 (11)
H11	0.0296 (5)	0.0255 (6)	0.0207 (5)	−0.0083 (4)	−0.0123 (4)	−0.0022 (4)
H12	0.0385 (6)	0.0216 (6)	0.0273 (6)	−0.0112 (5)	−0.0112 (5)	−0.0022 (4)
C1	0.01178 (18)	0.0128 (2)	0.00769 (19)	−0.00351 (15)	−0.00406 (14)	−0.00106 (15)
N2	0.02573 (18)	0.01371 (17)	0.01433 (17)	−0.00637 (13)	−0.01111 (14)	0.00085 (12)
H21	0.0391 (6)	0.0237 (6)	0.0238 (6)	−0.0096 (5)	−0.0153 (5)	0.0035 (4)
H22	0.0528 (8)	0.0282 (7)	0.0334 (7)	−0.0165 (6)	−0.0212 (6)	−0.0050 (5)
C2	0.01185 (18)	0.0113 (2)	0.00806 (19)	−0.00310 (14)	−0.00452 (14)	−0.00147 (14)
H2	0.0437 (7)	0.0168 (5)	0.0298 (6)	−0.0072 (4)	−0.0175 (5)	−0.0031 (4)
N3	0.01417 (14)	0.01249 (15)	0.00865 (14)	−0.00379 (11)	−0.00547 (11)	−0.00106 (11)
H3	0.0298 (5)	0.0189 (5)	0.0218 (5)	−0.0067 (4)	−0.0102 (4)	−0.0027 (4)
C3	0.00975 (17)	0.00965 (19)	0.00756 (18)	−0.00267 (14)	−0.00334 (13)	−0.00118 (14)
N4	0.01259 (13)	0.01247 (15)	0.00818 (13)	−0.00362 (11)	−0.00450 (10)	−0.00100 (10)
C4	0.01420 (19)	0.0090 (2)	0.0112 (2)	−0.00169 (15)	−0.00661 (15)	−0.00075 (15)
H4	0.0369 (6)	0.0216 (6)	0.0283 (6)	−0.0015 (4)	−0.0205 (5)	−0.0007 (4)
N5	0.01051 (12)	0.00880 (14)	0.00758 (13)	−0.00251 (10)	−0.00356 (10)	−0.00131 (10)
C5	0.0158 (2)	0.0086 (2)	0.0134 (2)	−0.00206 (15)	−0.00685 (16)	−0.00172 (15)
H5	0.0413 (7)	0.0149 (5)	0.0359 (7)	−0.0023 (4)	−0.0162 (5)	−0.0063 (4)
N6	0.01163 (13)	0.00919 (14)	0.00874 (13)	−0.00246 (10)	−0.00461 (10)	−0.00165 (10)
C6	0.01362 (19)	0.0086 (2)	0.0114 (2)	−0.00275 (14)	−0.00563 (15)	−0.00209 (15)
H6	0.0373 (6)	0.0235 (6)	0.0280 (6)	−0.0078 (5)	−0.0153 (5)	−0.0091 (4)
N7	0.01308 (13)	0.00891 (14)	0.00893 (13)	−0.00232 (10)	−0.00533 (10)	−0.00198 (10)
H7	0.0277 (5)	0.0152 (5)	0.0207 (5)	−0.0045 (4)	−0.0088 (4)	−0.0037 (4)
C7	0.00942 (17)	0.00889 (19)	0.00734 (18)	−0.00257 (13)	−0.00334 (13)	−0.00125 (14)
N8	0.01626 (14)	0.00942 (15)	0.01200 (15)	−0.00339 (11)	−0.00580 (11)	−0.00282 (11)
H81	0.0294 (5)	0.0224 (5)	0.0258 (5)	−0.0088 (4)	−0.0089 (4)	−0.0082 (4)
H82	0.0314 (5)	0.0194 (5)	0.0235 (5)	−0.0035 (4)	−0.0131 (4)	−0.0008 (4)
C8	0.01109 (18)	0.00951 (19)	0.00825 (18)	−0.00262 (14)	−0.00481 (14)	−0.00124 (14)
H8	0.0419 (7)	0.0167 (5)	0.0306 (6)	−0.0039 (4)	−0.0175 (5)	−0.0063 (4)
N9	0.01609 (14)	0.01160 (15)	0.01015 (14)	−0.00281 (11)	−0.00696 (11)	−0.00148 (11)
C9	0.01006 (17)	0.00860 (19)	0.00739 (18)	−0.00249 (13)	−0.00289 (13)	−0.00184 (14)
H91	0.0287 (5)	0.0170 (5)	0.0253 (5)	−0.0040 (4)	−0.0107 (4)	−0.0024 (4)
H92	0.0394 (6)	0.0281 (6)	0.0256 (6)	−0.0091 (5)	−0.0162 (5)	−0.0086 (4)
O4	0.0250 (3)	0.0304 (4)	0.0164 (3)	−0.0124 (3)	−0.0049 (2)	−0.0049 (3)
H41	0.0310 (6)	0.0396 (8)	0.0295 (6)	−0.0168 (5)	−0.0040 (5)	−0.0041 (5)
H42	0.0317 (6)	0.0366 (7)	0.0316 (7)	−0.0086 (5)	−0.0130 (5)	−0.0052 (5)
O5	0.0226 (3)	0.0131 (3)	0.0131 (3)	−0.0049 (2)	−0.0011 (2)	−0.00243 (19)
H51	0.0314 (6)	0.0294 (6)	0.0234 (6)	−0.0069 (5)	0.0020 (4)	−0.0095 (4)
H52	0.0357 (6)	0.0193 (5)	0.0265 (6)	−0.0063 (4)	−0.0062 (5)	−0.0061 (4)
O6	0.0137 (2)	0.0140 (3)	0.0171 (3)	−0.00256 (19)	−0.00400 (19)	−0.0016 (2)
H61	0.0204 (5)	0.0251 (6)	0.0349 (7)	−0.0017 (4)	−0.0070 (4)	−0.0066 (5)
H62	0.0277 (5)	0.0267 (6)	0.0245 (6)	−0.0086 (4)	−0.0025 (4)	−0.0028 (4)
O7	0.0156 (2)	0.0136 (3)	0.0144 (3)	−0.00590 (19)	−0.00466 (19)	−0.00061 (19)
H71	0.0275 (5)	0.0264 (6)	0.0276 (6)	−0.0108 (4)	−0.0091 (4)	−0.0049 (4)
H72	0.0328 (6)	0.0260 (6)	0.0214 (5)	−0.0104 (4)	−0.0074 (4)	−0.0001 (4)

Geometric parameters (Å, º)

O1—C10	1.2877 (7)	N6—N7	1.3585 (5)
O2—C10	1.2875 (8)	C6—H6	1.0864 (14)
O3—C10	1.2905 (7)	C6—C7	1.4024 (7)
N1—H12	1.0136 (14)	N7—H7	1.0578 (12)
N1—H11	1.0464 (11)	N7—C9	1.3587 (6)
N1—C1	1.3266 (6)	C7—C8	1.4720 (6)
C1—N2	1.3282 (7)	N8—H82	1.0138 (11)
C1—N3	1.3519 (6)	N8—H81	1.0280 (13)
N2—H22	0.9977 (15)	N8—C9	1.3218 (6)
N2—H21	1.0308 (12)	C8—H8	1.0964 (14)
C2—H2	1.0959 (14)	N9—H92	1.0041 (13)
C2—N4	1.2855 (5)	N9—H91	1.0371 (13)
C2—C3	1.4709 (7)	N9—C9	1.3279 (5)
N3—H3	1.0412 (13)	O4—H42	0.9713 (15)
N3—N4	1.3633 (5)	O4—H41	0.9795 (16)
C3—N5	1.3477 (5)	O5—H51	0.9682 (15)
C3—C4	1.3980 (7)	O5—H52	0.9739 (15)
C4—H4	1.0849 (12)	O6—H62	0.9797 (14)
C4—C5	1.3925 (8)	O6—H61	0.9808 (13)
N5—C7	1.3458 (6)	O7—H72	0.9736 (14)
C5—H5	1.0867 (14)	O7—H71	0.9846 (14)
C5—C6	1.3883 (6)	O8—H88	0.97 (5)
N6—C8	1.2833 (6)
O2—C10—O1	120.49 (5)	C6—C5—C4	118.86 (5)
O2—C10—O3	119.52 (5)	C8—N6—N7	118.56 (4)
O1—C10—O3	119.98 (5)	H6—C6—C5	121.56 (8)
H12—N1—H11	119.01 (11)	H6—C6—C7	119.60 (8)
H12—N1—C1	121.20 (8)	C5—C6—C7	118.83 (5)
H11—N1—C1	119.13 (8)	H7—N7—N6	121.71 (7)
N1—C1—N2	120.48 (4)	H7—N7—C9	122.14 (7)
N1—C1—N3	118.80 (4)	N6—N7—C9	116.11 (4)
N2—C1—N3	120.73 (5)	N5—C7—C6	122.83 (4)
H22—N2—H21	121.93 (12)	N5—C7—C8	116.31 (4)
H22—N2—C1	119.48 (9)	C6—C7—C8	120.86 (4)
H21—N2—C1	117.91 (9)	H82—N8—H81	116.94 (11)
H2—C2—N4	122.74 (8)	H82—N8—C9	118.65 (9)
H2—C2—C3	117.69 (8)	H81—N8—C9	123.65 (8)
N4—C2—C3	119.53 (4)	H8—C8—N6	123.73 (8)
H3—N3—C1	120.23 (7)	H8—C8—C7	118.56 (8)
H3—N3—N4	121.34 (7)	N6—C8—C7	117.71 (4)
C1—N3—N4	118.39 (4)	H92—N9—H91	119.22 (11)
N5—C3—C4	123.05 (4)	H92—N9—C9	122.37 (9)
N5—C3—C2	115.98 (4)	H91—N9—C9	118.29 (8)
C4—C3—C2	120.95 (4)	N8—C9—N9	123.27 (4)
C2—N4—N3	116.38 (4)	N8—C9—N7	119.67 (4)
H4—C4—C5	121.89 (9)	N9—C9—N7	117.06 (4)
H4—C4—C3	119.43 (9)	H42—O4—H41	104.84 (14)
C5—C4—C3	118.68 (4)	H51—O5—H52	106.36 (13)
C7—N5—C3	117.75 (4)	H62—O6—H61	111.78 (12)
H5—C5—C6	119.89 (9)	H72—O7—H71	107.04 (12)
H5—C5—C4	121.24 (8)

Hydrogen-bond geometry (Å, º)

D—H···A	D—H	H···A	D···A	D—H···A
N1—H11···O2	1.0464 (11)	1.6868 (12)	2.7309 (7)	175.17 (12)
N1—H12···O4^i	1.0136 (14)	2.0246 (16)	3.0128 (9)	164.28 (11)
N2—H21···O3	1.0308 (12)	1.8273 (14)	2.8408 (7)	166.86 (14)
N3—H3···O6^i	1.0412 (13)	1.7450 (14)	2.7855 (8)	177.38 (11)
C6—H6···O7^ii	1.0864 (14)	2.5721 (15)	3.6379 (9)	166.67 (10)
N7—H7···O1^iii	1.0578 (12)	1.6869 (13)	2.7400 (7)	173.13 (10)
N7—H7···O3^iii	1.0578 (12)	2.6284 (12)	3.3830 (7)	127.91 (8)
N8—H81···O2^iv	1.0280 (13)	1.7586 (14)	2.7690 (8)	166.64 (11)
N8—H82···O7^ii	1.0138 (11)	2.1202 (14)	2.9419 (7)	136.76 (11)
N9—H91···O3^iii	1.0371 (13)	1.7069 (14)	2.7424 (7)	175.88 (11)
N9—H92···O4^v	1.0041 (13)	2.1503 (14)	2.9706 (8)	137.67 (12)
O4—H41···O2^vi	0.9795 (16)	1.7994 (15)	2.7740 (10)	172.87 (15)
O4—H42···O6	0.9713 (15)	1.8762 (16)	2.8173 (10)	162.36 (16)
O5—H51···O1^vi	0.9682 (15)	1.8950 (15)	2.8623 (10)	176.99 (13)
O5—H52···N5	0.9739 (15)	1.9670 (14)	2.9397 (8)	176.52 (13)
O6—H61···O3^vii	0.9808 (13)	1.7320 (13)	2.6817 (8)	161.94 (13)
O6—H62···O7	0.9797 (14)	1.7737 (14)	2.7193 (9)	161.15 (12)
O7—H71···O1^vi	0.9846 (14)	1.7838 (13)	2.7677 (9)	177.26 (13)
O7—H72···O5^viii	0.9736 (14)	1.8137 (14)	2.7850 (9)	175.05 (12)

Symmetry codes: (i) x−1, y, z; (ii) −x+1, −y+1, −z+1; (iii) x+1, y−1, z+1; (iv) x+1, y, z+1; (v) x, y, z+1; (vi) −x, −y+1, −z; (vii) −x+1, −y+1, −z; (viii) −x+1, −y, −z+1.