1. Introduction
A pigment is a coloured organic or inorganic solid that is usually insoluble in its application media. Pigments are used when extreme resistance to light, weather or temperature is needed, e.g. in automotive paints and other outdoor applications. Due to its insolubility, a pigment will usually retain its crystal structure throughout the preparation and application process. Various physico-chemical properties of pigments, including colour, thermal stability, particle shape, light fastness and fluorescence, are dependent on the crystal packing of the pigment molecules.
Despite the commercial significance of non-azo organic pigments, only few full crystal structures of these materials have been reported, because single crystals are hard to obtain. Pigments are practically insoluble, thus solution growth methods are inappropriate. Pigments are produced by precipitation reactions leading to very fine powders. The crystallites are often so small that they lead to considerable line broadening in X-ray powder diffraction. Single-crystal growth by sublimation methods, on the other hand, rarely produces crystals of suitable quality for conventional X-ray single-crystal structure determination. As a result, the important crystallographic features, controlling the crystal chemistry of such materials, have remained unknown for many pigments.
The purpose of this research is to study the feasibility of obtaining the crystal structures of the polymorphs of pigments with the aid of crystal structure prediction (CSP) tools in combination with limited experimental data. A similar approach was reported by Erk (2001) and Schmidt et al. (2005). One can obtain a good impression of the state-of-the-art in crystal structure prediction by looking at the results of CSP tests organized in recent years by the Cambridge Crystallographic Data Centre (see Day et al., 2005, and references therein). Three CSP `blind tests' in 1999, 2001 and 2004 indicated that certain progress has been made in the ab initio prediction of structures of `small' (not more than 25 atoms) rigid molecules containing only C, H, N and O atoms.
As an example of an insoluble pigment with poor crystal quality, unsubstituted linear trans-quinacridone, often called quinacridone (QA), a parent compound of various related pigments (Fig. 1), was chosen. A review on various quinacridone-related compounds has been given by Lincke (2000). In patent reports, at least six polymorphic forms of quinacridone have been reported on the basis of X-ray powder data (Reidinger & Struve, 1955; Struve, 1955; Manger & Struve, 1957; Filho & Oliveira, 1992). Only three of them, α, β and γ, are widely accepted. The crystal structure of the γ polymorph (Mizuguchi et al., 2002) is available from the Cambridge Structural Database (CSD), while for the β polymorph only the cell parameters, space group and a powder pattern have been published (Paulus et al., 1989; Struve, 1955). The structure of the α polymorph, derived through calculations and reported in the CSD (Lincke & Finzel, 1996), is likely to be incorrect, as discussed later. Thus, this structural information was not used further in our work. Crystal structure prediction of quinacridone polymorphs has been mentioned previously (Leusen, 1996) by Accelrys (formerly Molecular Simulations Ltd) as an example of the use of the Cerius2 Polymorph Predictor, but the details of the prediction and the resulting structures remained unpublished.
2. Materials and methods
2.1. Experimental methods
The pure β and γ polymorphs of QA were provided by Agfa-Gevaert NV and used as received. The powder of α-QA was obtained by slow water dilution of a concentrated H2SO4 solution of QA. For that, a solution of γ-QA in concentrated sulfuric acid was placed in one test tube and water in a second one. The tubes were connected with each other and the water tube was heated to increase the water vapour pressure in the system. After the precipitation of a dark violet powder in the sulfuric acid solution, the mixture was filtered and the powder of α-QA was collected and dried. X-ray powder diffraction (XRPD) patterns of the samples were obtained using a Philips PW 1710 diffractometer by scanning the samples over 2θ angles ranging from 2 to 35°, using Cu Kα radiation with wavelength 1.5418 Å.
2.2. Computational methods
Crystal structure prediction was conducted using the Accelrys Polymorph Predictor (PP) module in the Cerius2 package (Accelrys, 2002). In short, this program generates randomly a number of crystal packings in a chosen space group. The energy of these structures is minimized using molecular mechanics force fields. These structures are subsequently ranked, e.g. according to the lattice energy. The structures that are highest in rank are expected to have the highest probability of being found experimentally. For more details on the theory and background of the prediction procedures implemented in the PP, refer to the work of Gdanitz et al. (1993) and Karfunkel et al. (1993, 1994).
The molecule of quinacridone was first drawn in the three-dimensional sketcher and then energetically minimized in the Dreiding 2.21 force field. Electrostatic potential derived (ESPD) point charges were assigned to the atoms in the molecule. These charges were obtained by ab initio Hartree–Fock calculations at the given molecular geometry using the quantum mechanics program GAUSSIAN94 (Frisch et al., 1998) with the 6-31G* basis set. Before starting the prediction procedure, the geometry of the molecule was optimized again in the force field, but this did not contribute any significant changes to the geometry.
The Dreiding 2.21 force field was also used to minimize the energy of the crystal structures in the final step of the PP procedure. The calculation of the van der Waals and Coulomb interactions was performed using the Ewald summation method; the dielectric constant ∊r was kept constant (unlike the default Cerius2 setting of distance proportionality).
Rietveld refinement was performed using the DBWS-9006 program (Wiles & Young, 1981), implemented in Cerius2. The refinement was performed in the range of 2θ from 2 to 35°. Only the scale factor, background, cell parameters and the peak profile parameters were varied; the fractional atomic coordinates were kept fixed at this stage. The atomic coordinates were optimized by lattice energy minimization with fixed unit-cell parameters between the Rietveld steps.
Ab initio structure prediction is usually conducted in the space groups most frequently found among organic solids. 18 space groups together account for more than 93% of the molecular crystals reported (Baur & Kassner, 1992). However, since the space groups of quinacridone polymorphs are already known from the literature [ for α-QA (although questionable) and P21/c for β and γ-QA], we used this knowledge to limit the choice of the space groups. The flat and rigid molecule of quinacridone has an inversion centre. As intramolecular symmetry is not recognized during the PP prediction process, the prediction should be done in a space group that is a subgroup of the desired space group; the inversion centre is expected to be part of the symmetry elements of the space group. For the present case, the subgroup should contain all symmetry elements of the desired group, excluding the inversion centre. Therefore, we conducted the search in the space groups P1, P21 and Pc only. It is not necessary to search in all three subgroups: with a sufficient number of trials all P21/c structures will be found among the structures of each of these subgroups. For practical reasons (to reduce the calculation time), the search was conducted in all three subgroups and the results were subsequently merged.
A single quinacridone molecule was placed in the asymmetric unit. A total of 2000 trial structures were generated in each space group. Then, the clustering algorithm was applied to remove duplicate structures which resulted in total of 497 structures. Minimization of the total potential energy was then performed on the remaining structures, now also relaxing the rigid-body constraints. After clustering of the minimized structures, 61 structures remained. After the prediction process, an additional step was performed to find the actual symmetry (with the Find Symmetry tool) followed by a further energy minimization of the crystal structures in the space group found if it differed from the initial space group. The latter step is necessary to reach the optimal structure in the space group found, as the atomic coordinates obtained by imposing the new symmetry may not correspond exactly to an energy minimum, in particular in the case of intramolecular symmetry, where small deviations in bond lengths lead to a substantial increase in energy. As the prediction procedure contains a Monte Carlo stage, it can give different results for different runs. Three runs were enough to obtain reproducible results.
In order to evaluate the polymorph prediction results, X-ray powder diffraction patterns of the predicted structures were calculated and compared with the experimental ones. Due to the limited amount of predicted structures, it was possible to conduct the comparison manually for every structure. The quality of the experimental data was not good enough for indexing the patterns and deducing the structure, but sufficient to evaluate the PP results. The structures with powder patterns closest to the experimental ones were subjected to the Rietveld refinement.
4. Discussion
We only consider the three widely accepted polymorphs (α, β, γ), which have distinct powder patterns. Additional crystal phases mentioned in patents were reported on the basis of powder patterns of poor quality, giving rise to quite some uncertainty concerning their existence. Moreover, one of these phases, designated as γ′ (Deuschel et al., 1963), was proven to be γ-QA. Variations in the crystallite size and shape caused, in this case, the differences in the powder patterns (Potts et al., 1994).
4.1. γ-Quinacridone
The structure of the γ polymorph, determined by single-crystal X-ray measurements, is known from the CSD (reference code QNACRD04). The prediction of this structure served as a test to evaluate the quality of the force field and polymorph prediction in general.
Structures 1 and 2 in Table 1 have powder patterns similar to the experimental one (Fig. 4). Both structures agree with the experimental crystal structure from the CSD. The difference between the structures is so small that they can be considered to represent the same structure. Thus, the structure of γ-QA is predicted correctly and as the global minimum. The molecular packing of γ-QA is shown in Fig. 3. Each quinacridone molecule forms hydrogen bonds with four neighbours and they are packed in a criss-cross manner in ribbons along [001].
| Figure 3 Two projections of the predicted crystal structure of γ-QA. Dashed yellow lines represent hydrogen bonds. In the view along the c axis, hydrogen bonds are perpendicular to the plane of the picture, connecting four neighbouring molecules to each quinacridone molecule at the crossing points. |
| Figure 4 Powder patterns of structures 1, 2 and 3 of Table 1 compared with the experimental powder pattern of γ-QA; structure 3 is added to show the large dissimilarity. |
Rietveld refinement performed on both structures resulted in R factors of 18.5 and 14.1% (structures 1 and 2, respectively) and nearly the same final structure. In Table 2, we cite data for the second structure, as it showed a better Rietveld fit. As a test, the experimental structure QNACRD04 was refined against our experimental powder data in the same way as it was for the predicted structure. The R factor of the fit was 10.8%, a value which is an indication for the order of the accuracy of fits for the other predicted structures, given the limited quality of the experimental powder diffraction pattern.
4.2. β-Quinacridone
In spite of the commercial significance of β-QA and many attempts to grow single crystals, its structure is still unknown in the open literature. At present, only the cell parameters, the space group and a sketchy powder pattern are available (Paulus et al., 1989; Struve, 1955). According to these data, β-QA has a monoclinic structure, space group P21/c and cell parameters as listed in Table 2.
The polymorph prediction resulted in several structures for β-QA; all fit the experimental powder pattern well (Fig. 5). However, none of the predicted structures reproduces the small peak at 14° in the experimental pattern, which can be attributed to an impurity.
| Figure 5 Powder patterns of structures 4, 5, 13 and 14 of Table 1 compared with the experimental powder pattern of β-QA. |
Structure 4, designated here as β1, is monoclinic, space group P21/c, and has a long c axis (∼31 Å). Each molecule forms hydrogen bonds with two neighbours, forming ribbons that are arranged in two directions, along [110] and [10] (Fig. 6). The cell parameters of this structure fit well to the data reported in the literature (Table 2).
| Figure 6 Two projections of the predicted crystal structure of β1-QA: views along the short and long axes. Dashed yellow lines represent hydrogen bonds. The projection along the c axis shows that each molecule is connected to two neighbouring molecules via hydrogen bonds. |
Structure 5, which is labelled as β2, has the same molecular packing and a similar unit cell as structure 4, though their cell lengths do not correspond. The packing is shown in Fig. 7. The distances between the planes of the molecules in the stacks in these two structures, β1 and β2, are 2.86 Å and 2.89 Å, and the angles between molecules in different stacks are 60.2° and 60.3°, respectively. This may give the impression that the structures are defined in different settings of the space group P21/c and that they can be adjusted to have the same cell parameters. However, the unique axis 3.988 Å in the β1 structure clearly differs from its value of 5.720 Å in the β2 structure. This shows that these two structures are indeed different. β1-QA has somewhat lower lattice energy and higher density as compared with β2-QA.
| Figure 7 Two projections of the predicted crystal structure of β2-QA: views along the short and long axes. Dashed yellow lines represent hydrogen bonds. |
Structures 13 and 14 are similar but defined in different space groups. Structure 13 is monoclinic, space group P21/c. Structure 14 has nearly the same arrangement of molecules and cell parameters as structure 13, but in space group P21, which is quite uncommon for a molecule with inversion symmetry. Both structures are designated here as β3. The unit-cell axes lengths of β3 are comparable with the two previously mentioned β structures, but the unique axes differ (Table 1). The crystal structures of β3 resemble that of β1 and β2 with the difference that the stacks of hydrogen-bonded molecules are arranged in a nearly parallel way, along [01] (Fig. 8). The hydrogen-bond patterns are similar for all the β structures; each molecule has four hydrogen bonds to two neighbouring molecules and the resulting ribbons are packed in stacks.
| Figure 8 Two projections of the predicted crystal structure of β3-QA: views along the short and long axes. Dashed yellow lines represent hydrogen bonds. |
Rietveld refinement of the β1 and β2 structures against the experimental powder pattern resulted in both cases in a reasonable fit with R values of 13.0 and 12.6%, respectively. The difference can be noted at 15° in Fig. 5, where β1 has two peaks and β2 only one. The corresponding peak in the experimental pattern is broad and makes it impossible to distinguish between these two structures. It is notable that both refined structures are very close in the energy density plot (Fig. 2b). Both structures of β3 undergo cell compression during the Rietveld refinement due to the shortening of the long axis. The energy minimization performed between the Rietveld steps forces molecules to leave the nearly parallel arrangement of stacks (see Fig. 9). The density is nevertheless hardly affected (Fig. 2b). The refinement converges at R values of 19.7% and 20.6% for the structures 13 and 14, respectively, and leads to nearly the same structure in both cases. Therefore, only the crystallographic data for the refined structure 13 are listed in Table 2.
As a result of the Rietveld refinement, β1 and β2 show better similarity with the experimental powder pattern, and the published data favour the β1 structure. However, considering that Rexp of the β structure reported by Paulus et al. (1989) is not known and that the energies of all three β structures mentioned here are within the error of the prediction (a few kcal mol−1), each of them remains a candidate for the experimental structure.
4.3. α-Quinacridone
The structure of α-QA is deposited in the CSD under the reference code QNACRD03. It has been derived from the γ-QA structure by systematic variation of atomic coordinates (considering chemical bonding constraints) in order to fit the experimental powder pattern of α-QA (Lincke & Finzel, 1996). The only difference with the γ structure is the angle between the molecules in the criss-cross arrangement. Using the Dreiding force field, we found that this structure converts during minimization readily to the structure of the γ polymorph; this makes the QNACRD03 structure questionable.
Several structures, with numbers 8–12 in Table 1, are predicted by the PP to have powder patterns similar to the experimental one (Fig. 11). Structures 8–11 are triclinic, space group . Each quinacridone molecule shares hydrogen bonds with two neighbours, forming ribbons similar to those in β-QA, along [100]. These are packed in such a manner that all molecules in the structure lie parallel to one another (Fig. 10). Structure 12 is monoclinic, space group P21, and it has a different packing motif. Crystallization of centrosymmetric molecules in a non-centrosymmetric space group is, however, very rare.
| Figure 10 Two projections of the predicted crystal structure of α-QA. Dashed yellow lines represent hydrogen bonds. |
| Figure 11 Powder patterns of structures 8, 9, 10, 11 and 12 of Table 1 compared with the experimental powder pattern of α-QA. |
The quality of the experimental powder pattern was not high enough to obtain an acceptable fit of the patterns. Structures 10 and 11, as being closest to the experimental powder pattern, were used further for the Rietveld refinement. The refinement converged at R values of 19.6% and 21.4%, respectively, and produced the structures at substantially much higher energy as compared with the initial values (Fig. 2). Rietveld refinement of the structure 12 resulted in a relatively low R value of 14.3%. Visual inspection of the resulting XRPD pattern showed, however, that the structure 12 fitted the experimental pattern badly, mainly due to the mismatch of the peaks in the area 25–30°. The low R value is misleading in this case.
Structure 10, showing the best Rietveld fit among all five structures for the α polymorph, can be regarded as the closest to the experimental structure. It is debatable if such a high R value is acceptable to confirm the similarity of the predicted structure and the experimental one. A higher quality powder pattern, at least, is necessary to arrive at a better fit.
What are the other structures in Table 1? Showing no similarity with the known experimental powder patterns, they have molecular packing and hydrogen-bond patterns that are similar to the polymorphs discussed above. Structure 3 has a criss-cross molecular motif, similar to that of γ-QA, but its powder pattern is very different from the experimental one (Fig. 4). Structures 6 and 7, as β1-QA, consist of hydrogen-bonded ribbons running in two directions. Structure 15 contains squeezed layers of molecular ribbons which resemble the packing of the β2 structure after the Rietveld refinement.
4.4. Quality of the force field
An important question is how certain one can be of the results of the polymorph prediction approach outlined in this paper. The confidence in the results heavily depends on the quality of the force field and the atomic charges used.
The Dreiding force field (Mayo et al., 1990), a generic force field, can be applied to a wide range of heteroatomic organic molecules with reasonable accuracy. It shows the best results when used with the charges derived from high-level quantum mechanics calculations (ESPD charges). It models hydrogen bonds in crystal structures reasonably well. In certain cases, for example hydrogen bonds in acids, limitations have been reported (Payne et al., 1998). The force-field description of the structures in the present case appears to be quite adequate, since all three polymorphs were found within 3 kcal mol−1. Moreover, γ-QA, usually mentioned in the literature as a stable polymorph at room temperature, is found to have the lowest lattice energy.
If calculated lattice energies are close to the sublimation enthalpies, this can serve as another validation of the force field used. Lattice energy for the structures listed in Table 1 is 50–53 kcal mol−1. The sublimation enthalpy of quinacridone, as for many other pigments, is hard to measure due to the difficulties in precise measurements at high sublimation temperature. Therefore, no experimental data were found in the literature.
A number of promising developments have been reported in the area of force-field development, such as a better description of the electrostatics by means of a multipole expansion (van Eijck et al., 2001), using a large number of charged cells to describe the charge density (Gavezzotti, 2002, 2003), or developing a tailor-made potential for the compound under investigation from quantum chemical calculations (Day & Price, 2003).
4.5. X-ray powder data quality and limitations
As current progress in CSP showed, the problem in predicting is not the generation of all the possible structures (for that, present algorithms are good enough), but rather the ranking and choosing the right structures from the list. In the case of pure ab initio prediction, the first structure in the ranking is not always the experimentally stable polymorph. Here we combine ab initio CSP with experimental XRPD patterns to increase the chance of arriving at the right structure. Comparing powder patterns of predicted structures and experimental patterns usually shows large deviations. The differences arise from the limited prediction power of the force field and from the bad quality of the experimental diffraction pattern. An absolute match of the powder patterns in such cases is mostly regarded as a pure coincidence and is rather rare. Thus, the effectiveness of powder pattern comparison, in order to choose the right structures from the list of predicted ones, varies from case to case.
To get an impression of the deviation from the experimental pattern, it is interesting to compare powder patterns for γ-QA. Fig. 12 shows the experimental powder pattern together with four calculated ones, namely calculated from the CSD structure QNACRD04 before and after energy minimization in the force field, the powder pattern of structure 2 in the list of predicted structures (minimized), and that of structure 2 after Rietveld refinement against the experimental data. In this figure it is remarkable that the patterns (c) and (d) are identical. This means that the structure predicted by the PP is the same as the experimental structure minimized using the same force field and, therefore, it is probable that these two structures are also identical.
| Figure 12 Comparison of powder patterns of γ-QA: (a) experimental powder pattern, (b) calculated from experimental structure QNACRD04, (c) calculated from QNACRD04 after minimization in the force field, (d) calculated from structure 2, (e) calculated from structure 2 after Rietveld refinement (peak profiles not fitted for clarity). |
Fig. 12 also demonstrates how the effect of the force field leading to the difference between patterns (b) and (c) is eliminated by Rietveld refinement. Three peaks at 12–15° are overlapping in the experimental pattern (a) and distinct in the one calculated from the CSD structure (b). Due to the inaccuracy of the force field, minimization of the energy results in changes in the structure (Fig. 13) and the three peaks merge into two in pattern (c). The predicted structure (d) also has two peaks in this area. After the Rietveld refinement three peaks are resolved again in the pattern (e). These five structures are assumed to be the same, but their powder patterns differ to a certain extent. This difference should be kept in mind when comparing other predicted structures with the experimental powder patterns.
| Figure 13 Projection of the γ-QA structure QNACRD04 before (magenta) and after (green) energy minimization. |
Rietveld refinement is often used to complete the prediction procedure and to refine the trial structures against powder data. Full Rietveld refinement is, however, only worthwhile if the experimental data are of excellent quality (e.g. collected using synchrotron radiation) and the input structure is close enough to the experimental one. In the case of small particle size or poor crystallinity, like for many pigments, the low quality of the powder diffraction patterns cannot be improved by using more powerful X-ray sources. Moreover, Rietveld refinement is not efficient if the candidate crystal structure is not close enough to the experimental one and too few peaks in the simulated and experimental powder patterns exactly match (which happens very often).
An example of the intrinsic ambiguity of powder data can be found in the structure of α-QA that was deposited in the CSD under the reference code QNACRD03 (Lincke & Finzel, 1996). This structure has a molecular packing very close to the criss-cross arrangement of molecules in γ-QA (and far from α-QA, as predicted by us); only the angle between the molecules is different. However, this criss-cross α-QA structure produces a powder pattern very close to the experimental pattern of α-QA, which leads to confusion. The reason to have more confidence in our variant of α-QA is that it represents a real local minimum in a reliable force field.
In this case and in general, experimental X-ray powder diffraction patterns are not always unique `fingerprints' of the crystal structures. Similarity of powder patterns of different polymorphs is not uncommon [for example, two polymorphs of terephthalic acid (Bailey & Brown, 1967; Brown, 1984)]. Thus, the methods that use the powder patterns or other types of experimental information solely to distinguish the polymorphs contain by default some degree of uncertainty. For more accuracy, a combination of different physico-chemical characterization methods should be applied.
4.6. Stability order
Another indication of the quality of the present prediction is the stability order of the polymorphs. According to the calculated lattice energy of the structures, the order of descending stability for the predicted polymorphs is γ–β–α or γ–α–β (Table 1), depending on which structure to accept for β-QA. As was mentioned before, experimental data about the stability of these polymorphs are not available. γ-QA is often mentioned in the literature as the most stable form and α-QA as the least stable form (Jaffe, 2001) (α-QA is therefore not commercially important).
Moreover, during our crystal growth experiments by sublimation (Panina et al., 2007), we observed that β-QA sublimed faster and at a lower temperature than γ-QA. This leads us to the conclusion that, at least at the temperature of sublimation (673 K), the β polymorph is less stable than the γ polymorph. Possibly, β-QA is less stable than γ-QA for the whole range of temperatures (monotropically related polymorphs), since no polymorphic transition could be observed till 773 K (Jones et al., 1975). Assuming such a monotropic relation, our experimental data confirm the predicted stability order for β- and γ-QA.
The stable γ polymorph was ranked highest in the prediction. Experimentally observed metastable polymorphs, in general, are not necessarily next in this ranking, even in the case of perfect prediction results, as kinetics determine the circumstances under which they are formed and kinetics are not part of the prediction method.
5. Conclusions
The structures of three polymorphs of quinacridone were predicted using Polymorph Predictor in combination with XRPD patterns of limited quality. The known γ structure was predicted correctly as first in the ranking. Three different structures were proposed for the unknown β-QA that closely match the experimental powder pattern. For α-QA, a better XRPD pattern is needed to obtain a reliable structure prediction.
The present work demonstrates a method to obtain crystal structures using Polymorph Predictor in combination with experimental XRPD patterns. Ab initio prediction (without any help of additional information, like XRPD) is still not very reliable, even for simple molecules, mainly due to the limited quality of the force field. The use of X-ray powder data helps to improve the reliability of the prediction. Rietveld refinement can be used to improve the result. The effectiveness of the latter step is limited by the flexibility of the molecule and/or the quality of the XRPD pattern.
In the present prediction method, the XRPD pattern need not be indexable. We expect this method to be useful in the case of rigid organic molecules, if a corresponding non-indexable powder pattern of reasonable quality is available. Pigments are a good example of the application of this method, since pigment molecules are often rigid and XRPD patterns of pigments are usually difficult to index because of the peak broadening as a result of small crystal sizes.
Note added in proof. While the paper was in press, another paper on the single-crystal X-ray structure determination of β-quinacridone was published (Nishimura et al., 2006). The structure reported in this paper fits very nicely our β1-structure, which is the highest in ranking for the predicted β-structures.