RESUMEN
CONTEXT: An accurate description of the molecular charge density is crucial for investigating intra- and inter-molecular properties. Among the different ways of describing and analyzing it, the widely used distributed multipole analysis (DMA) is an accurate method for decomposing the molecular charge density into atom-centered electric multipoles (monopole, dipole, quadrupole, and so on) that have a direct chemical interpretation. In this work, DMA was employed to decompose the molecular charge density of six chemically distinct molecules, namely, (2R)-2-amino-3-[(S)-prop-2-enylsulfinyl] propanoic acid (AAP), 4-amine-2-nitro-1,3,5 triazole (ANTA), (RS)-Propan-2-yl methylphosphonofluoridate (SARIN), chloromethane (CLMET), and 2-aminoacetic acid (GLY) into monopole, dipole, and quadrupole values. A hypothetical variation of ANTA built by exchanging all the nitrogen atoms with phosphorus that we named 4-phosphine-2-phosphite-1,3,5-phosphorine (ANTAP) was also studied. These molecules have different chemical structures bearing distinct carbon skeletons, electronegative atoms, and electron-withdrawing/donating groups. We found that although DFT multipole values can depend considerably on the exchange-correlation functional for specific atomic sites, the associated root-mean-square errors (RMSEs) compared to benchmark MP4 mainly were about [Formula: see text] The most significant variations were for monopoles and dipoles of sites highly polarized by adjacent atoms, and to a lesser degree, for the quadrupoles. The double hybrid B2PLYP and the hybrid meta M06-2X functionals, as expected in the framework of Jacob's ladder, overall give the most accurate results among the DFT methods. The MP2 DMA multipole values have an RMSE in relation to the MP4 benchmark mainly in the range [Formula: see text], thus representing a lower computational cost to obtain results with similar good accuracy without the ambiguity of choosing a DFT functional. The deviations of the HF multipoles from the benchmark in most cases were less than 20%, in agreement with the well-known fact that non-correlated charge densities have a slight dependence on the electronic correlation. We also confirmed that DMA values have a small dependence on the size of the basis set: deviations did not exceed 5% in most cases. However, the dependence of the DMA values on the size of the basis set increases with the rank of the electric multipole. To compute accurate values of DMA multipoles of an atom bonded to very electronegative atoms, especially dipoles and quadrupoles, a large basis set including diffuse functions is necessary. Despite that, for a given polarized basis set, the choice of the basis set to compute accurate DMA multipole values is not critical. METHODS: The molecular charge densities were computed using the electronic structure methods Hartree-Fock (HF), MP2, MP4, DFT/PBE, DFT/B3LYP, DFT/B3PW91, DFT/M06-2X, and DFT/B2PLYP implemented in the Gaussian 09 package. MP4 was the benchmark method. The DMA multipoles were obtained with the GDMA program of Stone. The 6-311G + + (d,p) basis set was used for the production calculations, and the augmented correlation-consistent Dunning's hierarchy of basis sets was employed to evaluate the dependence of the DMA multipoles on the basis set size.
RESUMEN
The quantification of bond strengths is a useful and general concept in chemistry. In this work, a Coulombic force model based on atomic electric charges computed using the accurate distributed multipole analysis (DMA) partition of the molecular charge density was employed to quantify the weakest N-NO2 and C-NO2 bond strengths of 19 nitrobenzene, 11 nitroazole, and 10 nitramine molecules. These bonds are known as trigger linkages because they are usually related to the initiation of an explosive. The three families of explosives combine different types of molecular properties and structures ranging from essentially aromatic molecules (nitrobenzenes) to others with moderate aromaticity (nitroazoles) and non-aromatic molecules with cyclic and acyclic skeletons (nitramines). We used the results to investigate the impact sensitivity of the corresponding explosives employing the trigger linkage concept. For this purpose, the computed Coulombic bond strength of the trigger linkages was used to build four sensitivity models that lead to an overall good agreement between the predicted values and available experimental sensitivity values even when the model included the three chemical families simultaneously. We discussed the role of the trigger linkages for determining the sensitivity of the explosives and rationalized eventual discrepancies in the models by examining alternative decomposition mechanisms and features of the molecular structures.