Role of Molecular, Crystal, and Surface Chemistry in Directing the Crystallization of Entacapone Polymorphs on the Au(111) Template Surface.

The pharmaceutical compound entacapone ((E)-2-cyano-3-(3,4-dihydroxy-5-nitrophenyl)-N,N-diethylprop-2-enamide) is important in the treatment of Parkinson's disease, exhibiting interesting polymorphic behavior upon crystallization from solution. It consistently produces its stable form A with a uniform crystal size distribution on the surface of an Au(111) template while concomitantly forming its metastable form D within the same bulk solution. Molecular modeling using empirical atomistic force-fields reveals more complex molecular and intermolecular structures for form D compared to form A, with the crystal chemistry of both polymorphs being dominated by van der Waals and π-π stacking interactions with lower contributions (ca. 20%) from hydrogen bonding and electrostatic interactions. Comparative lattice energies and convergence for the polymorphs are consistent with the observed concomitant polymorphic behavior. Synthon characterization reveals an elongated needle-like morphology for form D crystals in contrast to the more equant form A crystals with the surface chemistry of the latter exposing the molecules' cyano groups on its {010} and {011} habit faces. Density functional theory modeling of surface adsorption reveals preferential interactions between Au and the synthon GA interactions of form A on the Au surface. Molecular dynamics modeling of the entacapone/gold interface reveals the entacapone molecular structure within the first adsorbed layer to show nearly identical interaction distances, for both the molecules within form A or D with respect to the Au surface, while in the second and third layers when entacapone molecule-molecule interactions overtake the interactions between those of molecule-Au, the intermolecular structures are found to be closer to the form A structure than form D. In these layers, synthon GA (form A) could be reproduced with just two small azimuthal rotations (5° and 15°) whereas the closest alignment to a form D synthon requires larger azimuthal rotations (15° and 40°). The cyano functional group interactions with the Au template dominate interfacial interactions with these groups being aligned parallel to the Au surface and with nearest neighbor distances to Au atoms more closely matching those in form A than form D. The overall polymorph direction pathway thus encompasses consideration of molecular, crystal, and surface chemistry factors.

Evaluation of Force-Field Calculations of Lattice Energies on a Large Public Dataset, Assessment of Pharmaceutical Relevance, and Comparison to Density Functional Theory.

Crystal lattice energy is a key property affecting the ease of processing pharmaceutical materials during manufacturing, as well as product performance. We present an extensive comparison of 324 force-field protocols for calculating the lattice energies of single component, organic molecular crystals (further restricted to Z' less than or equal to one), corresponding to a wide variety of force-fields (DREIDING, Universal, CVFF, PCFF, COMPASS, COMPASSII), optimization routines, and other variations, which could be implemented as part of an automated workflow using the industry standard Materials Studio software. All calculations were validated using a large new dataset (SUB-BIG), which we make publicly available. This dataset comprises public domain sublimation data, from which estimated experimental lattice energies were derived, linked to 235 molecular crystals. Analysis of pharmaceutical relevance was performed according to two distinct methods based upon (A) public and (B) proprietary data. These identified overlapping subsets of SUB-BIG comprising (A) 172 and (B) 63 crystals, of putative pharmaceutical relevance, respectively. We recommend a protocol based on the COMPASSII force field for lattice energy calculations of general organic or pharmaceutically relevant molecular crystals. This protocol was the most highly ranked prior to subsetting and was either the top ranking or amongst the top 15 protocols (top 5%) following subsetting of the dataset according to putative pharmaceutical relevance. Further analysis identified scenarios where the lattice energies calculated using the recommended force-field protocol should either be disregarded (values greater than or equal to zero and/or the messages generated by the automated workflow indicate extraneous atoms were added to the unit cell) or treated cautiously (values less than or equal to -249 kJ/mol), as they are likely to be inaccurate. Application of the recommended force-field protocol, coupled with these heuristic filtering criteria, achieved an root mean-squared error (RMSE) around 17 kJ/mol (mean absolute deviation (MAD) around 11 kJ/mol, Spearman's rank correlation coefficient of 0.88) across all 226 SUB-BIG structures retained after removing calculation failures and applying the filtering criteria. Across these 226 structures, the estimated experimental lattice energies ranged from -60 to -269 kJ/mol, with a standard deviation around 29 kJ/mol. The performance of the recommended protocol on pharmaceutically relevant crystals could be somewhat reduced, with an RMSE around 20 kJ/mol (MAD around 13 kJ/mol, Spearman's rank correlation coefficient of 0.76) obtained on 62 structures retained following filtering according to pharmaceutical relevance method B, for which the distribution of experimental values was similar. For a diverse set of 17 SUB-BIG entries, deemed pharmaceutically relevant according to method B, this recommended force-field protocol was compared to dispersion corrected density functional theory (DFT) calculations (PBE + TS). These calculations suggest that the recommended force-field protocol (RMSE around 15 kJ/mol) outperforms PBE + TS (RMSE around 37 kJ/mol), although it may not outperform more sophisticated DFT protocols and future studies should investigate this. Finally, further work is required to compare our recommended protocol to other lattice energy calculation protocols reported in the literature, as comparisons based upon previously reported smaller datasets indicated this protocol was outperformed by a number of other methods. The SUB-BIG dataset provides a basis for these future studies and could support protocol refinement.

Off-the-shelf DFT-DISPersion methods: Are they now "on-trend" for organic molecular crystals?

Organic molecular crystals contain long-range dispersion interactions that can be challenging for solid-state methods such as density functional theory (DFT) to capture, and in some industrial sectors are overlooked in favor of classical methods to calculate atomistic properties. Hence, this publication addresses the critical question of whether dispersion corrected DFT calculations for organic crystals can reproduce the structural and energetic trends seen from experiment, i.e., whether the calculations can now be said to be truly "on-trend." In this work, we assess the performance of three of the latest dispersion-corrected DFT methods, in calculating the long-range, dispersion energy: the pairwise methods of D3(0) and D3(BJ) and the many-body dispersion method, MBD@rsSCS. We calculate the energetics and optimized structures of two homologous series of organic molecular crystals, namely, carboxylic acids and amino acids. We also use a classical force field method (using COMPASS II) and compare all results to experimental data where possible. The mean absolute error in lattice energies is 9.59 and 343.85 kJ/mol (COMPASS II), 10.17 and 16.23 kJ/mol (MBD@rsSCS), 10.57 and 18.76 kJ/mol [D3(0)], and 8.52 and 14.66 kJ/mol [D3(BJ)] for the carboxylic acids and amino acids, respectively. MBD@rsSCS produces structural and energetic trends that most closely match experimental trends, performing the most consistently across the two series and competing favorably with COMPASS II.

A perspective on using experiment and theory to identify design principles in dye-sensitized solar cells.

Dye-sensitized solar cells (DSCs) have been the subject of wide-ranging studies for many years because of their potential for large-scale manufacturing using roll-to-roll processing allied to their use of earth abundant raw materials. Two main challenges exist for DSC devices to achieve this goal; uplifting device efficiency from the 12 to 14% currently achieved for laboratory-scale 'hero' cells and replacement of the widely-used liquid electrolytes which can limit device lifetimes. To increase device efficiency requires optimized dye injection and regeneration, most likely from multiple dyes while replacement of liquid electrolytes requires solid charge transporters (most likely hole transport materials - HTMs). While theoretical and experimental work have both been widely applied to different aspects of DSC research, these approaches are most effective when working in tandem. In this context, this perspective paper considers the key parameters which influence electron transfer processes in DSC devices using one or more dye molecules and how modelling and experimental approaches can work together to optimize electron injection and dye regeneration.

Ab initio transition state searching in complex systems: fatty acid decarboxylation in minerals.

Because of the importance of mineral catalyzed decarboxylation reactions in both crude oil formation and, increasingly, biofuel production, we present a model study into the decarboxylation of the shortest fatty acid, propionic acid C(2)H(5)COOH, into an alkane and CO(2) catalyzed by a pyrophillite-like, phyllosilicate clay. To identify the decarboxylation pathway, we searched for a transition state between the reactant, comprised of the clay plus interlayer fatty acid, and the product, comprised of the clay plus interlayer alkane and carbon dioxide. Using linear and quadratic synchronous transit mechanisms we searched for a transition state followed by vibrational analysis to verify the intermediate found as a transition state. We employed a periodic cell, planewave, ab initio density functional theory computation to examine total energy differences, Mulliken charges, vibrational frequencies, and the frontier orbitals of the reactants, intermediates, and products. The results show that interpretation of vibrational data, Mulliken charges and Fermi-level orbital occupancies is necessary for the classification of a transition state in this type of mixed bulk surface plus interlayer species, clay-organic system.

Role of clay minerals in oil-forming reactions.

Mineral-catalyzed decarboxylation reactions are important in both crude oil formation and, increasingly, biofuel production. In this study we examined decarboxylation reactions of a model fatty acid, propionic acid, C(2)H(5)COOH, to an alkane, C(2)H(6), in a model of pyrophillite with an isomorphic substitution of aluminum in the tetrahedral layer. We model a postulated reaction mechanism (Almon, W. R.; Johns, W. D. 7th International Meeting on Organic Geochemistry 1975, Vol. 7) to ascertain the role of Al substitution and a counterion in decarboxylation reactions. We employ a periodic cell, planewave, ab initio DFT computation to examine the total energies and the frontier orbitals of different model sets, including the effects of charge on the reaction, the effect of Al substitution, and the role of Na counterions. The results show that an uncharged system with a sodium counterion is most feasible for catalyzing the decarboxylation reaction in an Al-substituted pyrophillite and, also, that analysis of the orbitals is a better indicator of a reaction than charge alone.