Synthesis and Characterization of Two Novel Naphthalenediimide/Zinc Phosphonate Crystalline Materials Precipitated from Different Solvents.

Hybrid naphthalenediimide/zinc phosphonate materials (NDI/Zn) were prepared by mixing solutions of N,N'-bis(2-phosphonoethyl)-1,4,5,8-naphthalenediimide (PNDI) and zinc nitrate, resulting in the precipitation of the desired compounds. Samples precipitated from water and N,N-dimethylformamide (DMF) were produced. The obtained samples had the expected elemental composition, and the presence of naphthalenediimides (NDI) was ascertained by infrared and UV-visible spectroscopy. All the samples were crystalline, according to powder X-ray diffraction. Nitrogen adsorption isotherms showed the presence of porosity in the NDI/Zn samples. Mesopores with a diameter = 4.1 nm were present in the sample from DMF, with total pore volume reaching 0.13 cm3/g.

Predicting crystal form stability under real-world conditions.

The physicochemical properties of molecular crystals, such as solubility, stability, compactability, melting behaviour and bioavailability, depend on their crystal form1. In silico crystal form selection has recently come much closer to realization because of the development of accurate and affordable free-energy calculations2-4. Here we redefine the state of the art, primarily by improving the accuracy of free-energy calculations, constructing a reliable experimental benchmark for solid-solid free-energy differences, quantifying statistical errors for the computed free energies and placing both hydrate crystal structures of different stoichiometries and anhydrate crystal structures on the same energy landscape, with defined error bars, as a function of temperature and relative humidity. The calculated free energies have standard errors of 1-2 kJ mol-1 for industrially relevant compounds, and the method to place crystal structures with different hydrate stoichiometries on the same energy landscape can be extended to other multi-component systems, including solvates. These contributions reduce the gap between the needs of the experimentalist and the capabilities of modern computational tools, transforming crystal structure prediction into a more reliable and actionable procedure that can be used in combination with experimental evidence to direct crystal form selection and establish control5.

Hydrogen bonding patterns and C-H...π interactions in the structure of the antiparkinsonian drug (R)-rasagiline mesylate determined using laboratory and synchrotron X-ray powder diffraction data.

The structure of (R)-rasagiline mesylate [(R)-RasH+·Mes-], an active pharmaceutical ingredient used to treat Parkinson's disease, is presented. The structure was determined from laboratory and synchrotron powder diffraction data, refined using the Rietveld method, and validated and optimized using dispersion-corrected DFT calculations. The unit-cell parameters obtained in both experiments are in good agreement and the refinement with both datasets converged to good agreement factors. The final parameters obtained from laboratory data were a = 5.4905 (8), b = 6.536 (2), c = 38.953 (3) Å, V = 1398.0 (4) Å3 and from synchrotron powder data were a = 5.487530 (10) Å, b = 6.528939 (12) Å, c = 38.94313 (9) Å, V = 1395.245 (5) Å3 with Z = 4 and space group P212121. Preferred orientation was properly accounted for using the synchrotron radiation data, leading to a March-Dollase parameter of 1.140 (1) instead of the 0.642 (1) value obtained from laboratory data. In the structure, (R)-RasH+ moieties form layers parallel to the ab plane connected by mesylate ions through N-H...O and C-H...O hydrogen bonds. These layers stack along the c axis and are further connected by C-H...π interactions. Hirshfeld surface analysis and fingerprint plot calculations indicate that the main interactions are: H...H (50.9%), H...C/C...H (27.1%) and H...O/O...H (21.1%).

Crystal structure from X-ray powder diffraction data, DFT-D calculation, Hirshfeld surface analysis, and energy frameworks of (RS)-trichlorme-thia-zide.

The structure of racemic (RS)-trichlorme-thia-zide [systematic name: (RS)-6-chloro-3-(di-chloro-meth-yl)-1,1-dioxo-3,4-di-hydro-2H-1λ6,2,4-benzo-thia-di-azine-7-sulfonamide], C8H8Cl3N3O4S2 (RS-TCMZ), a diuretic drug used in the treatment of oedema and hypertension, was determined from laboratory X-ray powder diffraction data using DASH [David et al. (2006 ▸). J. Appl. Cryst. 39, 910-915.], refined by the Rietveld method with TOPAS-Academic [Coelho (2018 ▸). J. Appl. Cryst. 51, 210-218], and optimized using DFT-D calculations. The extended structure consists of head-to-tail dimers connected by π-π inter-actions which, in turn, are connected by C-Cl⋯π inter-actions. They form chains propagating along [101], further connected by N-H⋯O hydrogen bonds to produce layers parallel to the ac plane that stack along the b-axis direction, connected by additional N-H⋯O hydrogen bonds. The Hirshfeld surface analysis indicates a major contribution of H⋯O and H⋯Cl inter-actions (32.2 and 21.7%, respectively). Energy framework calculations confirm the major contribution of electrostatic inter-actions (E elec) to the total energy (E tot). A comparison with the structure of S-TCMZ is also presented.

A jumping crystal predicted with molecular dynamics and analysed with TLS refinement against powder diffraction data.

By running a temperature series of molecular dynamics (MD) simulations starting from the known low-temperature phase, the experimentally observed phase transition in a 'jumping crystal' was captured, thereby providing a prediction of the unknown crystal structure of the high-temperature phase and clarifying the phase-transition mechanism. The phase transition is accompanied by a discontinuity in two of the unit-cell parameters. The structure of the high-temperature phase is very similar to that of the low-temperature phase. The anisotropic displacement parameters calculated from the MD simulations readily identified libration as the driving force behind the phase transition. Both the predicted crystal structure and the phase-transition mechanism were verified experimentally using TLS (translation, libration, screw) refinement against X-ray powder diffraction data.

How many ritonavir cases are there still out there?

Based on a thorough and critical analysis of the commercial crystal structure prediction studies of 41 pharmaceutical compounds, we conclude that for between 15 and 45% of all small-molecule drugs currently on the market the most stable experimentally observed polymorph is not the thermodynamically most stable crystal structure and that the appearance of the latter is kinetically hindered.

Crystal structure analysis of a star-shaped triazine compound: a combination of single-crystal three-dimensional electron diffraction and powder X-ray diffraction.

The solid-state structure of star-shaped 2,4,6-tris{(E)-2-[4-(dimethylamino)-phenyl]ethenyl}-1,3,5-triazine is determined from a powder sample by exploiting the respective strengths of single-crystal three-dimensional electron diffraction and powder X-ray diffraction data. The unit-cell parameters were determined from single crystal electron diffraction data. Using this information, the powder X-ray diffraction data were indexed, and the crystal structure was determined from the powder diffraction profile. The compound crystallizes in a noncentrosymmetric space group, P212121. The molecular conformation in the crystal structure was used to calculate the molecular dipole moment of 3.22 Debye, which enables the material to show nonlinear optical effects.

Validation of missed space-group symmetry in X-ray powder diffraction structures with dispersion-corrected density functional theory.

More than 600 molecular crystal structures with correct, incorrect and uncertain space-group symmetry were energy-minimized with dispersion-corrected density functional theory (DFT-D, PBE-D3). For the purpose of determining the correct space-group symmetry the required tolerance on the atomic coordinates of all non-H atoms is established to be 0.2 Å. For 98.5% of 200 molecular crystal structures published with missed symmetry, the correct space group is identified; there are no false positives. Very small, very symmetrical molecules can end up in artificially high space groups upon energy minimization, although this is easily detected through visual inspection. If the space group of a crystal structure determined from powder diffraction data is ambiguous, energy minimization with DFT-D provides a fast and reliable method to select the correct space group.

The application of tailor-made force fields and molecular dynamics for NMR crystallography: a case study of free base cocaine.

Motional averaging has been proven to be significant in predicting the chemical shifts in ab initio solid-state NMR calculations, and the applicability of motional averaging with molecular dynamics has been shown to depend on the accuracy of the molecular mechanical force field. The performance of a fully automatically generated tailor-made force field (TMFF) for the dynamic aspects of NMR crystallography is evaluated and compared with existing benchmarks, including static dispersion-corrected density functional theory calculations and the COMPASS force field. The crystal structure of free base cocaine is used as an example. The results reveal that, even though the TMFF outperforms the COMPASS force field for representing the energies and conformations of predicted structures, it does not give significant improvement in the accuracy of NMR calculations. Further studies should direct more attention to anisotropic chemical shifts and development of the method of solid-state NMR calculations.

Local structure in the disordered solid solution of cis- and trans-perinones.

The cis- and trans-isomers of the polycyclic aromatic compound perinone, C26H12N4O2, form a solid solution (Vat Red 14). This solid solution is isotypic to the crystal structures of cis-perinone (Pigment Red 194) and trans-perinone (Pigment Orange 34) and exhibits a combined positional and orientational disorder: In the crystal, each molecular position is occupied by either a cis- or trans-perinone molecule, both of which have two possible molecular orientations. The structure of cis-perinone exhibits a twofold orientational disorder, whereas the structure of trans-perinone is ordered. The crystal structure of the solid solution was determined by single-crystal X-ray analysis. Extensive lattice-energy minimizations with force-field and DFT-D methods were carried out on combinatorially complete sets of ordered models. For the disordered systems, local structures were calculated, including preferred local arrangements, ordering lengths, and probabilities for the arrangement of neighbouring molecules. The superposition of the atomic positions of all energetically favourable calculated models corresponds well with the experimentally determined crystal structures, explaining not only the atomic positions, but also the site occupancies and anisotropic displacement parameters.

Crystallographic and Dynamic Aspects of Solid-State NMR Calibration Compounds: Towards ab Initio NMR Crystallography.

The excellent results of dispersion-corrected density functional theory (DFT-D) calculations for static systems have been well established over the past decade. The introduction of dynamics into DFT-D calculations is a target, especially for the field of molecular NMR crystallography. Four (13) C ss-NMR calibration compounds are investigated by single-crystal X-ray diffraction, molecular dynamics and DFT-D calculations. The crystal structure of 3-methylglutaric acid is reported. The rotator phases of adamantane and hexamethylbenzene at room temperature are successfully reproduced in the molecular dynamics simulations. The calculated (13) C chemical shifts of these compounds are in excellent agreement with experiment, with a root-mean-square deviation of 2.0 ppm. It is confirmed that a combination of classical molecular dynamics and DFT-D chemical shift calculation improves the accuracy of calculated chemical shifts.

Structure of Pigment Yellow 181 dimethylsulfoxide N-methyl-2-pyrrolidone (1:1:1) solvate from XRPD + DFT-D.

With only a 2.6 Šresolution laboratory powder diffraction pattern of the θ phase of Pigment Yellow 181 (P.Y. 181) available, crystal-structure solution and Rietveld refinement proved challenging; especially when the crystal structure was shown to be a triclinic dimethylsulfoxide N-methyl-2-pyrrolidone (1:1:1) solvate. The crystal structure, which in principle has 28 possible degrees of freedom, was determined in three stages by a combination of simulated annealing, partial Rietveld refinement with dummy atoms replacing the solvent molecules and further simulated annealing. The θ phase not being of commercial interest, additional experiments were not economically feasible and additional dispersion-corrected density functional theory (DFT-D) calculations were employed to confirm the correctness of the crystal structure. After the correctness of the structure had been ascertained, the bond lengths and valence angles from the DFT-D minimized crystal structure were fed back into the Rietveld refinement as geometrical restraints (`polymorph-dependent restraints') to further improve the details of the crystal structure; the positions of the H atoms were also taken from the DFT-D calculations. The final crystal structure is a layered structure with an elaborate network of hydrogen bonds.

One hydrogen bond does not a separation make, or does it? Resolution of amines by diacetoneketogulonic acid.

Diacetoneketogulonic acid was used to separate primary amines from their racemic modifications and the selectivity of the acid was rationalized by lattice energy calculations and analyzing the weak interactions around the captured amines.

Validation of molecular crystal structures from powder diffraction data with dispersion-corrected density functional theory (DFT-D).

In 2010 we energy-minimized 225 high-quality single-crystal (SX) structures with dispersion-corrected density functional theory (DFT-D) to establish a quantitative benchmark. For the current paper, 215 organic crystal structures determined from X-ray powder diffraction (XRPD) data and published in an IUCr journal were energy-minimized with DFT-D and compared to the SX benchmark. The on average slightly less accurate atomic coordinates of XRPD structures do lead to systematically higher root mean square Cartesian displacement (RMSCD) values upon energy minimization than for SX structures, but the RMSCD value is still a good indicator for the detection of structures that deserve a closer look. The upper RMSCD limit for a correct structure must be increased from 0.25 Šfor SX structures to 0.35 Šfor XRPD structures; the grey area must be extended from 0.30 to 0.40 Å. Based on the energy minimizations, three structures are re-refined to give more precise atomic coordinates. For six structures our calculations provide the missing positions for the H atoms, for five structures they provide corrected positions for some H atoms. Seven crystal structures showed a minor error for a non-H atom. For five structures the energy minimizations suggest a higher space-group symmetry. For the 225 SX structures, the only deviations observed upon energy minimization were three minor H-atom related issues. Preferred orientation is the most important cause of problems. A preferred-orientation correction is the only correction where the experimental data are modified to fit the model. We conclude that molecular crystal structures determined from powder diffraction data that are published in IUCr journals are of high quality, with less than 4% containing an error in a non-H atom.

Distinguishing tautomerism in the crystal structure of (Z)-N-(5-ethyl-2,3-dihydro-1,3,4-thiadiazol-2-ylidene)-4-methylbenzenesulfonamide using DFT-D calculations and (13)C solid-state NMR.

The crystal structure of the title compound, C11H13N3O2S2, has been determined previously on the basis of refinement against laboratory powder X-ray diffraction (PXRD) data, supported by comparison of measured and calculated (13)C solid-state NMR spectra [Hangan et al. (2010). Acta Cryst. B66, 615-621]. The molecule is tautomeric, and was reported as an amine tautomer [systematic name: N-(5-ethyl-1,3,4-thiadiazol-2-yl)-p-toluenesulfonamide], rather than the correct imine tautomer. The protonation site on the molecule's 1,3,4-thiadiazole ring is indicated by the intermolecular contacts in the crystal structure: N-H...O hydrogen bonds are established at the correct site, while the alternative protonation site does not establish any notable intermolecular interactions. The two tautomers provide essentially identical Rietveld fits to laboratory PXRD data, and therefore they cannot be directly distinguished in this way. However, the correct tautomer can be distinguished from the incorrect one by previously reported quantitative criteria based on the extent of structural distortion on optimization of the crystal structure using dispersion-corrected density functional theory (DFT-D) calculations. Calculation of the (13)C SS-NMR spectrum based on the correct imine tautomer also provides considerably better agreement with the measured (13)C SS-NMR spectrum.

On the correlation between hydrogen bonding and melting points in the inositols.

Inositol, 1,2,3,4,5,6-hexahydroxycyclohexane, exists in nine stereoisomers with different crystal structures and melting points. In a previous paper on the relationship between the melting points of the inositols and the hydrogen-bonding patterns in their crystal structures [Simperler et al. (2006 ▶). CrystEngComm 8, 589], it was noted that although all inositol crystal structures known at that time contained 12 hydrogen bonds per molecule, their melting points span a large range of about 170 °C. Our preliminary investigations suggested that the highest melting point must be corrected for the effect of molecular symmetry, and that the three lowest melting points may need to be revised. This prompted a full investigation, with additional experiments on six of the nine inositols. Thirteen new phases were discovered; for all of these their crystal structures were examined. The crystal structures of eight ordered phases could be determined, of which seven were obtained from laboratory X-ray powder diffraction data. Five additional phases turned out to be rotator phases and only their unit cells could be determined. Two previously unknown melting points were measured, as well as most enthalpies of melting. Several previously reported melting points were shown to be solid-to-solid phase transitions or decomposition points. Our experiments have revealed a complex picture of phases, rotator phases and phase transitions, in which a simple correlation between melting points and hydrogen-bonding patterns is not feasible.

Structures of cefradine dihydrate and cefaclor dihydrate from DFT-D calculations.

The crystal structure of cefradine dihydrate, C16H19N3O4S·2H2O, is considered in the pharmaceutical sciences to be the epitome of an isolated-site hydrate. The structure from single-crystal X-ray data was described in 1976, but atomic coordinates were not published. The atomic coordinates are determined here by combining the information available from the published single-crystal data with a dispersion-corrected density functional theory (DFT-D) method that has been validated to reproduce molecular crystal structures very accurately. Additional proof for the correctness of the structure comes from comparison with cefaclor dihydrate, C15H14ClN3O4S·2H2O, which is isomorphous and for which more complete single-crystal data are available. H-atom positions have not previously been published for either compound. The DFT-D calculations confirm that both cefradine and cefaclor are present in the zwitterionic form in the two dihydrate structures. A potential ambiguity concerning the orientation of the cyclohexadienyl ring in cefradine dihydrate is also clarified, and on the basis of the calculated energies it is shown that disorder should not be expected at room temperature. The DFT-D methods can be applied to recover full structural data in cases where only partial information is available, and where it may not be possible or desirable to obtain new experimental data.

Reinterpretation of the monohydrate of clarithromycin from X-ray powder diffraction data as a trihydrate.

Noguchi, Fujiki, Iwao, Miura & Itai [Acta Cryst. (2012), E68, o667-o668] recently reported the crystal structure of clarithromycin monohydrate from synchrotron X-ray powder diffraction data. Voids in the crystal structure suggested the possible presence of two more water molecules. After successful location of the two additional water molecules, the Rietveld refinement still showed minor problems. These were resolved by noticing that one of the chiral centres in the molecule had been inverted. The corrected crystal structure of clarithromycin trihydrate, refined against the original data, is now reported. Dispersion-corrected density functional theory calculations were used to check the final crystal structure and to position the H atoms.