Derivation of a New Force Field for Crystal-Structure Prediction Using Global Optimization:  Nonbonded Potential Parameters for Hydrocarbons and Alcohols

2003 ◽  
Vol 107 (29) ◽  
pp. 7143-7154 ◽  
Author(s):  
Yelena A. Arnautova ◽  
Anna Jagielska ◽  
Jaroslaw Pillardy ◽  
Harold A. Scheraga
CrystEngComm ◽  
2021 ◽  
Author(s):  
Jianjun Hu ◽  
Wenhui Yang ◽  
Rongzhi Dong ◽  
Yuxin Li ◽  
Xiang Li ◽  
...  

Crystal structure prediction is now playing an increasingly important role in the discovery of new materials or crystal engineering.


2005 ◽  
Vol 61 (5) ◽  
pp. 528-535 ◽  
Author(s):  
Bouke P. van Eijck

In the third Cambridge blind test of crystal structure prediction, participants submitted extended lists of up to 100 hypothetical structures. In this paper these lists are analyzed for the two small semi-rigid molecules, hydantoin and azetidine, by performing a new energy minimization using an accurate force field, and grouping these newly minimized structures into clusters of equivalent structures. Many participants found the same low-energy structures, but no list appeared to be complete even for the structures with one independent molecule in the asymmetric unit. This may well be due to the fact that a cutoff at even 100 structures cannot ensure the presence of a structure that has a relatively high ranking in another force field. Moreover, some structures should have possibly been discarded because they correspond to transition states rather than true energy minima. The r.m.s. deviation between energies in corresponding clusters was calculated to compare the reported relative crystal energies for each pair of participants. Some groups of force fields show a reasonably good correspondence, yet the order of magnitude of their discrepancies is comparable to the energy differences between, say, the first ten structures of lowest energy. Therefore, even if we assume that energy is a sufficient criterion, it is not surprising that crystal structure predictions are still inconsistent and unreliable.


1998 ◽  
Vol 102 (17) ◽  
pp. 2904-2918 ◽  
Author(s):  
Ryszard J. Wawak ◽  
Jaroslaw Pillardy ◽  
Adam Liwo ◽  
Kenneth D. Gibson ◽  
Harold A. Scheraga

1999 ◽  
Vol 55 (4) ◽  
pp. 543-553 ◽  
Author(s):  
G. Filippini ◽  
A. Gavezzotti ◽  
J. J. Novoa

The crystal structures of two polymorphs of 4,4,5,5-tetramethyl-4,5-dihydro-1H-imidazol-1-oxyl 3-oxide (the 2-hydronitronylnitroxide radical, HNN) are analyzed by packing energy criteria. Other unobserved polymorphic crystal structures are generated using a polymorph predictor package and three different force fields, one of which is without explicit Coulomb-type terms. The relative importance of several structural motifs (hydrogen-bonded dimers, shape-interlocking dimers or extended hydrogen-bonded chains) is discussed. As usual, many crystal structures within a narrow energy range are generated by the polymorph predictor, confirming that ab initio crystal-structure prediction is still problematic. Comparisons of powder patterns generated from the atomic coordinates of the X-ray structure and from computational crystal structures confirm that although the energy ranking depends on the force field used, the X-ray structure of the \alpha polymorph was found to be among the most stable ones produced by the polymorph predictor, even using the chargeless force field.


2014 ◽  
Vol 70 (a1) ◽  
pp. C1541-C1541
Author(s):  
Jacco van de Streek ◽  
Kristoffer Johansson ◽  
Xiaozhou Li

The five Crystal-Structure Prediction (CSP) Blind Tests have shown that molecular-mechanics force fields are not accurate enough for crystal structure prediction[1]. The first--and only--method to successfully predict all four target crystal structures of one of the CSP Blind Tests was dispersion-corrected Density Functional Theory (DFT-D), and this is what we use for our work. However, quantum-mechanical methods (such as DFT-D), are too slow to allow simulations that include the effects of time and temperature, certainly for the size of molecules that are common in pharmaceutical industry. Including the effects of time and temperature therefore still requires molecular dynamics (MD) with less accurate force fields. In order to combine the accuracy of the successful DFT-D method with the speed of a force field to enable molecular dynamics, our group uses Tailor-Made Force Fields (TMFFs) as described by Neumann[2]. In Neumann's TMFF approach, the force field for each chemical compound of interest is parameterised from scratch against reference data from DFT-D calculations; in other words, the TMFF is fitted to mimic the DFT-D energy potential. Parameterising a dedicated force field for each individual compound requires an investment of several weeks, but has the advantage that the resulting force field is more accurate than a transferable force field. Combining crystal-structure prediction with DFT-D followed by molecular dynamics with a tailor-made force field allows us to calculate e.g. the temperature-dependent unit-cell expansion of each predicted polymorph, as well as possible temperature-dependent disorder. This is relevant for example when comparing the calculated X-ray powder diffraction patterns of the predicted crystal structures against experimental data.


Sign in / Sign up

Export Citation Format

Share Document