Search

Interpreting the initial decomposition mechanism is important for evaluating the thermal stability of explosives. In this study, we theoretically investigated the initial thermal decomposition reactions for two typical energetic materials, FOX-7 and RDX, in both the gas phase and crystal phase. Single molecular decomposition pathways in the gas phase are calculated using the density functional theory (DFT) method, and the crystal phase reactions are simulated through the MM/DFT-based ONIOM method. The calculation results indicate that the crystal environment has a significant influence on the initial thermal decomposition mechanism of FOX-7 and RDX. The cage effect induced by the crystal environment greatly confines molecular mobility and diffusion, rendering the generated small molecules to react with the remaining fragment and yield new decomposition channels compared with the gas phase condition. The crystal packing structures and intermolecular interactions (hydrogen bonds/π-π stacking) significantly increase the reaction barriers of FOX-7 and RDX, leading to the crystal phase reactions being more difficult to occur than in the gas phase. Since the practical application of explosives is mostly in the crystal state, it is important to consider the environmental effects on the initial decomposition reactions. The same insight can also be relevant for other energetic materials.

The deactivation process of salicylidene methylamine in the gas phase has been explored using static calculations (CASSCF, CASPT2, and CC2) and on-the-fly surface hopping dynamics simulations (CASSCF). Five minimum energy conical intersections (MECIs) were located upon the geometry optimization calculations. One corresponds to the excited state intramolecular proton transfer (ESIPT) process, and the remaining four arise from CN bond rotational motion. Our calculation results found that the molecule prefers to decay to the ground state through the four rotational motion related MECIs rather than the ESIPT related one. This mechanistic scenario is verified by the energy profiles connecting the Franck-Condon point and the MECIs at CASSCF, CASPT2, and CC2 levels. Our proposed new decay mechanism can explain the previous experimental findings of femtosecond pump-probe photoionization spectroscopy and can provide additional guidance to the rational design of photochemically switchable molecules.

The effect of a hydrogen bond on the photochemical synthesis of silver nanoparticles has been investigated via experimental and theoretical methods. In a benzophenone system, the photochemical synthesis process includes two steps, which are that hydrogen abstraction reaction and the following reduction reaction. We found that for the first step, an intermolecular hydrogen bond enhances the proton transfer. The efficiency of hydrogen abstraction increases with the hydrogen bond strength. For the second step, the hydrogen-bonded ketyl radical complex shows higher reducibility than the ketyl radical. The inductively coupled plasma-optical emission spectroscopy (ICP-OES) measurement exhibits a 2.49 times higher yield of silver nanoparticles in the hydrogen bond ketyl radical complex system than that for the ketyl radical system. Theoretical calculations show that the hydrogen bond accelerates electron transfer from the ketyl radical to the silver ion by raising the SOMO energy of the ketyl radical; thus, the SOMO-LUMO interaction is more favorable.

DFT and TDDFT methods have been performed to investigate the photoisomerization effect for dithiazolylarylene on solution. The weak S···N interaction and CH···N hydrogen bond restrain the rotation of the side-chain thiazolyl ring in open-isomer 1a, the higher stability of which prefers to show a high quantum yield of photoisomerization. The calculated UV-Vis spectrum at around 320 nm for open-isomer 1a is bathochromically shifted to 647 nm for closed-isomer 1b, in excellent agreement with the experimental photochromic phenomenon. The electron transition in ECD (electron circular dichroism) spectra for closed-isomer 1b with two chiral carbon atoms is dominated by ICT (intramolecular charge transition) and LE (local excitation) corresponding to one positive (440 nm) and one negative Cotton effect (650 nm), respectively, where the two chiral carbon atoms play a slight role in these transitions. The PES in the S(1) and S(0) states, respectively, indicates that the cyclization reaction from open-isomer 1a to closed-isomer 1b is allowed in the photoexcited state with high-conversion quantum efficiency, while it is forbidden in the thermodynamic process. In addition, the second-order nonlinear optical response for closed-isomer 1b is nearly six times larger than that for open-isomer 1a. It is also confirmed that the photoirradiation evokes the photoisomerization character to show dramatic difference in the second-order NLO response, which can be applied to designing photochromic materials and reversible NLO switches.

The rotational reorientation dynamics of oxazine 750 (OX750) in the first (with pump pulse at 660 nm) and a higher excited state (with pump pulse at 400 nm) in different polar solvents have been investigated using femtosecond time-resolved stimulated emission pumping fluorescence depletion (FS TR SEP FD) spectroscopy. In both excited states, three different anisotropy decay laws have been observed for OX750 in different solvents. Only in acetone and formamide could the anisotropy decays of OX750 be described by single-exponential functions, whereas the anisotropy decays have been found to exhibit biexponential behavior in other solvents. The slower anisotropy decay observed in all of the solvents has been assigned to the overall rotational relaxation of OX750 molecules, and a quantitative analysis of this time constant has been performed using the Stokes-Einstein-Debye hydrodynamic theory and the extended charge distribution model developed by Alavi and Waldeck. In both methanol and ethanol, a faster anisotropy decay on the order of picoseconds and a slower anisotropy decay on the hundreds of picoseconds time scale are observed. The most likely explanation for the faster anisotropy involves the rotation of the transition dipole moment in the excited state of OX750 resulting from the electron transfer (ET) reaction taking place from the alcoholic solvents to the OX750 chromophore. As a possible explanation, the wobbling-in-the-cone model has been used to analyze the biexponential anisotropy decays of OX750 in dimethylformamide (DMF) and dimethyl sulfoxide (DMSO). The observed faster anisotropy decays on the hundreds of femtoseconds time scale in DMF and DMSO are ascribed to the wobbling-in-the-cone motion of the ethyl group of OX750, which is sensitive to the strength of the hydrogen bond formed between the solvent and the protonation site of OX750.