Thermal unimolecular decomposition mechanism of 2,4,6-trinitrotoluene: a first-principles DFT study

Simple C–NO₂ homolysis, 4,6-dinitroanthranil (DNAt) production by dehydration, and the nitro-nitrite rearrangement–homolysis for gas-phase TNT decomposition were recently studied by Cohen et al. (J Phys Chem A 111:11074, 2007), based on DFT calculations. Apart from those three pathways, other possible initiation processes were suggested in this study, i.e., CH₃ removal, O elimination, H escape, OH removal, HONO elimination, and nitro oxidizing adjacent backbone carbon atom. The intermediate, 3,5-dinitro-2(or 4)-methyl phenoxy, is more favor to decompose into CO and 3,5-dinitro-2(or 4)-methyl-cyclopentadienyl than to loss NO following nitro-nitrite rearrangement. Below ~1,335 K, TNT condensing to DNAt by dehydration is kinetically the most favor process, and the formations of substituted phenoxy and following cyclopentadienyl include minor contribution. Above ~1,335 K, simple C–NO₂ homolysis kinetically dominates TNT decomposition; while the secondary process changes from DNAt production to CH₃ removal above ~2,112 K; DNAt condensed from TNT by dehydration yields to that by sequential losses of OH and H above ~1,481 K and to nitro-nitrite rearrangement–fragmentation above ~1,778 K; O elimination replaces DNAt production above ~2,491 K, playing the third role in TNT decomposition; H escaping directly from TNT thrives in higher temperature (above ~2,812 K), as the fourth largest process. The kinetic analysis indicates that CH₃ removal, O elimination, and H escape paths are accessible at the suggested TNT detonation time (~100–200 fs), besides C–NO₂ homolysis. HONO elimination and nitro oxidizing adjacent backbone carbon atom paths are negligible at all temperatures. The calculations also demonstrated that some important species observed by Rogers and Dacons et al. are thermodynamically the most favor products at all temperatures, possibly stemmed from the intermediates including 4,6-dinitro-2-nitroso-benzyl alcohol, 3,5-dinitroanline, 2,6-dinitroso-4-nitro-phenylaldehyde, 3,5-dinitro-1-nitrosobenzene, 3,5-dinitroso-1-nitrobenzene, and nitrobenzene. All transition states, intermediates, and products have been indentified, the structures, vibrational frequencies, and energies of them were verified at the uB3LYP/6-311++G(d,p) level. Our calculated energies have mean unsigned errors in barrier heights of 3.4–4.2 kcal/mol (Lynch and Truhlar in J Phys Chem A 105:2936, 2001), and frequencies have the recommended scaling factors for the B3-LYP/6-311+G(d,p) method (Andersson and Uvdal in J Phys Chem A 109:2937, 2005; Merrick et al. in J Phys Chem A 111:11683, 2007). All calculations corroborate highly with the previous experimental and theoretical results, clarifying some pertinent questions.     

<Emphasis Type="Italic">Ab initio</Emphasis> spectroscopic studies of non-rigid molecules: an application to acetic acid

The torsional levels of various isotopologues of acetic acid are determined from an ab initio potential energy surface using a flexible model depending on the OH-torsion and the methyl-torsion coordinates. Previous calculations for CH₃–COOH and CH₃–COOD are review and first theoretical energies of the one-deuterated species CH₂D–COOH are provided. The zero point vibrational energy correction and an exact definition for the methyl-torsional coordinate have been considered. The levels are compared with previous calculations (Senent in Mol Phys 99:1311, 2001) and experimental data (Havey et al. in J Mol Spectrosc 229:151, 2005). Isotopic effects on the torsional barriers and energies are discussed. For CH₂D–COOD, the deuteration splits by 25 cm⁻¹ the zero vibrational energy level.     

Separation of the vibrational,rotational, and translational motions of polyatomic molecules using curvilinear vibrational coordinates

A new method is suggested for separating the vibrational, rotational, and translational motions of polyatomic molecules using curvilinear vibrational coordinates that are linear with respect to the natural vibrational coordinates. It is shown that, in this case, Coriolis interactions between the vibrational and rotational motions are absent. The solutions of the anharmonic vibrational-rotational problems in the curvilinear and linear vibrational coordinates are compared. The absence of Coriolis interactions between the vibrational and rotational motions in the curvilinear vibrational coordinates is proved numerically. The same conclusion is additionally supported by calculations of the anharmonic vibrational energy levels for the H₂O, H₂S, NO₂, SO₂, and ClO₂ molecules in the linear and curvilinear vibrational coordinates using the Hamiltonian designed in the curvilinear vibrational coordinates with and without Coriolis vibrational-rotational interactions. Volgograd Pedagogical University. Translated fromZhurnal Strukturnoi Khimii, Vol. 36, No. 2, pp. 239–254, March–April, 1995. Translated by I. Izvekova     

Close-coupling study of rotational energy transfer of CO (upsilon=2) by collisions with He atoms

Quantum close-coupling scattering calculations of rotational energy transfer in the vibrationally excited CO due to collisions with He atom are presented for collision energies between 10(-5) and approximately 1000 cm-1 with CO being initially in the vibrational level upsilon=2 and rotational levels j=0,1,4, and 6. The He-CO interaction potential of Heijmen et al. [J. Chem. Phys. 107, 9921 (1997)] was adopted for the calculations. Cross sections for rovibrational transitions and state-to-state rotational energy transfer from selected initial rotational levels were computed and compared with recent measurements of Carty et al. [J. Chem. Phys. 121, 4671 (2004)] and available theoretical results. Comparison in all cases is found to be excellent, providing a stringent test for the scattering calculations as well as the reliability of the He-CO interaction potential by Heijmen et al.     

Pseudorotational dynamics of H3+ and Li3+

The origin of the pseudoprecession phenomenon is investigated through a computational study of the time evolution of H₃⁺ and Li₃⁺ by electron nuclear dynamics theory. In particular, the pseudorotation of both molecules is shown to induce a spatial rotation, which in turn leads to Coriolis coupling of the two orthogonal nuclear shape deformation modes. This effect is rooted in an anisotropy of the molecular ground state potential energy surface that is caused by the interaction between the D_3h ground state and a twofold degenerate first excited state. Computations are performed for a variety of vibrational energies. In addition, the impact of the anharmonicity of the ground state potential surface on the shape deformation modes and the coupling between them is discussed. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Experimental and theoretical study of the rovibrational relaxation of CH4 and CD4 by Ar

Experimental results are reported for the vibrational relaxation of the lowest bending modes of CH₄ and CD₄ br Ar in the temperature range of 140–376 K. Theoretical calculations are carried out in the framework of the semiclassical coupled-states approximation using asymptotic expressions of (3j) symbols and a first-order perturbation treatment. The confrontation of experimental and theoretical rate constants confirms the crucial role of rotational energy transfer upon the vibrational relaxation transfer.     

Photoelectron spectrum of NO2−: SAC-CI gradient study of vibrational-rotational structures

Three-dimensional accurate potential energy surfaces around the local minima of NO₂⁻ and NO₂ were calculated with the SAC/SAC-CI analytical energy gradient method. Therefrom, the ionization photoelectron spectra of NO₂⁻, the equilibrium geometries and adiabatic electron affinity of NO₂, and the vibrational frequencies including harmonicity and anharmonicity of NO₂⁻ and NO₂ were obtained. The calculated electron affinity was in reasonable agreement with the experimental value. The SAC-CI photoelectron spectra of NO₂⁻ at 350 K and 700 K including the rotational effects were calculated using the Franck–Condon approximation. The theoretical spectra reproduced well the fine experimental photoelectron spectra observed by Ervin et al. (J. Phys. Chem. 1988, 92, 5405). The results showed that the ionizations from many vibrational excited states as well as the vibrational ground state are included in the experimental photoelectron spectra especially at 700 K and that the rotational effects are important to reproduce the experimental photoelectron spectra of both temperatures. The SAC/SAC-CI theoretical results supported the analyses of the spectra by Ervin et al., except that we could show some small contributions from the asymmetric-stretching mode of NO₂⁻. © 2018 Wiley Periodicals, Inc.     

Vibrational population lifetimes of polyatomic molecules in liquids

CH-stretching modes were first excited by picosecond infrared pulses and the generated excess population was monitored by anti-Stokes scattering of subsequent ultrashort probe pulses. Experimental data are reported on five molecules: CHCl₃, CH₂Cl₂, CH₃CCl₃, CH₃CH₂OH, and CH₃I in the neat liquid and/or in solutions of CCl₄. The observed time constants vary between 1 and 100 ps depending upon the individual molecule and surrounding. Theoretical calculations show that rotational coupling, Fermi resonance, Coriolis coupling, and resonance energy transfer can strongly effect the vibrational population lifetime. The relevance of these processes is quite different for the various molecules investigated.     

State-to-state rotational transitions in H2+H2 collisions at low temperatures

We present quantum mechanical close-coupling calculations of collisions between two hydrogen molecules over a wide range of energies, extending from the ultracold limit to the superthermal region. The two most recently published potential energy surfaces for the H(2)-H(2) complex, the so-called Diep-Johnson (DJ) [J. Chem. Phys. 112, 4465 (2000); 113, 3480 (2000)] and Boothroyd-Martin-Keogh-Peterson (BMKP) [J. Chem. Phys. 116, 666 (2002)] surfaces, are quantitatively evaluated and compared through the investigation of rotational transitions in H(2)+H(2) collisions within rigid rotor approximation. The BMKP surface is expected to be an improvement, approaching chemical accuracy, over all conformations of the potential energy surface compared to previous calculations of H(2)-H(2) interaction. We found significant differences in rotational excitation/deexcitation cross sections computed on the two surfaces in collisions between two para-H(2) molecules. The discrepancy persists over a large range of energies from the ultracold regime to thermal energies and occurs for several low-lying initial rotational levels. Good agreement is found with experiment B. Mate et al., [J. Chem. Phys. 122, 064313 (2005)] for the lowest rotational excitation process, but only with the use of the DJ potential. Rate coefficients computed with the BMKP potential are an order of magnitude smaller.     

Potential energy functions of polyatomic molecules

A direct method is proposed for determining polyatomic potential energy functions, expressed in terms of normal coordinates, which yield a given set of vibrational excitation energies. The method is a modification of the semiclassical technique for computing vibrational energy levels of Percival and Pomphrey. The technique is used to derive potential functions for the NO₂, SO₂ and ClO₂ molecules. With these potentials twenty two higher vibrational excitations energies have been predicted for these molecules and these results differ from the experimental values by at most 3 cm^?1. The computed potential functions are not unique despite the apparent accuracy of the vibrational energy levels. Comparison with the RKR method indicates that the present method must be extended to include rotational perturbations.     

Exploring the reaction dynamics of O(3P)+ (X2 )→OH+(X3Σ−)+ H(2S) reaction with time‐dependent wave packet method

The O(³P)+ reaction has been investigated by employing time‐dependent quantum wave packet with split operator method on potential energy surface of the doublet ground‐state H₂O⁺(1²A″). The reaction probabilities and integral cross sections are calculated using centrifugal sudden approximation, which basically agree with the quasi‐classical results of Paniagua et al. [Phys. Chem. Chem. Phys. 2014, 16, 23594]. Moreover, the effect of vibrational and rotational excitation of reactant is investigated. The results show that the vibrational and rotational excitation effects on the integral cross section are not obvious. The little differences between Coriolis coupling results and centrifugal sudden approximation ones show that the cheaper centrifugal sudden calculations here reported are effective for this reaction.     

Large Amplitude Torsions in Nitrotoluene Isomers Studied by Rotational Spectroscopy and Quantum Chemistry Calculations

Rotational spectra of ortho-nitrotoluene (2-NT) and para-nitrotoluene (4-NT) have been recorded at low and room temperatures using a supersonic jet Fourier Transform microwave (MW) spectrometer and a millimeter-wave frequency multiplier chain, respectively. Supported by quantum chemistry calculations, the spectral analysis of pure rotation lines in the vibrational ground state has allowed to characterise the rotational energy, the hyperfine structure due to the ¹⁴N nucleus and the internal rotation splittings arising from the methyl group. For 2-NT, an anisotropic internal rotation of coupled −CH₃ and −NO₂ torsional motions was identified by quantum chemistry calculations and discussed from the results of the MW analysis. The study of the internal rotation splittings in the spectra of three NT isomers allowed to characterise the internal rotation potentials of the methyl group and to compare them with other mono-substituted toluene derivatives in order to study the isomeric influence on the internal rotation barrier.     

Electronic and Vibrational Coherence Dynamics in a Cyanine Dye Studied Using a Few‐Cycle Pulsed Laser

We report the relaxation times of electronic and vibrational coherence in the cyanine dye 1,1′,3,3,3′,3′‐hexamethyl‐4,4′,5,5′‐dibenzo‐2,2′‐indotricarbocyanine, measured using a 7.1 fs pulsed laser. The vibrational phase relaxation times are found to be between 380 and 680 fs in the ground and lowest excited singlet states. The vibrational dephasing times of the 294, 446, and 736 cm^?1 modes are relatively long among the six modes associated with excited‐state wave packets. The slower relaxations are explained in terms of a coupled triplet of vibrational modes, which preserves coherence by forming a tightly bound group to satisfy the condition of circa conservation of vibrational energy. Using data from the negative‐time range (i.e., when the probe pulse precedes the pump pulse), the electronic phase relaxation time is found to be 31±1 fs. The dynamic vibrational mode in the excited state (1171 cm^?1), detected in the positive‐time range, is also studied from the negative‐time traces under the same experimental conditions.     

Quantitative molecular thermochemistry based on path integrals

The calculation of thermochemical data requires accurate molecular energies and heat capacities. Traditional methods rely upon the standard harmonic normal-mode analysis to calculate the vibrational and rotational contributions. We utilize path-integral Monte Carlo for going beyond the harmonic analysis and to calculate the vibrational and rotational contributions to ab initio energies. This is an application and an extension of a method previously developed in our group [J. Chem. Phys. 118, 1596 (2003)].     

Dielectric relaxation studies of aromatic polyamides

Dielectric relaxation spectroscopy (DRS) is presented for a family of four aromatic polyamides trying to relate the structure of the lateral groups to the molecular mobility. A prominent sub-T_g absorption is always seen followed in some cases by remanent dielectric activity at room temperature and a subsequent increase of the loss permittivity. The low temperature relaxation is analyzed in terms of a Fuoss–Kirkwood equation to obtain the broadness and the strength of these relaxations as well as the activation energy (ranging from 10 to 11 Kcal/mol). The low frequency conductive peak shows in each case a half-width higher (1.30) than those corresponding to a single relaxation time peak (1.144). These values of the half-width are an indication of the complex character of these phenomena. A final discussion of the rotational barriers of the lateral chains rules out that such motions are the only molecular origin for the gamma relaxation. Instead, some kind of motion involving the main chain and where the interchain interactions play a significant role should be considered as responsible for that relaxation. © 1997 John Wiley & Sons, Inc. J Polym Sci B: Polym Phys 35 : 919–927, 1997     

Molecular and collective relaxations of ferroelectric side chain liquid crystalline polysiloxanes

The influences of mesogenic group chemical structures on dielectric relaxation behavior were investigated for ferroelectric side chain liquid crystalline polymers (FLCPs). The relaxation time and activation energies of the Goldstone mode, α‐, and β‐relaxations decrease with increasing spacer length because of the plasticizer effect of the spacer. Moreover, the relaxation intensity increases with increasing spacer length for FLCPs. An FLCP with a longer spacer length exhibits a higher mesogenic group mobility, and subsequently leads to easier reorientation toward the alternating electrical field. An increase in mesogenic core rigidity results in an increase in the relaxation time and activation energies, and a decrease in the relaxation intensities for the Goldstone mode, α‐, and β‐relaxations. Moreover, the β‐relaxation is suppressed and cannot be observed in the glassy state for FLCPs containing naphthyl biphenylcarboxylate as the mesogenic group. Shorter relaxation time, smaller activation energies, and higher intensity of the α‐, and β‐relaxations were obtained for FLCPs containing chiral moiety with a flexible heptyl alkyl chain. However, the relaxation intensity of the Goldstone mode for FLCPs containing this chiral moiety was smaller than that for FLCPs containing the chiral moiety with a butyl alkyl chain. For FLCPs containing a chiral moiety with two asymmetrical centers, their Goldstone mode relaxation showed larger amplitude. The α‐ and β‐relaxations are suppressed for these FLCPs because of the dense packing and memory effect of the smectic phase. The relationship between the chemical structure of the mesogenic group and dielectric relaxations is discussed in great detail. © 2006 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 44: 2035–2049, 2006     

Low-energy rotational inelastic collisions of H(+) + CO system

The quantum mechanical state-to-state rotational excitation cross sections have been computed using the ab initio ground electronic state potential energy surface of the system [M. Mladenovic and S. Schmatz, J. Chem. Phys. 109, 4456 (1998)] computed at coupled-cluster single and double and triple perturbative excitations method using correlation-consistent polarized valence quadruple zeta basis set where the asymptotic potential have been computed using the dipole moment, quadrupole moment, and the molecular polarizability components and fitted to this interaction potential. The anisotropy of the surface has been analyzed in terms of the multipolar expansion coefficients for the rigid-rotor surface. The integral cross sections for rotational excitations have been computed by solving close-coupled equations at very low collision energies (5-200 cm(-1)) and the corresponding rates have been obtained for a range of low temperatures (5-175 K). The j = 0 → j(') = 1 rotational excitation cross section (and rate) is found to be the dominant followed by the j = 0 → j(') = 2 in these collision energies. The close-coupling, coupled-state, and infinite-order sudden approximations coupling calculations have been performed in the energy range of 0.1-1.0 eV using vibrational ground potential. The rotational cross sections have been obtained by performing computationally accurate close-coupling calculations at 0.1 eV using vibrationally averaged potential (ν = 1) and compared with the results of vibrational ground potential.     

A theoretical study on CH<Subscript>2</Subscript>N<Subscript>2</Subscript> isomers: structure and energetics

Five CH₂N₂ isomers, namely cyanamide, carbodiimide, diazomethane, isocyanamide and nitrilimine, have been investigated at a high level of accuracy. The singles and doubles coupled-cluster method including a perturbational correction for connected triple excitations, CCSD(T), in conjunction with correlation-consistent basis sets ranging in size from triple to quintuple zeta have been employed. Extrapolation to the complete basis set limit has been used with treatments of core-valence correlation effects in order to accurately predict structures, relative energies as well as N–H and C–H bond dissociation energies. The latter required to also investigate the HNNC radical with the same methodology used for CH₂N₂ isomers, while HCNN and HNCN data are available in the literature by the same authors (Puzzarini and Gambi in J Chem Phys 122:064316, 2005). For all the species studied, harmonic vibrational frequencies have also been evaluated at the CCSD(T) level in order to obtain zero-point corrections to total energies.