The Gaussian-n, complete basis set, and Weizmann-1 quantum chemical models for heats of formation are applied to a set of molecules with relevance to the combustion or pyrolysis of chemical warfare materials. Most of these models generate standard deviations from experiment that are less than 2 kcal/mol. The structures and vibrational frequencies that are generated in the course of these calculations are in good agreement with experimental data. Detailed comparisons with respect to structural types indicate that the present computational models are likely to generate useful data for complex models of combustion and pyrolysis of chemical warfare materials.

The CBS-QB3 and G4 thermochemical models have been used to generate energetic, structural, and spectroscopic data on a set of molecules with trivalent or pentavalent phosphorus atoms that can serve as simulants of chemical warfare agents. Based on structural data, the conformational stabilities of these molecules are explained in terms of the anomeric interaction within the OPOC and OPSC fragments. For those cases where experimental data are available, comparisons have been made between calculated and previously reported vibrational frequencies. All varieties of bond dissociation energies have been examined except those for C-H and P═O bonds. In trivalent phosphorus molecules, the O-C and S-C bonds have the lowest dissociation energies. In the pentavalent phosphorus set, the S-C bonds, followed by P-S bonds, have the lowest dissociation energies. In the fluorinated simulant molecules, the P-F bond is strongest, and the P-C or O-C bonds are weakest.

A computational study of the kinetics of isomerization and elimination reactions of organophosphorus and organosulfur reactions is presented with a view to characterizing the predictive capabilities of widely applied techniques for processes that pertain to the destruction of chemical warfare agents. A set of 22 reactions has been studied, and the results have been compared to experimentally derived data. The BMK functional and the MG3S basis set have been used to compute minimum energy paths. Corrections have been added from CBS-QB3, CASSCF, and CASMP2 calculations. Thermal rate constants at experimental temperatures have been calculated with canonical variational transition state theory and small-curvature tunneling theory. The quality of these results may depend on recrossing of the variational transition state, the amount of radical or diradical character found in the minimum energy paths, or the accuracy of barrier heights.

Time-dependent density functional theory (TD-DFT) and electron propagator theory (EPT) are used to calculate the electronic transition energies and ionization energies, respectively, of species containing phosphorus or sulfur. The accuracy of TD-DFT and EPT, in conjunction with various basis sets, is assessed with data from gas-phase spectroscopy. TD-DFT is tested using 11 prominent exchange-correlation functionals on a set of 37 vertical and 19 adiabatic transitions. For vertical transitions, TD-CAM-B3LYP calculations performed with the MG3S basis set are lowest in overall error, having a mean absolute deviation from experiment of 0.22 eV, or 0.23 eV over valence transitions and 0.21 eV over Rydberg transitions. Using a larger basis set, aug-pc3, improves accuracy over the valence transitions via hybrid functionals, but improved accuracy over the Rydberg transitions is only obtained via the BMK functional. For adiabatic transitions, all hybrid functionals paired with the MG3S basis set perform well, and B98 is best, with a mean absolute deviation from experiment of 0.09 eV. The testing of EPT used the Outer Valence Green's Function (OVGF) approximation and the Partial Third Order (P3) approximation on 37 vertical first ionization energies. It is found that OVGF outperforms P3 when basis sets of at least triple-ζ quality in the polarization functions are used. The largest basis set used in this study, aug-pc3, obtained the best mean absolute error from both methods -0.08 eV for OVGF and 0.18 eV for P3. The OVGF/6-31+G(2df,p) level of theory is particularly cost-effective, yielding a mean absolute error of 0.11 eV.