Molecular dynamic simulation on comprehensive performance of ferrocene-based burning rate catalysts in a solid propellent

Abstract

To determine the compatibility and interaction between ferrocene (Fc) and other components in AP/HTPB/Fc-based BRCs propellant, molecular dynamics simulation (MD simulation) method was used to study the binding energy, cohesive energy density (CED), mean square displacement (MSD) and mechanical properties of mononuclear ferrocene (ethyl ferrocene, butyl ferrocene and octyl ferrocene), binuclear ferrocene (2,2-(ethyl ferrocene) propane, 2,2, -bis-(butylferrocene) propane and hydroxymethyl bis-ferrocenyl propane), ionic ferrocene (three ionic ferrocenes with ferrocene methyl dimethyl ammonium as cation, perchlorate, nitrate and picrate as anions respectively) and cyclic ferrocene (1- ferrocenyl methyl -2- ferrocenyl -3- methylbenzimidazolium, binuclear ferrocene bridged by bridged biphenylferrocene and dihydropyrazole).The binding energy results show that alkyl groups, polar groups (such as hydroxyl groups) and cyclic substituents can improve the interaction between Fc and AP/HTPB. Moreover, the electrostatic interaction between ionic ferrocene and AP is strong, but the van der Waals interaction between ionic ferrocene and HTPB is weak. The CED data show that the mixed models based on alkyl mononuclear ferrocene, catocene, cyclic ferrocene and ionic ferrocene all have high cohesive energy with good safety performance. The MSD results show that the order of anti-migration ability is cyclic ferrocene > binuclear alkyl ferrocene > mononuclear alkyl ferrocene > binuclear hydroxyl ferrocene > ionic ferrocene. Regarding to the value of mechanical parameters (K, G, E and K/G), the mixed model of AP/HTPB/ dihydropyrazole bridged binuclear ferrocene possess the best mechanical properties among the 12 models. Based on the comprehensive data, we conclude that dihydropyrazole bridged binuclear ferrocene is the most ideal ferrocene-based burning rate catalyst among 12 ferrocenes.