Lattice defects in quinacridone

Dedicated to Professor Dr Bernd O. Kolbesen

Various lattice defects in the α^I-phase of quinacridone (C₂₀H₁₂N₂O₂) were simulated using lattice-energy minimizations. α^I-Quinacridone forms a chain structure in P1, Z = 1. The molecules are connected by hydrogen bonds along [010], by π-stacking along [100] and by weak van der Waals interactions along [001]. α^I-Quinacridone is inherently nanocrystalline. Lattice defects were calculated in correspondingly large supercells with up to 4464 atoms, using a previously evaluated force field. Vacancies, vacancy aggregates and interstitial molecules are energetically very unfavourable. A misorientation of a single molecule (flip around [010] by 180°) causes an energy increase of 243.7 kJ mol⁻¹. Various edge and screw dislocations were investigated. A screw dislocation along [010] causes an energy increase of ΔE = 38.0 kJ mol⁻¹ per molecule, all other line dislocations are even worse. In contrast, the rotation of an entire chain around the chain axis [010] by 180° leads to only a very small energy increase (ΔE = 1.6 kJ mol⁻¹) and the real crystals probably contain a high number of such defects. Various planar defects were calculated, including different stacking disorders and misfit-layer structures with two different types of layers having different lateral periodicities. Stacking faults along [001] with herringbone stacking instead of parallel stacking are energetically quite favourable (ΔE = 2.2 kJ mol⁻¹); the same is true for domains with misoriented molecules in the [001] direction. As an example for a bulk defect, domains are calculated in which blocks of 4 × 4 chains are rotated by 180° around [010], which leads to an energy increase of only 1.1 kJ mol⁻¹. Twinning by mirroring at the (001) plane is energetically favourable (ΔE = 0.9 kJ mol⁻¹). This twinning was observed in an HRTEM image. It is probable that the crystallites also contain rotations of chains, layers or blocks around [010] by 180°, but these defects cause only a very slight modification of the molecular packing, which was not observable in the HRTEM image. The lattice defects in α^I-quinacridone investigated here provide an insight into lattice defects, their energies and local structures. Similar lattice defects are expected to occur also in other similar organic chain structures.

1. Introduction

Quinacridone (Fig. 1) is an organic semiconductor. Organic semiconductors are used for electronic devices such as organic light-emitting diodes, organic solar cells (organic photovoltaic) and organic field-effect transistors. The electronic properties of semiconductors strongly depend on the type and frequency of lattice defects, such as dislocation lines or impurities. The nature, local structure and frequency of lattice defects in inorganic semiconductors, e.g. silicon, has been widely studied for many decades (see e.g. Kolbesen & Cerva, 1999; Kolbesen, 2014). For organic semiconductors, much less work has been reported, hitherto (see e.g. Sherwood, 1969; Desvergne et al., 1974; Mokichev & Pakhomov, 1982; Ide et al., 1993; Cuppen et al., 2004; Kang et al., 2005; Walker et al., 2005; Chapman et al., 2006; Wu et al., 2017; Li et al., 2022).

Figure 1 Chemical diagram for quinacridone.

Typical lattice defects in inorganic solids are: point defects, e.g. vacancies, interstitial atoms, or dopant atoms; line defects, e.g. edge or screw dislocations; planar defects such as stacking faults or small-angle grain boundaries; and bulk defects such as voids. All these defects may occur in organic crystals, too (see e.g. Scheffen-Lauenroth et al., 1981; Olson et al., 2018). Additionally, there can be orientational defects (or orientational disorder) of side groups, or of entire molecules, e.g. head-to-tail disorder. The orientational defect may be located in single molecules, or can proliferate to chains or layers of molecules in an incorrect orientation. All these lattice defects were investigated in this work. Additionally, organic crystals can incorporate foreign molecules such as impurities, byproducts, decomposition products, water and solvent molecules which will be the subjects of further work.

As a model compound, we chose the α^I-phase of quinacridone. Quinacridone (Fig. 1) is an industrial organic pigment used for the colouration of paints, coatings, plastics and printing inks (Hunger & Schmidt, 2018). In the International Colour Index, it is registered as Pigment Violet 19 (Colour Index International, 2022; Abel, 1998). There are four polymorphic forms. The α^I- and α^II-phases result from the synthesis of the compound, depending on the synthetic conditions. Heating the α^I-phase in NaOH yields the β-phase, whereas heating in organic solvents, e.g. iso­butanol, results in the γ-phase. The crystal structure of the α^I-phase was determined by crystal structure prediction followed by Rietveld refinement (Leusen, 1994, 1996; Paulus et al., 2007). The structures of the β- and γ-phases were determined by single-crystal X-ray diffraction (Paulus et al., 1989; Potts et al., 1994; Mizuguchi, Sasaki & Tojo, 2002; Nishimura et al., 2006).

The α^II-phase (at that time there was no distinction between the phases α^I and α^II; both were called α-phase) was investigated by Lincke & Finzel (1996), using a limited-quality powder pattern of the α^II-phase. They constructed a structure in with two independent molecules per unit cell, both located on inversion centres. The molecules were assumed to form a criss-cross pattern like in γ-quinacridone. This packing was manually fitted to the powder pattern (Lincke & Finzel, 1996). However, the final fit was not convincing, and the structure remained questionable (Paulus et al., 2007). Finally, the structure of the α^II-phase was determined by electron diffraction (Gorelik et al., 2016), revealing a chain structure and, thereby, disproving the structural model of Lincke & Finzel (1996).

The α^I-phase crystallizes in with Z = 1, all other phases crystallize in P2₁/c, Z = 2.

The quinacridone molecules are always planar or close to planar. In all phases the molecules are situated on crystallographic inversion centres. In the α^I-, α^II- and β-phases, each molecule is connected to two neighbouring molecules by two hydrogen bonds each, resulting in chains. In the α^I-phase all chains are parallel (Fig. 2). The chains are not fully planar, but exhibit steps of approximately 0.6 Å between the molecules [Fig. 2(c)]. In the α^II-phase the chains form, in principle, a herringbone packing; but actually the α^II-phase exhibits a stacking disorder with a mixture of herringbone and parallel arrangements of neighbouring sheets, as revealed by electron diffraction and lattice-energy minimizations [Fig. 3(a); Gorelik et al., 2016]. In the β-phase the chains extend in two different directions, [110] and ], which form an angle of 69° [Fig. 3(b)]. In the three phases, the chains are stacked on top of each other, resulting in layers, but the stacking of the layers is different. In the α^I-phase, neighbouring layers are translationally equivalent, in the β-phase they are rotated by 69° and in the α^II-phase by 180°. Correspondingly, these three structures belong to the same polytype family. In the γ-phase, each molecule is connected to four neighbouring molecules, resulting in a criss-cross pattern [Fig. 3(c)]. In all phases, the molecules are densely packed. The densities are high (α^I = 1493 kg m⁻³, α^II = 1414 kg m⁻³, β = 1540 kg m⁻³, γ = 1470 kg m⁻³). The lattice energies are high, and all phases are fully insoluble in water and all solvents at ambient conditions. Quinacridone is thermally stable in air to temperatures of about 600–700 K, when it starts to sublime.

Figure 2 Crystal structure of α^I-quinacridone (without lattice defects). View directions (a) [100], (b) [010], (c) [001]. Colour code for all figures: C grey, H white, N blue, O red. Hydrogen bonds are shown as dashed lines. Note that in this work a different unit-cell setting is used than in earlier work (for details, see Section 2.1). The structure after unit-cell transformation and lattice-energy optimization is shown.

Figure 3 Crystal structures of other polymorphs of quinacridone: (a) α^II-phase (with stacking disorder), (b) β-phase and (c) γ-phase.

The crystal structures of quinacridone have a strong effect on the electronic structure, as is evident from the different colours of the individual polymorphs. The α^I-phase has a dull reddish-violet shade, the α^II-phase is red, the β-phase is reddish violet and the γ-phase is red, too (but slightly more yellowish than α^II-quinacridone). In solution, quinacridone is yellow to orange. A similar colour shift between solution and solid state is also observed for several other organic pigments used in car coatings, including diketo­pyrrolo­pyrrole or perylene pigments. In all these cases, the observed shift from a yellow to orange colour in solution to a red or violet in the solid state is caused by the stacking of the molecules, resulting in a classic blue-shifted H-aggregate. Thus, the red colour of most red cars is only a solid-state effect; without the intermolecular interactions of the molecules in the crystals, the cars would be yellow or orange.

The α^I-phase of quinacridone presents itself as an exemplary model system for the investigation of lattice defects, because the molecule is rigid and the crystal structure is quite simple. The unit cell contains only one molecule, which is situated on a crystallographic inversion centre. Similar crystal structures in , Z = 1, also consisting of planar molecules connected by double hydrogen bonds into parallel chains, have also been reported for other organic semiconductors, e.g. 2,9-di­methyl­quinacridone (Pigment Red 122; Mizuguchi, Senju & Sakai, 2002) and diketo­pyrrolo­pyrrole (Pigment Red 255; Mizuguchi et al., 1992).

There is a further advantage to using the α^I-phase as a model system. The molecules are connected in the a direction by π-stacking, in the b direction by hydrogen bonds and in the c direction by van der Waals interactions between C and H atoms. Hence, the effect of the different intermolecular interactions on the lattice defects can clearly be distinguished.

Industrially, the α^I-phase is obtained as follows. The starting material, 2,5-bis­(phenylamino)terephthalic acid is heated in molten polyphospho­ric acid to 100–150°C. The resulting dark-blue mixture of (probably protonated) quinacridone in hot polyphospho­ric acid is poured into ice-water. The α^I-quinacridone immediately precipitates as a nanocrystalline powder. The α^I-phase is always nanocrystalline. Recrystallization is not possible because of the insolubility in all media. Treatment with solvents in suspension at elevated temperature either does not improve the crystallinity, or results in a phase transformation into the more stable β- or γ-phases. The best powder which we ever obtained had a domain size of 20 nm, as determined by Rietveld refinement. Samples from industrial production may be close to amorphous. Correspondingly, the α^I-phase certainly contains a high number of lattice defects. However, nothing is known of the nature of these defects, their local structure or their frequency. Experimental investigations on the defects in this nanocrystalline powder would be difficult. In contrast, the simulation of the lattice defects by lattice-energy minimizations is feasible.

For the simulation of the lattice defects, correspondingly large superstructures are constructed, e.g. 5 × 5 × 5 unit cells for a point defect such as a vacancy. Such large supercells at present prevent the minimization of the lattice energy with quantum-mechanical methods. Hence, we used force-field methods, which were already successfully employed for the calculation of the crystal structures of α^I-quinacridone (Leusen, 1994, 1996) and 2-methyl-quinacridone (Schlesinger et al., 2020).

In this paper, the applied force field is evaluated first. Subsequently, various types of point defects, line defects, planar defects and bulk defects of α^I-quinacridone are investigated. Finally, we compare the simulation results with an experimental HRTEM image of α^I-quinacridone.

2. Calculation details and experimental details

2.1. Force fields

Two different force fields were tested: the Dreiding force field (Mayo et al., 1990) and a tailor-made force field fitted to lattice-energy calculations by dispersion-corrected density-functional theory (DFT-D), as described by Neumann et al. (2008). For the atomic charges, various approaches were tested: (1) the Gasteiger method (Gasteiger & Marsilli, 1980); (2) the charge-equilibrium method as implemented in Materials Studio (Accelrys, 2008); (3) charges calculated by the electrostatic potential (ESP) method, based on HF/6-31G** calculations and (4) charges associated with the tailor-made force-field. Different algorithms for the summation of the electrostatic energy were employed: atom-based, group-based and Ewald summation.

Finally, all lattice defects were calculated with the Dreiding force field, ESP charges and Ewald summation for Coulomb and van der Waals interactions.

DFT-D calculations of the four polymorphs of quinacridone without lattice defects were performed with the program CASTEP (Clark et al., 2005), with the PBE functional (Perdew et al., 1996) and a dispersion correction of Grimme (2006). The unit-cell parameters were optimized together with all atomic coordinates.

2.2. Programs

The best program for crystal modelling studies on organic compounds with force fields was, in our opinion, Cerius² (Accelrys, 2005). However, Cerius² requires an SGI workstation with IRIX operating system, which is outdated. Thus, we used its successor program, Materials Studio (Version 4.4, Accelrys, 2008), which runs on a Windows PC. Drawings were created with Mercury (Macrae et al., 2020) and SCHAKAL (Keller, 1999).

2.3. Model building and optimization

α^I-Quinacridone crystallizes in , Z = 1, with unit-cell parameters at room temperature of a = 3.8017 (15), b = 6.612 (3), c = 14.485 (6) Å, α = 100.68 (8), β = 94.40 (6), γ = 102.11 (5)°, V = 346.67 (11) ų (Paulus et al., 2007). For the simulation of lattice defects, the unit cell was transformed with a′ = −a, b′ = a + b, c′ = −b − c and the molecules were shifted by (½, 0, ½). In the resulting unit-cell setting, the a axis corresponds to the thickness of the molecule, the b axis to the width and the c axis to the length of the molecule. The molecule is situated on a crystallographic inversion centre at (½, ½, ½) in the centre of the unit cell. This setting facilitates the construction, evaluation and description of the various lattice defects.

For the construction of the lattice defects, correspondingly large supercells were set up. For example, for point defects such as vacancies we used a 3 × 3 × 3 or 5 × 5 × 5 supercell, for line defects such as edge or screw dislocations in the c direction a 9 × 6 × 1 supercell, and for planar defects in (100) direction a 8 × 1 × 1 supercell. Subsequently, the molecules were manually removed, rotated or shifted, to obtain a sensible starting point for the following lattice-energy minimizations. In complicated cases, such as screw or edge dislocations, several different models were set up and optimized, in order to find the energetically most favourable local structure for a given defect type. Screw and edge dislocations as well as twinning violate 3D periodicity of the crystal. To counteract this effect, we always employed them pairwise.

Exemplary manually constructed starting structures along with their optimized versions are shown in Sections 4.1, 5.2.3, 5.3.1 and 5.3.2.

For the optimization, the convergence criteria in Materials Studio were set to `ultrafine'. Calculation times were in the order of 2 min to 2 h on a standard PC.

For the calculation of point and line defects, the unit-cell parameters were kept fixed to the values of the optimized, undisturbed structure. This procedure reflects the assumption that the overall unit-cell parameters of a crystal do not change upon the occurrence of a point or line defect. Additional calculations were performed with variable unit-cell parameters. If the supercell is not too small, these calculations yielded very similar structures and energies as those with fixed unit-cell parameters. In contrast, calculations of planar defects such as stacking faults require the corresponding unit-cell parameters to be optimized. For example, a planar defect in the (100) plane requires a, β and γ to be optimized, whereas α, b and c can be fixed. Volume defects of limited size were treated like point and line defects, i.e. with fixed unit-cell parameters.

The optimized, undisturbed structure of α^I-quinacridone has a total force-field energy of −309.264 kcal mol⁻¹ = −1294.52 kJ mol⁻¹. Test calculations proved that the structure and energy are independent of the unit-cell setting. This energy value served as reference for all lattice defects. The relative energies of the various lattice defects described in this paper are given in kJ mol⁻¹ per molecule per lattice defect (1 kcal = 4.184 kJ), unless stated otherwise. For example, the energy of a vacancy in a 5 × 5 × 5 supercell is calculated as the difference between the total energy of the optimized supercell and 5³ − 1 times the lattice energy of the optimized initial cell with Z = 1. For line defects and planar defects, the energy is given per unit-cell lengths of the undisturbed structure in the corresponding directions, i.e. per molecule of the lattice defect.

For most lattice defects, the lattice energies of the optimized structures depend only slightly on the size of the unit cell. For example, the energy of a vacancy is 231.69 kJ mol⁻¹ in a 3 × 3 × 3 supercell and 233.37 kJ mol⁻¹ in a 5 × 5 × 5 supercell. The energy of the 3 × 3 × 3 supercell model is reduced to 230.66 kJ mol⁻¹, if the unit-cell parameters are optimized too.

2.4. Pair-distribution function analyses

X-ray powder diffraction data were recorded at 6ID-D MUCAT beamline at the Advanced Photon Source (APS) at Argonne National Laboratory, using a wavelength of 0.1428 Å. The sample was prepared on Kapton tapes and measured at 100 K. The pair distribution function (PDF) was calculated with PDFgetX2 (Qiu et al., 2004), using a Q_max of 21 Å⁻¹.

3. Force-field validation

The force field was evaluated by lattice-energy minimization of all four polymorphs of quinacridone, using different force fields, different atomic charges and different methods for the summation of the electrostatic energy. The best results were obtained for the Dreiding force field in combination with ab initio ESP charges and Ewald summation. Upon optimization, the crystal structures changed only slightly (Tables 1 and S2 in supporting information). The resulting lattice energies agree with the experimental observation that the β- and γ-phases are thermodynamically more stable than the α-phases (see Table S3). (Kinetically, all phases are stable at ambient conditions for decades.) This force field was used for the calculation of all lattice defects.

Table 1
Evaluation of the force field

Lattice-energy optimizations on α^I-quinacridone without lattice defects. `Dreiding' denotes the Dreiding force field with 6-31G**-ESP charges and Ewald summation.

	Experimental†	Dreiding	DFT-D
a (Å)	3.8017 (15)	3.9241	3.685
b (Å)	6.901 (3)	6.8938	6.386
c (Å)	14.766 (6)	14.9598	14.574
α (°)	99.30 (8)	98.999	101.486
β (°)	99.74 (6)	100.601	97.861
γ (°)	110.48 (6)	114.976	100.066
V (Å^3)	346.67 (11)	347.917	325.646

†After unit-cell transformation.

Interestingly, exactly this combination of a Dreiding force field, 6-31G**-ESP charges and Ewald summation was already successfully used in 1994 for the crystal structure prediction of quinacridone (Leusen, 1994, 1996), which led to the crystal structure determination of α^I-quinacridone.

The optimized crystal structure of α^I-quinacridone is shown in Fig. 2. The unit-cell parameters are given in Table 1. Henceforth, these unit-cell parameters are denoted as a₀, b₀, c₀, α₀, β₀ and γ₀. This structure served as the basis for the construction and lattice-energy minimizations of all lattice defects.

4. Results and discussion: lattice defects: point defects

4.1. Vacancies

Vacancies are frequently found in metals, e.g. silicon, or in ionic solids, e.g. metal halides, where they show up as missing atoms or ions. In molecular crystals, a vacancy corresponds to a missing entire molecule.

For α^I-quinacridone, vacancies were calculated with a supercell of 3 × 3 × 3 or 5 × 5 × 5 unit cells. The latter one contains 5³ − 1 = 124 molecules with 4464 atoms. The resulting optimized structure is shown in Figs. 4(a) and 4(b). Astonishingly, the molecules in the neighbourhood of the vacancy show no tendency to move towards the empty space. The molecular positions change by less than 0.1 Å. The reason might be that all neighbouring molecules are fixed to their neighbours by hydrogen bonds and π-stacking. The molecular packing of α^I-quinacridone is very space-efficient, and very favourable in terms of intermolecular energies. Hence, every molecular movement by more than 0.1 Å would cause these intermolecular interactions to deteriorate.

Figure 4 Vacancy. (a) Superstructure with 3 × 3 × 3 unit cells after optimization. Perspective view, approximately along [010]. The vacancy is marked by x. (b) Superstructure with 5 × 5 × 5 unit cells, containing 4464 atoms, after optimization. View direction approximately [010]. (c) Starting structure with a molecule in the centre of the vacancy, before optimization. In (c) the layer with the vacancy is shown in dark colours, the other layers are light grey.

The vacancy leads to a very high energy increase of 233.4 kJ mol⁻¹. This energy increase is caused by the loss of two hydrogen bonds (25.5 kJ mol⁻¹), but even more by the loss of the intermolecular van der Waals energy (154.2 kJ mol⁻¹) and the electrostatic energy (66.4 kJ mol⁻¹).

We also investigated the instance in which one of the neighbouring molecules is positioned in the centre of the vacancy, so that the H atoms of all NH groups are involved in hydrogen bonds, some of them being bifurcated [see Fig. 4(c)]. However, their geometry strongly deviates from usual hydrogen-bond geometries. Apparently, this arrangement is energetically unfavourable: upon optimization the central molecule moves back to its original position, and the refinement converges to the situation with a full vacancy in an otherwise undisturbed structure as shown in Fig. 4(a).

The high energy reveals that such an empty vacancy is very unlikely. Under real experimental conditions, the vacancies are probably filled with other molecules, such as water, solvents or byproduct molecules.

4.2. Vacancy aggregates

Vacancy aggregates, e.g. two of three missing neighbouring molecules, are energetically very unfavourable, too.

4.3. Interstitial molecules

α^I-Quinacridone has an efficient molecular packing with a density of 1493 kg m⁻³. Any attempt to squeeze an additional molecule into this structure leads to strong local distortions and a high energy increase. An example with an energy increase of about 670 kJ mol⁻¹ is shown in Fig. 5.

Figure 5 Attempt to calculate a structure with one interstitial molecule. This figure and all following figures show the structures after optimization. The additional molecule leads to a complete distortion of the whole structure in its vicinity.

4.4. Orientation faults

A rotation of a molecule by 180° around its long molecular axis or its medium molecular axis results in a mutual exchange of CO and NH groups. Such an orientational `head-to-tail disorder' of a single molecule is sterically quite possible, but results in energetically unfavourable C=O···O=C and N—H···H—N contacts. The misoriented molecule evades these contacts by rotating out of the plane by 11.80°, see Fig. 6. The energy is +243.7 kJ mol⁻¹.

Figure 6 Orientation fault of a single molecule.

4.5. Combination of orientation fault and vacancy

In order to provide more space for the misoriented molecule, we tested the combination of an orientation fault with a neighbouring vacancy. The misoriented molecule was placed in the centre of the available space and the structure was optimized. The resulting structure is shown in Fig. 7. The energy of ΔE = 366.4 kJ mol⁻¹ is lower than the sum of the structures with a vacancy and with a misoriented molecule, but still higher than for a misoriented molecule without a vacancy. Hence, a misoriented single molecule will be extremely rare in a real crystal.

Figure 7 Combination of orientation fault and vacancy. The central layer is shown as bold sticks, the other layers in wireframe style.

5. Line defects

Line defects can occur in all spatial directions. For α^I-quin­acridone the principal directions are:

(a) [100]: π-stacking of the molecules,

(b) [010]: hydrogen bonds,

(c) [001]: weak van der Waals interactions between the ends of the molecules.

Edge and screw dislocations cause a violation of the translational periodicity in their vicinity. Thus, a single edge or screw dislocation cannot be handled with a three-dimensionally periodic supercell. Therefore, all our models contain a pair of edge or screw dislocations so that the sum of the two Burgers vectors vanishes (see Fig. 8). The resulting pair of parallel dislocation lines can be regarded as a section of a dislocation loop. This approach spawns the problem that the local structure and the energy of edge and screw dislocations depend not only on the direction of the dislocation line and the Burgers vector, but also on the relative position and distance of the two dislocation lines in the supercell. However, a combination of two dislocation lines may well be present in a real crystal.

Figure 8 Edge dislocations along [100] with Burgers vectors (0,1,0) and . The two dislocation lines run in the view direction of [100] and are marked by S. The two `circuits around the dislocation line' are shown in black; the Burgers vectors are shown in blue.

Various models for edge and screw dislocations were tested. In the following we focus on those models which remained chemically sensible after optimization.

The total energy of a line defect depends on the length of the dislocation line in the crystal. In our calculations, the dislocation line is assumed to be infinite, and the models are translationally periodic in the direction of the dislocation line. In the following, the energies are given per dislocation line per translational period, i.e. per unit-cell parameter of the unit cell in the corresponding direction, e.g. per 6.386 Å in the [010] direction.

Apart from edge and screw dislocations, two other types of line dislocations were investigated: missing molecular chains and chains with incorrect molecular orientations.

5.1. Line vacancy along [010] (missing molecular chain)

In the [010] direction, the molecules are connected by hydrogen bonds. Removing one entire chain along [010] does not break any hydrogen bonds, but the van der Waals interactions to molecules in neighbouring chains are lost. The neighbouring chains only show a slight tendency to fill the empty space (see Fig. 9). The energy increases by 145.0 kJ mol⁻¹ per missing molecule. As for a vacancy and a point defect, in a real crystal, the empty space would be filled by other molecules.

Figure 9 Line vacancy after optimization with fixed unit-cell parameters.

Line vacancies along other spatial directions would cause the breaking of many hydrogen bonds and are energetically very unfavourable too.

5.2. Edge dislocations and dislocation loops

In edge dislocations, the Burgers vector is (nearly) perpendicular to the direction of the dislocation line. We investigated dislocations along [100], [010] and [001], with Burgers vectors (0,1,0), (0,0,1) and (1,0,0).

5.2.1. Dislocation line parallel to [100] and Burgers vector (0,1,0)

A model with two edge dislocations, amounting to a dislocation loop, was assembled by constructing a supercell of 1 × 3 × 8 unit cells and removing four consecutive molecules along [001]. The remaining molecules show no tendency to move into the vacancy, because they are held in place by hydrogen bonds (see Fig. 8). As such, a relatively large vacancy remains, which is energetically unfavourable at 295.0 kJ mol⁻¹ per dislocation line.

5.2.2. Dislocation line parallel to [100] and Burgers vector (0,0,1)

The dislocation line was generated in a similar fashion as the one before, by construction of a supercell of 1 × 6 × 3 unit cells and the removal of three molecules along [010] (see Fig. 10). As before, a large vacancy persists because the remaining molecules are fixed in place by hydrogen bonds, resulting in an energy of 104.1 kJ mol⁻¹.

Figure 10 Edge dislocation along [100] with the Burgers vector (0,0,1). View direction [100]. The positions of the dislocation lines are marked by S.

5.2.3. Dislocation line parallel to [010] and Burgers vector (1,0,0)

This model was constructed with a supercell of 8 × 1 × 8 unit cells and the removal of four neighbouring molecules parallel to [101] [Fig. 11(a)]. During the optimization, the molecules show a clear tendency to fill the vacancy by rotations of the entire molecules. The rotations start with the molecules in the layers close to the dislocation line, and continues with the layers in between. Fig. 11(b) shows the situation without reaching full convergence. The energy of this structure is by 243.8 kJ mol⁻¹ per dislocation line higher than the energy of the undisturbed quinacridone; hence this structure is energetically very unfavourable.

Figure 11 Edge dislocation along [010] and Burgers vector (1,0,0). (a) Starting structure and (b) structure after optimization. The vacancy is partially filled upon the rotation of the molecules. The positions of the dislocation lines are marked by S

5.2.4. Dislocation line parallel to [010] and Burgers vector (0,0,1)

The dislocation loop was constructed in a supercell of 10 × 1 × 3 unit cells by the removal of four neighbouring molecules in the [100] direction (see Fig. 12). The molecules barely move into the vacancy, making this model energetically unfavourable with an energy of 115.07 kJ mol⁻¹ per dislocation line.

Figure 12 Edge dislocation along [010], Burgers vector (0,0,1).

5.2.5. Dislocation line parallel to [001] and Burgers vector (1,0,0)

For this case, several models were constructed, which differ in the relative position of the two dislocation lines and the number of withdrawn molecules.

In Model 1, the dislocation lines are separated by 4a₀ − b₀, and one molecule was removed [see Fig. 13(a)]. In the optimized structure the hydrogen bonds are disrupted at the dislocation line, otherwise they are maintained through rotation and translation of the molecules. Although this structure looks chemically reasonable, the energy increase is high (210.8 kJ mol⁻¹). Apparently, the distortion of the structure causes a significant increase of the intermolecular energy. The removal of a second molecule increases the energy further by 245.5 kJ mol⁻¹.

Figure 13 Models for an edge dislocation along [001] with a Burgers vector (1,0,0): (a) Model 1, (b) Model 2 and (c) Model 3. The positions of the dislocation lines are marked by S.

In Model 2, the dislocation lines are separated by 2b₀, and two molecules were removed. In the optimization, this void is filled by the neighbouring molecules [Fig. 13(b)], which is in contrast to the maintenance of the structure observed with a single vacancy [Fig. 4(c)]. Some of the hydrogen bonds become distorted. The energy increase is 156.9 kJ mol⁻¹ per dislocation line and per missing molecule, which is slightly less than for a single vacancy (233.4 kJ mol⁻¹).

Model 3 exhibits a similar behaviour. Here, the dislocation lines are separated by 4b₀ with the removal of four molecules. The optimized structure [Fig. 13(c)] is a good impression of the actual local structure at an edge dislocation. The energy is 161.3 kJ mol⁻¹ per dislocation line, comparable to Model 2.

5.2.6. Dislocation line parallel to [001], Burgers vector (0,1,0)

Akin to a dislocation line parallel to [100] with a Burgers vector (0,1,0) (see Section 5.2.1), this lattice defect leads to a disruption of many hydrogen bonds and is energetically very unfavourable.

5.3. Screw dislocations

In screw dislocations, the Burgers vector is parallel to the dislocation line.

5.3.1. Dislocation line along [100]

The investigated structure is based on a 1 × 4 × 8 supercell with two dislocation lines with the Burgers vectors (1,0,0) and . The slip plane is (001), see Fig. 14. The starting model was constructed with chemically reliable chains having reliable hydrogen bond geometries [see Fig. 14(c) for chain No. 1]. During the optimization, the molecules shifted along [100], resulting in a structure with strongly distorted hydrogen bonds [see Fig. 14(b)]. Due to the distortion, the energy is 61.6 kJ mol⁻¹ per dislocation line, although no hydrogen bonds are broken.

Figure 14 Screw dislocation with the dislocation line along [100] with supercell 1 × 4 × 8. (a) View direction [100]. The positions of the dislocation lines are marked by . The slip planes are drawn in red. The numbers 1–8 denote the chain numbers. (b) Local structure at the chains 8/1 and 4/5. View direction (which corresponds to in the original cell). (c) Starting structure of chain 1.

5.3.2. Dislocation line along [010]

A screw dislocation along [010] does not break any hydrogen bonds, but entire chains of molecules are shifted in the chain direction. Our model has a slip plane parallel to (001) (see Fig. 15). The starting model consisted of an ensemble with various shifts in the y-direction (0.1, 0.3, 0.5, 0.7, 0.9) [see Fig. 15(a)]. In the optimized structures, all molecules moved to positions close to y = 0 or 0.5 [see Fig. 15(b)]. Apparently, these positions are energetically preferred.

Figure 15 Screw dislocation with the dislocation line along [010]. The numbers denote the y-coordinate of the molecular centres. The positions of the dislocation lines are marked by S. The slip planes are drawn in red. View direction [010]. (a) Starting structure, (b) optimized structure.

Apart from the longitudinal translation of the chains, the structure is barely distorted. Accordingly, the lattice energy is quite favourable with ΔE = 38.0 kJ mol⁻¹, because the hydrogen bonds are not broken and the van der Waals energy as well as the dense packing with π–π interactions are maintained.

5.3.3. Dislocation line along [001]

The [001]-direction corresponds to the long molecular axis. A screw dislocation along [001] causes a translation of the molecules along their long axis leading to the loss of hydrogen bonds, and unfavourable intermolecular contacts. Hence screw dislocations in this direction would be energetically very unfavourable.

5.4. Misorientation of an entire chain

5.4.1. Chain direction [010]

In [010] direction, the molecules are connected by hydrogen bonds into a chain. The rotation of all molecules in the chain by 180° around their medium molecular axis, which is almost parallel to [010], corresponds to a rotation of an entire chain around [010] by 180° (see Fig. 16). This rotation leads to a mutual exchange of CO and NH groups of all molecules in the chain. The same structure is formed by mirroring of the chain at a mirror plane perpendicular to [010]. These operations do not break any hydrogen bonds. Furthermore, the terminal benzene rings of the molecule do not move, i.e. the outer shape of the chain does not change. The chain itself has rod group (rod group No. 2; Kopský & Litvin, 2002). In contrast, the outer shape of the chain is close to having symmetry, which is a non-standard setting of (rod group No. 6).

Figure 16 Supercell with misoriented chain after optimization. View along b*. The misoriented chain is drawn in bold.

Correspondingly, the van der Waals interactions between the benzene rings at the end of the molecules is in the c direction do not change. However, the interaction with neighbouring molecules in the a direction (π stacking) is modified because the van der Waals interactions as well as the electrostatic interactions change. Surprisingly, this reorientation leads to an energy increase of only 1.57 kJ mol⁻¹. Apart from the misorientation of the chains, the crystal structure is not distorted at all. The low energy indicates that this lattice defect should frequently occur in the real crystal. Correspondingly, the crystals may contain different arrangements of the same rod. If these arrangements were periodic, they would be rod polytypes. However, the distribution of misoriented chains is not periodical, but almost statistical.

5.4.2. Other directions

The misorientation of a row of molecules in any other direction leads to unfavourable C=O⋯O=C and N—H⋯H—N contacts, which are energetically very unfavourable as shown in Section 4.4.

5.5. Misorientation of two neighbouring chains

As seen in Section 5.4.1, the misorientation of an entire chain in the [010]-direction causes an energy increase of only 1.57 kJ mol⁻¹. This prompted us to investigate if an ensemble of two neighbouring misoriented chains might be even more favourable. Each chain is surrounded by six neighbouring chains with mutual translation vectors of (1,0,0), (0,0,1), , , and . Correspondingly, there are three possibilities to group two neighbouring chains, which are described below.

5.5.1. Neighbours with translation of (1,0,0)

The simultaneous rotation of two neighbouring chains associated by a translation of (1,0,0) leads to an energy of 1.33 kJ mol⁻¹ per chain, which is even slightly more favourable than the rotation of two separate chains (Fig. 17).

Figure 17 Two misoriented chains associated by a translation of (1,0,0). The orange box highlights the misoriented chains. View along b*.

5.5.2. Neighbours with translation of (0,0,1)

Two chains with a mutual shift of (0,0,1) are connected only by a single van der Waals contact of two CH groups (see Fig. 18). Correspondingly, the energy (1.59 kJ mol⁻¹ per chain) is about the same as for separate chains.

Figure 18 Two misoriented chains associated by a translation of (0,0,1). The misoriented chains are framed in orange. View along b*.

5.5.3. Neighbours with translation of (1,0,1)

As in the previous model, the two chains are connected only by a single van der Waals contact of two CH groups (see Fig. 19). Surprisingly, this arrangement of two misoriented chains is even better than that of two chains linked by π–π stacking (Section 5.5.1). The energy is 1.28 kJ mol⁻¹.

Figure 19 Two misoriented chains associated by a translation of . The misoriented molecules are framed in orange. View along b*.

Apparently, the interactions between the terminal CH groups play a non-negligible role. The low energies of all sets of neighbouring chains indicate that a real crystal may also contain larger sets of misoriented chains, e.g. misoriented layers or volume regions, which will be described in Sections 6.2 and 7.2.

6. Planar defects

A simple chain structure can, in principle, exhibit various types of planar defects. We will focus on a few interesting ones.

6.1. Small-angle grain boundaries

A small-angle grain boundary interrupts the 3D translational periodicity of the crystal. For a calculation of a small-angle grain boundary with a 3D periodic model, the model must contain at least two small-angle grain boundaries, so that the sum of the Burgers vectors vanishes. A small-angle grain boundary can be regarded as a series of parallel edge dislocations. Actually, a periodic series of parallel edge dislocations have already been investigated in Section 5.2. The resulting energies were very high, which is an argument against the occurrence of small-angle grain boundaries.

6.2. Misorientation of a layer of molecules

As seen in Section 5.4.1, the rotation of an entire chain of molecules around [010] by 180° leads to an energy increase of only 1.57 kJ mol⁻¹. Similarly, all molecules of an entire layer can be rotated around [010] by 180°, leading to a layer with inverted molecular orientation. The resulting model can be regarded as a special case of a stacking fault, in which the molecules are rotated, but the molecular centres do not move. Hence the translational periodicity of the molecular centres is maintained.

It should be noted, that the rotation of molecules by 180° has a different effect than a twinning by reflection at a net plane, because the structure is triclinic, and a mirroring at a net plane would change the molecular orientations and positions in a different way, see Sections 6.3 and 8.

In this section, we describe the misorientation of molecules in a single layer. Models with a misorientation of a block of molecules in neighbouring layers are described in Section 7.2.

6.2.1. (100) layer

An inverted orientation of molecules in a layer parallel to (100) (see Fig. 20) leads to an energy increase of only 1.46 kJ mol⁻¹ per molecule, which is similar to the instance with a single misoriented chain.

Figure 20 Layer of misoriented molecules, parallel to (100). View along b*. The misoriented layer is shown in bold.

6.2.2. (101) layer

As seen in Section 5.5, the chains are not only linked in the [001] direction, but also in the direction, which corresponds to a layer parallel to (101). The model for a structure with a (101) layer with inverted molecular orientation was constructed in the following way: the unit cell was transformed with . The (101) layer of the original cell corresponds to the (100) layer of the transformed cell. Subsequently, a 5 × 1 × 1 supercell was assembled and one molecule was rotated (see Fig. 21). After optimization with free unit-cell parameters, the energy was only 0.78 kJ mol⁻¹ per molecule higher than that of an undisturbed crystal. The unit-cell parameters only changed by less than 0.17 Å and 1.1°.

Figure 21 Layer parallel to (101) with inverted orientation, calculated as a (100) layer in a unit cell with = c0 − a0. View along b*.

6.2.3. (001) layer

The rotation of molecules within a layer parallel to (001) (see Fig. 22) leads to an energy increase of only 1.34 kJ mol⁻¹ per molecule, similar to a rotation in a (100) layer.

Figure 22 Layer parallel to (001) with inverted molecular orientation. View direction [100]. The central layer consists of molecules with inverted orientation.

6.3. Other stacking faults

6.3.1. Herringbone instead of parallel stacking

A mixture of herringbone and parallel stacking of layers along (001) was experimentally observed in α^II-quinacridone by electron diffraction. The α^II-phase exhibits predominantly herringbone stacking with a minor contribution of parallel stacking [see Fig. 3(a)]. In contrast, the α^I-phase shows mainly parallel stacking. However, stacking faults with herringbone stacking cannot be ruled out. A corresponding model is shown in Fig. 23. The herringbone fragment corresponds to the structure of the α^II-phase, except for a misorientation of the central chain.

Figure 23 Structure of a stacking fault with local herringbone stacking. (a) View direction [010]. Herringbone contacts are marked with an H, parallel contacts with a P. (b) View direction [100].

The herringbone fragment in the α^II-phase corresponds to a twinning at (001) with a twin domain thickness of one layer. Twins at (001) with a larger domain thickness are described in Section 8.1.

The stacking fault in α^I-quinacridone leads to an energy increase of 4.36 kJ mol⁻¹, corresponding to 2.18 kJ mol⁻¹ per herringbone contact, which is slightly higher than for a pure herringbone stacking of the α^II-phase (ΔE = 1.53 kJ mol⁻¹).

6.3.2. Commensurate misfit-layer structures, i.e. stackings with different lateral periodicities

Fig. 24 shows a three-layer model in which the first and third layers contain eight molecules and the second one only seven. This model was obtained by chance. In our investigation of structures with a missing chain, we constructed a 8 × 1 × 3 supercell and removed the central molecule, which corresponds to the removal of an entire chain along [010]. In the subsequent optimization with free unit-cell parameters, the molecules of the central layer rotated by 13.73°, and the unit-cell parameter c shrank by 0.21 Å. Despite of the missing molecule, the density is almost as high as in the undisturbed structure (1472 versus 1493 kg m⁻³). The resulting structure is shown in Fig. 24. This structure contains eight molecules per supercell in the first and third layer but only seven molecules in the central layer, which corresponds to a commensurate misfit-layer structure with . In inorganic chemistry, such structures are usually modulated. In the case of molecular crystals, the modulation could affect the molecular position, orientation, and conformation (Wagner & Schönleber, 2009). Our model is modulated, too, especially concerning the shift of the individual molecules in the c direction, which is clearly visible in Fig. 24(a). The modulation curve is shown in Fig. 24(c). The energy of this structure is 76.7 kJ mol⁻¹ for the supercell with 8+7+8 molecules, corresponding to 11.0 kJ mol⁻¹ per molecule in the central layer.

Figure 24 Views of misfit-layer structure: (a) in direction [010] and (b) in direction [100]. (c) Modulation of the molecular position in the z direction. Top curves: molecules in the 8-layer, bottom: molecules in the 7-layer.

For inorganic compounds, such structures containing layers with atom deficiencies are known, for example the Magnéli phase Mo₈O₂₃ (Magnéli, 1948). Misfit-layer structures are known in inorganic chemistry, e.g. structures built from two chemically different layers with different lateral dimensions (Makovicky & Hyde, 1981). Similarly, organic chain-misfit structures are known for host–guest systems. In contrast, to the best of our knowledge, organic misfit-layer structures have not been observed experimentally for homomolecular organic compounds. There are several problems with such structures:

(1) Mixed oxidations states, such as Mo^VI/Mo^V in Mo₈O₂₃, are rare in organic compounds.

(2) The individual planes may very easily slip on each other. Consequently, the crystal is probably mechanically unstable.

(3) At the surface of a growing crystal, such a slip may occur very easily, and an additional molecule can be inserted in the deficient layer to achieve the energetically more favourable periodic stacking. This process is even facilitated by the typical high mobility of molecules on growing surfaces.

(4) According to the energy, such misfit layers are very rare. In the diffraction patterns, such a low number of misfit layers would cause only very minor effects (faint diffuse scattering), which would be too weak to be observed experimentally. Also in the pair-distribution function, the effect is too small to be recognised.

In spectroscopic methods (IR, Raman, SS-NMR, UV–vis etc.) the low concentration of misfit layers is probably below the detection limit. Furthermore, the layer misfit does not change the hydrogen bond system and causes only a small shift in the spectra. Only HRTEM, AFM or similar local methods would work.

7. Volume defects

7.1. Voids (three-dimensional vacancy aggregates)

Three-dimensional vacancy aggregates are energetically unfavourable, as already established for aggregates of line vacancies (see Section 5). In a real crystal, such voids will typically be filled with water or other molecules.

7.2. Domains of misoriented molecules

The crystal may contain entire regions with inverted orientation of. Two situations were considered: (i) large sheets of misoriented molecules with the sheets being parallel to (001), (ii) blocks of 4 × 4 misoriented chains parallel to [010].

7.2.1. Lamellar domains with misoriented molecules

Fig. 25 shows a crystal assembled from a periodic arrangement of domains, in which each domain consists of a sheet with a thickness of four molecules with inverted orientation of the molecules. Note that this model does not reflect a twinning, but only a rotation of the molecules. A twin would change the direction of the lattice vectors, but in our model, the vectors of the triclinic lattice were not modified. In the subsequent optimization, the inverted orientation did not cause a change in the molecular positions, and the general packing did not change, even the unit-cell parameters changed by less than 0.2 Å (see Table S4). The energy of the supercell increases by 4.50 kJ mol⁻¹, i.e. by 2.25 kJ mol⁻¹ per phase boundary.

Figure 25 Lamellar domains with rotated molecules. Blue and orange boxes mark the domains of the different molecular orientations. View direction [100].

7.2.2. Blocks of misoriented molecules

Fig. 26 shows a model which contains of blocks of 4 × 4 inverted chains parallel to [010]. As for all other models with inverted chains, neither the position nor the spatial orientation changes significantly. The resulting structure has an energy of 36.0 kJ mol⁻¹ per supercell. The supercell contains 32 modified molecule–molecule contacts at the domain boundaries; hence, the energy increase is as low as 1.1 kJ mol⁻¹ per modified molecule–molecule contact.

Figure 26 Blocks of misoriented molecules. The domains with a different orientation are marked with orange and blue. View direction b*.

8. Twinning

A simple triclinic compound such as α^I-quinacridone can exhibit twinning in various ways. Here we focus on those twins which might explain the twinning observed in the HRTEM image (see Section 10).

8.1. Twinning by mirroring at (001), Model 1

Twinning of the crystal structure of α^I-quinacridone at the (001) face corresponds to a mirroring of the structure at a plane parallel to (001), whereby the mirror plane is located between the molecules. The lattice-energy optimization reveals that a mirror plane would actually lead to unfavourable C—H⋯H—C contacts. An energetically more favourable local structure is formed by replacing the mirror plane with a local glide plane parallel to (001), having an intrinsic translation of about (−0.42, 0.24, 0) (see Fig. 27).

Figure 27 Twinning by mirroring at (001), Model 1. (a) View direction [100]. The local glide planes between the twin domains are drawn in red, the local screw axes in green. (b) View direction [010].

Since the individual twin domains have inversion symmetry, there is a second possibility to describe this local structure at the twin boundary: The combination of the local glide plane with the inversion centre generates a local screw axis parallel to c* with an intrinsic translation of 1/c₀* [Fig. 27(a)].

In the twin structure, the molecular planes in the different twin domains are inclined to each other by 44.9°. If the second twin domain has a thickness of one molecule only, the structure corresponds to the herringbone stacking fault described in Section 6.3.1, except for a misorientation of the central chain. If nearly all domains consist of only one layer of molecules, the structure of the α^II-phase is obtained.

The twin was modelled using domains with a thickness of three molecules – which should be sufficient. Optimization leads to an energy as low as 0.92 kJ mol⁻¹ per twin domain. Hence, this microscopic twinning is energetically very favourable.

8.2. Twinning by mirroring at (001), Model 2

In the previous section, the twinning mirror plane parallel to (001) was located between the molecules. Alternatively, the mirror plane can go through the molecular centres. Since the molecules themselves do not possess a corresponding mirror plane, they must be disordered. However, this disorder leads only to a mutual exchange of CO and NH groups within an entire layer of molecules, which has no great effect on the packing energy, as shown in Section 6.2.3. The adjacent layers follow an exact mirror image, see Fig. 28.

Figure 28 Twinning by mirroring at (001), Model 2. (a) View direction [100]. The local mirror planes are drawn in red. Note that the mirror planes do not act on the molecules, which they cut, but only on the neighbouring and next-neighbouring molecules. (b) View direction [010].

After optimization, the molecules form a wavy arrangement [see Fig. 28(b)]. This structure has an energy of 6.4 kJ mol⁻¹ per twin boundary.

8.3. Twinning by mirroring at (001), Model 3

A third attempt to construct a twin by mirroring the structure at the (001) face leads to the structure shown in Fig. 29. This model has an energy of only 0.77 kJ mol⁻¹ per twin boundary. On a first glance, the structure seems to be reasonable. However, the two twin domains have a different structure. The `left' twin domain does not correspond to the crystal structure of α^I-quinacridone, but to a different, hypothetical polymorph, which has not been observed experimentally.

Figure 29 Attempt to construct a structure with twinning by mirroring at (001), Model 3. (a) View direction [100]. (b) View direction [010]. The left domain, marked by an yellowish orange box, has a different structure than α^I-quinacridone.

9. Discussion: probability of lattice defects

If crystals are grown under thermodynamic control, the frequency of the individual lattice defects can be roughly estimated with the Boltzmann formula:

where n₁, n₂ are probabilities of two local structures, g is degeneracy and is energy difference between the two local structures. All investigated point defects (vacancy, interstitial molecule, misoriented molecule) are energetically so unfavourable that they will hardly occur. Point defects with foreign molecules (e.g. incorporation of water or byproducts) may occur, but were not investigated here.

Similarly, all investigated edge dislocations and screw dislocations are energetically very unfavourable. The only exception is a screw dislocation along [010], which is more favourable because it keeps the hydrogen bonded chains intact and consists only of a translation of the chains along the chain direction. This defect has a relative energy of 38.0 kJ mol⁻¹ per dislocation line, corresponding to a frequency of 5.3 × 10⁻⁸, which indicates that this lattice defect might very rarely occur in a real crystal.

In contrast, the rotation of an entire chain of molecules along the chain axis [010] leads to a surprisingly low energy increase of only 1.57 kJ mol⁻¹ per chain molecule. With the Boltzmann formula one may estimate that about one third of the chains have an inverted orientation. Calculations with two neighbouring chains reveal that this probability even slightly increases if two neighbouring chains are rotated simultaneously.

Similarly, the rotation of all molecules within an entire layer, or within a sheet consisting of several layers or within a block of chains parallel to [010], only requires a small energetic effort. Correspondingly, such defects can easily occur. The rotation of all molecules within a layer corresponds to a stacking fault in which the molecules are rotated, but the 3D-translational periodicity of the molecular positions is left unchanged.

Overall, the crystal should contain approximately one third of chains with inverted orientation. These misoriented chains are distributed almost statistically, with only a weak correlation between neighbouring chains. However, all these orientation defects do neither significantly change the unit-cell parameters, nor the molecular positions. This disorder does not lead to a periodic superstructure – even not to a local superstructure – and correspondingly, not to any superstructure reflections, but only to some diffuse scattering. Since the disorder affects only the exchange of CO by NH groups, the diffuse scattering should be quite weak. In single crystal data, the orientational disorder would surely be apparent from the electron density distribution at the NH and CO groups. In contrast, in the powder diffraction pattern of this nanocrystalline material, the disorder is not visible.

There are two additional lattice defects which are energetically accessible: herringbone stacking faults and twinning:

(1) The parallel stacking can be interrupted by a layer with herringbone stacking. Such a mixture of parallel and herringbone stacking of chains was already observed in α^II-quinacridone. There, the herringbone stacking dominates, and is interrupted by layers with parallel stacking. The energy differences between parallel and herringbone stacking are low, in the α^I-phase as well as in the α^II-phase. This suggests a continuous series of structures from a pure parallel packing to a pure herringbone packing, with the α^I and α^II phases being close to both ends of this series.

(2) The same local structure is observed in twinning along (001). The only difference is that in a twinning the domains are larger than in a single layer. This twinning requires an energy increase of 0.92 kJ mol⁻¹ only, hence it should frequently occur in a real crystal.

10. Experimental observation of lattice defects by HRTEM

A highly crystalline sample of α^I-quinacridone was obtained by Tetsuya Ogawa et al. by vacuum deposition of purified quinacridone on alkali halide single crystals at 140–170°C (Ogawa et al., 1999). Single crystals of sizes up to 700 × 100 × 20 nm were grown. The HRTEM image of one of the crystals was kindly provided by Tetsuya Ogawa. The image is shown in Fig. 30.

Figure 30 (a) HRTEM image of α^I-quinacridone, with kind permission of T. Ogawa. View direction [100], b axis vertical. The white box marks the section enlarged in (b). (b) Section of the HRTEM image, and the corresponding simulated twin Model 1. (c) Simulated structure containing a layer (100) with misoriented molecules. The layer is marked by an arrow.

At first glance, the molecules seem to be well ordered in this HRTEM image. However, the HRTEM image actually shows many lattice defects. For example, the region marked by a white box in Fig. 30(a) contains a twin boundary, which resembles the simulated twin Model 1 [see Fig. 30(b)].

The resolution of the image does not allow the observation of finer defects. For example, a misorientation of a chain of molecules has an effect too weak to be seen in the HRTEM. Even an entire layer of misoriented molecules could hardly be recognized [see Fig. 30(c)]. Correspondingly, the HRTEM image does not provide information on the frequency of misoriented chains. Point defects such as vacancies or small foreign molecules are too weak to be seen. Interstitial molecules would be recognizable through the distortion of the structure in their vicinity, but are absent in the HRTEM image. (A similar effect would be produced by radiation damage to the crystalline structure during TEM imaging, and, therefore would be difficult to interpret.) Other larger defects, such as edge or screw dislocations, grain boundaries, or aggregates of foreign molecules are either absent or not visible.

11. Lattice defects in industrial α^I-quinacridone

All the simulations in this paper imply a crystal growth near the thermodynamic equilibrium. The single crystals obtained by Ogawa et al. (1999) were grown by vacuum sublimation. Nevertheless, they contain a high number of lattice defects. In the industrial process, the crystallization of α^I-quinacridone is very far from thermodynamic equilibrium: a hot solution of probably protonated quinacridone in molten polyphospho­ric acid at 150°C is poured into ice water, whereby the quinacridone instantaneously precipitates as an insoluble, nanocrystalline, crude powder. Under these conditions, the crystals surely contain a high number of defects, including those which were calculated to be energetically unfavourable, e.g. various types of point defects, edge and screw dislocations. Even fragments of misfit-layer structures cannot be ruled out.

12. PDF investigations on α^I-quinacridone

We tried to elucidate the nature of lattice defects in industrial α^I-quinacridone using pair-distribution function (PDF) analysis (Schmidt, 2010).

The PDFs of the four quinacridone polymorphs are significantly different (see Fig. 31). However, the α^I-phase contains so many molecules at surfaces, grain boundaries, probably also in amorphous regions and at various lattice defects that it is impossible to disentangle all these effects. Furthermore, a misorientation of an entire chain, which is calculated to be very frequent, has only a very minor effect to the PDF, because most atoms remain on their original positions.

Figure 31 Pair-distribution functions of quinacridone polymorphs (Schmidt, 2010).

Hence, the PDF does not provide information on the individual lattice defects.

13. Conclusion

In this work, a large variety of lattice defects in α^I-quinacridone was investigated, showing energies and local structures in the vicinity of the defects. An overview is given in Fig. 32. Surprisingly, a rotational flip of an entire chain of molecules is a very facile lattice defect, having an energy of about 1.5 kJ mol⁻¹ only. Correspondingly, about one third of the chains is expected to be misoriented in a real crystal. The ease of the misorientation is caused by the symmetry of the molecular shape, which is close to having an additional mirror plane. Furthermore, twinning at the (001) plane requires only 0.92 kJ mol⁻¹. According to the HRTEM image, this twinning occurs frequently. Other lattice defects, such as vacancies, interstitial molecules, edge and screw dislocations, are energetically quite unfavourable. Nevertheless, they may occur in a real crystal, if the crystallization occurs far from thermodynamic equilibrium.

Figure 32 Overview of the major lattice defects described in this paper, sorted according to energy.

We expect that these results can be transferred to other organic structures with chains of molecules. It should be noted that the misorientation of a molecular chain is only possible if the chain has a nearly pseudosymmetric shape. All other lattice defects can be expected to be as rare as in α^I-quinacridone.

Supporting information

Computing details

(I)

Crystal data

Triclinic, P1	α = 99.30 (8)°
a = 3.8017 (15) Å	β = 99.74 (6)°
b = 6.901 (3) Å	γ = 110.48 (6)°
c = 14.766 (6) Å	V = 346.67 (11) Å^3

Data collection

h = →	l = →
k = →

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

	x	y	z	Uiso*/Ueq
C1	0.46611	0.28601	0.49611	0.0500*
C2	0.70403	0.60903	0.17059	0.0500*
C3	0.70147	0.48007	0.09875	0.0500*
C4	0.69451	0.27711	0.10453	0.0500*
C5	0.63696	0.20556	0.18514	0.0500*
C6	0.60349	0.33169	0.25514	0.0500*
C7	0.63043	0.53313	0.25095	0.0500*
C8	0.62651	0.67441	0.34183	0.0500*
C9	0.5687	0.5945	0.4233	0.0500*
C10	0.53175	0.37945	0.41936	0.0500*
H12	0.73379	0.74749	0.16799	0.0500*
H13	0.72778	0.52618	0.03838	0.0500*
H14	0.7037	0.1787	0.0477	0.0500*
H15	0.58786	0.06056	0.18336	0.0500*
H16	0.53638	0.12138	0.33818	0.0500*
H17	0.44601	0.14561	0.49211	0.0500*
N18	0.56514	0.25474	0.33842	0.0500*
O19	0.66269	0.87869	0.33989	0.0500*

(commensurate-structure)

Crystal data

Triclinic, P1	α = 99.4129°
a = 31.0621 Å	β = 101.5896°
b = 6.9301 Å	γ = 116.5790°
c = 45.0783 Å

Data collection

h = →	l = →
k = →

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

	x	y	z	Uiso*/Ueq
C1	0.08188	0.34219	0.17175	0.00000*
C2	0.06822	0.43396	0.28134	0.00000*
C3	0.06775	0.57262	0.30763	0.00000*
C4	0.07570	0.78903	0.30788	0.00000*
C5	0.08336	0.86461	0.28134	0.00000*
C6	0.08269	0.72509	0.25490	0.00000*
C7	0.07633	0.51117	0.25478	0.00000*
C8	0.07799	0.38655	0.22750	0.00000*
C9	0.08263	0.46460	0.20062	0.00000*
C10	0.08711	0.43463	0.14639	0.00000*
H11	0.06205	0.27687	0.28164	0.00000*
H12	0.06135	0.51485	0.32660	0.00000*
H13	0.07550	0.88973	0.32712	0.00000*
H14	0.08911	1.02054	0.28117	0.00000*
H15	0.08942	0.94390	0.23013	0.00000*
H16	0.07690	0.18421	0.16890	0.00000*
N17	0.08752	0.79978	0.22959	0.00000*
O18	0.07631	0.20170	0.22749	0.00000*
C19	0.09372	0.77380	0.17790	0.00000*
C20	0.10249	0.67255	0.06711	0.00000*
C21	0.10065	0.52970	0.04029	0.00000*
C22	0.09275	0.31347	0.04013	0.00000*
C23	0.08759	0.24267	0.06728	0.00000*
C24	0.09019	0.38572	0.09414	0.00000*
C25	0.09650	0.59949	0.09422	0.00000*
C26	0.09647	0.72725	0.12188	0.00000*
C27	0.09264	0.65078	0.14895	0.00000*
C28	0.08814	0.68069	0.20317	0.00000*
H29	0.10845	0.82916	0.06673	0.00000*
H30	0.10522	0.58408	0.02085	0.00000*
H31	0.09116	0.20931	0.02051	0.00000*
H32	0.08214	0.08737	0.06757	0.00000*
H33	0.08502	0.16987	0.11924	0.00000*
H34	0.09879	0.93194	0.18076	0.00000*
N35	0.08686	0.31390	0.11978	0.00000*
O36	0.09841	0.91268	0.12199	0.00000*
C37	0.20541	0.33790	0.17838	0.00000*
C38	0.19687	0.44075	0.28919	0.00000*
C39	0.19842	0.58347	0.31589	0.00000*
C40	0.20606	0.79920	0.31595	0.00000*
C41	0.21157	0.87031	0.28886	0.00000*
C42	0.20930	0.72745	0.26212	0.00000*
C43	0.20308	0.51391	0.26211	0.00000*
C44	0.20327	0.38615	0.23452	0.00000*
C45	0.20691	0.46197	0.20739	0.00000*
C46	0.21030	0.42945	0.15292	0.00000*
H47	0.19082	0.28403	0.28961	0.00000*
H48	0.19384	0.52926	0.33533	0.00000*
H49	0.20730	0.90292	0.33547	0.00000*
H50	0.21712	1.02573	0.28854	0.00000*
H51	0.21463	0.94326	0.23699	0.00000*
H52	0.20032	0.17976	0.17556	0.00000*
N53	0.21279	0.79927	0.23649	0.00000*
O54	0.20148	0.20108	0.23449	0.00000*
C55	0.21741	0.76976	0.18445	0.00000*
C56	0.22551	0.66640	0.07356	0.00000*
C57	0.22350	0.52303	0.04674	0.00000*
C58	0.21571	0.30701	0.04664	0.00000*
C59	0.21066	0.23663	0.07383	0.00000*
C60	0.21319	0.37978	0.10064	0.00000*
C61	0.21952	0.59352	0.10068	0.00000*
C62	0.21955	0.72158	0.12833	0.00000*
C63	0.21591	0.64570	0.15545	0.00000*
C64	0.21247	0.67811	0.20989	0.00000*
H65	0.23160	0.82320	0.07317	0.00000*
H66	0.22786	0.57696	0.02726	0.00000*
H67	0.21407	0.20257	0.02703	0.00000*
H68	0.20529	0.08154	0.07416	0.00000*
H69	0.20797	0.16406	0.12576	0.00000*
H70	0.22250	0.92792	0.18727	0.00000*
N71	0.20986	0.30817	0.12629	0.00000*
O72	0.22146	0.90692	0.12841	0.00000*
C73	0.32822	0.33319	0.18229	0.00000*
C74	0.32400	0.44549	0.29408	0.00000*
C75	0.32739	0.59184	0.32114	0.00000*
C76	0.33492	0.80723	0.32105	0.00000*
C77	0.33865	0.87464	0.29353	0.00000*
C78	0.33487	0.72875	0.26653	0.00000*
C79	0.32862	0.51525	0.26660	0.00000*
C80	0.32745	0.38463	0.23873	0.00000*
C81	0.33033	0.45864	0.21141	0.00000*
C82	0.33305	0.42424	0.15677	0.00000*
H83	0.31794	0.28889	0.29459	0.00000*
H84	0.32431	0.54070	0.34095	0.00000*
H85	0.33742	0.91351	0.34080	0.00000*
H86	0.34412	1.02983	0.29313	0.00000*
H87	0.33900	0.94195	0.24110	0.00000*
H88	0.32293	0.17470	0.17947	0.00000*
N89	0.33721	0.79809	0.24064	0.00000*
O90	0.32545	0.19915	0.23868	0.00000*
C91	0.34051	0.76544	0.18831	0.00000*
C92	0.35043	0.66508	0.07788	0.00000*
C93	0.34943	0.52347	0.05129	0.00000*
C94	0.34126	0.30678	0.05107	0.00000*
C95	0.33488	0.23389	0.07795	0.00000*
C96	0.33672	0.37577	0.10463	0.00000*
C97	0.34330	0.59000	0.10473	0.00000*
C98	0.34276	0.71707	0.13222	0.00000*
C99	0.33883	0.64077	0.15929	0.00000*
C100	0.33602	0.67498	0.21388	0.00000*
H101	0.35686	0.82256	0.07759	0.00000*
H102	0.35493	0.57941	0.03207	0.00000*
H103	0.34042	0.20392	0.03162	0.00000*
H104	0.32919	0.07813	0.07819	0.00000*
H105	0.33075	0.15879	0.12963	0.00000*
H106	0.34568	0.92372	0.19109	0.00000*
N107	0.33272	0.30302	0.13016	0.00000*
O108	0.34467	0.90241	0.13228	0.00000*
C109	0.45180	0.33052	0.18136	0.00000*
C110	0.44883	0.44837	0.29348	0.00000*
C111	0.45342	0.59758	0.32068	0.00000*
C112	0.46079	0.81245	0.32037	0.00000*
C113	0.46353	0.87710	0.29257	0.00000*
C114	0.45872	0.72864	0.26547	0.00000*
C115	0.45247	0.51534	0.26570	0.00000*
C116	0.45062	0.38261	0.23775	0.00000*
C117	0.45352	0.45584	0.21036	0.00000*
C118	0.45699	0.42136	0.15583	0.00000*
H119	0.44274	0.29189	0.29413	0.00000*
H120	0.45126	0.54891	0.34077	0.00000*
H121	0.46388	0.92030	0.34017	0.00000*
H122	0.46903	1.03222	0.29205	0.00000*
H123	0.46146	0.93844	0.23965	0.00000*
H124	0.44655	0.17226	0.17863	0.00000*
N125	0.46030	0.79588	0.23937	0.00000*
O126	0.44797	0.19594	0.23771	0.00000*
C127	0.46327	0.76078	0.18696	0.00000*
C128	0.47654	0.66498	0.07730	0.00000*
C129	0.47628	0.52438	0.05090	0.00000*
C130	0.46818	0.30790	0.05075	0.00000*
C131	0.46110	0.23417	0.07750	0.00000*
C132	0.46216	0.37493	0.10400	0.00000*
C133	0.46862	0.58884	0.10398	0.00000*
C134	0.46697	0.71410	0.13119	0.00000*
C135	0.46240	0.63709	0.15816	0.00000*
C136	0.45895	0.67160	0.21264	0.00000*
H137	0.48311	0.82264	0.07698	0.00000*
H138	0.48237	0.58104	0.03179	0.00000*
H139	0.46797	0.20605	0.03145	0.00000*
H140	0.45541	0.07842	0.07777	0.00000*
H141	0.45559	0.15726	0.12895	0.00000*
H142	0.46805	0.91819	0.18954	0.00000*
N143	0.45740	0.30112	0.12941	0.00000*
O144	0.46853	0.89855	0.13105	0.00000*
C145	0.57977	0.32426	0.15872	0.00000*
C146	0.56931	0.42233	0.26906	0.00000*
C147	0.57251	0.56739	0.29617	0.00000*
C148	0.58068	0.78421	0.29643	0.00000*
C149	0.58507	0.85395	0.26912	0.00000*
C150	0.58189	0.71008	0.24211	0.00000*
C151	0.57530	0.49575	0.24196	0.00000*
C152	0.57558	0.36843	0.21438	0.00000*
C153	0.58005	0.44608	0.18746	0.00000*
C154	0.58571	0.41746	0.13349	0.00000*
H155	0.56252	0.26420	0.26926	0.00000*
H156	0.56857	0.51355	0.31570	0.00000*
H157	0.58294	0.88942	0.31623	0.00000*
H158	0.59064	1.00961	0.26891	0.00000*
H159	0.58702	0.92587	0.21695	0.00000*
H160	0.57482	0.16647	0.15597	0.00000*
N161	0.58537	0.78228	0.21648	0.00000*
O162	0.57357	0.18288	0.21430	0.00000*
C163	0.59107	0.75458	0.16463	0.00000*
C164	0.60756	0.66653	0.05572	0.00000*
C165	0.60826	0.52780	0.02956	0.00000*
C166	0.59945	0.31003	0.02913	0.00000*
C167	0.59102	0.23334	0.05550	0.00000*
C168	0.59184	0.37347	0.08194	0.00000*
C169	0.59862	0.58812	0.08210	0.00000*
C170	0.59656	0.71255	0.10922	0.00000*
C171	0.59110	0.63328	0.13595	0.00000*
C172	0.58544	0.66191	0.18994	0.00000*
H173	0.61427	0.82458	0.05551	0.00000*
H174	0.61538	0.58639	0.01073	0.00000*
H175	0.59978	0.20958	0.00991	0.00000*
H176	0.58474	0.07640	0.05556	0.00000*
H177	0.58535	0.15545	0.10683	0.00000*
H178	0.59606	0.91249	0.16741	0.00000*
N179	0.58667	0.29839	0.10720	0.00000*
O180	0.59881	0.89850	0.10924	0.00000*
C181	0.70646	0.33123	0.15457	0.00000*
C182	0.69038	0.41364	0.26349	0.00000*
C183	0.69140	0.55342	0.29029	0.00000*
C184	0.70041	0.77206	0.29108	0.00000*
C185	0.70736	0.84826	0.26449	0.00000*
C186	0.70581	0.70867	0.23765	0.00000*
C187	0.69873	0.49311	0.23712	0.00000*
C188	0.70042	0.36952	0.20978	0.00000*
C189	0.70618	0.45097	0.18327	0.00000*
C190	0.71290	0.42676	0.12958	0.00000*
H191	0.68335	0.25467	0.26333	0.00000*
H192	0.68512	0.49425	0.30916	0.00000*
H193	0.70131	0.87372	0.31072	0.00000*
H194	0.71335	1.00499	0.26462	0.00000*
H195	0.71277	0.92929	0.21309	0.00000*
H196	0.70136	0.17293	0.15160	0.00000*
N197	0.71090	0.78511	0.21248	0.00000*
O198	0.69815	0.18323	0.20948	0.00000*
C199	0.71850	0.76342	0.16108	0.00000*
C200	0.73236	0.67486	0.05148	0.00000*
C201	0.73085	0.53334	0.02465	0.00000*
C202	0.72230	0.31575	0.02416	0.00000*
C203	0.71599	0.24163	0.05096	0.00000*
C204	0.71812	0.38313	0.07781	0.00000*
C205	0.72498	0.59815	0.07818	0.00000*
C206	0.72392	0.72343	0.10564	0.00000*
C207	0.71858	0.64329	0.13231	0.00000*
C208	0.71193	0.66766	0.18606	0.00000*
H209	0.73902	0.83296	0.05140	0.00000*
H210	0.73615	0.58974	0.00543	0.00000*
H211	0.72108	0.21301	0.00453	0.00000*
H212	0.70993	0.08495	0.05095	0.00000*
H213	0.71198	0.16461	0.10263	0.00000*
H214	0.72371	0.92192	0.16409	0.00000*
N215	0.71364	0.30832	0.10314	0.00000*
O216	0.72626	0.90979	0.10589	0.00000*
C217	0.83278	0.34108	0.15694	0.00000*
C218	0.81312	0.41728	0.26489	0.00000*
C219	0.81100	0.55186	0.29085	0.00000*
C220	0.82045	0.77118	0.29182	0.00000*
C221	0.83025	0.85235	0.26597	0.00000*
C222	0.83054	0.71580	0.23960	0.00000*
C223	0.82348	0.50028	0.23907	0.00000*
C224	0.82646	0.37901	0.21208	0.00000*
C225	0.83252	0.46091	0.18563	0.00000*
C226	0.83886	0.43602	0.13188	0.00000*
H227	0.80632	0.25868	0.26477	0.00000*
H228	0.80184	0.48828	0.30896	0.00000*
H229	0.81934	0.86910	0.31096	0.00000*
H230	0.83652	1.00957	0.26617	0.00000*
H231	0.83879	0.93865	0.21535	0.00000*
H232	0.82774	0.18288	0.15398	0.00000*
N233	0.83675	0.79415	0.21469	0.00000*
O234	0.82455	0.19354	0.21190	0.00000*
C235	0.84480	0.77335	0.16343	0.00000*
C236	0.85596	0.67925	0.05320	0.00000*
C237	0.85419	0.53709	0.02634	0.00000*
C238	0.84603	0.32024	0.02599	0.00000*
C239	0.84052	0.24791	0.05299	0.00000*
C240	0.84300	0.39022	0.07988	0.00000*
C241	0.84951	0.60456	0.08014	0.00000*
C242	0.84909	0.73112	0.10775	0.00000*
C243	0.84448	0.65244	0.13459	0.00000*
C244	0.83822	0.67751	0.18840	0.00000*
H245	0.86212	0.83638	0.05296	0.00000*
H246	0.85894	0.59242	0.00700	0.00000*
H247	0.84447	0.21662	0.00632	0.00000*
H248	0.83481	0.09199	0.05311	0.00000*
H249	0.83741	0.17287	0.10482	0.00000*
H250	0.84996	0.93176	0.16644	0.00000*
N251	0.83918	0.31680	0.10536	0.00000*
O252	0.85126	0.91712	0.10798	0.00000*
C253	0.95747	0.34499	0.16579	0.00000*
C254	0.93951	0.42760	0.27430	0.00000*
C255	0.93733	0.56285	0.30018	0.00000*
C256	0.94584	0.78031	0.30068	0.00000*
C257	0.95538	0.85982	0.27468	0.00000*
C258	0.95582	0.72265	0.24845	0.00000*
C259	0.94919	0.50813	0.24819	0.00000*
C260	0.95193	0.38577	0.22117	0.00000*
C261	0.95752	0.46590	0.19453	0.00000*
C262	0.96318	0.43860	0.14057	0.00000*
H263	0.93327	0.27025	0.27449	0.00000*
H264	0.92913	0.50174	0.31864	0.00000*
H265	0.94449	0.87853	0.31967	0.00000*
H266	0.96145	1.01645	0.27467	0.00000*
H267	0.96363	0.94379	0.22397	0.00000*
H268	0.95254	0.18705	0.16291	0.00000*
N269	0.96173	0.79964	0.22341	0.00000*
O270	0.95016	0.20069	0.22112	0.00000*
C271	0.96928	0.77660	0.17203	0.00000*
C272	0.97924	0.67842	0.06149	0.00000*
C273	0.97755	0.53605	0.03468	0.00000*
C274	0.96976	0.31994	0.03451	0.00000*
C275	0.96451	0.24879	0.06163	0.00000*
C276	0.96685	0.39118	0.08844	0.00000*
C277	0.97309	0.60486	0.08855	0.00000*
C278	0.97280	0.73193	0.11617	0.00000*
C279	0.96866	0.65467	0.14315	0.00000*
C280	0.96308	0.68213	0.19714	0.00000*
H281	0.98519	0.83505	0.06113	0.00000*
H282	0.98215	0.59065	0.01526	0.00000*
H283	0.96829	0.21609	0.01491	0.00000*
H284	0.95914	0.09363	0.06190	0.00000*
H285	0.96142	0.17465	0.11346	0.00000*
H286	0.97436	0.93478	0.17493	0.00000*
N287	0.96329	0.31874	0.11402	0.00000*
O288	0.97469	0.91728	0.11628	0.00000*
C289	0.05801	0.28657	0.49655	0.00000*
C290	−0.00194	0.27761	0.59165	0.00000*
C291	−0.01182	0.39585	0.61531	0.00000*
C292	−0.00297	0.61493	0.61671	0.00000*
C293	0.01508	0.71391	0.59381	0.00000*
C294	0.02413	0.59533	0.56994	0.00000*
C295	0.01711	0.37939	0.56900	0.00000*
C296	0.02978	0.27863	0.54507	0.00000*
C297	0.04619	0.38236	0.52192	0.00000*
C298	0.07454	0.40311	0.47484	0.00000*
H299	−0.00910	0.11779	0.59099	0.00000*
H300	−0.02570	0.32118	0.63156	0.00000*
H301	−0.01007	0.70057	0.63407	0.00000*
H302	0.02144	0.87189	0.59445	0.00000*
H303	0.04072	0.83615	0.54845	0.00000*
H304	0.05399	0.13030	0.49353	0.00000*
N305	0.03931	0.69294	0.54794	0.00000*
O306	0.02753	0.09248	0.54484	0.00000*
C307	0.06735	0.71361	0.50297	0.00000*
C308	0.12730	0.72258	0.40787	0.00000*
C309	0.13718	0.60434	0.38421	0.00000*
C310	0.12832	0.38525	0.38282	0.00000*
C311	0.11028	0.28627	0.40571	0.00000*
C312	0.10123	0.40486	0.42959	0.00000*
C313	0.10825	0.62080	0.43052	0.00000*
C314	0.09558	0.72155	0.45446	0.00000*
C315	0.07917	0.61782	0.47760	0.00000*
C316	0.05082	0.59707	0.52468	0.00000*
H317	0.13446	0.88240	0.40854	0.00000*
H318	0.15105	0.67901	0.36796	0.00000*
H319	0.13543	0.29962	0.36546	0.00000*
H320	0.10392	0.12830	0.40507	0.00000*
H321	0.08463	0.16404	0.45107	0.00000*
H322	0.07137	0.86989	0.50599	0.00000*
N323	0.08604	0.30725	0.45159	0.00000*
O324	0.09783	0.90770	0.45468	0.00000*
C325	0.19667	0.31364	0.50640	0.00000*
C326	0.13458	0.29904	0.60064	0.00000*
C327	0.12342	0.41451	0.62383	0.00000*
C328	0.13249	0.63407	0.62541	0.00000*
C329	0.15179	0.73591	0.60305	0.00000*
C330	0.16169	0.61932	0.57946	0.00000*
C331	0.15466	0.40327	0.57845	0.00000*
C332	0.16810	0.30437	0.55478	0.00000*
C333	0.18476	0.40897	0.53176	0.00000*
C334	0.21280	0.42966	0.48456	0.00000*
H335	0.12741	0.13909	0.59989	0.00000*
H336	0.10840	0.33751	0.63960	0.00000*
H337	0.12453	0.71765	0.64248	0.00000*
H338	0.15828	0.89423	0.60382	0.00000*
H339	0.17877	0.86147	0.55820	0.00000*
H340	0.19278	0.15761	0.50339	0.00000*
N341	0.17746	0.71845	0.55770	0.00000*
O342	0.16600	0.11853	0.55459	0.00000*
C343	0.20566	0.74001	0.51279	0.00000*
C344	0.26335	0.74512	0.41675	0.00000*
C345	0.27252	0.62562	0.39281	0.00000*
C346	0.26392	0.40699	0.39153	0.00000*
C347	0.24679	0.30976	0.41478	0.00000*
C348	0.23835	0.42947	0.43888	0.00000*
C349	0.24514	0.64498	0.43975	0.00000*
C350	0.23306	0.74673	0.46394	0.00000*
C351	0.21725	0.64404	0.48730	0.00000*
C352	0.18925	0.62343	0.53453	0.00000*
H353	0.27038	0.90467	0.41736	0.00000*
H354	0.28570	0.69900	0.37629	0.00000*
H355	0.27051	0.32034	0.37398	0.00000*
H356	0.24058	0.15210	0.41421	0.00000*
H357	0.22249	0.18995	0.46060	0.00000*
H358	0.20958	0.89606	0.51580	0.00000*
N359	0.22387	0.33312	0.46113	0.00000*
O360	0.23516	0.93263	0.46416	0.00000*
C361	0.33733	0.33351	0.51380	0.00000*
C362	0.27439	0.31199	0.60771	0.00000*
C363	0.26151	0.42317	0.63032	0.00000*
C364	0.27079	0.64339	0.63226	0.00000*
C365	0.29137	0.74899	0.61054	0.00000*
C366	0.30217	0.63512	0.58717	0.00000*
C367	0.29549	0.41947	0.58605	0.00000*
C368	0.30950	0.32280	0.56250	0.00000*
C369	0.32568	0.42781	0.53940	0.00000*
C370	0.35244	0.44905	0.49172	0.00000*
H371	0.26728	0.15195	0.60682	0.00000*
H372	0.24500	0.34255	0.64538	0.00000*
H373	0.26196	0.72453	0.64909	0.00000*
H374	0.29801	0.90778	0.61155	0.00000*
H375	0.31908	0.87838	0.56612	0.00000*
H376	0.33388	0.17833	0.51075	0.00000*
N377	0.31803	0.73571	0.56555	0.00000*
O378	0.30779	0.13775	0.56236	0.00000*
C379	0.34505	0.75754	0.52028	0.00000*
C380	0.39975	0.76234	0.42302	0.00000*
C381	0.40796	0.64220	0.39856	0.00000*
C382	0.39972	0.42401	0.39721	0.00000*
C383	0.38389	0.32793	0.42085	0.00000*
C384	0.37635	0.44821	0.44542	0.00000*
C385	0.38273	0.66321	0.44638	0.00000*
C386	0.37114	0.76467	0.47093	0.00000*
C387	0.35631	0.66238	0.49453	0.00000*
C388	0.32960	0.64124	0.54220	0.00000*
H389	0.40658	0.92164	0.42369	0.00000*
H390	0.42025	0.71496	0.38176	0.00000*
H391	0.40561	0.33684	0.37930	0.00000*
H392	0.37785	0.17046	0.42019	0.00000*
H393	0.36156	0.20936	0.46727	0.00000*
H394	0.34837	0.91246	0.52326	0.00000*
N395	0.36284	0.35238	0.46792	0.00000*
O396	0.37257	0.94937	0.47116	0.00000*
C397	0.48858	0.31235	0.49788	0.00000*
C398	0.43165	0.28864	0.59418	0.00000*
C399	0.42125	0.40114	0.61821	0.00000*
C400	0.42934	0.61986	0.62026	0.00000*
C401	0.44671	0.72297	0.59750	0.00000*
C402	0.45561	0.60870	0.57303	0.00000*
C403	0.44991	0.39430	0.57163	0.00000*
C404	0.46273	0.29847	0.54727	0.00000*
C405	0.47722	0.40380	0.52387	0.00000*
C406	0.50220	0.42863	0.47566	0.00000*
H407	0.42515	0.12934	0.59309	0.00000*
H408	0.40729	0.32186	0.63408	0.00000*
H409	0.42217	0.70213	0.63800	0.00000*
H410	0.45269	0.88088	0.59869	0.00000*
H411	0.47025	0.85090	0.55170	0.00000*
H412	0.48629	0.15936	0.49477	0.00000*
N413	0.46937	0.70846	0.55090	0.00000*
O414	0.46196	0.11488	0.54686	0.00000*
C415	0.49360	0.73137	0.50466	0.00000*
C416	0.54811	0.74888	0.40736	0.00000*
C417	0.55781	0.63467	0.38312	0.00000*
C418	0.55068	0.41788	0.38155	0.00000*
C419	0.53439	0.31755	0.40478	0.00000*
C420	0.52572	0.43274	0.42925	0.00000*
C421	0.53099	0.64610	0.43040	0.00000*
C422	0.51888	0.74374	0.45501	0.00000*
C423	0.50488	0.63980	0.47863	0.00000*
C424	0.47995	0.61508	0.52689	0.00000*
H425	0.55423	0.90723	0.40820	0.00000*
H426	0.57035	0.71096	0.36667	0.00000*
H427	0.55764	0.33470	0.36380	0.00000*
H428	0.52897	0.16092	0.40389	0.00000*
H429	0.51190	0.19257	0.45093	0.00000*
H430	0.49577	0.88410	0.50773	0.00000*
N431	0.51257	0.33450	0.45161	0.00000*
O432	0.51985	0.92768	0.45542	0.00000*
C433	0.63176	0.26880	0.49486	0.00000*
C434	0.57725	0.25129	0.59216	0.00000*
C435	0.56755	0.36551	0.61639	0.00000*
C436	0.57468	0.58230	0.61796	0.00000*
C437	0.59097	0.68263	0.59474	0.00000*
C438	0.59964	0.56743	0.57026	0.00000*
C439	0.59438	0.35407	0.56911	0.00000*
C440	0.60648	0.25643	0.54450	0.00000*
C441	0.62048	0.36037	0.52089	0.00000*
C442	0.64541	0.38509	0.47262	0.00000*
H443	0.57113	0.09294	0.59131	0.00000*
H444	0.55501	0.28921	0.63284	0.00000*
H445	0.56772	0.66548	0.63572	0.00000*
H446	0.59639	0.83925	0.59562	0.00000*
H447	0.61346	0.80760	0.54859	0.00000*
H448	0.62959	0.11607	0.49178	0.00000*
N449	0.61279	0.66567	0.54791	0.00000*
O450	0.60551	0.07249	0.54409	0.00000*
C451	0.63678	0.68782	0.50164	0.00000*
C452	0.69372	0.71154	0.40533	0.00000*
C453	0.70412	0.59904	0.38131	0.00000*
C454	0.69603	0.38032	0.37925	0.00000*
C455	0.67866	0.27720	0.40202	0.00000*
C456	0.66975	0.39147	0.42648	0.00000*
C457	0.67546	0.60587	0.42788	0.00000*
C458	0.66263	0.70170	0.45224	0.00000*
C459	0.64814	0.59638	0.47565	0.00000*
C460	0.62316	0.57154	0.52386	0.00000*
H461	0.70022	0.87084	0.40643	0.00000*
H462	0.71807	0.67832	0.36544	0.00000*
H463	0.70320	0.29805	0.36151	0.00000*
H464	0.67267	0.11929	0.40083	0.00000*
H465	0.65512	0.14927	0.44781	0.00000*
H466	0.63907	0.84081	0.50475	0.00000*
N467	0.65600	0.29171	0.44861	0.00000*
O468	0.66340	0.88529	0.45266	0.00000*
C469	0.78031	0.24264	0.47924	0.00000*
C470	0.72560	0.23785	0.57650	0.00000*
C471	0.71739	0.35798	0.60096	0.00000*
C472	0.72564	0.57617	0.60231	0.00000*
C473	0.74147	0.67226	0.57867	0.00000*
C474	0.74901	0.55197	0.55410	0.00000*
C475	0.74263	0.33697	0.55314	0.00000*
C476	0.75422	0.23552	0.52859	0.00000*
C477	0.76905	0.33780	0.50499	0.00000*
C478	0.79576	0.35894	0.45732	0.00000*
H479	0.71878	0.07855	0.57583	0.00000*
H480	0.70511	0.28522	0.61776	0.00000*
H481	0.71975	0.66334	0.62022	0.00000*
H482	0.74751	0.82972	0.57933	0.00000*
H483	0.76380	0.79082	0.53225	0.00000*
H484	0.77699	0.08772	0.47626	0.00000*
N485	0.76252	0.64780	0.53160	0.00000*
O486	0.75278	0.05081	0.52836	0.00000*
C487	0.78803	0.66668	0.48572	0.00000*
C488	0.85097	0.68819	0.39181	0.00000*
C489	0.86385	0.57701	0.36920	0.00000*
C490	0.85457	0.35679	0.36726	0.00000*
C491	0.83399	0.25119	0.38898	0.00000*
C492	0.82320	0.36506	0.41235	0.00000*
C493	0.82987	0.58071	0.41347	0.00000*
C494	0.81586	0.67738	0.43702	0.00000*
C495	0.79968	0.57237	0.46012	0.00000*
C496	0.77292	0.55114	0.50780	0.00000*
H497	0.85808	0.84823	0.39270	0.00000*
H498	0.88035	0.65763	0.35414	0.00000*
H499	0.86340	0.27565	0.35043	0.00000*
H500	0.82735	0.09240	0.38797	0.00000*
H501	0.80628	0.12181	0.43340	0.00000*
H502	0.79148	0.82185	0.48877	0.00000*
N503	0.80733	0.26447	0.43397	0.00000*
O504	0.81757	0.86244	0.43716	0.00000*
C505	0.91969	0.26018	0.48673	0.00000*
C506	0.86201	0.25507	0.58277	0.00000*
C507	0.85283	0.37457	0.60671	0.00000*
C508	0.86144	0.59320	0.60799	0.00000*
C509	0.87857	0.69043	0.58474	0.00000*
C510	0.88701	0.57072	0.56064	0.00000*
C511	0.88022	0.35521	0.55978	0.00000*
C512	0.89230	0.25345	0.53558	0.00000*
C513	0.90810	0.35615	0.51222	0.00000*
C514	0.93610	0.37676	0.46500	0.00000*
H515	0.85498	0.09551	0.58216	0.00000*
H516	0.83966	0.30118	0.62324	0.00000*
H517	0.85484	0.67984	0.62554	0.00000*
H518	0.88478	0.84809	0.58532	0.00000*
H519	0.90287	0.81024	0.53892	0.00000*
H520	0.91578	0.10413	0.48372	0.00000*
N521	0.90149	0.66707	0.53839	0.00000*
O522	0.89020	0.06756	0.53536	0.00000*
C523	0.92869	0.68655	0.49312	0.00000*
C524	0.99078	0.70114	0.39888	0.00000*
C525	1.00194	0.58567	0.37569	0.00000*
C526	0.99287	0.36612	0.37411	0.00000*
C527	0.97357	0.26427	0.39647	0.00000*
C528	0.96366	0.38086	0.42007	0.00000*
C529	0.97069	0.59691	0.42108	0.00000*
C530	0.95726	0.69581	0.44474	0.00000*
C531	0.94059	0.59122	0.46776	0.00000*
C532	0.91256	0.57052	0.51496	0.00000*
H533	0.99795	0.86110	0.39963	0.00000*
H534	1.01696	0.66267	0.35992	0.00000*
H535	1.00083	0.28253	0.35704	0.00000*
H536	0.96707	0.10596	0.39570	0.00000*
H537	0.94659	0.13872	0.44132	0.00000*
H538	0.93258	0.84258	0.49613	0.00000*
N539	0.94790	0.28174	0.44183	0.00000*
O540	0.95936	0.88165	0.44493	0.00000*
C541	0.03163	0.22638	0.82163	0.00000*
C542	0.02287	0.32763	0.93241	0.00000*
C543	0.02471	0.47048	0.95923	0.00000*
C544	0.03261	0.68670	0.95939	0.00000*
C545	0.03776	0.75750	0.93224	0.00000*
C546	0.03517	0.61446	0.90538	0.00000*
C547	0.02886	0.40069	0.90530	0.00000*
C548	0.02888	0.27293	0.87764	0.00000*
C549	0.03271	0.34940	0.85058	0.00000*
C550	0.03722	0.31949	0.79635	0.00000*
H551	0.01691	0.17101	0.93280	0.00000*
H552	0.02013	0.41609	0.97867	0.00000*
H553	0.03419	0.79086	0.97901	0.00000*
H554	0.04322	0.91280	0.93195	0.00000*
H555	0.04034	0.83031	0.88028	0.00000*
H556	0.02656	0.06824	0.81876	0.00000*
N557	0.03849	0.68628	0.87974	0.00000*
O558	0.02695	0.08749	0.87753	0.00000*
C559	0.04347	0.65799	0.82778	0.00000*
C560	0.05713	0.56622	0.71818	0.00000*
C561	0.05761	0.42756	0.69189	0.00000*
C562	0.04966	0.21116	0.69164	0.00000*
C563	0.04200	0.13557	0.71818	0.00000*
C564	0.04267	0.27509	0.74462	0.00000*
C565	0.04902	0.48901	0.74474	0.00000*
C566	0.04736	0.61362	0.77202	0.00000*
C567	0.04272	0.53558	0.79890	0.00000*
C568	0.03825	0.56554	0.85313	0.00000*
H569	0.06330	0.72331	0.71788	0.00000*
H570	0.06400	0.48533	0.67293	0.00000*
H571	0.04986	0.11045	0.67241	0.00000*
H572	0.03625	−0.02036	0.71835	0.00000*
H573	0.03594	0.05628	0.76939	0.00000*
H574	0.04845	0.81597	0.83062	0.00000*
N575	0.03783	0.20039	0.76993	0.00000*
O576	0.04904	0.79848	0.77203	0.00000*
C577	0.15607	0.22358	0.82750	0.00000*
C578	0.14612	0.32175	0.93803	0.00000*
C579	0.14780	0.46413	0.96484	0.00000*
C580	0.15560	0.68024	0.96501	0.00000*
C581	0.16084	0.75139	0.93789	0.00000*
C582	0.15851	0.60899	0.91108	0.00000*
C583	0.15226	0.39531	0.91097	0.00000*
C584	0.15256	0.26825	0.88335	0.00000*
C585	0.15669	0.34550	0.85637	0.00000*
C586	0.16227	0.31804	0.80238	0.00000*
H587	0.14017	0.16513	0.93839	0.00000*
H588	0.14320	0.40953	0.98426	0.00000*
H589	0.15706	0.78408	0.98461	0.00000*
H590	0.16621	0.90655	0.93762	0.00000*
H591	0.16393	0.82552	0.88606	0.00000*
H592	0.15099	0.06540	0.82460	0.00000*
N593	0.16206	0.68143	0.88550	0.00000*
O594	0.15066	0.08289	0.88324	0.00000*
C595	0.16788	0.65519	0.83373	0.00000*
C596	0.18585	0.57258	0.72522	0.00000*
C597	0.18803	0.43733	0.69934	0.00000*
C598	0.17952	0.21988	0.69884	0.00000*
C599	0.16997	0.14036	0.72484	0.00000*
C600	0.16953	0.27753	0.75107	0.00000*
C601	0.17616	0.49205	0.75133	0.00000*
C602	0.17342	0.61441	0.77835	0.00000*
C603	0.16783	0.53428	0.80499	0.00000*
C604	0.16218	0.56158	0.85895	0.00000*
H605	0.19208	0.72993	0.72503	0.00000*
H606	0.19623	0.49845	0.68088	0.00000*
H607	0.18086	0.12166	0.67985	0.00000*
H608	0.16390	−0.01627	0.72485	0.00000*
H609	0.16172	0.05639	0.77555	0.00000*
H610	0.17282	0.81313	0.83661	0.00000*
N611	0.16362	0.20054	0.77611	0.00000*
O612	0.17520	0.79949	0.77840	0.00000*
C613	0.28056	0.22683	0.83609	0.00000*
C614	0.26940	0.32092	0.94632	0.00000*
C615	0.27117	0.46309	0.97318	0.00000*
C616	0.27932	0.67994	0.97353	0.00000*
C617	0.28484	0.75227	0.94653	0.00000*
C618	0.28236	0.60996	0.91964	0.00000*
C619	0.27584	0.39562	0.91938	0.00000*
C620	0.27626	0.26906	0.89177	0.00000*
C621	0.28088	0.34774	0.86493	0.00000*
C622	0.28713	0.32267	0.81112	0.00000*
H623	0.26323	0.16379	0.94657	0.00000*
H624	0.26641	0.40775	0.99252	0.00000*
H625	0.28088	0.78355	0.99320	0.00000*
H626	0.29055	0.90819	0.94641	0.00000*
H627	0.28794	0.82730	0.89470	0.00000*
H628	0.27539	0.06841	0.83308	0.00000*
N629	0.28617	0.68337	0.89416	0.00000*
O630	0.27409	0.08306	0.89154	0.00000*
C631	0.29257	0.65910	0.84258	0.00000*
C632	0.31223	0.58290	0.73464	0.00000*
C633	0.31435	0.44831	0.70867	0.00000*
C634	0.30491	0.22900	0.70771	0.00000*
C635	0.29510	0.14783	0.73355	0.00000*
C636	0.29481	0.28438	0.75992	0.00000*
C637	0.30187	0.49990	0.76045	0.00000*
C638	0.29890	0.62117	0.78744	0.00000*
C639	0.29284	0.53927	0.81389	0.00000*
C640	0.28649	0.56415	0.86764	0.00000*
H641	0.31904	0.74150	0.73475	0.00000*
H642	0.32351	0.51189	0.69056	0.00000*
H643	0.30601	0.13108	0.68856	0.00000*
H644	0.28883	−0.00939	0.73335	0.00000*
H645	0.28657	0.06153	0.78417	0.00000*
H646	0.29761	0.81729	0.84554	0.00000*
N647	0.28861	0.20602	0.78483	0.00000*
O648	0.30080	0.80663	0.78762	0.00000*
C649	0.40685	0.23676	0.83844	0.00000*
C650	0.39299	0.32532	0.94804	0.00000*
C651	0.39450	0.46684	0.97487	0.00000*
C652	0.40305	0.68443	0.97536	0.00000*
C653	0.40937	0.75855	0.94856	0.00000*
C654	0.40723	0.61705	0.92171	0.00000*
C655	0.40037	0.40203	0.92134	0.00000*
C656	0.40144	0.27675	0.89388	0.00000*
C657	0.40678	0.35688	0.86721	0.00000*
C658	0.41343	0.33252	0.81346	0.00000*
H659	0.38634	0.16722	0.94812	0.00000*
H660	0.38920	0.41044	0.99409	0.00000*
H661	0.40427	0.78717	0.99499	0.00000*
H662	0.41542	0.91523	0.94857	0.00000*
H663	0.41338	0.83557	0.89689	0.00000*
H664	0.40164	0.07826	0.83543	0.00000*
N665	0.41172	0.69186	0.89637	0.00000*
O666	0.39910	0.09039	0.89363	0.00000*
C667	0.41889	0.66894	0.84495	0.00000*
C668	0.43497	0.58653	0.73603	0.00000*
C669	0.43395	0.44675	0.70922	0.00000*
C670	0.42494	0.22811	0.70844	0.00000*
C671	0.41799	0.15192	0.73503	0.00000*
C672	0.41954	0.29151	0.76186	0.00000*
C673	0.42662	0.50706	0.76240	0.00000*
C674	0.42493	0.63065	0.78974	0.00000*
C675	0.41917	0.54921	0.81625	0.00000*
C676	0.41245	0.57342	0.86994	0.00000*
H677	0.44201	0.74550	0.73619	0.00000*
H678	0.44023	0.50592	0.69035	0.00000*
H679	0.42404	0.12644	0.68880	0.00000*
H680	0.41200	−0.00482	0.73490	0.00000*
H681	0.41258	0.07089	0.78643	0.00000*
H682	0.42400	0.82725	0.84792	0.00000*
N683	0.41445	0.21506	0.78704	0.00000*
O684	0.42720	0.81695	0.79004	0.00000*
C685	0.53428	0.24560	0.83489	0.00000*
C686	0.51780	0.33365	0.94379	0.00000*
C687	0.51710	0.47238	0.96996	0.00000*
C688	0.52590	0.69015	0.97038	0.00000*
C689	0.53433	0.76684	0.94401	0.00000*
C690	0.53352	0.62671	0.91757	0.00000*
C691	0.52674	0.41206	0.91742	0.00000*
C692	0.52879	0.28763	0.89029	0.00000*
C693	0.53426	0.36689	0.86357	0.00000*
C694	0.53992	0.33827	0.80957	0.00000*
H695	0.51109	0.17560	0.94400	0.00000*
H696	0.50997	0.41379	0.98879	0.00000*
H697	0.52557	0.79060	0.98960	0.00000*
H698	0.54061	0.92378	0.94395	0.00000*
H699	0.54001	0.84473	0.89268	0.00000*
H700	0.52929	0.08769	0.83210	0.00000*
N701	0.53868	0.70179	0.89231	0.00000*
O702	0.52654	0.10168	0.89027	0.00000*
C703	0.54559	0.67592	0.84080	0.00000*
C704	0.55604	0.57784	0.73046	0.00000*
C705	0.55284	0.43277	0.70334	0.00000*
C706	0.54467	0.21596	0.70308	0.00000*
C707	0.54028	0.14622	0.73040	0.00000*
C708	0.54346	0.29009	0.75741	0.00000*
C709	0.55005	0.50442	0.75755	0.00000*
C710	0.54978	0.63174	0.78514	0.00000*
C711	0.54531	0.55409	0.81206	0.00000*
C712	0.53964	0.58272	0.86602	0.00000*
H713	0.56283	0.73596	0.73026	0.00000*
H714	0.55678	0.48662	0.68382	0.00000*
H715	0.54241	0.11074	0.68328	0.00000*
H716	0.53471	−0.00944	0.73061	0.00000*
H717	0.53834	0.07431	0.78256	0.00000*
H718	0.55053	0.83371	0.84354	0.00000*
N719	0.53998	0.21790	0.78303	0.00000*
O720	0.55179	0.81729	0.78521	0.00000*
C721	0.66209	0.23940	0.81255	0.00000*
C722	0.64881	0.33520	0.92222	0.00000*
C723	0.64907	0.47580	0.94862	0.00000*
C724	0.65718	0.69228	0.94876	0.00000*
C725	0.66426	0.76601	0.92201	0.00000*
C726	0.66319	0.62525	0.89552	0.00000*
C727	0.65673	0.41134	0.89553	0.00000*
C728	0.65838	0.28608	0.86832	0.00000*
C729	0.66295	0.36309	0.84135	0.00000*
C730	0.66641	0.32858	0.78687	0.00000*
H731	0.64224	0.17754	0.92253	0.00000*
H732	0.64299	0.41914	0.96773	0.00000*
H733	0.65738	0.79413	0.96807	0.00000*
H734	0.66995	0.92176	0.92174	0.00000*
H735	0.66976	0.84292	0.87057	0.00000*
H736	0.65730	0.08199	0.80997	0.00000*
N737	0.66795	0.69906	0.87011	0.00000*
O738	0.65683	0.10163	0.86847	0.00000*
C739	0.67356	0.66966	0.81816	0.00000*
C740	0.67652	0.55181	0.70604	0.00000*
C741	0.67194	0.40260	0.67884	0.00000*
C742	0.66456	0.18773	0.67915	0.00000*
C743	0.66183	0.12308	0.70695	0.00000*
C744	0.66663	0.27154	0.73405	0.00000*
C745	0.67288	0.48484	0.73382	0.00000*
C746	0.67474	0.61757	0.76176	0.00000*
C747	0.67184	0.54434	0.78916	0.00000*
C748	0.66837	0.57882	0.84369	0.00000*
H749	0.68262	0.70829	0.70539	0.00000*
H750	0.67409	0.45127	0.65875	0.00000*
H751	0.66148	0.07988	0.65935	0.00000*
H752	0.65633	−0.03204	0.70746	0.00000*
H753	0.66389	0.06174	0.75987	0.00000*
H754	0.67880	0.82792	0.82088	0.00000*
N755	0.66506	0.20430	0.76015	0.00000*
O756	0.67738	0.80424	0.76181	0.00000*
C757	0.78484	0.23474	0.81121	0.00000*
C758	0.77492	0.33510	0.92163	0.00000*
C759	0.77593	0.47671	0.94823	0.00000*
C760	0.78410	0.69340	0.94845	0.00000*
C761	0.79048	0.76629	0.92157	0.00000*
C762	0.78864	0.62440	0.89489	0.00000*
C763	0.78205	0.41018	0.89479	0.00000*
C764	0.78260	0.28310	0.86730	0.00000*
C765	0.78653	0.35941	0.84023	0.00000*
C766	0.78934	0.32520	0.78564	0.00000*
H767	0.76849	0.17761	0.92193	0.00000*
H768	0.77042	0.42077	0.96745	0.00000*
H769	0.78494	0.79625	0.96790	0.00000*
H770	0.79616	0.92205	0.92133	0.00000*
H771	0.79460	0.84139	0.86989	0.00000*
H772	0.77967	0.07646	0.80843	0.00000*
N773	0.79264	0.69716	0.86936	0.00000*
O774	0.78068	0.09776	0.86724	0.00000*
C775	0.79713	0.66699	0.81723	0.00000*
C776	0.80136	0.55469	0.70544	0.00000*
C777	0.79796	0.40834	0.67838	0.00000*
C778	0.79044	0.19295	0.67847	0.00000*
C779	0.78671	0.12554	0.70599	0.00000*
C780	0.79049	0.27143	0.73299	0.00000*
C781	0.79673	0.48492	0.73292	0.00000*
C782	0.79790	0.61554	0.76079	0.00000*
C783	0.79502	0.54154	0.78811	0.00000*
C784	0.79230	0.57593	0.84275	0.00000*
H785	0.80741	0.71129	0.70493	0.00000*
H786	0.80105	0.45948	0.65857	0.00000*
H787	0.78794	0.08667	0.65872	0.00000*
H788	0.78124	−0.02965	0.70639	0.00000*
H789	0.78635	0.05823	0.75842	0.00000*
H790	0.80242	0.82547	0.82005	0.00000*
N791	0.78814	0.20209	0.75888	0.00000*
O792	0.79990	0.80103	0.76084	0.00000*
C793	0.90794	0.23041	0.81507	0.00000*
C794	0.89984	0.33377	0.92596	0.00000*
C795	0.90185	0.47715	0.95278	0.00000*
C796	0.90965	0.69317	0.95288	0.00000*
C797	0.91470	0.76354	0.92570	0.00000*
C798	0.91216	0.62040	0.89888	0.00000*
C799	0.90584	0.40666	0.89884	0.00000*
C800	0.90580	0.27859	0.87119	0.00000*
C801	0.90944	0.35448	0.84407	0.00000*
C802	0.91289	0.32206	0.78963	0.00000*
H803	0.89376	0.17698	0.92635	0.00000*
H804	0.89749	0.42321	0.97227	0.00000*
H805	0.91128	0.79761	0.97249	0.00000*
H806	0.92006	0.91864	0.92536	0.00000*
H807	0.91739	0.83612	0.87376	0.00000*
H808	0.90285	0.07226	0.81225	0.00000*
N809	0.91549	0.69201	0.87323	0.00000*
O810	0.90390	0.09325	0.87111	0.00000*
C811	0.91995	0.66228	0.82114	0.00000*
C812	0.92849	0.55944	0.71033	0.00000*
C813	0.92694	0.41672	0.68363	0.00000*
C814	0.91930	0.20099	0.68358	0.00000*
C815	0.91378	0.12988	0.71067	0.00000*
C816	0.91606	0.27273	0.73741	0.00000*
C817	0.92227	0.48627	0.73741	0.00000*
C818	0.92208	0.61402	0.76500	0.00000*
C819	0.91844	0.53821	0.79213	0.00000*
C820	0.91505	0.57072	0.84660	0.00000*
H821	0.93453	0.71615	0.70991	0.00000*
H822	0.93152	0.47092	0.66419	0.00000*
H823	0.91805	0.09726	0.66405	0.00000*
H824	0.90823	−0.02555	0.71098	0.00000*
H825	0.91072	0.05691	0.76253	0.00000*
H826	0.92503	0.82042	0.82397	0.00000*
N827	0.91257	0.20090	0.76303	0.00000*
O828	0.92388	0.79910	0.76503	0.00000*

(Twin-Model_1)

Crystal data

Triclinic, P1	α = 91.1075°
a = 3.9129 Å	β = 91.1768°
b = 6.9202 Å	γ = 116.3521°
c = 86.6554 Å

Data collection

h = →	l = →
k = →

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

	x	y	z	Uiso*/Ueq
C1	0.82040	0.09909	0.07274	0.00000*
C2	0.30811	−0.00352	0.12595	0.00000*
C3	0.21630	0.09026	0.13875	0.00000*
C4	0.28351	0.30722	0.13874	0.00000*
C5	0.43822	0.42923	0.12568	0.00000*
C6	0.52388	0.33542	0.11285	0.00000*
C7	0.46895	0.12066	0.11293	0.00000*
C8	0.57865	0.04363	0.09973	0.00000*
C9	0.71880	0.17038	0.08663	0.00000*
C10	0.96344	0.23775	0.06041	0.00000*
H11	0.25412	−0.16195	0.12619	0.00000*
H12	0.09836	−0.00063	0.14808	0.00000*
H13	0.21721	0.37502	0.14810	0.00000*
H14	0.48678	0.58602	0.12549	0.00000*
H15	0.66968	0.59871	0.10062	0.00000*
H16	0.78653	−0.05477	0.07147	0.00000*
N17	0.65605	0.45531	0.10045	0.00000*
O18	0.55970	−0.14268	0.09981	0.00000*
C19	0.90085	0.52164	0.07536	0.00000*
C20	1.40813	0.62111	0.02195	0.00000*
C21	1.49764	0.52584	0.00916	0.00000*
C22	1.43146	0.30908	0.00935	0.00000*
C23	1.27937	0.18871	0.02248	0.00000*
C24	1.19531	0.28396	0.03526	0.00000*
C25	1.24985	0.49859	0.03508	0.00000*
C26	1.14167	0.57639	0.04830	0.00000*
C27	1.00227	0.45033	0.06145	0.00000*
C28	0.75764	0.38296	0.08770	0.00000*
H29	1.46154	0.77938	0.02163	0.00000*
H30	1.61282	0.61488	−0.00031	0.00000*
H31	1.49542	0.23956	0.00004	0.00000*
H32	1.23131	0.03204	0.02276	0.00000*
H33	1.05083	0.02143	0.04753	0.00000*
H34	0.93467	0.67547	0.07662	0.00000*
N35	1.06442	0.16482	0.04768	0.00000*
O36	1.16093	0.76273	0.04819	0.00000*
C37	0.88942	0.57239	0.23896	0.00000*
C38	0.37280	0.46213	0.29214	0.00000*
C39	0.27979	0.55405	0.30500	0.00000*
C40	0.34666	0.77100	0.30514	0.00000*
C41	0.50185	0.89478	0.29217	0.00000*
C42	0.58834	0.80271	0.27926	0.00000*
C43	0.53445	0.58811	0.27920	0.00000*
C44	0.64585	0.51300	0.26595	0.00000*
C45	0.78595	0.64134	0.25291	0.00000*
C46	1.03177	0.71274	0.22671	0.00000*
H47	0.31900	0.30368	0.29225	0.00000*
H48	0.16126	0.46182	0.31427	0.00000*
H49	0.27976	0.83752	0.31456	0.00000*
H50	0.55020	1.05157	0.29209	0.00000*
H51	0.73312	1.06728	0.26715	0.00000*
H52	0.85761	0.41904	0.23759	0.00000*
N53	0.72026	0.92406	0.26693	0.00000*
O54	0.62856	0.32704	0.26594	0.00000*
C55	0.96508	0.99375	0.24187	0.00000*
C56	1.48165	1.10401	0.18870	0.00000*
C57	1.57464	1.01208	0.17583	0.00000*
C58	1.50778	0.79513	0.17569	0.00000*
C59	1.35262	0.67135	0.18866	0.00000*
C60	1.26614	0.76343	0.20158	0.00000*
C61	1.32002	0.97803	0.20164	0.00000*
C62	1.20864	1.05314	0.21489	0.00000*
C63	1.06855	0.92481	0.22792	0.00000*
C64	0.82273	0.85340	0.25413	0.00000*
H65	1.53545	1.26245	0.18858	0.00000*
H66	1.69314	1.10431	0.16656	0.00000*
H67	1.57467	0.72861	0.16627	0.00000*
H68	1.30427	0.51457	0.18875	0.00000*
H69	1.12137	0.49886	0.21368	0.00000*
H70	0.99690	1.14711	0.24324	0.00000*
N71	1.13424	0.64208	0.21390	0.00000*
O72	1.22592	1.23909	0.21490	0.00000*
C73	0.95333	0.04461	0.40547	0.00000*
C74	0.44591	−0.05483	0.45888	0.00000*
C75	0.35637	0.04045	0.47167	0.00000*
C76	0.42253	0.25721	0.47148	0.00000*
C77	0.57466	0.37757	0.45835	0.00000*
C78	0.65877	0.28231	0.44557	0.00000*
C79	0.60422	0.06768	0.44575	0.00000*
C80	0.71242	−0.01012	0.43253	0.00000*
C81	0.85187	0.11593	0.41938	0.00000*
C82	1.09658	0.18328	0.39313	0.00000*
H83	0.39249	−0.21311	0.45920	0.00000*
H84	0.24117	−0.04859	0.48113	0.00000*
H85	0.35854	0.32673	0.48078	0.00000*
H86	0.62273	0.53424	0.45806	0.00000*
H87	0.80331	0.54484	0.43330	0.00000*
H88	0.91950	−0.10922	0.40421	0.00000*
N89	0.78970	0.40145	0.43315	0.00000*
O90	0.69315	−0.19647	0.43264	0.00000*
C91	1.03381	0.46717	0.40809	0.00000*
C92	1.54624	0.56975	0.35488	0.00000*
C93	1.63807	0.47596	0.34208	0.00000*
C94	1.57085	0.25900	0.34210	0.00000*
C95	1.41610	0.13700	0.35515	0.00000*
C96	1.33042	0.23082	0.36799	0.00000*
C97	1.38536	0.44557	0.36790	0.00000*
C98	1.27564	0.52262	0.38110	0.00000*
C99	1.13545	0.39587	0.39420	0.00000*
C100	0.89072	0.32851	0.42042	0.00000*
H101	1.60024	0.72818	0.35464	0.00000*
H102	1.75605	0.56685	0.33275	0.00000*
H103	1.63716	0.19120	0.33273	0.00000*
H104	1.36753	−0.01979	0.35535	0.00000*
H105	1.18456	−0.03247	0.38021	0.00000*
H106	1.06770	0.62103	0.40936	0.00000*
N107	1.19820	0.11093	0.38038	0.00000*
O108	1.29461	0.70893	0.38102	0.00000*
C109	0.75774	0.80706	0.57270	0.00000*
C110	1.35972	0.96523	0.62588	0.00000*
C111	1.47349	0.88492	0.63875	0.00000*
C112	1.40706	0.66815	0.63891	0.00000*
C113	1.23106	0.53268	0.62594	0.00000*
C114	1.12374	0.61310	0.61302	0.00000*
C115	1.17734	0.82761	0.61295	0.00000*
C116	1.04456	0.89076	0.59969	0.00000*
C117	0.88368	0.75072	0.58666	0.00000*
C118	0.59601	0.65570	0.56045	0.00000*
H119	1.41361	1.12376	0.62598	0.00000*
H120	1.60683	0.98549	0.64802	0.00000*
H121	1.48925	0.61017	0.64833	0.00000*
H122	1.18268	0.37584	0.62587	0.00000*
H123	0.95974	0.33768	0.60093	0.00000*
H124	0.78712	0.95913	0.57132	0.00000*
N125	0.97215	0.48068	0.60070	0.00000*
O126	1.06162	1.07665	0.59967	0.00000*
C127	0.68732	0.38844	0.57565	0.00000*
C128	0.08530	0.23006	0.52248	0.00000*
C129	−0.02820	0.31030	0.50959	0.00000*
C130	0.03898	0.52722	0.50941	0.00000*
C131	0.21473	0.66273	0.52238	0.00000*
C132	0.32161	0.58235	0.53531	0.00000*
C133	0.26782	0.36781	0.53540	0.00000*
C134	0.40053	0.30470	0.54866	0.00000*
C135	0.56140	0.44478	0.56168	0.00000*
C136	0.84910	0.53980	0.58789	0.00000*
H137	0.03112	0.07148	0.52239	0.00000*
H138	−0.16196	0.20950	0.50034	0.00000*
H139	−0.04225	0.58536	0.49996	0.00000*
H140	0.26340	0.81962	0.52243	0.00000*
H141	0.48552	0.85780	0.54740	0.00000*
H142	0.65792	0.23636	0.57703	0.00000*
N143	0.47307	0.71480	0.54763	0.00000*
O144	0.38347	0.11880	0.54868	0.00000*
C145	0.96254	0.49322	0.73895	0.00000*
C146	1.56455	0.65067	0.79213	0.00000*
C147	1.67825	0.57017	0.80499	0.00000*
C148	1.61175	0.35337	0.80513	0.00000*
C149	1.43570	0.21807	0.79216	0.00000*
C150	1.32840	0.29867	0.77925	0.00000*
C151	1.38215	0.51322	0.77919	0.00000*
C152	1.24946	0.57656	0.76594	0.00000*
C153	1.08844	0.43666	0.75290	0.00000*
C154	0.80051	0.34200	0.72670	0.00000*
H155	1.61849	0.80921	0.79224	0.00000*
H156	1.81161	0.67061	0.81426	0.00000*
H157	1.69390	0.29525	0.81455	0.00000*
H158	1.38726	0.06122	0.79208	0.00000*
H159	1.16418	0.02336	0.76714	0.00000*
H160	0.99214	0.64536	0.73758	0.00000*
N161	1.17668	0.16639	0.76692	0.00000*
O162	1.26672	0.76250	0.76593	0.00000*
C163	0.89159	0.07445	0.74186	0.00000*
C164	0.28960	−0.08300	0.68869	0.00000*
C165	0.17590	−0.00249	0.67582	0.00000*
C166	0.24240	0.21430	0.67568	0.00000*
C167	0.41844	0.34960	0.68866	0.00000*
C168	0.52573	0.26900	0.70157	0.00000*
C169	0.47199	0.05446	0.70163	0.00000*
C170	0.60467	−0.00889	0.71488	0.00000*
C171	0.76569	0.13101	0.72791	0.00000*
C172	1.05361	0.22567	0.75411	0.00000*
H173	0.23566	−0.24153	0.68857	0.00000*
H174	0.04256	−0.10293	0.66655	0.00000*
H175	0.16024	0.27242	0.66627	0.00000*
H176	0.46688	0.50645	0.68874	0.00000*
H177	0.68995	0.54431	0.71367	0.00000*
H178	0.86198	−0.07768	0.74323	0.00000*
N179	0.67745	0.40129	0.71389	0.00000*
O180	0.58742	−0.19483	0.71489	0.00000*
C181	0.16682	1.17899	0.90517	0.00000*
C182	0.76881	1.33725	0.95834	0.00000*
C183	0.88229	1.25697	0.97123	0.00000*
C184	0.81508	1.04005	0.97141	0.00000*
C185	0.63934	0.90457	0.95844	0.00000*
C186	0.53247	0.98498	0.94550	0.00000*
C187	0.58629	1.19953	0.94542	0.00000*
C188	0.45360	1.26268	0.93216	0.00000*
C189	0.29272	1.12262	0.91914	0.00000*
C190	0.00501	1.02765	0.89293	0.00000*
H191	0.82301	1.49584	0.95844	0.00000*
H192	1.01604	1.35775	0.98048	0.00000*
H193	0.89629	0.98187	0.98085	0.00000*
H194	0.59066	0.74768	0.95838	0.00000*
H195	0.36854	0.70955	0.93341	0.00000*
H196	0.19624	1.33108	0.90379	0.00000*
N197	0.38100	0.85256	0.93318	0.00000*
O198	0.47069	1.44858	0.93214	0.00000*
C199	0.09631	0.76035	0.90812	0.00000*
C200	−0.50560	0.60229	0.85494	0.00000*
C201	−0.61935	0.68263	0.84206	0.00000*
C202	−0.55290	0.89941	0.84191	0.00000*
C203	−0.37692	1.03485	0.85488	0.00000*
C204	−0.26962	0.95440	0.86780	0.00000*
C205	−0.32324	0.73989	0.86787	0.00000*
C206	−0.19049	0.67671	0.88113	0.00000*
C207	−0.02961	0.81672	0.89415	0.00000*
C208	0.25807	0.91168	0.92036	0.00000*
H209	−0.55950	0.44376	0.85483	0.00000*
H210	−0.75268	0.58209	0.83279	0.00000*
H211	−0.63506	0.95742	0.83249	0.00000*
H212	−0.32852	1.19169	0.85495	0.00000*
H213	−0.10560	1.22980	0.87989	0.00000*
H214	0.06690	0.60827	0.90950	0.00000*
N215	−0.11803	1.08680	0.88012	0.00000*
O216	−0.20757	0.49080	0.88114	0.00000*