The results of the sixth blind test of organic crystal structure prediction methods are presented and discussed, highlighting progress for salts, hydrates and bulky flexible molecules, as well as on-going challenges.

Reliable prediction of molecular crystal energetics is a vital goal for computational chemistry. Here we show that accurate results can be obtained from a monomer-based many-body expansion truncated at the two-body level, with the monomer and dimer calculations suitably embedded in a model of the crystalline environment. By including the two dominant effects--electrostatics and exchange-repulsion--we are able to capture the important nonadditive terms in the energy, and approach very closely results from full periodic second-order Møller-Plesset calculations. The advantage of the current scheme is that extension to coupled-cluster and explicitly correlated F12 methods is completely straightforward. We demonstrate the approach through calculations on carbon dioxide, hydrogen fluoride, and ice XIh and XIc. In accord with previous studies, we find these two ice polymorphs to be very close in energy, with our periodic coupled-cluster single double triple-F12 calculation giving the hexagonal structure more stable by around 0.3 kJ mol(-1).

Generating sets of trial structures that sample the configurational space of crystal packing possibilities is an essential step in the process of ab initio crystal structure prediction (CSP). One effective methodology for performing such a search relies on low-discrepancy, quasi-random sampling, and our implementation of such a search for molecular crystals is described in this paper. Herein we restrict ourselves to rigid organic molecules and, by considering their geometric properties, build trial crystal packings as starting points for local lattice energy minimization. We also describe a method to match instances of the same structure, which we use to measure the convergence of our packing search toward completeness. The use of these tools is demonstrated for a set of molecules with diverse molecular characteristics and as representative of areas of application where CSP has been applied. An important finding is that the lowest energy crystal structures are typically located early and frequently during a quasi-random search of phase space. It is usually the complete sampling of higher energy structures that requires extended sampling. We show how the procedure can first be refined, through targetting the volume of the generated crystal structures, and then extended across a range of space groups to make a full CSP search and locate experimentally observed and lists of hypothetical polymorphs. As the described method has also been created to lie at the base of more involved approaches to CSP, which are being developed within the Global Lattice Energy Explorer (Glee) software, a few of these extensions are briefly discussed.

Structures of the α and β phases of resorcinol, a major commodity chemical in the pharmaceutical, agrichemical, and polymer industries, were the first polymorphic pair of molecular crystals solved by X-ray analysis. It was recently stated that "no additional phases can be found under atmospheric conditions" (Druzbicki, K. et al. J. Phys. Chem. B 2015, 119, 1681). Herein is described the growth and structure of a new ambient pressure phase, ε, through a combination of optical and X-ray crystallography and by computational crystal structure prediction algorithms. α-Resorcinol has long been a model for mechanistic crystal growth studies from both solution and vapor because prisms extended along the polar axis grow much faster in one direction than in the opposite direction. Research has focused on identifying the absolute sense of the fast direction-the so-called "resorcinol riddle"-with the aim of identifying how solvent controls crystal growth. Here, the growth velocity dissymmetry in the melt is analyzed for the β phase. The ε phase only grows from the melt, concomitant with the β phase, as polycrystalline, radially growing spherulites. If the radii are polar, then the sense of the polar axis is an essential feature of the form. Here, this determination is made for spherulites of β resorcinol (ε, point symmetry 222, does not have a polar axis) with additives that stereoselectively modify growth velocities. Both β and ε have the additional feature that individual radial lamellae may adopt helicoidal morphologies. We correlate the appearance of twisting in β and ε with the symmetry of twist-inducing additives.

We show how an embedded many-body expansion (EMBE) can be used to calculate accurate ab initio energies of water clusters and ice structures using wavefunction-based methods. We use the EMBE described recently by Bygrave et al. [J. Chem. Phys. 137, 164102 (2012)], in which the terms in the expansion are obtained from calculations on monomers, dimers, etc., acted on by an approximate representation of the embedding field due to all other molecules in the system, this field being a sum of Coulomb and exchange-repulsion fields. Our strategy is to separate the total energy of the system into Hartree-Fock and correlation parts, using the EMBE only for the correlation energy, with the Hartree-Fock energy calculated using standard molecular quantum chemistry for clusters and plane-wave methods for crystals. Our tests on a range of different water clusters up to the 16-mer show that for the second-order Møller-Plesset (MP2) method the EMBE truncated at 2-body level reproduces to better than 0.1 mE(h)/monomer the correlation energy from standard methods. The use of EMBE for computing coupled-cluster energies of clusters is also discussed. For the ice structures Ih, II, and VIII, we find that MP2 energies near the complete basis-set limit reproduce very well the experimental values of the absolute and relative binding energies, but that the use of coupled-cluster methods for many-body correlation (non-additive dispersion) is essential for a full description. Possible future applications of the EMBE approach are suggested.