Screened dipolar interactions in some molecular crystals

Energy of charged states in the RDX crystal: Trapping of charge-transfer pairs as a possible mechanism for initiating detonation

Calculations for the crystalline energetic material RDX (1,3,5-trinitro-1,3,5-triazacyclohexane) yield the effective polarizability (17.2 angstroms3), local electric field tensor, effective dipole moment (9.40 D), and dipole-dipole energy (-27.2 kJ/mol). Fourier-transform techniques give the polarization energy P for a single charge in the perfect crystal as -1.14 eV; the charge-dipole energy W(D) is zero if the crystal carries no bulk dipole moment. Polarization energies for charge-transfer (CT) pairs combine with the Coulomb energy E(C) to give the screened Coulomb energy E(scr); screening is nearly isotropic with E(scr) approximately = E(C)2.6. For CT pairs W(D) reduces to a term deltaW(D) arising from the interaction of the charge on each ion with the change in dipole moment on the other ion relative to the neutral molecule. The dipole moments are calculated as 7.40 D for the neutral molecule and 6.84 D and 7.44 D for the anion and cation, giving the lowest two CT pairs at -1.34 eV and -0.94 eV. The changes in P and W(D) near a molecular vacancy yield traps with depths that reach 400 meV for single charges and 185 meV for the nearest-neighbor CT pair. Divacancies yield traps with depths nearly equal to the sum of those produced by the separate vacancies. These results are consistent with a mechanism in which detonation of RDX is initiated by mechanical generation of CT pairs that localize at vacancies, recombine, and release energy sufficient to break bonds; crystals of molecules with lower dipole moments should be less sensitive.

Energy of charged states in the acetanilide crystal: trapping of charge-transfer states at vacancies as a possible mechanism for optical damage

Calculations for the acetanilide crystal yield the effective polarizability (16.6 A(3)), local electric field tensor, effective dipole moment (5.41 D), and dipole-dipole energy (-12.8 kJ/mol). Fourier-transform techniques are used to calculate the polarization energy P for a single charge in the perfect crystal (-1.16 eV); the charge-dipole energy W(D) is zero if the crystal carries no bulk dipole moment. Polarization energies for charge-transfer (CT) pairs combine with the Coulomb energy E(C) to give the screened Coulomb energy E(scr); screening is nearly isotropic, with E(scr) approximately E(C)/2.7. For CT pairs W(D) reduces to a term deltaW(D) arising from the interaction of the charge on each ion with the change in dipole moment on the other ion relative to the neutral molecule. The dipole moments calculated by density-functional theory methods with the B3LYP functional at the 6-311++G(**) level are 3.62 D for the neutral molecule, changing to 7.13 D and 4.38 D for the anion and cation, relative to the center of mass. Because of the large change in the anion, deltaW(D) reaches -0.9 eV and modifies the sequence of CT energies markedly from that of E(scr), giving the lowest two CT pairs at -1.98 eV and -1.41 eV. The changes in P and W(D) near a vacancy are calculated; W(D) changes for the individual charges because the vacancy removes a dipole moment and modifies the crystal dielectric response, but deltaW(D) and E(C) do not change. A vacancy yields a positive change DeltaP that scatters a charge or CT pair, but the change DeltaW(D) can be negative and large enough to outweigh DeltaP, yielding traps with depths that can exceed 150 meV for single charges and for CT pairs. Divacancies yield traps with depths nearly equal to the sum of those produced by the separate vacancies and so they can exceed 300 meV. These results are consistent with a mechanism of optical damage in which vacancies trap optically generated CT pairs that recombine and release energy; this can disrupt the lattice around the vacancy, thereby favoring trapping and recombination of CT pairs generated by subsequent photon absorption, leading to further lattice disruption. Revisions to previous calculations on trapping of CT pairs in anthracene are reported.