 Research Article
 Open Access
Crystal structures of increasingly large molecules: meeting the challenges with CRYSTALS software
 Pascal Parois^{1}Email author,
 Richard I Cooper^{1} and
 Amber L Thompson^{1}
https://doi.org/10.1186/s1306501501054
© Parois et al.; licensee Springer. 2015
 Received: 13 October 2014
 Accepted: 13 May 2015
 Published: 24 May 2015
Abstract
Background
The size and complexity of molecules being studied by single crystal diffraction is growing year by year, resulting in an increase in the difficulties encountered during structure determination. From the crystallisation itself and sample handling, to structure solution and refinement, specific problems due to larger molecules are discussed.
Results
During refinement, several methods are available to deal with the problems encountered with large structures within the software Crystals. Hydrogens atoms can neither be found easily nor refined freely, but restraints can be applied automatically. Special scattering factors can be used to model complex disorder. Finally chemical information can be included in the form of restraints in order to help the determination of a good model.
Multicollinearity problems are more likely in the refinement of large structures; to some extent more precise and accurate algorithms can help. Also, if the global minimum is less well defined, faster refinement enables more cycles to be carried out, a necessity for good convergence. The efficiency of the algorithms in Crystals have been increased to help address these issues.
Conclusions
Thus, crystal structures are getting larger and their complexity is increasing. Recent developments in precision and speed during the least squares in Crystals is helping the structural scientist to deal with larger structures more efficiently.
Keywords
 Crystals software
 Large molecule crystal structures
 Leastsquares refinement
Background
In the century since William Henry and William Lawrence Bragg reported their Nobel prizewinning studies of single crystal Xray diffraction [1], the determination of chemical structure using the method has become a mature analytical technique. Xray diffractometers have become increasingly available within university research facilities and data measurement, structure solution and refinement are all now carried out by nonexpert users. Whereas fifty years ago it took weeks to collect the data and months to process it, nowadays under the right conditions, an entire dataset can be collected in a matter of minutes and solution and refinement to publication quality can be completed in a similar amount of time. Thus what typically took one year, can now be completed in less than an hour. However, while this is true for the trivial case, there are many examples that are much more complex and scientists continue to push the limits with ever more challenging materials being studied.
Results and discussion
The challenges associated with determining the structure of larger molecules
The challenge begins with obtaining viable crystals: large molecules are often more difficult to crystallise. There are several reasons for this. Firstly, larger molecules generally form a more viscous solution and it is harder for them to organise themselves in a regular way. The archetypal “large molecules” are proteins, however, these contain many small but significant interactions leading to the formation of secondary structure domains, which in turn organise into the tertiary structure of the protein. In contrast, synthetic molecules generally do not have such a high density of strong directing interactions. Secondly, the search for appropriate crystallization conditions can be much more extensive. Proteins are typically crystallised from aqueous solutions; the structures of proteins in the absence of water may be of specialist interest [5], but conditions which reasonably closely resemble their invivo environment are usually desirable. The chemical diversity of synthetic molecules is much larger than for proteins, which means that the ideal solvent for crystallisation of synthetic molecules could come from a list of tens or hundreds, or a mixture thereof, and the search space becomes extensive.
Sample handling
Larger molecules often have multiple possible conformations of similar energy. Sometimes these conformations involve a slight alteration in the orientation of small part of the molecule e.g. a terminal tertiary butyl group, or a bridging phenyl group. This leads to one of the most time consuming aspects of small molecule crystal structure determination: dealing with disorder [6].
Trifluoromethyl substituted phenyl groups often exhibit disorder due to the torsional flexibility about the C—C bond which enables the CF _{3} group to rotate, or librate. The difference in energy between the situation where there are two fluorine atoms above the plane of the aromatic ring and one below, and vice versa, is negligible. Not only is the energy difference between the two states small, but the barrier to rotation between the states is also small, of the order of k _{ B } T (where k _{ B } is the Boltzmann constant and T is the temperature). Thus, at high temperatures it may possible to move between them, with the result that the CF _{3} group is spinning. This is generally referred to as “dynamic disorder”. At low temperature, there may be no kinetic pathway between the two states giving rise to “static disorder”. This means, that whereas dynamic disorder can be reduced by lowering the temperature of the data collection, generally, changing the measurement conditions will not change the degree of static disorder. Sudden cooling can also quench dynamic disorder to static disorder, as the sudden temperature change traps chemical groups in a local thermodynamic minimum. This has implications for the modern practice of using the oildrop flash freezing technique to mount samples at low temperature as this could exacerbate problems with disorder.
It may be possible to use thermal cycling to reduce this problem, however this not necessarily advisable as it can introduce other problems – quench cooling may trap a high temperature phase in a metastable form, whereas thermal cycling will allow it to form a lowtemperature pseudosymmetric phase leading to subsequent difficulties in the analysis. Crystals of large molecules can often be difficult to handle and it can be hard work to find diffracting samples at all. Therefore it would be a courageous crystallographer who is prepared to risk subjecting the sample to further environmental stresses once a crystal is on the diffractometer. Nevertheless, there has been success reported using the socalled “credit card” method with protein crystals: a sample is annealed by warming briefly to room temperature by momentarily blocking the coldstream of the cryodevice using a small piece of plastic or card. Devices have been developed to facilitate this, called “flippers” [7–9].
Restraints
In CRYSTALS, restraints are stored in LIST 16 and are generally of the form
thus a simple distance restraint might be written
The chance of repeated chemical motifs is much higher in large molecules, which means that rather than restraining a given distance to a specific value, sets of distances and angles can be restrained to be equivalent with a defined tolerance. In CRYSTALS these can be applied very efficiently using the command SAME. In this type of restraint, the first group of atoms named is the “target” and all following groups are mapped onto it in order specified. CRYSTALS then uses the connectivity of the first group to work out which distances and angles should be restrained and effectively decomposes the SAME into a series of DISTANCE and ANGLE restraints. For example the restraints below might be used to restrain two pyridine rings:
In this example, these four lines are decomposed into fortyeight separate DISTANCE and ANGLE restraints which map the geometry of the second ring onto the first and additionally uses the inherent mirror symmetry of the pyridine ring to map the geometry of mirrored versions of the first and second ring back onto the first. Thermal similarity restraints and vibrational restraints (which restrain given displacement parameters to be similar and obey the Hirshfeld test [10]) can similarly be simplified using SIMU and DELU between chemically identical groups.
Hydrogen treatment
Hydrogen atoms cannot be refined freely in structures of large molecules due to their small Xray scattering factors, relatively large atomic displacements (resulting from their terminal positions), and overall poor resolution data. However, hydrogen atoms are often visible in the difference Fourier map. In CRYSTALS, hydrogen atoms are added geometrically to sp^{1}, sp^{2} and sp^{3} carbon atoms, with the option of adding those bonded to heteroatoms from the difference map [13]. By default CRYSTALS then refines the hydrogen atom positions with soft restraints in a separate leastsquares cycle prior to inclusion in the final refinement using a riding model. This approach ensures the best possible fit to the measured data while maintaining a sensible geometry. In the extreme case, even with restraints, refining hydrogen atoms may be unstable, in which case, they can be added geometrically or to satisfy hydrogen bonds (using options within the software).
Special shapes
During a refinement, atoms are initially modelled using isotropic displacement parameters and later using a 3×3 tensor describing an anisotropic probability density function (PDF). However, these are often poor approximations of the real displacements of atoms within molecule. For example, the cyclopentadiene ring in a metallocene complex usually librates about an axis described by a vector from the metal to the centre of the ring. Thus the displacement of each individual atom is better described by a curved PDF [14].
In CRYSTALS, three nonstandard parameterisations of atom positions are implemented: the spherical shell; the line; and the torus. The spherical shell is the simplest and is described using the position of its centroid, a standard displacement parameter, U _{iso}, and an absolute magnitude which corresponds to the radius of the shell. The line and torus require additional descriptors for the orientation and, in the case of the line, the radius parameter is replaced by a length. These are the declination (the angle between the line axis or torus normal and the z axis of orthogonal coordinate system used in CRYSTALS) and the azimuth (the angle between the projection of the line axis or torus normal onto the xy plane and the x axis of the orthogonal coordinate).
These sets of parameters are typically used to model disordered anions like \(\mathrm{PF}_{6}^{}\) (where the spherical shell is used to model the six fluorine atoms), disordered solvent in channels (line), and most frequently, the torus which can model librating CF _{3} groups, methylcyclopentadiene ligands or benzene. In most cases, a combined model is used where the “special shape” is partially occupied and has conventional partially occupied atoms embedded in it to reflect the fact that the electron density is nonuniform. The high degree of parameter correlation in these combined models means that restraints are required to ensure convergence and to maintain a physically reasonable model.
Multicollinearity
Multicollinearity can occur in a leastsquares crystallographic model when two or more of the physical parameters are approximately related by a linear function. A simple example would be attempting to refine the relative occupancies of a site containing an unknown mixture of Ga^{+} and Ga^{3+} ions – the Xray data do not contain enough information (or rather it contains too much noise) to distinguish the two species [15]. As a consequence the relative occupancies of the two ions do not have much influence on the fit of the model to the data, however small changes in the data can cause wild fluctuations in the occupancies of the model, and they can easily take on physically meaningless values, provided that the sum of the two species remains at about unity.
In the case of perfectly multicollinear parameters, or combinations of parameters, the normal matrix A A ^{ t } is rank deficient and has no inverse, leading to computational problems carrying out a refinement. Noise in the experimental data makes perfect multicollinearity unlikely in practice and usually an inverse can be computed, though it may be inaccurate and highly sensitive to small errors in the data and small changes in the model.
As the complexity of a crystallographic model increases, the chance of approximate multicollinearity increases: atoms in disordered regions which occupy similar positions in the average structure will have very similar contributions to the structure factor equations; pseudosymmetry can also result in similar scattering contributions from related parts of the model. At the same time, multicollinearity becomes more of a problem as fewer data are included in the model and this effect can be caused by disorder or solvent or poor crystallinity limiting the Xray data resolution.
The increased tendency to disorder and associated resolution reduction in observable Xray data in structures of large molecules can introduce multicollinearity problems. In such refinements, the minimum may be poorly defined and restraints are added to ensure that the model is chemically sensible. The nature of these chemical data is usually to give very local (atom and bondcentred) information related to geometry and displacement parameters rather than giving information that defines a global minimum. As a result it can take several cycles of leastsquares refinement to find the global minimum with respect to all of the restraints. Consider a single alkyl chain: If all the bonds are 0.1 Å too short, restraining each to the correct length would require several cycles of shifting positions until all the bond length restraints are optimally satisfied. This is because the refinement can only attempt to satisfy the local equation for each bond at every step and the refinement has no way of knowing that the terminal atom will ultimately have to move a total of 0.1 Å for every bond in the chain.
This effect, together with the necessarily poorly defined minimum means the structure will generally be slow to converge, indeed, it may never converge and the atoms will tend to shift by small amounts with each additional cycle of leastsquares. The temptation is then to apply more restraints to encourage the refinement to obey preconceived ideas or to apply shiftlimiting restraints to bring the refinement to convergence. In fact it is often better to reduce the number of restraints and apply constraints instead, which remove correlated parameters where possible.
Solvent
Perhaps one of the biggest problems associated with large structures is the inclusion of disordered solvent molecules within the structure. It is often the case that solvent close to an ordered molecule is reasonably well ordered, but becomes increasingly disordered further away [20]. In extreme cases, the lack of longrange order is so severe that individual atoms cannot be identified. This phenomenon is well known in smallmolecule structures, but the problem increases as molecules get larger.
Since invariant and semiinvariant phase relationships of direct methods are derived using the assumption that the electron density is atomic, extreme disorder of solvent may lead to difficulties at the structure solution stage. The advent of chargeflipping methods [21, 22] provides one solution to this problem and SuperFlip [23] has been integrated seamlessly into CRYSTALS [24].
The solvent contribution must be included in the model to avoid systematic errors in Y _{ c } which will lead to a systematic error in the other parameters because the function minimised during refinement (Equation 1 [25]) contains the observed quantity Y _{ o } which has contributions from all atoms, including solvent, however disordered it may be.
The developer Response to Larger Structures
One of the key problems associated with the refinement of large molecules is the amount of time it takes: in addition to the time spent building the model, dealing with restraints, checking the validity of the parameters etc., there is the time spent by the computer actually processing the data. It may initially appear that computer processing is not the major bottleneck compared to the scientist actually modelling and refining a complex structure. However, while these computations of seconds or minutes are running the crystallographer is not free to do another task. It has been reported that users of software will lose attention and want to perform other tasks if forced to wait for more than ten seconds [32]. Studies consistently show that switching focus in this manner significantly impacts performance [33, 34] (regardless of a subject’s gender or degree of control over switching tasks [35]) and that the ability to switch tasks efficiently is actually negatively correlated with the subject’s own belief about their ability to multitask [36].
Improving the speed of the calculation during cycles of leastsquares refinement therefore helps to improve the user experience with bigger structures, reducing frustration and time penalties associated with frequent task switching.

calculation of the derivatives (xsflsx subroutine)

formation of the normal matrix A A ^{ t } (adlhsblock subroutine)

inversion of the normal matrix (xchols subroutine)
In versions of CRYSTALS 14.4x and earlier, the calculation of the derivatives (19 %) is relatively modest, the formation of the normal matrix (72 %) is the most time consuming part and the inversion (5 %) is negligible. However, these results were obtained on a very small structure of 253 parameters. Using a larger structure with 5,173 parameters, the inversion takes 17 % of the time while the calculation of the derivative is negligible at 1.5 %. The formation of the normal matrix remains the most consuming part at 78 %.^{2} The matrix inversion does not scale very well with the number of parameters (Fig. 13) hence the change in proportion. The time of the calculation depends simulatenously on the number of reflections and the number of parameters. The complexity increases linearly with the reflections, but with the parameters there is a much steeper increase. A small increase in the number of parameters results in a large increase in the calculation time, mainly due to the Cholesky inversion [37]. It should be noted that in general the number of parameters is proportional to the number of reflections, so that as larger molecules are studied both will tend to increase. We therefore chose to look for ways to improve the speed of two key parts of the algorithm: the formation of the normal matrix; and the matrix inversion.
Derivatives and normal matrix formation
Reflections are retrieved from a file one by one and the structure factor and its partial derivatives with respect to all parameters being refined are calculated. The derivatives, representing one row of the design matrix A, are used to update the normal matrix N. Updating the normal matrix reflection by reflection is a historical artefact, dating from an era when memory was a scarce resource. Nowadays, this approach is not optimal as it involves a large number of small i/o and memory transfers. Matrix multiplication is heavily bounded by memory bandwidth and less dependent on pure computing power [38]. The main difficulty is to feed data quickly enough to the computing units.
To address this in the latest version of CRYSTALS, reflections are loaded in batches. The derivatives are still calculated sequentially, but they are stored in a temporary partial design matrix block. The accumulation is done once the batch has been processed. This approach has the advantage that it groups similar operations together (i/o, computation of derivatives, matrix multiplications...), while the use of batches instead of the whole dataset maintains the flexibility to work with large sets of data that potentially would not fit in memory. In addition, the matrix multiplication to form the normal matrix is expensive enough to justify a call to an external dynamic library.
Further speed gains are made by optimising the calculation of the normal matrix (A A ^{ t }), a special case of a general matrix multiplication where only half of the values need to be computed. In this case, much lower level modifications involving techniques such as blocking, loop unrolling, cache optimisation, vectorisation, etc are required. Optimisation is challenging as it decreases the readability and flexibility of the source code and the complexity of modern microprocessors makes it difficult for the developers to obtain and maintain skills in these areas. When writing their own software crystallographers can optimise their code to a certain extent, but it is very unlikely that they would match results of high performance libraries given the complexity of modern CPUs and the time they can afford for programming. For example, the simple task of multiplying two matrices efficiently requires a very complicated algorithm [39, 40] which has been refined over many years.
Matrix inversion
The second bottleneck spotted in CRYSTALS is the matrix inversion for solving the linearised problem. This operation is standard in the LAPACK [41], an extension of the Level3 BLAS library. Several methods are available from Cholesky decomposition, L D L ^{ t } factorisation and eigen value decomposition. The home made Cholesky inversion used in CRYSTALS is replaced with an L D L ^{ t } factorisation via the calls to the subroutine SSYTRF and SSYTRI. SSYTRF uses the BunchKaufman diagonal pivoting method for the factorisation. The new method is not only faster, but also scales much better with the number of parameters. Both the MKL and OpenBlas libraries are multithreaded and will exploit multiple CPU cores when present.
The gains in performance gave us the opportunity to implement two features which will improve precision while the time penalties involved are masked by the overall gains.
Firstly, small numerical errors can occur when using single precision in the derivative and structure factor calculations. Therefore the accumulation is done in double precision^{3}.
Secondly, the normal matrix can be preconditioned before inversion. Preconditioning improves the accuracy of the inversion by reducing the loss of precision during floating point operations due to truncation or roundoff errors [42]. Preconditioning is applied as a premultiplication and postmultiplication by the diagonal matrix C on the matrix to invert, N.
Global optimisation and comparisons
An attempt to parallelise the calculation of the derivatives and the normal matrix is in progress, but due to historical constraints in the source code, for now, only a minor part surrounding the normal matrix accumulation has been improved. Parallelisation of the matrix inversion is automatically dealt with within the LAPACK library. Overall, CRYSTALS benefits from having a small number of additional CPU cores available, but scalability is far from ideal.
The old Fortran code in crystals is slowly being converted to modern Fortran to ensure a smooth transition to the future. This will enable further improvements and new optimisations in the future.
Conclusions
We have shown that the size of structures being published is increasing, and that this is directly correlated with the occurrence of disorder. Disorder leads to more manual effort constructing suitable structural models while the larger size of the optimisation problem leads to ever slower refinements. The developments within CRYSTALS, described above, are helping to keep pace with these changes. There are still some areas in CRYSTALS where some improvements are needed. The scalability using multithreads is probably going to give a significant improvement in the near future. Reorganisation of the code to allow a better vectorisation should also give good results on modern CPUs. This is also where the current microprocessors manufacturers improve their CPUs.
Methods
Benchmarks have been realised on a 8 cores Intel®; Xeon®; CPU E52665 with the HyperThreading Technology and Turbo Boost Technology deactivated. This CPU includes the latest 256 bits vector instructions AVX which have been used in crystals during the benchmarks (SHELXL does not seem to use these instructions). The operating system was Centos Linux 6.5 in 64 bits. CRYSTALS is open source and free to use, as such the source code can be downloaded from the website [43]. The new implementation of the software described in this publication is present in revision 5473 and earlier under subversion. CRYSTALS was compiled using the Gnu compiler collection (gcc) 4.9 and the library OpenBlas 0.2.11 [44]. At present, the Windows build of CRYSTALS is using Intel®; Math Kernel Library (Intel®; MKL) which offers similar performance to OpenBlas. SHELXL2014 64 bits was used as provided by the author’s website. The dataset used for the refinements is from Kondratuk D. V. et al. [30] with 49,666 reflections. The number of parameters refined varies between 1,261 to 6,913 depending on the test.
Profiling of the code used Valgrind, an instrumentation framework for building dynamic analysis tools [45] and kcachegrind as a visualisation software. Due to a huge overhead, only small systems up to a few hundreds parameters can be used. For bigger molecules, perf, a Linux profiling with hardware performance counters is used with virtually no impact on the execution time.
Time measurements made on CRYSTALS were made using timers in the source code around the area of interest. For SHELXL, the refinement cycle duration was taken as the difference between two timesamps printed during the refinement.
Endnotes
^{1} An assembly, as defined by the International Union of Crystallography’s CIF dictionary is “a cluster of atoms that show longrange positional disorder but are locally ordered. Within each such cluster of atoms [a group] is used to identify the sites that are simultaneously occupied”.
^{2} These values were not obtained with Valgrind but using perf which has very little overhead performance penalty.
^{3} Multiplication and addition have the same cost either for single precision or double precision. However, vector units can hold half the number of double precision numbers hence the performance penalty.
Declarations
Acknowledgements
The authors gratefully acknowledge support from EPSRC grant EP/K013009/1 (RIC and PP) and Diamond Light Source for a continuing award of beamtime (MT9981) without which none of the refinements underpinning this work would have been possible.
Authors’ Affiliations
References
 Bragg WH, Bragg WL. Xrays and Crystal Structure. London: G. Bell and Sons; 1915.Google Scholar
 Allen FH. The cambridge structural database: a quarter of a million crystal structures and rising. Acta Crystallographica. 2002; B58:380–8.View ArticleGoogle Scholar
 Bruno IJ, Cole JC, Edgington PR, Kessler M, Macrae CF, McCabe P, et al. New software for searching the Cambridge structural database and visualizing crystal structures Acta Crystallographica. 2002; B58:389–97.Google Scholar
 Betteridge PW, Carruthers JR, Cooper RI, Prout K, Watkin DJ. CRYSTALS version 12: software for guided crystal structure analysis. J Appl Crystallography. 2003; 36:1487.Google Scholar
 Fitzpatrick PA, Steinmetz ACU, Ringe D, Kilbanov AM. Enzyme crystal structure in a neat organic solvent. Proc Natl Acad Sci. 1993; 90:8653–657.View ArticleGoogle Scholar
 Müller P. Practical suggestions for better crystal structures. Crystallogr Rev. 2009; 15(1):57–83.View ArticleGoogle Scholar
 Harp JM, Timm DE, Bunick GJ. Macromolecular crystal annealing: overcoming increased mosaicity associated with cryocrystallography. Acta Crystallogr Section D. 1998; 54(4):622–8.View ArticleGoogle Scholar
 Harp JM, Hanson BL, Timm DE, Bunick GJ. Macromolecular crystal annealing: evaluation of techniques and variables. Acta Crystallogr Section D. 1999; 55(7):1329–1334.View ArticleGoogle Scholar
 Hanson BL, Harp JM, Bunick GJ. The welltempered protein crystal: annealing macromolecular crystals. Methods Enzymol. 2003; 368:217–35.Google Scholar
 Hirshfeld FL. Can xray data distinguish bonding effects from vibrational smearing?Acta Crystallogr Section A. 1976; 32:239–44.View ArticleGoogle Scholar
 Davis NKS, Thompson AL, Anderson HL. Bisanthracene fused porphyrins: Synthesis, crystal structure, and nearir absorption. Org Lett. 2010; 12(9):2124–127.View ArticleGoogle Scholar
 Cooper RI, Thorn A, Watkin DJ. Crystals enhancements: asymmetric restraints. J Appl Crystallogr. 2012; 45:1057–60.View ArticleGoogle Scholar
 Cooper RI, Thompson AL, Watkin DJ. CRYSTALS enhancements: dealing with hydrogen atoms in refinement. J Appl Crystallogr. 2010; 43(5 Part 1):1100–7. doi:10.1107/S0021889810025598.View ArticleGoogle Scholar
 Reilly AM, McLean KR, Morrison CA, Rankin DWH. Using moleculardynamics simulations to understand and improve the treatment of anharmonic vibrations. I. Study of positional parameters. Acta Crystallogr. 2011; A67:336–45.View ArticleGoogle Scholar
 Sutton KJ, Barnett SA, Christensen KE, Nowell H, Thompson AL, Allan DR, et al. DetOx: a program for determining anomalous scattering factors of mixedoxidationstate species. J Synchrotron Radiat. 2013; 20(1):200–4. doi:10.1107/S0909049512044007.View ArticleGoogle Scholar
 Fattah J, Twyman JM, Heyes SJ, Watkin DJ, Edwards AJ, Prout K, et al. Combination of cross polarization/magic angle spinning NMR and xray crystallography: structure and dynamics in a lowsymmetry molecular crystal, potassium penicillin V. J Am Chem Soc. 1993; 115:5636–650.View ArticleGoogle Scholar
 Wendeler M, Fattah J, Twyman JM, Edwards AJ, Dobson CM, Heyes SJ, et al. Combination of cp/mas nmr spectroscopy and xray crystallography: Structure and dynamics in molecular crystals of hydrogen, lithium, sodium, rubidium, and cesium penicillin v. J Am Chem Soc. 1997; 119(41):9793–803.View ArticleGoogle Scholar
 Anderson KM, Afarinkia K, Yu HW, Goeta AE, Steed JW. When Z ^{′}=2 Is Better than Z ^{′}=1 – Supramolecular Centrosymmetric HydrogenBonded Dimers in Chiral Systems. Cryst Growth Des. 2006; 6(9):2109–113.View ArticleGoogle Scholar
 Kilah NL, Wise MD, Serpell CJ, Thompson AL, White NG, Christensen KE, et al. Enhancement of anion recognition exhibited by a halogenbonding rotaxane host system. J Am Chem Soc. 2010; 132(34):11893–5.View ArticleGoogle Scholar
 Davis NKS, Thompson AL, Anderson HL. A porphyrin fused to four anthracenes. J Am Chem Soc. 2011; 133(1):30–1.View ArticleGoogle Scholar
 Oszlányi G, Süto A. Practical suggestions for better crystal structures. Acta Crystallogr. 2004; A60:134–41.View ArticleGoogle Scholar
 Palatinus L. The chargeflipping algorithm in crystallography. Acta Crystallogr. 2013; B69(1):1–16.View ArticleGoogle Scholar
 Palatinus L, Chapuis G. SUPERFLIP – a computer program for the solution of crystal structures by charge flipping in arbitrary dimensions. J Appl Crystallogr. 2007; 40:786–90.View ArticleGoogle Scholar
 Palatinus L, van der Lee A. Symmetry determination following structure solution in P1. J Appl Crystallogr. 2008; 41:975–84.View ArticleGoogle Scholar
 Watkin DJ. Structure refinement: some background theory and practical strategies. J Appl Crystallogr. 2008; 41:491–522.View ArticleGoogle Scholar
 van der Sluis P, Spek AL. BYPASS: an effective method for the refinement of crystal structures containing disordered solvent regions. Acta Crystallogr. 1990; A46:194–201.View ArticleGoogle Scholar
 Spek A. Singlecrystal structure validation with the program PLATON. J Appl Crystallogr. 2003; 36:7–13.View ArticleGoogle Scholar
 Sheldrick GM. A short history of SHELX. Acta Crystallogr Section A. 2008; 64(1):112–22. doi:10.1107/S0108767307043930.View ArticleGoogle Scholar
 Spek AL. PLATON SQUEEZE: a tool for the calculation of the disordered solvent contribution to the calculated structure factors. Acta Crystallogr Section C. 2015; 71(1):9–18. doi:10.1107/S2053229614024929.View ArticleGoogle Scholar
 Kondratuk DV, Sprafke JK, O’Sullivan MC, Perdigao LMA, Saywell A, Malfois M, et al.Verniertemplated synthesis, crystal structure, and supramolecular chemistry of a 12porphyrin nanoring. Chem A European J. 2014. doi:10.1002/chem.201403714.Google Scholar
 Sprafke JK, Kondratuk DV, Wykes M, Thompson AL, Hoffmann M, Drevinskas R, et al. Beltshaped πsystems: Relating geometry to electronic structure in a sixporphyrin nanoring. J Am Chem Soc. 2011; 133(43):17262–73. doi:10.1021/ja2045919.View ArticleGoogle Scholar
 Nielsen J. Usability Engineering. Interactive Technologies. Oval road, London NW1 7DX: Academic Press; 1994. http://books.google.co.uk/books?id=DBOowF7LqIQC.Google Scholar
 Monsell S. Task switching. Trends Cogn Sci. 2003; 7(3):134–40. doi:10.1016/S13646613(03)000287.View ArticleGoogle Scholar
 Gopher D, Armony L, Greenshpan Y. Switching tasks and attention policies. J Exp Psychol: General. 2000; 129(3):308–39.View ArticleGoogle Scholar
 Buser T, Peter N. Multitasking. Exp Econ. 2012; 15(4):641–55. doi:10.1007/s1068301293188.View ArticleGoogle Scholar
 Sanbonmatsu DM, Strayer DL, MedeirosWard N, Watson JM. Who multitasks and why? Multitasking ability, perceived multitasking ability, impulsivity, and sensation seeking. PLoS ONE. 2013; 8(1):54402.Google Scholar
 Rollett JS. Computing Methods in Crystallography, 1st edn. Headington Hill Hall, Oxford: Pergamon Press; 1965.Google Scholar
 Gunnels J, Gustavson F, Pingali K, Yotov K. Applied Parallel Computing. State of the Art in Scientific Computing In: Kågsträm B, Elmroth E, Dongarra J, Waśniewski J, editors. Lecture Notes in Computer Science. Berlin Heidelberg: Springer: 2007. p. 919–28.Google Scholar
 Goto K, van de Geijn R. Highperformance implementation of the level3 blas. ACM Trans Math Softw. 2008; 35(1):4–1414. doi:10.1145/1377603.1377607.View ArticleGoogle Scholar
 Goto K, van de Geijn RA. Anatomy of highperformance matrix multiplication. ACM Trans Math Softw. 2008; 34(3):12–11225. doi:10.1145/1356052.1356053.View ArticleGoogle Scholar
 Anderson E, Bai Z, Bischof C, Blackford S, Demmel J, Dongarra J, et al. LAPACK Users’ Guide, 3rd edn. Philadelphia, PA: Society for Industrial and Applied Mathematics; 1999.View ArticleGoogle Scholar
 Prince E, Boggs PT. International Tables for Crystallography, Volume C: Mathematical, physical and chemical tables, 1st edn. Dordrecht, Boston, London: Kluwer Academic Publishers; 1995, p. 598. Chap. 8.1.3.3 Conditioning.Google Scholar
 Chemical Crystallography, Department of chemistry, University of Oxford. http://www.xtl.ox.ac.uk.
 Wang Q, Zhang X, Zhang Y, Yi Q. Augem: Automatically generate high performance dense linear algebra kernels on x86 cpus. In: Proceedings of the International Conference on High Performance Computing, Networking, Storage and Analysis. SC ’13. New York, NY, USA: ACM: 2013. p. 25–12512. doi:10.1145/2503210.2503219.Google Scholar
 Nethercote N, Seward J. Valgrind: A framework for heavyweight dynamic binary instrumentation. In: Proceedings of the 2007 ACM SIGPLAN Conference on Programming Language Design and Implementation. PLDI ’07. New York, NY, USA: ACM: 2007. p. 89–100. doi:10.1145/1250734.1250746.Google Scholar
Copyright
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.