Harmonic force field for nitro compounds.

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 116°, excluding the 90° 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).