research papers\(\def\hfill{\hskip 5em}\def\hfil{\hskip 3em}\def\eqno#1{\hfil {#1}}\)

Journal logoJOURNAL OF
APPLIED
CRYSTALLOGRAPHY
ISSN: 1600-5767

Crystal structure prediction of organic pigments: quinacridone as an example

CROSSMARK_Color_square_no_text.svg

aIMM Dept. for Solid State Chemistry, Radboud University Nijmegen, Toernooiveld 1, 6525 ED Nijmegen, The Netherlands, bIPI, University of Bradford, Bradford, West Yorkshire, BD7 1DP, UK, cAkzo Nobel Chemicals BV, Dept. CPT, PO Box 9300, 6800 SB Arnhem, The Netherlands, and dAgfa-Gevaert NV, Septestraat 27, B-2640 Mortsel, Belgium
*Correspondence e-mail: hugo.meekes@science.ru.nl

(Received 30 June 2006; accepted 20 October 2006)

The structures of the α, β and γ polymorphs of quinacridone (Pigment Violet 19) were predicted using Polymorph Predictor software in combination with X-ray powder diffraction patterns of limited quality. After generation and energy minimization of the possible structures, their powder patterns were compared with the experimental ones. On this basis, candidate structures for the polymorphs were chosen from the list of all structures. Rietveld refinement was used to validate the choice of structures. The predicted structure of the γ polymorph is in accordance with the experimental structure published previously. Three possible structures for the β polymorph are proposed on the basis of X-ray powder patterns comparison. It is shown that the α structure in the Cambridge Structural Database is likely to be in error, and a new α structure is proposed. The present work demonstrates a method to obtain crystal structures of industrially important pigments when only a low-quality X-ray powder diffraction pattern is available.

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[Erk, P. (2001). High Performance Pigments, edited by H. M. Smith, ch. 8, pp. 103-123. Weinheim: Wiley-VCH.]) and Schmidt et al. (2005[Schmidt, M. U., Ermrich, M. & Dinnebier, R. E. (2005). Acta Cryst. B61, 37-45.]). 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[Day, G. M. et al. (2005). Acta Cryst. B61, 511-527.], 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[link]), was chosen. A review on various quinacridone-related compounds has been given by Lincke (2000[Lincke, G. (2000). Dyes Pigments, 44, 101-122.]). In patent reports, at least six polymorphic forms of quinacridone have been reported on the basis of X-ray powder data (Reidinger & Struve, 1955[Reidinger, A. D. & Struve, W. S. (1955). US Patent 2844484.]; Struve, 1955[Struve, W. S. (1955). US Patent 2844485.]; Manger & Struve, 1957[Manger, C. W. & Struve, W. S. (1957). US Patent 2844581.]; Filho & Oliveira, 1992[Filho, D. S. & Oliveira, C. M. F. (1992). J. Mater. Sci. 27, 5101-5107.]). Only three of them, α, β and γ, are widely accepted. The crystal structure of the γ polymorph (Mizuguchi et al., 2002[Mizuguchi, J., Sasaki, T. & Tojo, K. (2002). Z. Kristallogr. 217, 249-250.]) 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[Paulus, E. F., Dietz, E., Kroh, F., Prokschy, F. & Lincke, G. (1989). 12th European Crystallographic Meeting, published in Z. Kristallogr. (1990, Suppl. 2), pp. 23-24.]; Struve, 1955[Struve, W. S. (1955). US Patent 2844485.]). The structure of the α polymorph, derived through calculations and reported in the CSD (Lincke & Finzel, 1996[Lincke, G. & Finzel, H. U. (1996). Cryst. Res. Technol. 31, 441-452.]), 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[Leusen, F. J. J. (1996). J. Cryst. Growth, 166, 900-903.]) 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.

[Figure 1]
Figure 1
Linear trans-quinacridone.

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[Accelrys (2002). Cerius2, version 4.71, User Guide. Accelrys Inc., San Diego, USA.]). 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[Gdanitz, R. J., Karfunkel, H. R. & Leusen, F. J. J. (1993). J. Mol. Graph. 11, 275-276.]) and Karfunkel et al. (1993[Karfunkel, H. R., Rohde, B., Leusen, F. J. J., Gdanitz, R. J. & Rihs, G. (1993). J. Comput. Chem. 14, 1125-1135.], 1994[Karfunkel, H. R., Leusen, F. J. J. & Gdanitz, R. J. (1994). J. Comput. Aided Mater. Des. 1, 177-185.]).

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[Frisch, M. J. et al. (1998). GAUSSIAN94 (Revision D.4), Gaussian Inc., Pittsburgh, USA.]) 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[Wiles, D. B. & Young, R. A. (1981). J. Appl. Cryst. 14, 149-151.]), 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.

2.3. Crystal structure prediction procedure

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[Baur, W. H. & Kassner, D. (1992). Acta Cryst. B48, 356-369.]). However, since the space groups of quinacridone polymorphs are already known from the literature [[P\bar{1}] 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.

3. Results

The prediction, as a result of three consecutive runs, ended with 103 structures ranging in lattice energy from −53 to −28 kcal mol−1 (1 cal = 4.184 J). Fig. 2[link] shows the predicted structures in terms of their lattice energy and density. We used density as a parameter to present the results in a clearer way and to distinguish structures with the same energy. A higher density might suggest a more stable structure, but this indicator is not as strong as low lattice energies.

[Figure 2]
Figure 2
(a) Results of the polymorph prediction as an energy-density distribution plot. Each star represents one predicted structure. The structures assigned to known polymorphs on the basis of powder patterns are marked with red symbols. Rietveld-refined structures of red-marked structures are represented by corresponding blue symbols. (b) Detailed view of the top-left area of (a).

Table 1[link] gives the information of the first 15 low-energy structures sorted with respect to the lattice energy. The last column shows the attribution of the structures to known polymorphs on the basis of the XRPD pattern. Structures in bold are considered to be close enough to the experimental structures to serve as input for the Rietveld refinement. The resulting structures after the Rietveld refinement are shown in Fig. 2[link]. Table 2[link] gives the crystallographic data of the refined structures and published experimental data on β- and γ-QA for comparison.

Table 1
Fifteen low-energy predicted structures ranked with respect to the lattice energy

Besides crystallographic details, the density and the lattice energies are given. The column `Crystal form' gives the assignment of the structure to known polymorphs on the basis of available XRPD patterns; structures that are not assigned have powder patterns which show no resemblance with the experimental ones. The structures in bold were subsequently subjected to Rietveld refinement.

    Cell axes Cell angles      
No. Space group a (Å) b (Å) c (Å) α (°) β (°) γ (°) d (g cm−3) Elattice (kcal mol−1) Crystal form
1 P21/c 14.130 4.008 13.424 90 70.77 90 1.445 −53.098 [\boldgamma]
2 P21/c 14.255 3.824 13.671 90 77.31 90 1.427 −52.930 [\boldgamma]
3 P21/c 11.348 6.280 13.031 90 130.48 90 1.469 −51.986
4 P21/c 5.734 3.966 31.096 90 104.51 90 1.515 −51.522 [\boldbeta_1]
5 P21/c 3.988 5.720 30.288 90 94.04 90 1.505 −51.352 [\boldbeta_2]
6 P21 28.184 4.933 4.977 90 81.94 90 1.514 −51.278
7 P21/c 5.019 4.891 28.088 90 94.05 90 1.508 −51.015
8 [P\bar{1}] 5.924 16.784 3.809 109.91 91.51 100.58 1.489 −50.634 α
9 [P\bar{1}] 15.767 3.874 5.813 89.67 95.71 97.23 1.480 −50.633 α
10 [{\bi P}\bar{\bf 1}] 6.922 3.873 14.770 85.27 85.57 63.07 1.476 −50.563 [\boldalpha]
11 [{\bi P}\bar{\bf 1}] 6.271 3.935 15.019 80.49 73.91 81.94 1.484 −50.555 [\boldalpha]
12 P21 3.887 29.214 6.920 90 116.52 90 1.475 −50.451 α
13 P21/c 5.745 31.288 3.931 90 90.29 90 1.468 −50.440 [\boldbeta_3]
14 P21 3.837 31.228 5.874 90 89.18 90 1.474 −50.373 [\boldbeta_3]
15 P21/c 15.887 7.013 31.143 90 11.87 90 1.454 −49.773

Table 2
Crystallographic data of predicted quinacridone polymorphs after Rietveld refinement and known experimental data for γ (Mizuguchi et al., 2002[Mizuguchi, J., Sasaki, T. & Tojo, K. (2002). Z. Kristallogr. 217, 249-250.]) and β (Paulus et al., 1989[Paulus, E. F., Dietz, E., Kroh, F., Prokschy, F. & Lincke, G. (1989). 12th European Crystallographic Meeting, published in Z. Kristallogr. (1990, Suppl. 2), pp. 23-24.]) polymorphs

The numbers labelling the predicted structures refer to the entries of Table 1[link].

  γ (No. 2) γ-exp β1 (No. 4) β2 (No. 5) β3 (No. 13) β-exp α (No. 10)
Space group P21/c P21/c P21/c P21/c P21/c P21/c [P\bar{1}]
a (Å) 13.998 13.700 5.621 4.057 5.557 5.692 6.748
b (Å) 3.890 3.840 4.059 5.622 30.674 3.975 3.776
c (Å) 13.383 13.350 31.177 30.938 4.149 30.02 14.400
α (°) 90 90 90 90 90 90 85.34
β (°) 79.04 100.09 100.94 98.17 89.00 96.76 81.76
γ (°) 90 90 90 90 90 90 70.20
Z 2 2 2 2 2 2 1
V3) 716 691 698 696 707 675 341
d (g cm−1) 1.450 1.485 1.485 1.486 1.467 1.447 1.519

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[Deuschel, W., Gundel, F., Wuest, H. & Daubach, E. (1963). US Patent 3074950.]), 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[Potts, G. D., Jones, W., Bullock, J. F., Andrews, S. J. & Maginn, S. J. (1994). 22, 2565-2566.]).

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[link] 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[link]. Each quinacridone molecule forms hydrogen bonds with four neighbours and they are packed in a criss-cross manner in ribbons along [001].[link]

[Figure 3]
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]
Figure 4
Powder patterns of structures 1, 2 and 3 of Table 1[link] 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[link], 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[Paulus, E. F., Dietz, E., Kroh, F., Prokschy, F. & Lincke, G. (1989). 12th European Crystallographic Meeting, published in Z. Kristallogr. (1990, Suppl. 2), pp. 23-24.]; Struve, 1955[Struve, W. S. (1955). US Patent 2844485.]). According to these data, β-QA has a monoclinic structure, space group P21/c and cell parameters as listed in Table 2[link].

The polymorph prediction resulted in several structures for β-QA; all fit the experimental powder pattern well (Fig. 5[link]). 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]
Figure 5
Powder patterns of structures 4, 5, 13 and 14 of Table 1[link] 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 [1[\bar{1}]0] (Fig. 6[link]). The cell parameters of this structure fit well to the data reported in the literature (Table 2[link]).

[Figure 6]
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[link]. 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]
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[link]). 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 [[\bar{1}]01] (Fig. 8[link]). 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]
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[link], 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[link]). 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[link]). The density is nevertheless hardly affected (Fig. 2b[link]). 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[link].

[Figure 9]
Figure 9
Structure of β3-QA minimized in the force field (red) and after the Rietveld refinement (blue).

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[Paulus, E. F., Dietz, E., Kroh, F., Prokschy, F. & Lincke, G. (1989). 12th European Crystallographic Meeting, published in Z. Kristallogr. (1990, Suppl. 2), pp. 23-24.]) 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[Lincke, G. & Finzel, H. U. (1996). Cryst. Res. Technol. 31, 441-452.]). 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[link], are predicted by the PP to have powder patterns similar to the experimental one (Fig. 11). Structures 8–11 are triclinic, space group [P\bar{1}]. 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[link][link]). 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]
Figure 10
Two projections of the predicted crystal structure of α-QA. Dashed yellow lines represent hydrogen bonds.
[Figure 11]
Figure 11
Powder patterns of structures 8, 9, 10, 11 and 12 of Table 1[link] 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[link]). 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[link]? 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[link]). 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[Mayo, S. L., Olafson, B. D. & Goddard, W. A. III (1990). J. Phys. Chem. 94, 8897-8909.]), 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[Payne, R. S., Roberts, R. J., Rowe, R. C. & Docherty, R. (1998). J. Comput. Chem. 19, 1-20.]). 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[link] 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[Eijck, B. P. van, Mooij, W. T. M. & Kroon, J. (2001). J. Comput. Chem. 22, 805-815.]), using a large number of charged cells to describe the charge density (Gavezzotti, 2002[Gavezzotti, A. (2002). J. Phys. Chem. B, 106, 4145-4154.], 2003[Gavezzotti, A. (2003). J. Phys. Chem. B, 107, 2344-2353.]), or developing a tailor-made potential for the compound under investigation from quantum chemical calculations (Day & Price, 2003[Day, G. M. & Price, S. L. (2003). J. Am. Chem. Soc. 125, 16434-16443.]).

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[link] 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]
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[link] 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[link]) 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]
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[Lincke, G. & Finzel, H. U. (1996). Cryst. Res. Technol. 31, 441-452.]). 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[Bailey, M. & Brown, C. J. (1967). Acta Cryst. 22, 387-391.]; Brown, 1984[Brown, C. J. (1984). Acta Cryst. C40, 1762.])]. 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[link]), 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[Jaffe, E. E. (2001). High Performance Pigments, edited by H. M. Smith, ch. 18, pp. 279-306. Weinheim: Wiley-VCH.]) (α-QA is therefore not commercially important).

Moreover, during our crystal growth experiments by sublimation (Panina et al., 2007[Panina, N., Janssen, F. F. B. J., Deroover, G., Meekes, H. & Vlieg, E. (2007). In preparation.]), 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[Jones, F., Okui, N. & Patterson, D. (1975). J. Soc. Dyers Col. 91, 361-365.]). 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[Nishimura, N., Senju, T. & Mizuguchi, J. (2006). Acta Cryst. E62, o4683-o4685.]). The structure reported in this paper fits very nicely our β1-structure, which is the highest in ranking for the predicted β-structures.

Supporting information


Computing details top

Figures top
[Figure 1]
[Figure 2]
[Figure 3]
[Figure 4]
[Figure 5]
[Figure 6]
[Figure 7]
[Figure 8]
[Figure 9]
[Figure 10]
[Figure 11]
[Figure 12]
[Figure 13]
(quin01) top
Crystal data top
?β = 70.7739°
Mr = ?V = ? Å3
Monoclinic, P21/CZ = ?
a = 14.1296 Å? radiation, λ = ? Å
b = 4.0083 Å × × mm
c = 13.4242 Å
Data collection top
h = ??l = ??
k = ??
Refinement top
Crystal data top
?β = 70.7739°
Mr = ?V = ? Å3
Monoclinic, P21/CZ = ?
a = 14.1296 Å? radiation, λ = ? Å
b = 4.0083 Å × × mm
c = 13.4242 Å
Data collection top
Refinement top
Fractional atomic coordinates and isotropic or equivalent isotropic displacement parameters (Å2) top
xyzUiso*/Ueq
C140.418980.115990.970430.00000*
C150.420330.178131.073420.00000*
C160.502630.059951.102430.00000*
H220.505590.101971.176570.00000*
N230.342600.223780.940610.00000*
C240.264740.394451.003830.00000*
C250.257930.467581.108220.00000*
C260.337630.356691.140500.00000*
C270.190420.497200.963050.00000*
C280.107530.678111.026690.00000*
C290.098990.748971.131780.00000*
C300.173540.644341.173130.00000*
O310.334570.420851.233000.00000*
H320.342580.176070.870180.00000*
H330.196890.442080.886970.00000*
H340.053850.758440.996580.00000*
H350.038470.879761.178140.00000*
H360.166020.698781.249910.00000*
Bond lengths (Å) top
C14—C16i1.403C25—C261.405
C14—C151.411C25—C301.414
C14—N231.340C26—O311.255
C15—C161.424C27—C281.403
C15—C261.413C27—H331.019
C16—C14i1.403C28—C291.405
C16—H221.024C28—H341.022
N23—C241.335C29—C301.407
N23—H320.965C29—H351.021
C24—C251.404C30—H361.024
C24—C271.397
Symmetry code: (i) x1, y, z2.
(quin10) top
Crystal data top
?β = 85.5673°
Mr = ?γ = 63.0734°
Triclinic, P1V = ? Å3
a = 6.9225 ÅZ = ?
b = 3.8729 Å? radiation, λ = ? Å
c = 14.7698 Å × × mm
α = 85.2710°
Data collection top
h = ??l = ??
k = ??
Refinement top
Crystal data top
?β = 85.5673°
Mr = ?γ = 63.0734°
Triclinic, P1V = ? Å3
a = 6.9225 ÅZ = ?
b = 3.8729 Å? radiation, λ = ? Å
c = 14.7698 Å × × mm
α = 85.2710°
Data collection top
Refinement top
Fractional atomic coordinates and isotropic or equivalent isotropic displacement parameters (Å2) top
xyzUiso*/Ueq
C300.289280.621390.312030.00000*
C290.188810.728850.387930.00000*
C280.027950.662200.385840.00000*
C270.142910.491450.306990.00000*
C160.212190.040250.006980.00000*
N230.156610.240390.156750.00000*
C260.255500.318650.154440.00000*
C250.172160.444740.233090.00000*
C240.042480.389300.231200.00000*
C150.134780.160380.075630.00000*
C140.078520.121460.081060.00000*
H220.366810.073620.013840.00000*
O310.441880.338650.155430.00000*
H320.300220.226210.157440.00000*
H330.299430.443020.304550.00000*
H340.101380.739470.441020.00000*
H350.274290.857590.444760.00000*
H360.447400.675080.314720.00000*
Bond lengths (Å) top
C30—C291.406C16—C151.422
C30—C251.413N23—C241.335
C30—H361.021N23—C141.342
C29—C281.405N23—H320.970
C29—H351.022C26—C251.407
C28—C271.402C26—C151.413
C28—H341.023C26—O311.258
C27—C241.396C25—C241.404
C27—H331.016C15—C141.411
C16—C14i1.402C14—C16i1.402
C16—H221.016
Symmetry code: (i) x, y, z.
(quin11) top
Crystal data top
?β = 73.9054°
Mr = ?γ = 81.9402°
Triclinic, P1V = ? Å3
a = 6.2707 ÅZ = ?
b = 3.9352 Å? radiation, λ = ? Å
c = 15.0191 Å × × mm
α = 80.4854°
Data collection top
h = ??l = ??
k = ??
Refinement top
Crystal data top
?β = 73.9054°
Mr = ?γ = 81.9402°
Triclinic, P1V = ? Å3
a = 6.2707 ÅZ = ?
b = 3.9352 Å? radiation, λ = ? Å
c = 15.0191 Å × × mm
α = 80.4854°
Data collection top
Refinement top
Fractional atomic coordinates and isotropic or equivalent isotropic displacement parameters (Å2) top
xyzUiso*/Ueq
C300.153790.110870.311050.00000*
C290.020280.265590.388460.00000*
C280.197110.415590.389440.00000*
C270.279000.414760.311470.00000*
C160.215300.169750.008820.00000*
N230.226700.285280.159580.00000*
C260.187440.049490.153600.00000*
C250.069590.105810.233300.00000*
C240.145190.266580.233740.00000*
C150.101260.030280.075190.00000*
C140.114650.142450.082520.00000*
H220.373040.288240.017120.00000*
O310.373390.218270.153580.00000*
H320.370570.415370.160890.00000*
H330.435540.523460.311060.00000*
H340.293980.527030.446110.00000*
H350.082540.269950.444170.00000*
H360.311900.005910.311630.00000*
Bond lengths (Å) top
C30—C291.405C16—C151.422
C30—C251.412N23—C241.336
C30—H361.021N23—C141.343
C29—C281.405N23—H320.970
C29—H351.020C26—C251.406
C28—C271.404C26—C151.413
C28—H341.024C26—O311.258
C27—C241.397C25—C241.404
C27—H331.016C15—C141.411
C16—C14i1.402C14—C16i1.402
C16—H221.016
Symmetry code: (i) x, y, z.
(quin12) top
Crystal data top
?β = 116.5216°
Mr = ?V = ? Å3
Monoclinic, P21Z = ?
a = 3.8873 Å? radiation, λ = ? Å
b = 29.2141 Å × × mm
c = 6.9195 Å
Data collection top
h = ??l = ??
k = ??
Refinement top
Crystal data top
?β = 116.5216°
Mr = ?V = ? Å3
Monoclinic, P21Z = ?
a = 3.8873 Å? radiation, λ = ? Å
b = 29.2141 Å × × mm
c = 6.9195 Å
Data collection top
Refinement top
Fractional atomic coordinates and isotropic or equivalent isotropic displacement parameters (Å2) top
xyzUiso*/Ueq
C10.461180.174960.133700.00000*
C20.366150.136780.043250.00000*
C30.201930.097540.170790.00000*
C40.147130.097000.385550.00000*
N50.010830.059930.509380.00000*
C60.236620.135130.476020.00000*
H70.188290.134250.632830.00000*
C80.394450.174360.350680.00000*
H90.462170.202140.416770.00000*
C100.087920.058070.098000.00000*
C110.057300.018820.229050.00000*
C120.095830.021930.441290.00000*
C130.164850.022690.163060.00000*
C140.313700.059530.306560.00000*
C150.352240.056330.518810.00000*
C160.244540.014860.584800.00000*
O170.106980.058360.088350.00000*
H180.002840.060390.652770.00000*
H190.580750.203240.040980.00000*
H200.419820.137780.115190.00000*
H210.131520.026450.009650.00000*
H220.277920.011090.738230.00000*
N230.420920.097620.239110.00000*
C240.558410.134610.363700.00000*
C250.613640.134860.578570.00000*
C260.498180.095430.650570.00000*
C270.649240.172940.274730.00000*
C280.809430.211960.401700.00000*
C290.877550.211990.618900.00000*
C300.780420.173770.707810.00000*
O310.516930.095620.836880.00000*
H320.407200.098140.095710.00000*
H330.600300.172340.117800.00000*
H340.878780.239960.337200.00000*
H351.000070.239860.713610.00000*
H360.834810.174480.866420.00000*
Bond lengths (Å) top
C1—C21.407C13—H211.016
C1—C81.405C14—C151.411
C1—H191.022C14—N231.343
C2—C31.413C15—C161.422
C2—H201.021C15—C261.413
C3—C41.404C16—H221.016
C3—C101.407N23—C241.336
C4—N51.335N23—H320.970
C4—C61.396C24—C251.404
N5—C121.342C24—C271.397
N5—H180.970C25—C261.406
C6—H71.016C25—C301.412
C6—C81.402C26—O311.258
C8—H91.022C27—C281.404
C10—C111.414C27—H331.016
C10—O171.258C28—C291.405
C11—C121.411C28—H341.024
C11—C131.423C29—C301.405
C12—C161.401C29—H351.020
C13—C141.402C30—H361.021
(quin13) top
Crystal data top
?β = 90.2924°
Mr = ?V = ? Å3
Monoclinic, P21/CZ = ?
a = 5.7454 Å? radiation, λ = ? Å
b = 31.2880 Å × × mm
c = 3.9308 Å
Data collection top
h = ??l = ??
k = ??
Refinement top
Crystal data top
?β = 90.2924°
Mr = ?V = ? Å3
Monoclinic, P21/CZ = ?
a = 5.7454 Å? radiation, λ = ? Å
b = 31.2880 Å × × mm
c = 3.9308 Å
Data collection top
Refinement top
Fractional atomic coordinates and isotropic or equivalent isotropic displacement parameters (Å2) top
xyzUiso*/Ueq
C140.092460.040010.093620.00000*
C150.118350.038390.088660.00000*
C160.209830.002270.179080.00000*
H220.362460.004570.306760.00000*
N230.184350.077500.189930.00000*
C240.086870.115040.113330.00000*
C250.124350.117230.063970.00000*
C260.221550.078140.169010.00000*
C270.201190.152140.217380.00000*
C280.103910.192430.152960.00000*
C290.110300.195040.016260.00000*
C300.223990.157750.125860.00000*
O310.403030.078910.347440.00000*
H320.326810.077690.321420.00000*
H330.355210.149720.340630.00000*
H340.187720.219590.230120.00000*
H350.184000.224150.061070.00000*
H360.378920.160480.249070.00000*
Bond lengths (Å) top
C14—C16i1.401C25—C261.407
C14—C151.411C25—C301.413
C14—N231.342C26—O311.260
C15—C161.422C27—C281.403
C15—C261.414C27—H331.014
C16—C14i1.401C28—C291.404
C16—H221.015C28—H341.023
N23—C241.336C29—C301.406
N23—H320.970C29—H351.020
C24—C251.404C30—H361.019
C24—C271.396
Symmetry code: (i) x, y, z.
(quin14) top
Crystal data top
?β = 89.1833°
Mr = ?V = ? Å3
Monoclinic, P21Z = ?
a = 3.8370 Å? radiation, λ = ? Å
b = 31.2284 Å × × mm
c = 5.8740 Å
Data collection top
h = ??l = ??
k = ??
Refinement top
Crystal data top
?β = 89.1833°
Mr = ?V = ? Å3
Monoclinic, P21Z = ?
a = 3.8370 Å? radiation, λ = ? Å
b = 31.2284 Å × × mm
c = 5.8740 Å
Data collection top
Refinement top
Fractional atomic coordinates and isotropic or equivalent isotropic displacement parameters (Å2) top
xyzUiso*/Ueq
C10.537630.194810.785530.00000*
C20.447170.157820.661510.00000*
C30.524330.116850.751100.00000*
C40.698320.113990.962520.00000*
N50.790190.076151.050880.00000*
C60.782970.150711.086960.00000*
H70.903330.147731.241080.00000*
C80.702090.191370.999880.00000*
H90.763490.218221.091610.00000*
C100.437570.077890.643980.00000*
C110.533270.037800.737820.00000*
C120.711400.038830.949600.00000*
C130.461750.002670.636280.00000*
C140.560900.040820.744900.00000*
C150.739160.039760.956640.00000*
C160.810860.000701.058130.00000*
O170.260780.078990.461770.00000*
H180.920220.076141.193790.00000*
H190.481500.224310.719930.00000*
H200.326260.161110.506730.00000*
H210.338050.004490.482680.00000*
H220.934640.002521.211710.00000*
N230.481460.078190.644150.00000*
C240.572380.116030.733160.00000*
C250.747080.118780.944380.00000*
C260.834690.079801.050910.00000*
C270.485650.152860.610150.00000*
C280.565000.193500.698420.00000*
C290.732070.196720.911950.00000*
C300.824090.159681.034760.00000*
O311.011710.080861.233040.00000*
H320.351320.078200.501280.00000*
H330.364420.149990.456340.00000*
H340.500320.220510.608690.00000*
H350.788730.226110.978220.00000*
H360.945820.162871.189280.00000*
Bond lengths (Å) top
C1—C21.406C13—H211.015
C1—C81.404C14—C151.411
C1—H191.021C14—N231.342
C2—C31.413C15—C161.422
C2—H201.020C15—C261.414
C3—C41.404C16—H221.015
C3—C101.408N23—C241.336
C4—N51.336N23—H320.970
C4—C61.396C24—C251.404
N5—C121.341C24—C271.396
N5—H180.970C25—C261.407
C6—H71.015C25—C301.413
C6—C81.402C26—O311.259
C8—H91.022C27—C281.403
C10—C111.414C27—H331.014
C10—O171.259C28—C291.404
C11—C121.411C28—H341.023
C11—C131.422C29—C301.406
C12—C161.401C29—H351.019
C13—C141.402C30—H361.020
(quin15) top
Crystal data top
?β = 11.8368°
Mr = ?V = ? Å3
Monoclinic, P21/CZ = ?
a = 15.9300 Å? radiation, λ = ? Å
b = 7.0129 Å × × mm
c = 31.1400 Å
Data collection top
h = ??l = ??
k = ??
Refinement top
Crystal data top
?β = 11.8368°
Mr = ?V = ? Å3
Monoclinic, P21/CZ = ?
a = 15.9300 Å? radiation, λ = ? Å
b = 7.0129 Å × × mm
c = 31.1400 Å
Data collection top
Refinement top
Fractional atomic coordinates and isotropic or equivalent isotropic displacement parameters (Å2) top
xyzUiso*/Ueq
C141.543260.095221.561320.00000*
C151.533900.105891.555820.00000*
C161.489570.199351.493360.00000*
H221.480000.343721.486960.00000*
N231.584100.190501.618890.00000*
C241.619050.102021.674230.00000*
C251.609050.097701.672340.00000*
C261.569440.199041.613720.00000*
C271.664630.215081.733620.00000*
C281.698290.131101.791170.00000*
C291.684330.068261.787830.00000*
C301.639820.182751.728770.00000*
O311.567200.378751.613820.00000*
H321.587120.328891.620310.00000*
H331.674570.358931.735500.00000*
H341.733060.214501.835270.00000*
H351.707600.130451.828970.00000*
H361.630010.327201.727360.00000*
Bond lengths (Å) top
C14—C16i1.401C25—C261.407
C14—C151.411C25—C301.413
C14—N231.342C26—O311.261
C15—C161.422C27—C281.402
C15—C261.414C27—H331.014
C16—C14i1.401C28—C291.403
C16—H221.014C28—H341.023
N23—C241.336C29—C301.406
N23—H320.971C29—H351.020
C24—C251.404C30—H361.019
C24—C271.395
Symmetry code: (i) x+3, y, z3.
(quin02) top
Crystal data top
?β = 77.3096°
Mr = ?V = ? Å3
Monoclinic, P21/CZ = ?
a = 14.2551 Å? radiation, λ = ? Å
b = 3.8236 Å × × mm
c = 13.6709 Å
Data collection top
h = ??l = ??
k = ??
Refinement top
Crystal data top
?β = 77.3096°
Mr = ?V = ? Å3
Monoclinic, P21/CZ = ?
a = 14.2551 Å? radiation, λ = ? Å
b = 3.8236 Å × × mm
c = 13.6709 Å
Data collection top
Refinement top
Fractional atomic coordinates and isotropic or equivalent isotropic displacement parameters (Å2) top
xyzUiso*/Ueq
C140.081240.130120.931840.00000*
C150.079000.136791.035500.00000*
C160.003490.004721.103740.00000*
H220.007070.006751.179360.00000*
N230.157810.252620.865640.00000*
C240.235190.387350.891770.00000*
C250.240830.405340.992870.00000*
C260.161060.277571.063150.00000*
C270.310120.509300.815690.00000*
C280.392670.652470.840110.00000*
C290.399750.668570.940900.00000*
C300.324320.546461.017630.00000*
O310.163550.289481.154230.00000*
H320.158010.243130.795110.00000*
H330.304420.495980.742710.00000*
H340.447200.745460.784860.00000*
H350.459980.771980.958420.00000*
H360.330630.561961.090790.00000*
Bond lengths (Å) top
C14—C16i1.403C25—C261.405
C14—C151.411C25—C301.414
C14—N231.340C26—O311.254
C15—C161.423C27—C281.403
C15—C261.413C27—H331.020
C16—C14i1.403C28—C291.405
C16—H221.024C28—H341.021
N23—C241.335C29—C301.407
N23—H320.964C29—H351.020
C24—C251.404C30—H361.025
C24—C271.397
Symmetry code: (i) x, y, z2.
(quin03) top
Crystal data top
?β = 130.4828°
Mr = ?V = ? Å3
Monoclinic, P21/CZ = ?
a = 11.3480 Å? radiation, λ = ? Å
b = 6.2795 Å × × mm
c = 13.0313 Å
Data collection top
h = ??l = ??
k = ??
Refinement top
Crystal data top
?β = 130.4828°
Mr = ?V = ? Å3
Monoclinic, P21/CZ = ?
a = 11.3480 Å? radiation, λ = ? Å
b = 6.2795 Å × × mm
c = 13.0313 Å
Data collection top
Refinement top
Fractional atomic coordinates and isotropic or equivalent isotropic displacement parameters (Å2) top
xyzUiso*/Ueq
C140.418080.127700.616540.00000*
C150.418150.182780.511550.00000*
C160.501810.052740.394000.00000*
H220.505830.089820.315290.00000*
N230.341040.246570.728580.00000*
C240.266340.425880.746320.00000*
C250.258340.494940.648500.00000*
C260.333630.367300.532730.00000*
C270.197040.542310.864390.00000*
C280.118090.732140.886570.00000*
C290.108190.802010.789980.00000*
C300.178080.685290.670910.00000*
O310.325430.419800.444590.00000*
H320.339410.200240.800090.00000*
H330.203630.489920.934570.00000*
H340.067840.818810.972540.00000*
H350.051050.940100.805990.00000*
H360.170860.740690.601220.00000*
Bond lengths (Å) top
C14—C16i1.403C25—C261.405
C14—C151.411C25—C301.414
C14—N231.340C26—O311.255
C15—C161.423C27—C281.404
C15—C261.412C27—H331.019
C16—C14i1.403C28—C291.405
C16—H221.024C28—H341.023
N23—C241.336C29—C301.406
N23—H320.965C29—H351.020
C24—C251.404C30—H361.024
C24—C271.396
Symmetry code: (i) x+1, y, z1.
(quin04) top
Crystal data top
?β = 104.5126°
Mr = ?V = ? Å3
Monoclinic, P21/CZ = ?
a = 5.7342 Å? radiation, λ = ? Å
b = 3.9658 Å × × mm
c = 31.0958 Å
Data collection top
h = ??l = ??
k = ??
Refinement top
Crystal data top
?β = 104.5126°
Mr = ?V = ? Å3
Monoclinic, P21/CZ = ?
a = 5.7342 Å? radiation, λ = ? Å
b = 3.9658 Å × × mm
c = 31.0958 Å
Data collection top
Refinement top
Fractional atomic coordinates and isotropic or equivalent isotropic displacement parameters (Å2) top
xyzUiso*/Ueq
C500.938550.556280.540560.00000*
C510.150050.374710.540470.00000*
C520.210040.322490.499290.00000*
H580.362520.196780.498230.00000*
N590.876140.617460.578620.00000*
C600.004820.505320.617950.00000*
C610.219900.324790.621890.00000*
C620.286090.256430.582140.00000*
C630.918760.575720.655160.00000*
C640.048890.472230.697480.00000*
C650.266910.298100.702110.00000*
C660.352750.223440.664600.00000*
O670.469310.076930.583760.00000*
H680.733070.749170.578270.00000*
H690.760980.701950.651350.00000*
H700.985450.523710.724670.00000*
H710.362740.223540.732890.00000*
H720.511390.095980.668800.00000*
Bond lengths (Å) top
C50—C52i1.401C61—C621.407
C50—C51ii1.411C61—C661.414
C50—N591.341C62—O671.260
C51—C50iii1.411C63—C60ii1.395
C51—C521.422C63—C64ii1.402
C51—C621.414C63—H691.015
C52—C50i1.401C64—C63iii1.402
C52—H581.014C64—C651.404
N59—C60ii1.336C64—H70iii1.021
N59—H680.970C65—C661.406
C60—N59iii1.336C65—H711.021
C60—C611.405C66—H721.020
C60—C63iii1.395H70—C64ii1.021
Symmetry codes: (i) x+1, y+1, z+1; (ii) x+1, y, z; (iii) x1, y, z.
(quin05) top
Crystal data top
?β = 94.0361°
Mr = ?V = ? Å3
Monoclinic, P21/CZ = ?
a = 3.9882 Å? radiation, λ = ? Å
b = 5.7197 Å × × mm
c = 30.2878 Å
Data collection top
h = ??l = ??
k = ??
Refinement top
Crystal data top
?β = 94.0361°
Mr = ?V = ? Å3
Monoclinic, P21/CZ = ?
a = 3.9882 Å? radiation, λ = ? Å
b = 5.7197 Å × × mm
c = 30.2878 Å
Data collection top
Refinement top
Fractional atomic coordinates and isotropic or equivalent isotropic displacement parameters (Å2) top
xyzUiso*/Ueq
C500.963530.614750.540580.00000*
C510.145840.403400.540230.00000*
C520.177970.289290.499030.00000*
H580.303780.135610.497770.00000*
N590.920870.727230.578640.00000*
C600.052650.650490.617770.00000*
C610.236230.440980.621430.00000*
C620.285100.322410.581670.00000*
C630.000270.785100.655080.00000*
C640.124320.710670.697180.00000*
C650.303020.499580.701450.00000*
C660.359610.364640.663890.00000*
O670.465570.141450.583140.00000*
H680.788660.869690.578380.00000*
H690.870600.937460.651480.00000*
H700.084810.809030.724500.00000*
H710.393950.444620.732010.00000*
H720.490120.211910.667880.00000*
Bond lengths (Å) top
C50—C52i1.401C60—C631.396
C50—C51ii1.411C61—C621.407
C50—N591.342C61—C661.413
C51—C50iii1.411C62—O671.260
C51—C521.422C63—C641.402
C51—C621.414C63—H69iii1.014
C52—C50i1.401C64—C651.403
C52—H581.014C64—H701.022
N59—C60ii1.336C65—C661.406
N59—H680.970C65—H711.020
C60—N59iii1.336C66—H721.020
C60—C611.405H69—C63ii1.014
Symmetry codes: (i) x+1, y+1, z+1; (ii) x+1, y, z; (iii) x1, y, z.
(quin06) top
Crystal data top
?β = 81.9358°
Mr = ?V = ? Å3
Monoclinic, P21Z = ?
a = 28.1840 Å? radiation, λ = ? Å
b = 4.9326 Å × × mm
c = 4.9767 Å
Data collection top
h = ??l = ??
k = ??
Refinement top
Crystal data top
?β = 81.9358°
Mr = ?V = ? Å3
Monoclinic, P21Z = ?
a = 28.1840 Å? radiation, λ = ? Å
b = 4.9326 Å × × mm
c = 4.9767 Å
Data collection top
Refinement top
Fractional atomic coordinates and isotropic or equivalent isotropic displacement parameters (Å2) top
xyzUiso*/Ueq
C10.453870.249510.198500.00000*
C20.415190.075520.174360.00000*
C30.372680.092550.361420.00000*
C40.369960.286880.569290.00000*
N50.330580.314510.750900.00000*
C60.408550.458000.593890.00000*
H70.405840.595200.746760.00000*
C80.450700.441570.408140.00000*
H90.478860.568140.425000.00000*
C100.332090.072590.355620.00000*
C110.290420.044070.548270.00000*
C120.291590.159190.746840.00000*
C130.248190.203120.555240.00000*
C140.208410.159530.753160.00000*
C150.209580.043730.951730.00000*
C160.251810.202780.944760.00000*
O170.333140.251030.173420.00000*
H180.330610.453120.889430.00000*
H190.484440.236920.062400.00000*
H200.418370.061390.019790.00000*
H210.246270.351340.416070.00000*
H220.253730.351001.083930.00000*
N230.169420.314850.749100.00000*
C240.130040.287230.930710.00000*
C250.127320.092891.138580.00000*
C260.167910.072251.144380.00000*
C270.091450.458350.906110.00000*
C280.049300.441911.091860.00000*
C290.046130.249851.301500.00000*
C300.084810.075871.325640.00000*
O310.166860.250691.326580.00000*
H320.169390.453460.610570.00000*
H330.094160.595540.753240.00000*
H340.021140.568481.075000.00000*
H350.015560.237271.437600.00000*
H360.081630.061051.480210.00000*
Bond lengths (Å) top
C1—C21.406C13—H211.014
C1—C81.403C14—C151.411
C1—H191.020C14—N231.342
C2—C31.413C15—C161.422
C2—H201.018C15—C261.415
C3—C41.405C16—H221.014
C3—C101.408N23—C241.336
C4—N51.336N23—H320.971
C4—C61.396C24—C251.405
N5—C121.342C24—C271.396
N5—H180.971C25—C261.408
C6—H71.013C25—C301.413
C6—C81.402C26—O311.261
C8—H91.022C27—C281.402
C10—C111.415C27—H331.013
C10—O171.261C28—C291.403
C11—C121.411C28—H341.022
C11—C131.422C29—C301.406
C12—C161.401C29—H351.020
C13—C141.401C30—H361.018
(quin07) top
Crystal data top
?β = 94.0515°
Mr = ?V = ? Å3
Monoclinic, P21/CZ = ?
a = 5.0189 Å? radiation, λ = ? Å
b = 4.8913 Å × × mm
c = 28.0879 Å
Data collection top
h = ??l = ??
k = ??
Refinement top
Crystal data top
?β = 94.0515°
Mr = ?V = ? Å3
Monoclinic, P21/CZ = ?
a = 5.0189 Å? radiation, λ = ? Å
b = 4.8913 Å × × mm
c = 28.0879 Å
Data collection top
Refinement top
Fractional atomic coordinates and isotropic or equivalent isotropic displacement parameters (Å2) top
xyzUiso*/Ueq
C140.640440.026510.541470.00000*
C150.436000.171020.540280.00000*
C160.296010.194050.498180.00000*
H220.147120.332930.496280.00000*
N230.777760.058810.580390.00000*
C240.730600.092810.619670.00000*
C250.533040.296730.622340.00000*
C260.387130.332370.581790.00000*
C270.883900.040580.658240.00000*
C280.845320.193550.700370.00000*
C290.649800.398870.703420.00000*
C300.493890.450890.664790.00000*
O310.208080.513570.582730.00000*
H320.916930.196700.580440.00000*
H331.024070.108620.655530.00000*
H340.958070.156260.728600.00000*
H350.620770.511120.733940.00000*
H360.353640.601320.667990.00000*
Bond lengths (Å) top
C14—C16i1.401C25—C261.408
C14—C151.411C25—C301.413
C14—N231.342C26—O311.261
C15—C161.422C27—C281.402
C15—C261.415C27—H331.013
C16—C14i1.401C28—C291.402
C16—H221.014C28—H341.022
N23—C241.337C29—C301.405
N23—H320.971C29—H351.020
C24—C251.405C30—H361.018
C24—C271.396
Symmetry code: (i) x1, y, z+1.
(quin08) top
Crystal data top
?β = 91.5084°
Mr = ?γ = 100.5753°
Triclinic, P1V = ? Å3
a = 5.9241 ÅZ = ?
b = 16.7841 Å? radiation, λ = ? Å
c = 3.8090 Å × × mm
α = 109.9087°
Data collection top
h = ??l = ??
k = ??
Refinement top
Crystal data top
?β = 91.5084°
Mr = ?γ = 100.5753°
Triclinic, P1V = ? Å3
a = 5.9241 ÅZ = ?
b = 16.7841 Å? radiation, λ = ? Å
c = 3.8090 Å × × mm
α = 109.9087°
Data collection top
Refinement top
Fractional atomic coordinates and isotropic or equivalent isotropic displacement parameters (Å2) top
xyzUiso*/Ueq
C300.020110.321450.248470.00000*
C290.191890.395170.437960.00000*
C280.401670.386630.587400.00000*
C270.436200.304090.555020.00000*
C160.212950.001940.174770.00000*
N230.303840.154250.360100.00000*
C260.100610.160200.019480.00000*
C250.056760.238170.211080.00000*
C240.264390.231000.372890.00000*
C150.058710.078880.004490.00000*
C140.154280.079570.181430.00000*
H220.368790.004450.299750.00000*
O310.281620.163160.152260.00000*
H320.446830.153760.488830.00000*
H330.586550.296990.663260.00000*
H340.528190.439980.719780.00000*
H350.164470.455180.464970.00000*
H360.130380.329170.140060.00000*
Bond lengths (Å) top
C30—C291.407C16—C151.422
C30—C251.414N23—C241.336
C30—H361.021N23—C141.341
C29—C281.405N23—H320.970
C29—H351.021C26—C251.407
C28—C271.402C26—C151.414
C28—H341.022C26—O311.259
C27—C241.395C25—C241.404
C27—H331.015C15—C141.411
C16—C14i1.401C14—C16i1.401
C16—H221.015
Symmetry code: (i) x, y, z.
(quin09) top
Crystal data top
?β = 95.7095°
Mr = ?γ = 97.2323°
Triclinic, P1V = ? Å3
a = 15.7670 ÅZ = ?
b = 3.8742 Å? radiation, λ = ? Å
c = 5.8128 Å × × mm
α = 89.6689°
Data collection top
h = ??l = ??
k = ??
Refinement top
Crystal data top
?β = 95.7095°
Mr = ?γ = 97.2323°
Triclinic, P1V = ? Å3
a = 15.7670 ÅZ = ?
b = 3.8742 Å? radiation, λ = ? Å
c = 5.8128 Å × × mm
α = 89.6689°
Data collection top
Refinement top
Fractional atomic coordinates and isotropic or equivalent isotropic displacement parameters (Å2) top
xyzUiso*/Ueq
C300.319330.033090.260050.00000*
C290.393070.095800.153470.00000*
C280.385300.258410.062890.00000*
C270.303380.298680.169240.00000*
C160.002500.176940.210440.00000*
N230.154000.230700.168190.00000*
C260.158990.122070.239880.00000*
C250.236770.003760.151060.00000*
C240.230210.176140.061930.00000*
C150.078240.065340.127330.00000*
C140.079470.114700.084250.00000*
H220.005490.302970.363440.00000*
O310.161680.298540.422500.00000*
H320.153680.361130.311200.00000*
H330.296810.417220.324970.00000*
H340.438730.348270.142260.00000*
H350.452520.069440.233500.00000*
H360.326530.152000.416560.00000*
Bond lengths (Å) top
C30—C291.407C16—C151.422
C30—C251.413N23—C241.336
C30—H361.020N23—C141.341
C29—C281.405N23—H320.970
C29—H351.021C26—C251.407
C28—C271.402C26—C151.414
C28—H341.022C26—O311.259
C27—C241.395C25—C241.404
C27—H331.015C15—C141.411
C16—C14i1.401C14—C16i1.401
C16—H221.014
Symmetry code: (i) x, y, z.

Experimental details

(quin01)(quin10)(quin11)(quin12)
Crystal data
Chemical formula????
Mr????
Crystal system, space groupMonoclinic, P21/CTriclinic, P1Triclinic, P1Monoclinic, P21
Temperature (K)????
a, b, c (Å)14.1296, 4.0083, 13.42426.9225, 3.8729, 14.76986.2707, 3.9352, 15.01913.8873, 29.2141, 6.9195
α, β, γ (°)90, 70.7739, 9085.2710, 85.5673, 63.073480.4854, 73.9054, 81.940290, 116.5216, 90
V3)????
Z????
Radiation type?, λ = ? Å?, λ = ? Å?, λ = ? Å?, λ = ? Å
µ (mm1)????
Crystal size (mm) × × × × × × × ×
Data collection
Diffractometer????
Absorption correction????
No. of measured, independent and
observed (?) reflections
?, ?, ? ?, ?, ? ?, ?, ? ?, ?, ?
Rint????
Refinement
R[F2 > 2σ(F2)], wR(F2), S ?, ?, ? ?, ?, ? ?, ?, ? ?, ?, ?
No. of reflections????
No. of parameters????
No. of restraints????
Δρmax, Δρmin (e Å3)?, ??, ??, ??, ?


(quin13)(quin14)(quin15)(quin02)
Crystal data
Chemical formula????
Mr????
Crystal system, space groupMonoclinic, P21/CMonoclinic, P21Monoclinic, P21/CMonoclinic, P21/C
Temperature (K)????
a, b, c (Å)5.7454, 31.2880, 3.93083.8370, 31.2284, 5.874015.9300, 7.0129, 31.140014.2551, 3.8236, 13.6709
α, β, γ (°)90, 90.2924, 9090, 89.1833, 9090, 11.8368, 9090, 77.3096, 90
V3)????
Z????
Radiation type?, λ = ? Å?, λ = ? Å?, λ = ? Å?, λ = ? Å
µ (mm1)????
Crystal size (mm) × × × × × × × ×
Data collection
Diffractometer????
Absorption correction????
No. of measured, independent and
observed (?) reflections
?, ?, ? ?, ?, ? ?, ?, ? ?, ?, ?
Rint????
Refinement
R[F2 > 2σ(F2)], wR(F2), S ?, ?, ? ?, ?, ? ?, ?, ? ?, ?, ?
No. of reflections????
No. of parameters????
No. of restraints????
Δρmax, Δρmin (e Å3)?, ??, ??, ??, ?


(quin03)(quin04)(quin05)(quin06)
Crystal data
Chemical formula????
Mr????
Crystal system, space groupMonoclinic, P21/CMonoclinic, P21/CMonoclinic, P21/CMonoclinic, P21
Temperature (K)????
a, b, c (Å)11.3480, 6.2795, 13.03135.7342, 3.9658, 31.09583.9882, 5.7197, 30.287828.1840, 4.9326, 4.9767
α, β, γ (°)90, 130.4828, 9090, 104.5126, 9090, 94.0361, 9090, 81.9358, 90
V3)????
Z????
Radiation type?, λ = ? Å?, λ = ? Å?, λ = ? Å?, λ = ? Å
µ (mm1)????
Crystal size (mm) × × × × × × × ×
Data collection
Diffractometer????
Absorption correction????
No. of measured, independent and
observed (?) reflections
?, ?, ? ?, ?, ? ?, ?, ? ?, ?, ?
Rint????
Refinement
R[F2 > 2σ(F2)], wR(F2), S ?, ?, ? ?, ?, ? ?, ?, ? ?, ?, ?
No. of reflections????
No. of parameters????
No. of restraints????
Δρmax, Δρmin (e Å3)?, ??, ??, ??, ?


(quin07)(quin08)(quin09)
Crystal data
Chemical formula???
Mr???
Crystal system, space groupMonoclinic, P21/CTriclinic, P1Triclinic, P1
Temperature (K)???
a, b, c (Å)5.0189, 4.8913, 28.08795.9241, 16.7841, 3.809015.7670, 3.8742, 5.8128
α, β, γ (°)90, 94.0515, 90109.9087, 91.5084, 100.575389.6689, 95.7095, 97.2323
V3)???
Z???
Radiation type?, λ = ? Å?, λ = ? Å?, λ = ? Å
µ (mm1)???
Crystal size (mm) × × × × × ×
Data collection
Diffractometer???
Absorption correction???
No. of measured, independent and
observed (?) reflections
?, ?, ? ?, ?, ? ?, ?, ?
Rint???
Refinement
R[F2 > 2σ(F2)], wR(F2), S ?, ?, ? ?, ?, ? ?, ?, ?
No. of reflections???
No. of parameters???
No. of restraints???
Δρmax, Δρmin (e Å3)?, ??, ??, ?

 

Acknowledgements

This work is supported by the Institute for the Promotion of Innovation by Science and Technology in Flanders (IWT).

References

First citationAccelrys (2002). Cerius2, version 4.71, User Guide. Accelrys Inc., San Diego, USA.  Google Scholar
First citationBailey, M. & Brown, C. J. (1967). Acta Cryst. 22, 387–391.  CSD CrossRef CAS IUCr Journals Web of Science Google Scholar
First citationBaur, W. H. & Kassner, D. (1992). Acta Cryst. B48, 356–369.  CrossRef CAS Web of Science IUCr Journals Google Scholar
First citationBrown, C. J. (1984). Acta Cryst. C40, 1762.  CrossRef Web of Science IUCr Journals Google Scholar
First citationDay, G. M. et al. (2005). Acta Cryst. B61, 511–527.  Web of Science CSD CrossRef CAS IUCr Journals Google Scholar
First citationDay, G. M. & Price, S. L. (2003). J. Am. Chem. Soc. 125, 16434–16443.  Web of Science CrossRef PubMed CAS Google Scholar
First citationDeuschel, W., Gundel, F., Wuest, H. & Daubach, E. (1963). US Patent 3074950.  Google Scholar
First citationEijck, B. P. van, Mooij, W. T. M. & Kroon, J. (2001). J. Comput. Chem. 22, 805–815.  Web of Science CrossRef Google Scholar
First citationErk, P. (2001). High Performance Pigments, edited by H. M. Smith, ch. 8, pp. 103–123. Weinheim: Wiley-VCH.  Google Scholar
First citationFilho, D. S. & Oliveira, C. M. F. (1992). J. Mater. Sci. 27, 5101–5107.  CrossRef CAS Web of Science Google Scholar
First citationFrisch, M. J. et al. (1998). GAUSSIAN94 (Revision D.4), Gaussian Inc., Pittsburgh, USA.  Google Scholar
First citationGavezzotti, A. (2002). J. Phys. Chem. B, 106, 4145–4154.  Web of Science CrossRef CAS Google Scholar
First citationGavezzotti, A. (2003). J. Phys. Chem. B, 107, 2344–2353.  Web of Science CrossRef CAS Google Scholar
First citationGdanitz, R. J., Karfunkel, H. R. & Leusen, F. J. J. (1993). J. Mol. Graph. 11, 275–276.  CrossRef Google Scholar
First citationJaffe, E. E. (2001). High Performance Pigments, edited by H. M. Smith, ch. 18, pp. 279–306. Weinheim: Wiley-VCH.  Google Scholar
First citationJones, F., Okui, N. & Patterson, D. (1975). J. Soc. Dyers Col. 91, 361–365.  CrossRef CAS Google Scholar
First citationKarfunkel, H. R., Leusen, F. J. J. & Gdanitz, R. J. (1994). J. Comput. Aided Mater. Des. 1, 177–185.  CrossRef CAS Google Scholar
First citationKarfunkel, H. R., Rohde, B., Leusen, F. J. J., Gdanitz, R. J. & Rihs, G. (1993). J. Comput. Chem. 14, 1125–1135.  CrossRef CAS Web of Science Google Scholar
First citationLeusen, F. J. J. (1996). J. Cryst. Growth, 166, 900–903.  CrossRef CAS Web of Science Google Scholar
First citationLincke, G. (2000). Dyes Pigments, 44, 101–122.  Web of Science CrossRef CAS Google Scholar
First citationLincke, G. & Finzel, H. U. (1996). Cryst. Res. Technol. 31, 441–452.  CrossRef CAS Web of Science Google Scholar
First citationManger, C. W. & Struve, W. S. (1957). US Patent 2844581.  Google Scholar
First citationMayo, S. L., Olafson, B. D. & Goddard, W. A. III (1990). J. Phys. Chem. 94, 8897–8909.  CrossRef CAS Web of Science Google Scholar
First citationMizuguchi, J., Sasaki, T. & Tojo, K. (2002). Z. Kristallogr. 217, 249–250.  Web of Science CrossRef CAS Google Scholar
First citationNishimura, N., Senju, T. & Mizuguchi, J. (2006). Acta Cryst. E62, o4683–o4685.  CSD CrossRef IUCr Journals Google Scholar
First citationPanina, N., Janssen, F. F. B. J., Deroover, G., Meekes, H. & Vlieg, E. (2007). In preparation.  Google Scholar
First citationPaulus, E. F., Dietz, E., Kroh, F., Prokschy, F. & Lincke, G. (1989). 12th European Crystallographic Meeting, published in Z. Kristallogr. (1990, Suppl. 2), pp. 23–24.  Google Scholar
First citationPayne, R. S., Roberts, R. J., Rowe, R. C. & Docherty, R. (1998). J. Comput. Chem. 19, 1–20.  CrossRef CAS Google Scholar
First citationPotts, G. D., Jones, W., Bullock, J. F., Andrews, S. J. & Maginn, S. J. (1994). 22, 2565–2566.  Google Scholar
First citationReidinger, A. D. & Struve, W. S. (1955). US Patent 2844484.  Google Scholar
First citationSchmidt, M. U., Ermrich, M. & Dinnebier, R. E. (2005). Acta Cryst. B61, 37–45.  Web of Science CSD CrossRef CAS IUCr Journals Google Scholar
First citationStruve, W. S. (1955). US Patent 2844485.  Google Scholar
First citationWiles, D. B. & Young, R. A. (1981). J. Appl. Cryst. 14, 149–151.  CrossRef CAS Web of Science IUCr Journals Google Scholar

© International Union of Crystallography. Prior permission is not required to reproduce short quotations, tables and figures from this article, provided the original authors and source are cited. For more information, click here.

Journal logoJOURNAL OF
APPLIED
CRYSTALLOGRAPHY
ISSN: 1600-5767
Follow J. Appl. Cryst.
Sign up for e-alerts
Follow J. Appl. Cryst. on Twitter
Follow us on facebook
Sign up for RSS feeds