Harmonic force field for nitro compounds

EP Bellido and JM Seminario, JOURNAL OF MOLECULAR MODELING, 18, 2805-2811 (2012).

DOI: 10.1007/s00894-011-1273-3

Molecular simulations leading to sensors for the detection of explosive compounds require force field parameters that can reproduce the mechanical and vibrational properties of energetic materials. We developed precise harmonic force fields for alanine polypeptides and glycine oligopeptides using the FUERZA procedure that uses the Hessian tensor (obtained from ab initio calculations) to calculate precise parameters. In this work, we used the same procedure to calculate generalized force field parameters of several nitro compounds. We found a linear relationship between force constant and bond distance. The average angle in the nitro compounds was 116A degrees, excluding the 90A degrees angle of the carbon atoms in the octanitrocubane. The calculated parameters permitted the accurate molecular modeling of nitro compounds containing many functional groups. Results were acceptable when compared with others obtained using methods that are specific for one type of molecule, and much better than others obtained using methods that are too general (these ignore the chemical effects of surrounding atoms on the bonding and therefore the bond strength, which affects the mechanical and vibrational properties of the whole molecule).

Return to Publications page