Anharmonic Vibrational Properties of Explosives from Temperature-Dependent Raman

peri-Fused polyaromatic molecular contacts for perovskite solar cells

Molecule-based selective contacts have become a crucial component to ensure high-efficiency inverted perovskite solar cells1-5. These molecules always consist of a conjugated core with heteroatom substitution to render the desirable carrier-transport capability6-9. So far, the design of successful conjugation cores has been limited to two N-substituted π-conjugated structures, carbazole and triphenylamine, with molecular optimization evolving around their derivatives2,5,10-12. However, further improvement of the device longevity has been hampered by the concomitant limitations of the molecular stability induced by such heteroatom-substituted structures13,14. A more robust molecular contact without sacrificing the electronic properties is in urgent demand, but remains a challenge. Here we report a peri-fused polyaromatic core structure without heteroatom substitution that yields superior carrier transport and selectivity over conventional heteroatom-substituted core structures. This core structure produced a relatively chemically inert and structurally rigid molecular contact, which considerably improved the performance of perovskite solar cells in terms of both efficiency and durability. The champion device showed an efficiency up to 26.1% with greatly improved longevity under different accelerated-ageing tests.

Selective Phonon Stimulation Mechanism to Tune Thermal Transport

In this paper, we determine the degree to which changes can be induced in the equilibrium thermal diffusivity and conductivity of a material via a selective nonequilibrium infrared stimulation mechanism for phonons. Using the molecular crystal RDX, we use detailed momentum-dependent coupling information across the entire Brillouin zone and the phonon gas model to show that stimulating selected modes in the spectrum of a target material can induce substantial changes in the overall thermal transport properties. Specifically in the case of RDX, stimulating modes at ∼22.74 cm^-1 over a linewidth of 1 cm^-1 can lead to enhanced scattering rates that reduce the overall thermal diffusivity and conductivity by 15.58 and 12.46%, respectively, from their equilibrium values. Due to the rich spectral content in the materials, however, stimulating modes near ∼1140.67 cm^-1 over a similar bandwidth can produce an increase in the thermal diffusivity and conductivity by 55.73 and 144.07%, respectively. The large changes suggest a mechanism to evoke substantially modulated thermal transport properties through light-matter interaction.

Mode-Selective Vibrational Energy Transfer Dynamics in 1,3,5-Trinitroperhydro-1,3,5-triazine (RDX) Thin Films

The coupling of inter- and intramolecular vibrations plays a critical role in initiating chemistry during the shock-to-detonation transition in energetic materials. Herein, we report on the subpicosecond to subnanosecond vibrational energy transfer (VET) dynamics of the solid energetic material 1,3,5-trinitroperhydro-1,3,5-triazine (RDX) by using broadband, ultrafast infrared transient absorption spectroscopy. Experiments reveal VET occurring on three distinct time scales: subpicosecond, 5 ps, and 200 ps. The ultrafast appearance of signal at all probed modes in the mid-infrared suggests strong anharmonic coupling of all vibrations in the solid, whereas the long-lived evolution demonstrates that VET is incomplete, and thus thermal equilibrium is not attained, even on the 100 ps time scale. Density functional theory and classical molecular dynamics simulations provide valuable insights into the experimental observations, revealing compression-insensitive time scales for the initial VET dynamics of high-frequency vibrations and drastically extended relaxation times for low-frequency phonon modes under lattice compression. Mode selectivity of the longest dynamics suggests coupling of the N-N and axial NO₂ stretching modes with the long-lived, excited phonon bath.

3-Phonon Scattering Pathways for Vibrational Energy Transfer in Crystalline RDX

A long-held belief is that shock energy induces initiation of an energetic material through an energy up-pumping mechanism involving phonon scattering through doorway modes. In this paper, a Fermi's golden rule-based 3-phonon theoretical analysis of energy up-pumping in RDX is presented that considers possible doorway pathways through which energy transfer occurs. On average, modes with frequencies up to 102 cm^-1 scatter quickly and transfer over 99% of the vibrational energy to other low-frequency modes up to 102 cm^-1 within 0.16 ps. These low-frequency modes scatter less than 0.5% of the vibrational energy directly to modes with significant nitrogen-nitrogen (NN) activity. The midfrequency modes from 102 to 1331 cm^-1 further up-pump the energy to these modes within 5.6 ps. The highest-frequency modes scatter and redistribute a small fraction of the vibrational energy to all other modes, which last over 2000 ps. The midfrequency modes between 457 and 462 cm^-1 and between 831 and 1331 cm^-1 are the most critical for vibrational heating of the NN modes and phenomena, leading to initiation in energetics. In contrast, modes stimulated by the shock with frequencies up to 102 cm^-1 dominate vibrational cooling of the NN modes.

Predicting the impact sensitivities of energetic materials through zone-center phonon up-pumping

The development of new energetic materials (EMs) is accompanied by significant hazards, prompting interest in their computational design. Before reliable in silico design strategies can be realized, however, approaches to understand and predict EM response to mechanical impact must be developed. We present here a fully ab initio model based on phonon up-pumping that successfully ranks the relative impact sensitivity of a series of organic EMs. The methodology depends only on the crystallographic unit cell and Brillouin zone center vibrational frequencies. We, therefore, expect this approach to become an integral tool in the large-scale screening of potential EMs.

Sub-picosecond to Sub-nanosecond Vibrational Energy Transfer Dynamics in Pentaerythritol Tetranitrate

The time scale associated with shock-induced detonation is a key property of energetic materials that remains poorly understood. Herein, we test aspects of one potential mechanism, the phonon up-pumping mechanism, where shock compression excites lattice phonon modes, transferring energy to intramolecular vibrations leading to chemical bond cleavage and reaction. Using ultrafast infrared pump-probe spectroscopy on pentaerythritol tetranitrate (PETN), we reveal sub-picosecond vibrational energy transfer (VET) from the photoexcited band at 1660 cm^-1 into every other infrared-active mode in the probed frequency range 800-1800 cm^-1. Energy transfer processes remain incomplete at 150 ps. Computational predictions from density functional theory are used in tandem to elucidate VET pathways in PETN.

Intermolecular Energy Transfer Dynamics at a Hot-Spot Interface in RDX Crystals

The phonon mediated vibrational up-pumping mechanisms assume an intact lattice and climbing of a vibrational ladder using strongly correlated multiphonon dynamics under equilibrium or near-equilibrium conditions. Important dynamic processes far from-equilibrium in regions of large temperature gradient after the onset of decomposition reactions in energetic solids are relatively unknown. In this work, we present a classical molecular dynamics (MD) simulation-based study of such processes using a nonreactive and a reactive potential to study a fully reacted and unreacted zone in RDX (1,3,5-trinitro-1,3,5-triazocyclohexane) crystal under nonequilibrium conditions. The energy transfer rate is evaluated as a function of temperature difference between the reacted and unreacted regions, and for different widths and cross-sectional area of unreacted RDX layers. Vibrational up-pumping processes probed using velocity autocorrelation functions indicate that the mechanisms at high-temperature interfaces are quite different from the standard phonon-based models proposed in current literature. In particular, the up-pumping of high-frequency vibrations are seen in the presence of small molecule collisions at the hot-spot interface with strong contributions from bending modes. It also explains some major difference in the order of decomposition of C-N and N-N bonds as seen in recent literature on initiation chemistry.

Temperature dependence of thermal expansion tensors of energetic materials

Unit-cell values as well as thermal expansion tensors for 13 energetic materials are calculated from variable-temperature X-ray diffraction data. The thermal expansion tensors and their temperature dependence are reported numerically, algebraically and graphically.

Low temperature Raman study of a liquid crystalline system 4-decyloxy benzoic acid (4DBA)

The Raman spectra of a liquid crystalline system, 4-Decyloxy benzoic acid (4DBA) have been recorded at different temperatures within the interval 300-78 K in order to identify the structural changes in crystalline state of a nematogen and to understand the molecular alignment therein. The earlier predicted dimer structure of 4DBA was optimized with DFT method and the theoretical Raman spectra of dimer as well as monomer have been calculated for comparison with the experimental spectra. The mode specific quartic coupling coefficient; A(i,ω) and phonon frequency; ωi have been calculated using temperature dependent anharmonic perturbation theory. The precise band shape analysis of Raman bands at ~807, ~881, ~1255, ~1282, ~1436, ~1576, ~1604, ~2881 and ~3081cm(-1) gives signature of temperature induced slow crystal modification. The structural changes leading to crystal modification have been discussed.

Deep ultraviolet resonance Raman excitation enables explosives detection

We measured the 229 nm absolute ultraviolet (UV) Raman cross-sections of the explosives trinitrotoluene (TNT), pentaerythritol tetranitrate (PETN), cyclotrimethylene-trinitramine (RDX), the chemically related nitroamine explosive HMX, and ammonium nitrate in solution. The 229 nm Raman cross-sections are 1000-fold greater than those excited in the near-infrared and visible spectral regions. Deep UV resonance Raman spectroscopy enables detection of explosives at parts-per-billion (ppb) concentrations and may prove useful for stand-off spectroscopic detection of explosives.

Raman and infrared fingerprint spectroscopy of peroxide-based explosives

A comparative study of the vibrational spectroscopy of peroxide-based explosives is presented. Triacetone triperoxide (TATP) and hexamethyl-enetriperoxide-diamine (HMTD), now commonly used by terrorists, are examined as well as other peroxide-ring structures: DADP (diacetone diperoxide); TPTP [3,3,6,6,9,9-Hexaethyl-1,2,4,5,7,8-hexaoxo-nonane (tripentanone triperoxide)]; DCypDp {6,7,13,14-Tetraoxadispiro [4.2.4.2]tetradecane (dicyclopentanone diperoxide)}; TCypDp {6,7,15,16,22,23-Hexaoxatrispiro[4.2.4.2.4.2] henicosane (tricyclopentanone triperoxide)}; DCyhDp {7,8,15,16-tetraoxadispiro [5.2.5.2] hexadecane (dicyclohexanone diperoxide)}; and TCyhTp {7,8,14,15,21,22-hexaoxatrispiro [5.2.5.2.5.2] tetracosane (tricyclohexanone triperoxide)}. Both Raman and infrared (IR) spectra were measured and compared to theoretical calculations. The calculated spectra were obtained by calculation of the harmonic frequencies of the studied compounds, at the density functional theory (DFT) B3LYP/cc-pVDZ level of theory, and by the use of scaling factors. It is found that the vibrational features related to the peroxide bonds are strongly mixed. As a result, the spectrum is congested and highly sensitive to minor changes in the molecule.

Comparison of the numerical stability of methods for anharmonic calculations of vibrational molecular energies

On model examples, we compare the performance of the vibrational self-consistent field, variational, and four perturbational schemes used for computations of vibrational energies of semi-rigid molecules, with emphasis on the numerical stability. Although the accuracy of the energies is primarily dependent on the quality of the potential energy surface, approximate approaches to the anharmonic vibrational problem often do not converge to the same results due to the approximations involved. For furan, the sensitivity to variations of the anharmonic potential was systematically investigated by adding random noise to the cubic and quartic constants. The self-consistent field methods proved to be the most resistant to the potential variations. The second order perturbational techniques are sensitive to random degeneracies and provided the least stable results. However, their stability could be significantly improved by a simple generalization of the perturbational formula. The variational configuration interaction is practically limited by the size of the matrix that can be diagonalized for larger molecules; however, relatively fewer states need to be involved than for smaller ones, in favor of the computing.

Theoretical Studies of Energy Transfer Rates of Secondary Explosives

Understanding the mechanism of shock-induced chemical reaction in secondary explosives is necessary to pursue the development and the safe use of new explosives having high performance and low sensitivity. In an effort to understand the mechanism, the energy transfer rates of such secondary explosives as PETN(I), PETN(II), delta-HMX, alpha-HMX, beta-HMX, RDX, ANTA, DMN, and NM have been evaluated based on the formula derived by Fried and Ruggiero [Fried, L. E.; Ruggiero, A. J. J. Phys. Chem. 1994, 98, 9786]. The energy transfer rates were determined in terms of the density of vibrational states and the unharmonic vibron-phonon coupling term, which were calculated by using a flexible potential containing both intra- and intermolecular terms. For the secondary explosives, a good correlation was found between the energy transfer rates and the impact sensitivity. The energy transfer rates are several times faster for the explosives with higher sensitivity such as PETN, HMX, and RDX than those with lower sensitivity such as ANTA, DMN, and NM. The calculations presented suggest the energy transfer rate in secondary explosive crystals is a significant factor in their sensitivity and introduction of double bond, or hydrogen bonds, or caged structure into secondary explosives is expected to decrease the sensitivity.