Leucopterin, the white pigment in butterfly wings: structural analysis by PDF fit, FIDEL fit, Rietveld refinement, solid-state NMR and DFT-D.

Leucopterin (C6H5N5O3) is the white pigment in the wings of Pieris brassicae butterflies, and other butterflies; it can also be found in wasps and other insects. Its crystal structure and its tautomeric form in the solid state were hitherto unknown. Leucopterin turned out to be a variable hydrate, with 0.5 to about 0.1 molecules of water per leucopterin molecule. Under ambient conditions, the preferred state is the hemihydrate. Initially, all attempts to grow single crystals suitable for X-ray diffraction were to no avail. Attempts to determine the crystal structure by powder diffraction using the direct-space method failed, because the trials did not include the correct, but rare, space group P2/c. Attempts were made to solve the crystal structure by a global fit to the pair distribution function (PDF-Global-Fit), as described by Prill and co-workers [Schlesinger et al. (2021). J. Appl. Cryst. 54, 776-786]. The approach worked well, but the correct structure was not found, because again the correct space group was not included. Finally, tiny single crystals of the hemihydrate could be obtained, which allowed at least the determination of the crystal symmetry and the positions of the C, N and O atoms. The tautomeric state of the hemihydrate was assessed by multinuclear solid-state NMR spectroscopy. 15N CPMAS spectra showed the presence of one NH2 and three NH groups, and one unprotonated N atom, which agreed with the 1H MAS and 13C CPMAS spectra. Independently, the tautomeric state was investigated by lattice-energy minimizations with dispersion-corrected density functional theory (DFT-D) on 17 different possible tautomers, which also included the prediction of the corresponding 1H, 13C and 15N chemical shifts in the solid. All methods showed the presence of the 2-amino-3,5,8-H tautomer. The DFT-D calculations also confirmed the crystal structure. Heating of the hemihydrate results in a slow release of water between 130 and 250 °C, as shown by differential thermal analysis and thermogravimetry (DTA-TG). Temperature-dependent powder X-ray diffraction (PXRD) showed an irreversible continuous shift of the reflections upon heating, which reveals that leucopterin is a variable hydrate. This observation was also confirmed by PXRD of samples obtained under various synthetic and drying conditions. The crystal structure of a sample with about 0.2 molecules of water per leucopterin was solved by a fit with deviating lattice parameters (FIDEL), as described by Habermehl et al. [Acta Cryst. (2022), B78, 195-213]. A local fit, starting from the structure of the hemihydrate, as well as a global fit, starting from random structures, were performed, followed by Rietveld refinements. Despite dehydration, the space group remains P2/c. In both structures (hemihydrate and variable hydrate), the leucopterin molecules are connected by 2-4 hydrogen bonds into chains, which are connected by further hydrogen bonds to neighbouring chains. The molecular packing is very efficient. The density of leucopterin hemihydrate is as high as 1.909 kg dm-3, which is one of the highest densities for organic compounds consisting of C, H, N and O only. The high density might explain the good light-scattering and opacity properties of the wings of Pieris brassicae and other butterflies.

Ambiguous structure determination from powder data: four different structural models of 4,11-di-fluoro-quinacridone with similar X-ray powder patterns, fit to the PDF, SSNMR and DFT-D.

Four different structural models, which all fit the same X-ray powder pattern, were obtained in the structure determination of 4,11-di-fluoro-quinacridone (C20H10N2O2F2) from unindexed X-ray powder data by a global fit. The models differ in their lattice parameters, space groups, Z, Z', molecular packing and hydrogen bond patterns. The molecules form a criss-cross pattern in models A and B, a layer structure built from chains in model C and a criss-cross arrangement of dimers in model D. Nevertheless, all models give a good Rietveld fit to the experimental powder pattern with acceptable R-values. All molecular geometries are reliable, except for model D, which is slightly distorted. All structures are crystallochemically plausible, concerning density, hydrogen bonds, intermolecular distances etc. All models passed the checkCIF test without major problems; only in model A a missed symmetry was detected. All structures could have probably been published, although 3 of the 4 structures were wrong. The investigation, which of the four structures is actually the correct one, was challenging. Six methods were used: (1) Rietveld refinements, (2) fit of the crystal structures to the pair distribution function (PDF) including the refinement of lattice parameters and atomic coordinates, (3) evaluation of the colour, (4) lattice-energy minimizations with force fields, (5) lattice-energy minimizations by two dispersion-corrected density functional theory methods, and (6) multinuclear CPMAS solid-state NMR spectroscopy (1H, 13C, 19F) including the comparison of calculated and experimental chemical shifts. All in all, model B (perhaps with some disorder) can probably be considered to be the correct one. This work shows that a structure determination from limited-quality powder data may result in totally different structural models, which all may be correct or wrong, even if they are chemically sensible and give a good Rietveld refinement. Additionally, the work is an excellent example that the refinement of an organic crystal structure can be successfully performed by a fit to the PDF, and the combination of computed and experimental solid-state NMR chemical shifts can provide further information for the selection of the most reliable structure among several possibilities.

Structure determination from unindexed powder data from scratch by a global optimization approach using pattern comparison based on cross-correlation functions.

A method of ab initio crystal structure determination from powder diffraction data for organic and metal-organic compounds, which does not require prior indexing of the powder pattern, has been developed. Only a reasonable molecular geometry is required, needing knowledge of neither unit-cell parameters nor space group. The structures are solved from scratch by a global fit to the powder data using the new program FIDEL-GO (`FIt with DEviating Lattice parameters - Global Optimization'). FIDEL-GO uses a similarity measure based on cross-correlation functions, which allows the comparison of simulated and experimental powder data even if the unit-cell parameters deviate strongly. The optimization starts from large sets of random structures in various space groups. The unit-cell parameters, molecular position and orientation, and selected internal degrees of freedom are fitted simultaneously to the powder pattern. The optimization proceeds in an elaborate multi-step procedure with built-in clustering of duplicate structures and iterative adaptation of parameter ranges. The best structures are selected for an automatic Rietveld refinement. Finally, a user-controlled Rietveld refinement is performed. The procedure aims for the analysis of a wide range of `problematic' powder patterns, in particular powders of low crystallinity. The method can also be used for the clustering and screening of a large number of possible structure candidates and other application scenarios. Examples are presented for structure determination from unindexed powder data of the previously unknown structures of the nanocrystalline phases of 4,11-difluoro-, 2,9-dichloro- and 2,9-dichloro-6,13-dihydro-quinacridone, which were solved from powder patterns with 14-20 peaks only, and of the coordination polymer dichloro-bis(pyridine-N)copper(II).

Structure determination of organic compounds by a fit to the pair distribution function from scratch without prior indexing.

A method for the ab initio crystal structure determination of organic compounds by a fit to the pair distribution function (PDF), without prior knowledge of lattice parameters and space group, has been developed. The method is called 'PDF-Global-Fit' and is implemented by extension of the program FIDEL (fit with deviating lattice parameters). The structure solution is based on a global optimization approach starting from random structural models in selected space groups. No prior indexing of the powder data is needed. The new method requires only the molecular geometry and a carefully determined PDF. The generated random structures are compared with the experimental PDF and ranked by a similarity measure based on cross-correlation functions. The most promising structure candidates are fitted to the experimental PDF data using a restricted simulated annealing structure solution approach within the program TOPAS, followed by a structure refinement against the PDF to identify the correct crystal structure. With the PDF-Global-Fit it is possible to determine the local structure of crystalline and disordered organic materials, as well as to determine the local structure of unindexable powder patterns, such as nanocrystalline samples, by a fit to the PDF. The success of the method is demonstrated using barbituric acid as an example. The crystal structure of barbituric acid form IV solved and refined by the PDF-Global-Fit is in excellent agreement with the published crystal structure data.

Comparison and evaluation of pair distribution functions, using a similarity measure based on cross-correlation functions.

An approach for the comparison of pair distribution functions (PDFs) has been developed using a similarity measure based on cross-correlation functions. The PDF is very sensitive to changes in the local structure, i.e. small deviations in the structure can cause large signal shifts and significant discrepancies between the PDFs. Therefore, a comparison based on pointwise differences (e.g. R values and difference curves) may lead to the assumption that the investigated PDFs as well as the corresponding structural models are not in agreement at all, whereas a careful visual inspection of the investigated structural models and corresponding PDFs may reveal a relatively good match. To quantify the agreement of different PDFs for those cases an alternative approach is introduced: the similarity measure based on cross-correlation functions. In this paper, the power of this application of the similarity measure to the analysis of PDFs is highlighted. The similarity measure is compared with the classical R wp values as representative of the comparison based on pointwise differences as well as with the Pearson product-moment correlation coefficient, using polymorph IV of barbituric acid as an example.

Structure determination from powder data without prior indexing, using a similarity measure based on cross-correlation functions.

A method to refine organic crystal structures from powder diffraction data with incorrect lattice parameters has been developed. The method is based on the similarity measure developed by de Gelder et al. [J. Comput. Chem. (2001), 22, 273-289], using the cross- and auto-correlation functions of a simulated and an experimental powder pattern. The lattice parameters, molecular position, molecular orientation and selected intramolecular degrees of freedom are optimized until the similarity measure reaches a maximum; subsequently, a Rietveld refinement is carried out. The program FIDEL (FIt with DEviating Lattice parameters) implements this method. The procedure is also suitable for unindexed powder data, powder diagrams of very low quality and powder diagrams of non-phase-pure samples. Various applications are shown, including structure determinations from powder data using crystal structure predictions by standard force-field methods. Other useful applications include the automatic structure determination from powder data starting from the crystal structures of isostructural compounds (e.g. a solvate, hydrate or chemical derivative), or from crystal data measured at a different temperature or pressure.