research papers
LCDiXRay: a userfriendly program for powder diffraction indexing of columnar liquid crystals
^{a}Centro di Eccelenza CEMIF.CAL, LASCAMM CRINSTM della Calabria, Dipartimento di Chimica e Tecnologie Chimiche, Università della Calabria, 87036 Arcavacata di Rende (CS), Italy, ^{b}Dipartimento di Ingegneria per l'Ambiente ed il Territorio e Ingegneria Chimica, Università della Calabria, 87036 Arcavacata di Rende (CS), Italy, ^{c}Dipartimento di Fisica, Università della Calabria, 87036 Arcavacata di Rende (CS), Italy, and ^{d}Dipartimento di Informatica, Modellistica, Elettronica e Sistemistica, Università della Calabria, 87036 Arcavacata di Rende (CS), Italy
^{*}Correspondence email: nicolas.godbert@unical.it
The formulation of a standard computerized procedure for the indexing of powder Xray diffraction (PXRD) patterns of columnar liquid crystals, with the determination of all structural information extracted from a properly indexed PXRD spectrum and the attribution of the columnar in situ determination of their structural parameters such as type, cell parameters, area, intermolecular stacking distance between consecutive discoids and, in the case of ordered mesophases, the estimation of the number of molecules constituting each discoid.
symmetry, is presented. In particular, the proposed program notably accelerates the identification of columnar mesophases together with the1. Introduction
Liquid crystals (LCs) are often characterized by powder Xray diffraction (PXRD) analysis since the molecular organization within the
gives rise to well defined PXRD patterns. Surprisingly, to the best of our knowledge, no specific computerized tools have been developed for the resolution of PXRD patterns. To date, the analysis of PXRD patterns has been performed manually by scientists through semi `housemade' rules which are time consuming and often originate indexing mistakes, resulting in misleading determinations of the symmetries and/or cell parameters.Columnar LCs (or discotics) can be divided according to their cell geometry into three different categories, columnar hexagonal (Col_{h}), columnar rectangular (Col_{r}) and columnar oblique (Col_{o}), which were first classified by Levelut in 1983 (Levelut, 1983). In addition, columnar tetrahedral structures (Col_{t}, also denoted Col_{squ}) have been observed for phthalocyanin liquid crystals, probably as a result of their peculiar structure (Ohta et al., 1991; Komatsu et al., 1994), and have been also encountered for nonconventional Tshaped polyphilic triblock molecules based on rodlike biphenyl core systems (Chen et al., 2005). However, a Col_{t} phase can be considered as a Col_{r} of equal parameters; in this regard, a detailed discussion will be presented below. Moreover, a columnar organization showing characteristics of both columnar and smectic phases is known as a lamellocolumnar phase (Col_{L}). However, the indexing of PXRD patterns related to smectic mesophases will not be considered in the present work. While discotics may also show lowordered nematic mesophases [discotic nematic (N_{D}) and columar nematic (N_{C})], these LC phases will not be discussed in this report because of the lack of singularity of their PXRD patterns, due to the high intrinsic disorder. A comprehensive description of these mesophases is given in the excellent review written by Laschat et al. (2007), to which interested readers are referred. Moreover, highly ordered columnar mesophases presenting threedimensional ordered structures such as H phase or plastic phases (Chandrasekhar et al., 2002) will not be considered in this report.
The present article reports on standardized guidelines for the indexing of columnar LCDiXRay. This protocol is based on necessary initial hypotheses coupled with mathematical expressions specific to each type. In this way the generation of various sets of possible indexing values is allowed. Comparing the obtained sets with experimental data, the identification of the most likely solution is achieved, reducing the probability of error. Furthermore, from the generated data set reproducing the experimental PXRD pattern, it is straightforward to access all the structural parameters of the such as geometry and dimension, area, intermolecular stacking distances, and number of molecules within the discoid constituting the columns for ordered mesophases. All the classical columnar mesophases included in this report are summarized in Table 1 with related schemes and structural parameters.
PXRD patterns and the implementation of these guidelines in the program

1.1. General definition and properties of columnar liquid crystals
Since their first identification by Xray diffraction data in 1977 by Chandrasekhar et al. (1977), columnar (or discotic) liquid crystals have received much attention and have attracted growing interest, especially in the past decade, for their applicative role in optoelectronic devices such as light emitting diodes, photovoltaic cells, transistors, displays and photoconductors (SchmidtMende et al., 2002; Hesse et al., 2010; Zheng et al., 2011; Shimizu et al., 2007; Kaafarani, 2011; Sergeyev et al., 2007). The high interest shown towards both organic and inorganic columnar liquid crystals is testified by the high number of dedicated peer reviews, regarding their structure, properties and applications, published over the past few years (Laschat et al., 2007; Kaafarani, 2011; Sergeyev et al., 2007; Kumar, 2006; Kato et al., 2006; Tschierske, 2007). The typical disclike structure of discotic molecules comprises a flat and rigid aromatic core surrounded by numerous long flexible alkyl chains. The columnar organization is directed by π–π stacking interactions of the aromatic cores of the discshaped molecules (discoidal or ellipsoidal). However, discoids can be obtained by selfassembly of several molecules, held together via secondary interactions such as π–π stacking or hydrogenbonding networks. Rodlike molecules have been also reported to selfassemble into ellipsoidal discs, originally observed through association of three single molecules (Guillon et al., 1987) and, as recently reported, association of 15–19 single molecules (Shimogaki et al., 2011). The presence of secondary interactions often contributes to increasing the order and stability of the resulting widening the temperature range of the columnar organization. In conventional discotics, this can be achieved, for instance, by both increasing the dimension of the aromatic core and introducing, within the alkyl chains, some chemical functions favouring hydrogen bonds (Kumar, 2006; Gehringer et al., 2005).
The concept of order in discotics is often ambiguous. There is a difference between the intrinsic order of columnar mesophases (within/between columns) and induced order of the bulk material obtained through alignment. Indeed, discotics can be macroscopically aligned (in a homogenous manner) in two different ways: (i) in a homeotropic manner (i.e. all column directors are placed orthogonally to the substrate plane) or (ii) in a planar orientation (i.e. all column directors are parallel to the substrate). Alignment, required depending on the application sought, can be achieved by appropriate techniques or by using suitable substrates (Li et al., 2010). Aligned samples could be described as a highly ordered discotic However, the terminology Col_{ho}, Col_{ro}, Col_{oo}, indicating ordered hexagonal, rectangular and oblique columnar phases, respectively, is not referred to such aligned mesophases but only used for intrinsic ordered phases, i.e. when the intracolumnar degree of order is high enough to observe in the PXRD pattern the reflection referred to as h_{0}, due to the diffraction generated by stacking of discoids. In the case of ordered columnar phases generated by stacking of molecules within the discoids, supplementary reflections (h_{i}, i > 0) can be observed in the wideangle region. The difference between an ordered and a disordered columnar phase is illustrated in Fig. 1, together with their corresponding expected PXRD patterns. Accordingly, Col_{hd}, Col_{rd} and Col_{od} are referred to as disordered hexagonal, rectangular and oblique columnar phases, respectively.
2. Indexing of a PXRD pattern of columnar liquid crystals
For practical reasons, all the diffraction angles (2θ) will be converted into the distances through the Bragg law,
where θ is the diffraction angle, λ is the wavelength of the monochromatic Xray beam, h, k, l are the of the associated reflection and n is an integer. Hereafter, all diffractions will be considered as their distance equivalents instead of their diffraction angles.
As clearly seen from the schematized PXRD pattern of a columnar LC presented in Fig. 2, two different angle regions can be defined: the smallangle area (from 2θ ≃ 0 to ca 2θ = 12–15°) and the highangle area (for 2θ > ca 12–15°), the latter characterized by the broad halo generated by the slow motion of the flexible molten alkyl chains. The centre of this broad peak is conventionally referred to as the h_{ch} reflection halo.
Within the smallangle area, only reflection peaks relative to the intercolumnar d_{hk}_{0} distances are observable. In the wideangle region characterized by the broad halo h_{ch}, only intracolumnar d_{00l} reflection peaks (note d_{001} = h_{0}), if any, can be present, together with all possible reflection peaks due to intradiscoidal stacking (h_{i}). The indexing of a columnar PXRD pattern will be performed by comparison between the experimental distances d_{hk0} observed in the smallangle region of the spectrum and the calculated data sets of d_{hk}_{0} distances obtained through initial assumptions and mathematical expressions of the specific columnar relative to its twodimensional geometry. For columnar mesophases Col_{h} with hexagonal geometry,
For columnar mesophases Col_{r} with rectangular geometry,
For columnar mesophases Col_{o} with planeparallel geometry
Here a_{h}, a_{r} and b_{r}, and a_{o}, b_{o} and γ are the cell parameters of the columnar hexagonal, rectangular and oblique mesophases, respectively (see Table 1 for cell drawings).
2.1. Indexing of the columnar hexagonal (Col_{h})
As already reported in various studies (Laschat et al., 2007; Chandrasekhar et al., 2002), the indexing of a columnar hexagonal is rather straightforward, because of the high symmetry of the p6mm The calculation of a data set of distances, as showed by equation (2), requires the determination of the unique unknown cell parameter a_{h}, implying the need of just one initial hypothesis. Taking into consideration the first reflection peak of the PXRD pattern (d_{h1}_{k1}_{0}; Fig. 2) and assuming its identity as d_{100}, the following relations are easily deduced from equation (2):
and
From the resulting equation (6), it is now evident that all distances and therefore all the reflection peaks observed in the PXRD pattern will be in typical ratio values with respect to the first indexed peak. These ratio values calculated from equation (6) are reported in Table 2.

However, not all the d_{hk}_{0} peaks will be necessarily observed, since the number and nature of the observed peaks depend highly on the intercolumnar degree of order. Note that the first observed peak d_{h1}_{k1}_{0}, usually indexed as d_{100}, could be d_{110} or even of higher in the case of very large diameter discoids. In this case, new correlation ratios between peaks have to be determined.
2.2. Indexing of columnar rectangular mesophases (Col_{r})
The indexing of the PXRD pattern of a Col_{r} is more challenging because of the absence of a mathematical sequence as previously shown for the Col_{h} Equation (3), which rules out distances of a rectangular geometry, clearly shows the dependency on two unknown unitcell constants, a_{r} and b_{r}. Therefore, in the case of a Col_{r} PXRD pattern, an initial hypothesis for the first two reflection peaks of the spectrum (d_{h1}_{k1}_{0} and d_{h2}_{k2}_{0}; Fig. 2), which are often observed rather close to each other, has to be imposed. The reason of the splitting of the first peak when passing from a Col_{h} to a Col_{r} has already been described (Laschat et al., 2007; Chandrasekhar et al., 2002). For each hypothetical value attributed to the couple h_{1}k_{1}/h_{2}k_{2} used for indexing the first two peaks, a corresponding set of d_{hk}_{0} values can be calculated and compared with the experimental d_{hk}_{0} values. The calculation of these sets is based on the determination of the unknown constants A (A = 1/a_{r}^{2}) and B (B = 1/b_{r}^{2}) derived from equation (3), which can be rewritten for the chosen couple h_{1}k_{1}/h_{2}k_{2} as
By combining these two equations, both A and B can be expressed and calculated as a function of the hypothesized first two reflection peaks.
Note that the most frequently encountered h_{1}k_{1}/h_{2}k_{2} couples for Col_{r} mesophases are 11/20 and 20/11, but other combinations have also been reported, such as 11/02, 01/11, 11/31, 20/02 and 02/13 (Kilian et al., 2000; Morale et al., 2003; Pucci et al., 2005; Venkatesan et al., 2008; Maringa et al., 2008; Kaller et al., 2009; Camerel et al., 2006; Seo et al., 2007; Wuckert et al., 2009; Kaller et al., 2010; Amaranatha Reddy et al., 2005). Once the ideal calculated set has been identified, the unitcell parameters can be easily extracted, being equal to d_{100} (a_{r}) and d_{010} (b_{r}), and the area can be calculated (see Table 1).
A direct comparison between the observed d_{hk}_{0} distances of the PXRD pattern and the calculated sets generated through these equations gives rise to the best possible indexing of the studied spectrum.
The final analysis step to be performed is the eventual determination of the _{r} Four different space groups (c2mm, p2mm, p2gg and p2mg) can be encountered and the accurate determination can be rather tricky. The spacegroup attribution is based on the presence and/or absence of reflection peaks, which implies the availability of a rather high number of reflections in the PXRD pattern, not so often occurring. Extinction rules relative to all space groups and the corresponding illustrations are reported in Table 1.
belonging to the studied Col2.3. Indexing of the columnar tetragonal (Col_{t})
The Col_{t} can be considered as a Col_{r} with a unique cell parameter (a_{r} = b_{r}). Hence, equation (3) can be written as follows:
with a_{t} the parameter of the Col_{t} (see Table 1).
As for the Col_{h} adopting a similar analysis, the distances are dependent on a sole unknown constant (a_{t}), and the resulting ratio values with respect to the first reflection, assumed to be d_{100}, can be extrapolated from equation (9); the results are collected in Table 3.

Again, not all reflections will be necessarily present in the PXRD pattern, though for the p4mm of the tetragonal all reflections are theoretically allowed. Moreover, for very large discoids, the first observed peak d_{h1}_{k1}_{0} could differ from d_{100}, and hence new correlation ratios must be determined through equation (9). The Col_{t} is, however, less frequently encountered and often is observed in the case of highly specific shaped LCs (Ohta et al., 1991; Komatsu et al., 1994; Chen et al., 2005).
2.4. Indexing of the columnar oblique (Col_{o})
The columnar oblique _{o}) represents the most complicated case in the indexing procedure. The presence of three consecutive rather close first reflection peaks in the smallangle region of the PXRD pattern is, however, the first indication of a possible Col_{o} As clearly shown by equation (4), the introduction of the three unknown unitcell parameters (a_{o}, b_{o} and γ) within the rather complex mathematical expression characterizes the Col_{o} Fortunately, the Col_{o} is more rarely encountered because strong core–core interactions between molecules are required to develop this phase (Laschat et al., 2007). Similarly to the case of Col_{r} mesophases previously described, initial hypotheses for the resolution of the Col_{o} PXRD patterns have to be formulated. In this case, the indexing of the first three peaks of the spectrum (d_{h1}_{k1}_{0}, d_{h2}_{k22}0 and d_{h}_{3}_{k3}_{0}; Fig. 2) has to be hypothesized. The calculation of the distance sets will be performed through the determination of three unknown constants A [A = 1/(a_{o}^{2}sin^{2}γ)], B [B = 1/(b_{o}^{2}sin^{2}γ)] and C (C = −2cosγ/a_{o}b_{o}sin^{2}γ) derived from equation (4).
(ColThe PXRD pattern of a Col_{o} will be indexed by comparison between the experimental data and the sets of d_{hk}0 values generated via the procedure summarized in the following scheme:
Note that the symmetry of the Col_{o} corresponds to the p1 allowing all the hk0 reflections to be present. As already mentioned, Col_{o} mesophases are rarely observed and therefore it is difficult to define which initial triad values h_{1}k_{1}, h_{2}k_{2}, h_{3}k_{3} are more frequently encountered for the indexing of the first three peaks. Examples reported for initial attribution are 20/11/1, 10/11/2, 11/20/1 and 10/01/11 (Morale et al., 2003; Trzaska et al., 1999; Choi et al., 2011; Pucci et al., 2011).
2.5. Indexing refinement
Once the most likely nature of the columnar d_{hk}_{0} values with a more accurate fitting with respect to all the observed experimental ones, i.e. reducing the discrepancy between observed and calculated data, it is necessary to proceed to a reevaluation (or refinement) of the values of the initial peaks, taking into account the maximum number of data available. Following this calculation, all the d_{hk0} peaks are redetermined using the refined values, obviously through the mathematical expression of the identified mesophase.
has been identified, in order to obtain a data set of calculatedFor the Col_{h} the value of d_{100} can be reevaluated from equation (10) (Zelcer et al., 2007):
where N_{hk0} is the number of hk0 observed reflections. For Col_{r} mesophases, the reevaluation of the first two interplanar distances (d_{100} and d_{010}) can be performed. In order to take into account the maximum number of available data, it is recommended to determine first the mean value of d_{100} through equation (11) when the number of observed h00 reflections is higher than the number of 0k0 observed reflection peaks,
(where N_{h00} is the number of h00 observed reflections), and use this reevaluated value of d_{100} in equation (12) for reevaluating d_{010}:
(where N_{hk0} is the number of hk0 observed reflections). In the opposite case (when N_{0k0} > N_{h}_{00}), the reevaluation of d_{010} should be performed first. This can be achieved by rewriting equation (11) for d_{010}. Then, d_{100} will be obtained through the appropriate transcription of equation (12).
The Col_{t} procedure can be performed through an analogous methodology. Since in this particular rectangular d_{100} is equal to d_{010}, only the reevaluation of the first interplanar distance value is necessary. This can be achieved via
where N_{hk0} is the number of hk0 observed reflections. Finally, the Col_{o} procedure requires the redetermination of the two first interplanar distance values d_{100} and d_{010} as well as the reevaluation of the γ angle, which can be achieved by the following equations:
where N_{h00} is the number of h00 observed reflections,
where N_{0k0} is the number of 0k0 observed reflections, and
where N_{hk0} is the number of hk0 observed reflections.
2.6. Number of molecules in discoids
Once the indexing of a PXRD pattern of a columnar i.e. when the h_{0} = d_{001} reflection peak is present). The number of molecules per (or cross section) (z) can be estimated according to (Lehmann et al., 2006)
has been performed, the number of molecules within the discoid at the origin of the columnar stacking can be determined, but only for ordered phases (where ρ is the density of the phase, N_{A} is Avogadro's constant, S is the columnar area, h_{0} is the height of the columnar slice and M is the molecular weight of the constitutive molecule. S can be easily calculated from the cell parameters, keeping in mind that its expression depends on the geometry of the cell (see Table 1). The density can be estimated unless measured experimentally. Density values ranging from 0.9 g cm^{−3} up to 1.2 g cm^{−3} have been reported depending on the nature of the studied liquid crystals (full organic molecules, metallomesogens, neutral or ionic mesogens) (Gunyakov et al., 2003; Kaller et al., 2009, 2010; Ionescu et al., 2012). The number of discoids present in the (Z_{disc}) is characteristic of the geometry and the proposed model (see Table 1). Consequently, the number of molecules per discoid (N) is easily accessible by dividing the number of molecules in the cross section (z) by the corresponding Z_{disc} value. Therefore, a hypothesis on how mesogen molecules are eventually organized to form the discoidal shape can be formulated. This finding can be further supported by experimental data when intradiscoidal stacking gives rise to sufficiently intense reflections in the wideangle region of the PXRD pattern (see Fig. 2). In the absence of such information, only theoretical modelling can shed light on the discoidal molecular assembly. However, it has to be mentioned that this calculation is based on the initial hypothesis that columns are organized in the with a straight vertical stacking order. In reality the intracolumnar stacking distance h might be higher or lower than h_{0}. This is the case for tilted (h > h_{0}) and undulating (h < h_{0}) columnar organizations (Weber et al., 1991; Donnio et al., 1997). Neither case will be considered by LCDiXRay in its current version, but they will be considered in future upgrades.
3. Software description and examples
3.1. Algorithm and software implementation
A userfriendly program has been implemented in a Java objectoriented framework in order to identify correctly the via indexing of PXRD patterns and, as a consequence, determine the optimal parameters best fitting the observed data. It is worth underlining the advances in terms of computational speed of the proposed approach. The framework has been designed on the basis of the previously described indexing procedures for Col_{h}, Col_{r}, Col_{t} and Col_{o} mesophases.
The core idea is to rewrite the experimental data d_{i} as a function of the peak indices h_{i}, k_{i} and unknown parameters { p_{k}}^{n}_{k = 1} for all types of according to the following equation:
types, a relative interface has been initialized with appropriate methods as shown in the class diagram reported in Fig. 3The mathematical equations of the interplanar distances [equations (2)–(4)] can be described in the form of equation (18) with { p_{k}}^{n}_{k = 1} representing the unknown structural parameters of the hypothesized Given an nlength userselected input data set { d_{k}}^{n}_{k = 1}, coupled with n different indices ( h_{i}, k_{i}), the related n parameters are computed by solving the system formed by the n equations derived from equation (18). Once such parameters are obtained, an estimate d_{gen} of the entire experimental data set is generated by varying the indices h and k in an appropriate userselected discrete grid. An optimization search is performed with the indexing couples ( h_{i}, k_{i}) as control parameters, by considering a rootmeansquare difference (RMSD) optimization process:
The output of the above procedure represents the set of optimal parameters providing the minimum RMSD index value within the range of userselected discrete indices, according to the chosen generation model.
The generation model, i.e. hexagonal, rectangular or oblique, and the data set extracted from the experimental PXRD pattern are needed as initial input. The Java program is characterized by dynamic panels that are updated according to the chosen generation model and provides a test on the consistency of the data set for the indexing problem. The class diagram of the designed user interface is depicted in Fig. 4.
The `InteractiveTableModel' provides useful tools to choose an appropriate number of observed peaks representing the initial hypothesis for the indexing algorithm from which all the h_{max} and k_{max}) represent the maximum admissible values to generate all the observed diffraction peaks. An appropriate choice of parameters h_{max} and k_{max} is required to avoid a set of d_{hk}_{0} indexing that does not have any realistic significance but which could mathematically result in a lower RMSD index.
parameter(s) will be determined. It is also possible to select the Miller index range defining the admissible indexing values. A similar range definition is required in the generation process of the estimated data set. To this end, two (Finally, once a convergent result has been obtained which suits the user's expectations, a popup menu allows them, as a first option, to proceed to the
of the indexing, following the methodologies previously described. A second option allowing the determination of the number of molecules within the discoid is also accessible only for ordered mesophases.Ultimately, all the generated data, the indexing set and the optimal parameters can be exported for further manipulations.
3.2. Practical examples
3.2.1. Ordered hexagonal columnar mesophase
Although the indexing of a Col_{h} is rather straightforward to perform because of the characteristic ratios between the observed peaks of a PXRD pattern (see §2.1), the following example of a Col_{ho} discotic exhibited by discotic (I) (see Fig. 5) (McKenna et al., 2005), a poly(propylene imine) dendrimer based on a triphenylene, is a perfect example to illustrate at first glance the performance of the LCDiXRay program. The experimental and calculated literature data for (I) are presented in Fig. 6, together with the initial window frame relative to the loaded data.
The hexagonal model has been selected on the LCDiXRay program window, as well as the first observed peak (d_{1} = 58.1 Å), which will be taken into account to determine the parameter. On the far right of the program window, only the first Miller index subset is active, corresponding to the userselected discrete grid of indices allowed for the selected initial peak. Having chosen as a first subset the values 0 and 2 for h_{1} and k_{1} minimum and maximum, respectively, the only possible values to be considered for d_{1} are d_{hk}_{0} = d_{100}, d_{200}, d_{110}, d_{210} and d_{220}. Finally, on the bottom right of the window, the maximum values of h and k that will be considered to index all observed peaks (d_{2}–d_{9}) have both been fixed equal to 5. Note that the 11th observed peak d_{11} = 3.55 Å has been assigned as h_{0} and the tenth peak d_{10} = 4.5 Å has been assigned as h_{CH}; both values will therefore be removed from the fitting procedure. The active window obtained after a validation check is presented in Fig. 7, clearly showing the expected fitting for a hexagonal with an RMSD index of 0.4770 for a parameter of a_{h} = 67.09 Å. To reproduce exactly the literature data reported in Fig. 6 (McKenna et al., 2005), the initial selected peak (from which all calculations are performed) must be d_{3} = 28.9 Å (see supplementary information^{1}), the initial choice made by the authors. In these conditions, an RMSD index of 0.3933 and a parameter a_{h} = 66.74 Å are obtained.
An active popup window in LCDiXRay allows the user to proceed to the data to reduce eventually the RMSD between experimental and calculated interplanar distances, according to the method previously described. The refined data reported in Fig. 8 and obtained from the optimization results based on the selected observed peak at 58.1 Å show the lowering of the RMSD index to 0.2324 with a final parameter a_{h} = 66.65 Å. Note that this RMSD index is even lower than that obtained reproducing the literature indexing (i.e. considering d_{3} = 28.9 Å as initial selected peak).
Through this example, we can clearly see the genuine performance of LCDiXRay, allowing in a few clicks the generation of the optimized indexing despite the initial assumption made.
3.2.2. Ordered rectangular columnar with determination of the number of molecules within a discoid
A relevant example of the indexing of a Col_{r} is represented by the discotic exhibited by the cyclopalladated photoconductive Nile red complex (II) (Fig. 9) (Ionescu et al., 2012).
The active LCDiXRay window of the loaded interplanar distances observed on the PXRD pattern of discotic (II) with its initial assignment, together with the LCDiXRay result window obtained after checking the hypothesis of a Col_{r} model and of the obtained preliminary results, are shown in Fig. 10.
Starting from the loaded data (active window of LCDiXRay; Fig. 10a) a twoclick procedure (check and supplementary Figs. S2 and S3 ) allows the visualization of the most probable indexing of the of (II) through a Col_{r} model (Fig. 10b). We obtain a final RMSD value of 0.1901 with parameters a_{r} = 37.09 Å and b_{r} = 65.04 Å. In this case, the first two observed peaks were initially chosen to perform all calculations, with a selected discrete grid of ranging from 0 to 2 (h_{imin} = k_{imin} = 0; h_{imax} = k_{imax} = 2, with i = 1 or 2).
A second activated window allows the determination of the number of molecules within the (a). Following equation (2), the number of molecules present within the cross section has been determined for four different standard density values. As expected for metallomesogens, a density of 1.2 g cm^{−3} can reasonably be attributed for discotic (II), resulting in four molecules being present in the cross section. According to the obtained indexing together with the extinction rules (see Table 1), the of the for (II) can be deduced as p2gg. Unequivocally, since in the p2gg two discoids are present in the (Z_{disc} = 2; Fig. 11b), each disc contains two molecules, forming a dimerlike discoid. Theoretical calculations have been performed confirming the existence of a hydrogenbond network established between two sidebyside molecules placed in a headtotail arrangement as illustrated in Fig. 11(c) (Ionescu et al., 2012).
area of the as illustrated in Fig. 113.2.3. Disordered oblique columnar mesophase
Finally, to illustrate the performance and the swiftness of the LCDiXRay program, the identification of a columnar oblique is presented through the indexing of the PXRD pattern recorded at 393 K for the cobalt complex (III) (see Fig. 12) (Morale et al., 2003).
As shown in Fig. 13(a) for the Col_{o} model, all subsets of the LCDiXRay program window are now activated, since three interplanar distances have to be initially selected to index as a Col_{o} The literaturereported interplanar distances have been loaded and, again, in only one validation click, the results window allows us to check the feasibility of the Col_{o} hypothesis. The results window (Fig. 13b) gives the indexing with the smallest RMSD between experimental and calculated distances, based as always on the userselected choice of the initial The obtained results are in agreement with reported literature data.
Note that in this example the Figs. S4 and S5 ). This is due to the small number of observed reflection peaks on which obviously all the calculations are based. In particular, the lack of 0k0 reflection peaks means that the d_{010} reflection value has to be fixed, excluding it from the Since for Col_{o} mesophases three initial peaks are required, the tool can only achieve realistic results when a large number of reflection peaks are observed, which is unfortunately rarely encountered for this type of mesophase.
of the data performed following the previously described procedure will lead to a slight increase of the RMSD value, which varies from 0.2038 (unrefined) to 0.2218 (refined) (supplementary3.3. Experimental errors
LCDiXRay performs the procedure from the loaded interplanar distances values, although the scattering angles 2θ are the experimental data collected over an Xray powder The 2θ experimental error is highly dependent on the experimental technique used (reflection versus transmission). Furthermore, the interplanar distance accuracy is intrinsically generally much lower for large d values than for low d values. Not surprisingly, throughout the literature on PXRD of mesophases, data accuracy is most of the time left undiscussed and experimental data are reported in d values in ångström with one (in most cases) or at most two digits. In some cases, to take into account such error in the accurate determination of distances, several indexings of possible are proposed. For all these reasons, the current version of LCDiXRay will not take into account this unpredictable experimental error. Consequently, the RMSD index given by LCDiXRay is mostly indicative and only relevant for comparisons, and the accuracy of the constant values obtained through the program must be critically evaluated by the user. Similarly, in its current version only one set of data is proposed by the software, which corresponds to the lowest RMSD value obtained from the loaded values and initial restrictions. The proposal of several solutions, with restricted deviations reflecting the experimental errors obtained during the recording of the PXRD pattern, will be implemented in the near future.
4. Conclusions
LCDiXRay is a userfriendly program for powder diffraction indexing of columnar liquid crystals. The main objective of this powerful tool is to accelerate the determination of the exact nature of the presented by columnar LCs. Furthermore, the program contains all the mathematical expressions required for data and determination of structural parameters of the identified mesophases (unitcell geometry and dimensions, crossarea section, number of molecules within the discoid for ordered phases). The determination of columnar mesophases, in particular Col_{r} and Col_{o}, is often performed manually; LCDiXRay is able in a few minutes to perform all the timeconsuming operations and provide the most likely indexing of a recorded PXRD pattern. Further implementations of LCDiXray are planned in due course and will concern also the determination of the most likely for Col_{r} mesophases on the basis of the extinction rules characterizing them, as well as the possibility of performing a of Col_{r} data as a Col_{t} when the cell parameters are compatible (a_{r} ≃ b_{r}). Finally, it is expected that a similar procedure will be introduced for the indexing of calamitic mesophases, allowing the inclusion in LCDiXray of lamellar columnar and nonclassical mesophases, such as plastic columnar phases or the threedimensionally highly ordered mesophases (H phase).
The LCDIXRay program is freely available from the authors on request. We would like to suggest that users send feedback for future improvements. A request form is available in the supporting information .
Supporting information
Additional figures. DOI: 10.1107/S1600576714003240/nb5093sup1.pdf
Request form . DOI: 10.1107/S1600576714003240/nb5093sup2.pdf
Acknowledgements
This work was supported by the European Community's Seventh Framework Programme (FP7 2007–2013), through the MATERIA project (PONa3_00370).
References
Camerel, F., Donnio, B., Bourgogne, C., Schmutz, M., Guillon, D., Davidson, P. & Ziessel, R. (2006). Chem. Eur. J. 16, 4261–4274. Web of Science CSD CrossRef
Chandrasekhar, S., Krishna Presad, S., Shankar Rao, D. S. & Balagurusamy, V. S. K. (2002). Proc. Indian Natl Sci. Acad. 68A, 175–191.
Chandrasekhar, S., Sadashiva, B. K. & Suresh, K. A. (1977). Pramana, 9, 471–480. CrossRef CAS Web of Science
Chen, B., Baumeister, U., Pelzl, G., Das, M. K., Zeng, X., Ungar, G. & Tschierske, C. (2005). J. Am. Chem. Soc. 127, 16578–16591. Web of Science CrossRef PubMed CAS
Choi, J., Han, J., Ryu, M. & Cho, B. (2011). Bull. Korean Chem. Soc. 32, 781–782. Web of Science CrossRef CAS
Donnio, B., Heinrich, B., GulikGrzywicki, Th., Delacroix, H., Guillon, D. & Bruce, D. W. (1997). Chem. Mater. 9, 2951–2965. Web of Science CrossRef CAS
Gehringer, L., Bourgogne, C., Guillon, D. & Donnio, B. (2005). J. Mater. Chem. 15, 1696–1703. Web of Science CrossRef CAS
Grolik, J., Dudek, L., Eilmes, J., Eilmes, A., Górecki, M., Frelek, J., Heinrich, B. & Donnio, B. (2012). Tetrahedron, 68, 3875–3884. Web of Science CrossRef CAS
Guillon, D., Skoulios, A. & Malthete, J. (1987). Europhys. Lett. 3, 67–72. CrossRef CAS Web of Science
Gunyakov, V. A., Shestakov, N. P. & Shibli, S. M. (2003). Liq. Cryst. 30, 871–875. Web of Science CrossRef CAS
Hesse, H. C., Weickert, J., AlHussein, M., Dössel, L., Feng, X., Müllen, K. & SchmidtMende, L. (2010). Sol. Energ. Mater. Sol. Cells, 94, 560–567. Web of Science CrossRef CAS
Ionescu, A., Godbert, N., Crispini, A., Termine, R., Golemme, A. & Ghedini, M. (2012). J. Mater. Chem. 22, 23617–23626. Web of Science CrossRef CAS
Kaafarani, B. R. (2011). Chem. Mater. 23, 378–396. Web of Science CrossRef CAS
Kaller, M., Deck, C., Meister, A., Hause, G., Baro, A. & Laschat, S. (2010). Chem. Eur. J. 16, 6326–6337. Web of Science CrossRef CAS PubMed
Kaller, M., Tussetschlager, S., Fischer, P., Deck, C., Baro, A., Giesselmann, F. & Laschat, S. (2009). Chem. Eur. J. Chem. 15, 9530–9542. Web of Science CrossRef CAS
Kato, T., Mizoshita, N. & Kishimoto, K. (2006). Angew. Chem. Int. Ed. 45, 38–68. Web of Science CrossRef CAS
Kilian, D., Knawby, D., Athanassopoulou, M. A., Trzaska, S. T., Swager, T. M., Wrobel, S. & Haase, W. (2000). Liq. Cryst. 27, 509–521. CrossRef CAS
Komatsu, T., Ohta, K., Watanabe, T., Ikemoto, H., Fujimoto, T. & Yamamoto, I. (1994). J. Mater. Chem. 4, 537–540. CrossRef CAS Web of Science
Kumar, S. (2006). Chem. Soc. Rev. 35, 83–109. Web of Science CrossRef PubMed CAS
Laschat, S., Baro, A., Steinke, N., Giesselmann, F., Hägele, C., Scalia, G., Judele, R., Kapatsina, E., Sauer, S., Schreivogel, A. & Tosoni, M. (2007). Angew. Chem. Int. Ed. 46, 4832–4887. Web of Science CrossRef CAS
Lehmann, M., Köhn, C., Meier, H., Renker, S. & Oehlhof, A. (2006). J. Mater. Chem. 16, 441–451. Web of Science CrossRef CAS
Levelut, A. M. (1983). J. Chim. Phys. 80, 149–161. CAS
Li, J., He, Z., Zhao, H., Gopee, H., Kong, X., Xu, M., An, X., Jing, X. & Cammidge, A. N. (2010). Pure Appl. Chem. 82, 1993–2003. CAS
Maringa, N., Lenoble, J., Donnio, B., Guillon, D. & Deschenaux, R. (2008). J. Mater. Chem. 18, 1524–1534. Web of Science CrossRef CAS
McKenna, M. D., Barberá, J., Marcos, M. & Serrano, J. L. (2005). J. Am. Chem. Soc. 127, 619–625. Web of Science CrossRef PubMed CAS
Morale, F., Date, R. W., Guillon, D., Bruce, D. W., Finn, R. L., Wilson, C., Blake, A. J., Schröder, M. & Donnio, B. (2003). Chem. Eur. J. 9, 2484–2501. Web of Science CSD CrossRef PubMed CAS
Ohta, K., Watanabe, T., Hasebe, H., Morizumi, Y., Fujimoto, T., Yamamoto, I., Lelièvre, D. & Simon, J. (1991). Mol. Cryst. Liq. Cryst. 196, 13–26. CrossRef CAS Web of Science
Pucci, D., Barberio, G., Bellusci, A., Crispini, A., La Deda, M., Ghedini, M. & Szerb, E. I. (2005). Eur. J. Inorg. Chem. 12, 2457–2463. Web of Science CrossRef
Pucci, D., Crispini, A., Ghedini, M., Szerb, E. I. & La Deda, M. (2011). Dalton Trans. 40, 4614–4622. Web of Science CSD CrossRef CAS PubMed
Reddy, R. A., Raghunathan, V. A. & Sadashiva, B. K. (2005). Chem. Mater. 17, 274–283. Web of Science CrossRef CAS
SchmidtMende, L., Fechtenkötter, A., Müllen, K., Friend, R. & MacKenzie, J. (2002). Physica E, 14, 263–267. CAS
Seo, J., Kim, S., Gihm, S. H., Park, C. R. & Park, S. Y. (2007). J. Mater. Chem. 17, 5052–5057. Web of Science CrossRef CAS
Sergeyev, S., Pisula, W. & Geerts, Y. H. (2007). Chem. Soc. Rev. 36, 1902–1929. Web of Science CrossRef PubMed CAS
Shimizu, Y., Oikawa, K., Nakayama, K. & Guillon, D. (2007). J. Mater. Chem. 17, 4223–4229. Web of Science CrossRef CAS
Shimogaki, T., Dei, S., Ohta, K. & Matsumoto, A. (2011). J. Mater. Chem. 21, 10730–10737. Web of Science CrossRef CAS
Trzaska, S. T., Zheng, H. & Swager, T. M. (1999). Chem. Mater. 11, 130–134. Web of Science CrossRef CAS
Tschierske, C. (2007). Chem. Soc. Rev. 36, 1930–1970. Web of Science CrossRef PubMed CAS
Venkatesan, K., Kouwer, P. H. J., Yagi, S., Müller, P. & Swager, T. M. (2008). J. Mater. Chem. 18, 400–407. Web of Science CSD CrossRef CAS
Weber, P., Guillon, D. & Skoulios, A. (1991). Liq. Cryst. 9, 369–382. CrossRef CAS Web of Science
Wuckert, E., Hägele, C., Giesselmann, F., Baro, A. & Laschat, S. (2009). Beilstein J. Org. Chem. 5, 57. Web of Science CrossRef PubMed
Zelcer, A., Donnio, B., Bourgogne, C., Cukiernik, F. D. & Guillon, D. (2007). Chem. Mater. 19, 1992–2006. Web of Science CrossRef CAS
Zheng, Q., Fang, G., Sun, N., Qin, P., Fan, X., Cheng, F., Yuan, L. & Zhao, X. (2011). Sol. Energ. Mater. Sol. Cells, 95, 2200–2205. Web of Science CrossRef CAS
© International Union of Crystallography. Prior permission is not required to reproduce short quotations, tables and figures from this article, provided the original authors and source are cited. For more information, click here.