Ab initio MRD-CI calculations for breaking a chemical bond in a molecule in a crystal or other solid environment. III. Me2N?NO2 decomposition of dimethylnitramine in a large crystalline environment

Recently we extended our strategy for MRD-CI (multireference double excitation-configuration interaction) calculations, based on localized/local orbitals and an “effective” CI Hamiltonian, for molecular decompositions of large molecules to breaking a chemical bond in a molecule in a crystalline or other solid environment. Our technique begins with an explicit quantum chemical SCF calculation for a reference molecule surrounded by a number of other molecules in the multipole environment of more distant neighbors. The resulting canonical molecular orbitals are then localized, and the localized occupied and virtual orbitals in the region of interest are included explicitly in the MRD-CI with the remainder of the occupied localized orbitals being folded into an “effective” CI Hamiltonian. The MRD-CI calculations are then carried out for breaking a bond in the reference molecule. This method is completely general in that the space treated explicitly, as well as the surrounding space, may contain voids, defects, deformations, dislocations, impurities, dopants, edges and surfaces, boundaries, etc. Dimethylnitramine is the smallest prototype of the energetic R₂N—NO₂ nitramines, such as the 6-member ring RDX or the 8-member ring HMX. Decomposition of energetic compounds is initiated in the solid by a breaking of the target bond. Thus, it is crucial to know the difference in energy between breaking a bond in an isolated energetic molecule versus in the molecule in a solid. In the present study, we have carried out MRD-CI calculations for the Me₂N—NO₂ dissociation of dimethylnitramine in a dimethylnitramine crystal. The cases we investigated were one dimethylnitramine molecule (surrounded by 53 and 685 neighboring dimethylnitramine molecules represented by multipoles), three dimethylnitramine molecules, and three dimethylnitramine molecules (surrounded by 683 neighbors). All multipoles were cumulative atomic multipoles up through quadrupoles. The MRD-CI calculations on dimethylnitramine required large numbers of reference configurations from which were allowed all single and double excitations.