Prediction of densities for solid energetic molecules with molecular surface electrostatic potentials

The densities of high energetic molecules in the solid state were calculated with a simplified scheme based on molecular surface electrostatic potentials (MSEP). The MSEP scheme for density estimation, originally developed by Politzer et al., was further modified to calculate electrostatic potential on a simpler van der Waals surface. Forty‐one energetic molecules containing at least one nitro group were selected from among a variety of molecular types and density values, and were used to test the suitability of the MSEP scheme for predicting the densities of solid energetic molecules. For comparison purposes, we utilized the group additivity method (GAM) incorporating the parameter sets developed by Stine (Stine‐81) and by Ammon (Ammon‐98 and ‐00). The absolute average error in densities from our MSEP scheme was 0.039 g/cc. The results based on our MSEP scheme were slightly better than the GAM results. In addition, the errors in densities generated by the MSEP scheme were almost the same for various molecule types, while those predicted by GAM were somewhat dependent upon the molecule types. © 2008 Wiley Periodicals, Inc. J Comput Chem, 2008

[1]  M. Parrinello,et al.  Exploring polymorphism: the case of benzene. , 2005, Angewandte Chemie.

[2]  William Jones,et al.  Beyond the isotropic atom model in crystal structure prediction of rigid molecules: atomic multipoles versus point charges , 2005 .

[3]  Chang Kon Kim,et al.  Prediction of physicochemical properties of organic molecules using van der Waals surface electrostatic potentials , 2004, J. Comput. Chem..

[4]  C. Pickard,et al.  Systematic prediction of crystal structures: An application to sp 3 -hybridized carbon polymorphs , 2004 .

[5]  Scott M. Woodley,et al.  Prediction of crystal structures using evolutionary algorithms and related techniques , 2004 .

[6]  J. Dunitz.  Are crystal structures predictable? , 2003, Chemical communications.

[7]  Herman L. Ammon,et al.  New Atom/Functional Group Volume Additivity Data Bases for the Calculation of the Crystal Densities of C-, H-, N-, O-, F-, S-, P-, Cl-, and Br-Containing Compounds , 2001 .

[8]  B. Rice,et al.  Theoretical Studies of Solid Nitromethane , 2000 .

[9]  H. Scheraga,et al.  Crystal Structure Prediction by Global Optimization as a Tool for Evaluating Potentials: Role of the Dipole Moment Correction Term in Successful Predictions† , 2000 .

[10]  H. Ammon,et al.  A New Atom/Functional Group Volume Additivity Data Base for the Calculation of the Crystal Densities of C, H, N, O and F‐Containing Compounds , 1998 .

[11]  Betsy M. Rice,et al.  ISOTHERMAL-ISOBARIC MOLECULAR DYNAMICS SIMULATIONS OF 1,3,5,7-TETRANITRO-1,3,5,7-TETRAAZACYCLOOCTANE (HMX) CRYSTALS , 1998 .

[12]  Donald L. Thompson,et al.  Intermolecular Potential for the Hexahydro-1,3,5-trinitro-1,3,5-s-triazine Crystal (RDX): A Crystal Packing, Monte Carlo, and Molecular Dynamics Study , 1997 .

[13]  J. Murray,et al.  Relationships of molecular surface electrostatic potentials to some macroscopic properties , 1996 .

[14]  A. Gavezzotti,et al.  Are Crystal Structures Predictable , 1994 .

[15]  Per-Anders Persson,et al.  Rock Blasting and Explosives Engineering , 1993 .

[16]  L. E. Wangen,et al.  A multiblock partial least squares algorithm for investigating complex chemical systems , 1989 .

[17]  A. Höskuldsson.  PLS regression methods , 1988 .

[18]  Michael J. Frisch,et al.  Self‐consistent molecular orbital methods 25. Supplementary functions for Gaussian basis sets , 1984 .

[19]  J. Stine.  Prediction of crystal densities of organic explosives by group additivity , 1981 .

[20]  J. Agrawal.  Recent trends in high-energy materials , 1998 .

[21]  A. Gavezzotti,et al.  Self-organization of small organic molecules in liquids, solutions and crystals: static and dynamic calculations , 1998 .

[22]  Peter Politzer,et al.  Quantitative treatments of solute/solvent interactions , 1994 .