Terahertz spectrum and normal-mode relaxation in pentaerythritol tetranitrate: effect of changes in bond-stretching force-field terms

Terahertz (THz) active normal-mode relaxation in crystalline pentaerythritol tetranitrate (PETN) was studied using classical molecular dynamics simulations for energy and density conditions corresponding to room temperature and atmospheric pressure. Two modifications to the fully flexible non-reactive force field due to Borodin et al. [J. Phys. Chem. B 112, 734 (2008)] used in a previous study of THz-active normal-mode relaxation in PETN [J. Chem. Phys. 134, 014513 (2011)] were considered to assess the sensitivity of the earlier predictions to details of the covalent bond-stretching terms in the force field. In the first modification the harmonic bond-stretching potential was replaced with the Morse potential to study the effect of bond anharmonicity on the THz-region mode relaxation. In the second modification the C-H and nitro-group N-O bond lengths were constrained to constant values to mimic lower quantum occupation numbers for those high-frequency modes. The results for relaxation times of the initially excited modes were found to be insensitive to either force-field modification. Overall time scales for energy transfer to other modes in the system were essentially unaffected by the force-field modifications, whereas the detailed pathways by which the energy transfer occurs are more complicated for the Morse potential than for the harmonic-bond and fixed-bond cases. Terahertz infrared absorption spectra constructed using calculated normal-mode frequencies, transition dipoles, and relaxation times for THz-active modes were compared to the spectra obtained from the Fourier transform of the dipole-dipole time autocorrelation function (DDACF). Results from the two approaches are in near agreement with each other and with experimental results in terms of main peak positions. Both theoretical methods yield narrower peaks than observed experimentally and in addition predict a weaker peak at ω ～ 50 cm(-1) that is weak or absent experimentally. Peaks obtained using the DDACF approach are broader than those obtained from the normal-mode relaxation method.     

The nonaqueous microtitration of pentaerythritol tetranitrate (PETN)

Summary PETN has been determined by nonaqueous titration as an acid in pyridine and l, l,3,3-tetramethylguanidine. The titrant was tetrabutylammonium hydroxide. Endpoints were determined potentiometrically by a glass/modified calomel electrode system, or by two polarized platinum electrodes. The method may be applicable to other very weakly acidic compounds for which no other titrimetric methods are currently available.

---

Zusammenfassung Pentaeerythrit-Tetranitrat wurde als Säure gelöst in Pyridin und 1,1,3,3-Tetramethylguanidin mit Tetrabutylammoniumhydroxid titriert. Der Endpunkt wurde potentiometrisch mit dem System Glaselektrode/modifizierte Kalomelelektrode oder mit zwei polarisierten Platinelektroden bestimmt. Das Verfahren dürfte sich auch für andere sehr schwache Säuren eignen, für die andere maßanalytische Methoden zur Zeit nicht verfügbar sind.

---

---

Work performed under the auspices of the U. S. Department of Energy by the Lawrence Livermore Laboratory under contract number W-7405-ENG-48.     

季戊四醇四硝酸酯的分子结构和热解机理

运用SCF-AM1-MO方法，在RHF和UHF水平上分别全优化计算了季戊四醇四硝酸酯(太安)的分子结构和热解机理。计算所得三种构象的能量相差很小、且相互转化的能垒也较小。其均裂O-NO~2键(生成双自由基)的热解引发反应活化能较小，而通过α-H转移产生RCHO和HONO的环消除反应需要较高活化能。     

Atoms-in-molecules study of intra- and intermolecular bonding in the pentaerythritol tetranitrate crystal

Chemical bonding in the pentaerythritol tetranitrate crystal based on the experimental electron density obtained from X-ray diffraction data at 100 K and theoretical calculations at the experimental molecular geometry have been analyzed in terms of the Quantum Theory of Atoms in Molecules. Features of the intra- and intermolecular bond critical points and the oxygen atom lone-pair locations are discussed. Numerous intermolecular bonding interactions, including O...H and O...O, have been found and characterized. Atomic charges and atomic energies were integrated and compared with those for similar compounds. The N-O topological bond orders have been calculated for the first time, and the PETN atomic valences have been estimated.     

Ab initio all-electron periodic Hartree-Fock study of hydrostatic compression of pentaerythritol tetranitrate

We present a computational study of hydrostatic compression effects on the pentaerythritol tetranitrate (PETN) energetic material up to 22.7 GPa by means of the ab initio all-electron periodic Hartree-Fock quantum mechanical method with the STO-3G Gaussian basis set. We fitted the calculated volume-energy relation to the energy SJEOS polynomial function from which we obtained the compression dependence of the pressure (P), the bulk modulus (B), and its pressure derivative (B'). We also fitted the experimental volume-pressure relation to the pressure SJEOS polynomial function, which allowed us to calculate the experimental bulk modulus (B(exp)) and its pressure derivative (). Our calculated values, B = 6.73 GPa and B' = 24.63, are in reasonable agreement with the values B(exp) = 8.48 GPa and = 14.42 from our fit to the experimental X-ray data and with the value B(exp) = 9.8 GPa that was derived from the experimental elastic constants. In addition, we present a discussion on how the lattice vectors and the internal coordinates (i.e., bond lengths, bond angles, and torsion angles) of the C(CH(2)ONO(2))(4) molecules in the PETN lattice change during hydrostatic compression of the crystal. Our calculated results suggest that the C(CH(2)ONO(2))(4) molecules cannot be considered as being rigid but are in fact flexible, accommodating lattice compression through torsions, bendings in their bond angles, and contractions in their bond lengths. At pressures higher than about 8 GPa, however, both the C(CH(2)ONO(2))(4) molecules and the c lattice vector seem to stiffen somewhat. The a lattice vector does not exhibit this stiffening. As a consequence, the pressure dependence of the c/a ratio shows a minimum at about 8 GPa.     

Raman spectra of shock compressed pentaerythritol tetranitrate single crystals: anisotropic response

To gain insight into the anisotropic sensitivity of shocked pentaerythritol tetranitrate (PETN) single crystals, single-pulse Raman spectroscopy was used to examine the response of crystals shocked along the [100] (insensitive) and [110] (sensitive) orientations. High-resolution Raman spectra revealed several orientation-dependent features under shock compression: (i) substantially different stress dependence of the Raman shift for the CH(2) and NO(2) stretching modes for the two orientations, (ii) discontinuity in the stress dependence of the Raman shift for the CH(2) stretching modes above 4 GPa for the [110] orientation, and (iii) large broadening for the CH(2) and NO(2) asymmetric stretching modes for stresses above 4 GPa for the [110] orientation. The present data in combination with previous static pressure results provide support for conformational changes in PETN molecules for shock compression along the [110] (sensitive) orientation. Implications of the present results for the anisotropic sensitivity of shocked PETN are discussed.     

Melting of pentaerythritol tetranitrate (PETN) nanoconfined in controlled pore glasses (CPG)

The melting temperature depression of pentaerythritol tetranitrate, nanoconfined in controlled pore glasses (CPG), was systematically studied by differential scanning calorimetry (DSC). The solid–liquid interfacial energy σ _sl was obtained from the Gibbs–Thomson equation fit to melting temperature vs. reciprocal pore diameter. The pore size distribution of the CPG pores was also estimated from the DSC data. Pore sizes obtained from the manufacturer by BET are compared with those determined from the DSC curves using either the curves directly or by assuming spherical shaped confining cavities. The thermal mass vales are in better agreement with the BET estimation.     

Thermodynamic analysis of pure and impurity doped pentaerythritol tetranitrate crystals grown at room temperature

Pentaerythritol tetranitrate (PETN) powders are used to initiate other explosives. During long-term storage, changes in powder properties can cause changes in the initiation performance. Changes in the morphology and surface area of aging powders are observed due to sublimation and growth of PETN crystals through coarsening mechanisms, (e.g. Ostwald ripening, sintering, etc.). In order to alleviate the sublimation of PETN crystals under service conditions, stabilization methods such as thermal cycling and doping with certain impurities during or after the crystallization of PETN have been proposed. In this report we present our work on the effect of impurities on the morphology and activation energy of the PETN crystals. The pure and impurity doped crystals of PETN were grown from supersaturated acetone solution by solvent evaporation technique at room temperature. The difference in the morphology of the impurity-doped PETN crystal compared to pure crystal was examined by optical microscopy. The changes in the activation energies and the evaporation rates are determined by thermogravimetry (TG). Our activation energies of evaporation agree with earlier reported enthalpies of vaporization. The morphology and activation energy of PETN crystals doped with Ca, Na, and Fe cations are similar to that for pure PETN crystal, whereas the Zn-ion-doped PETN crystals have different morphology and decreased activation energy.     

Spectrophotometric determination of pentaerythritol tetranitrate (PETN) in waste water from lead styphnate primer plants

An accurate spectrophotometric method is proposed for the determination of pentaerythritol tetranitrate (PETN) in waste water from lead styphnate primer plants by use of phenoldisulphonic acid. The waste water is filtered through a sintered glass crucible and the PETN is determined in the filtrate and the residue. In the determination of PETN in the filtrate, sodium hydroxide is added and the PETN is extracted with methylene chloride (in alkaline solution, styphnate and TNT are not extracted). The methylene chloride solution is then evaporated to dryness, the residue is treated with a solution of phenoldisulphonic acid in sulphuric acid, water and ammonia are added, and the yellow colour is measured. In the determination of PETN in the residue, the PETN is dissolved in acetone, an aliquot of the acetone solution is treated with water and sodium hydroxide, the PETN is extracted with methylene chloride and the colour is developed as above. Various factors affecting the determination were investigated. The solubility of PETN in water was studied.     

Orientation dependent far-infrared terahertz absorptions in single crystal pentaerythritol tetranitrate (PETN) using terahertz time-domain spectroscopy

Terahertz time-domain spectroscopy (THZ-TDS) has been used to measure the absorption spectra in the range 7-100 cm(-1) (0.2-3 THz) of single crystal pentaerythritol tetranitrate (PETN). Absorption was measured in transmission mode as a function of incident polarization with the incident and transmitted wave vectors oriented along the crystallographic directions [100], <10(a/c)(2)>, and <110>. Samples were rotated with respect to the incident polarization while absorption was measured at both 300 and 20 K. Comparatively minor differences were observed among the three orientations. Two broad absorptions at 72 and >90 cm(-1), and several weaker absorptions at 36, 55, 80, and 82 cm(-1), have been observed at cryogenic temperatures.     

Decomposition of pentaerythritol tetranitrate [C(CH2ONO2)4] following electronic excitation

We report the experimental and theoretical study of the decomposition of gas phase pentaerythritol tetranitrate (PETN) [C(CH(2)ONO(2))(4)] following electronic state excitation. PETN has received major attention as an insensitive, high energy explosive; however, the mechanism and dynamics of the decomposition of this material are not clear yet. The initial decomposition mechanism of PETN is explored with nanosecond energy resolved spectroscopy and quantum chemical theory employing the ONIOM algorithm at the complete active space self-consistent field (CASSCF) level. The nitric oxide (NO) molecule is observed as an initial decomposition product from PETN at three UV excitation wavelengths (226, 236, and 248 nm) with a pulse duration of 8 ns. Energies of the three excitation wavelengths coincide with the (0-0), (0-1), and (0-2) vibronic bands of the NO A (2)Σ(+) ← X (2)Π electronic transition, respectively. A unique excitation wavelength independent dissociation channel is observed for PETN, which generates the NO product with a rotationally cold (~20 K) and a vibrationally hot (~1300 K) distribution. Potential energy surface calculations at the ONIOM(CASSCF:UFF) level of theory illustrate that conical intersections play an important role in the decomposition mechanism. Electronically excited S(1) PETN returns to the ground state through the (S(1)/S(0))(CI) conical intersection, and undergoes a nitro-nitrite isomerization to generate the NO product.     

Local mode relaxation in polychlorotrifluoroethylene

The effects of temperature and pressure on the dielectric relaxation time τ, relaxation strength Δε, distribution of relaxation time β , and shape of absorption curve of the local mode relaxation in polychlorotrifluoroethylene are discussed from the viewpoint of intermolecular interaction. Changes of τ and Δε with temperature and pressure are explainable in terms of changes in intermolecular distance. The temperature and pressure dependences of the dielectric retardation spectrum for the local mode relaxation are different from each other.     

Raman scattering studies of the high-pressure stability of pentaerythritol tetranitrate, C(CH2ONO2)4

High-pressure Raman scattering studies have been performed on a crystalline energetic material, pentaerythritol tetranitrate C(CH2ONO2)4 (PETN), an important secondary explosive. In situ, ambient-temperature investigations employed diamond anvil cell techniques and nitrogen as a quasi-hydrostatic-pressure-transmitting medium. The pressure-induced alterations in the profiles of the Raman lines, including positions, bandwidths, and intensities, were studied in a compression sequence up to about 31.3 GPa and in a subsequent decompression to ambient conditions. The observed changes of the Raman spectra implied that PETN gradually densified and compressed smoothly up to the highest investigated pressures. Compression below 12 GPa gradually shifted all Raman peaks to higher frequencies without significantly changing their relative intensities or bandwidths. At higher pressures, the peak intensities of the Raman spectra decreased considerably and the bands broadened significantly. The Raman spectrum of the material quenched from 31.3 GPa to ambient conditions indicated that no pressure-driven permanent reconstructive modification or decomposition of the PETN structure occurred. That is, the spectral changes were completely reversible upon compression and subsequent decompression to ambient conditions.     

Analysis of 2,4,6-trinitrotoluene, pentaerythritol tetranitrate and cyclo-1,3,5-trimethylene-2,4,6-trinitramine using negative corona discharge ion mobility spectrometry

In this study, the capability of negative corona discharge ion mobility spectrometry (IMS) for quantitative magnitude of several explosives including 2,4,6-trinitrotoluene (TNT), pentaerythritol tetranitrate (PETN) and cyclo-1,3,5-trimethylene-2,4,6-trinitramine (RDX) has been evaluated for the first time. The total current obtained with the negative corona discharge was about 100 times larger than that of IMS based on ⁶³Ni, which results in a lower detection limit and a wider linear dynamic range. The detection limits for PETN, TNT and RDX were 8×10⁻¹¹, 7×10⁻¹¹ and 3×10⁻¹⁰ g, respectively. The calibration plots for these explosives showed linear dynamic ranges of about four orders of magnitude.     

基于太赫兹透射光谱的乙醇水溶液Debye模型研究

利用太赫兹透射光谱测量了液态乙醇和液态水以及不同浓度乙醇水溶液在22℃的介电常数.并利用Levenberg-Marquardt算法拟合得到了它们的Debye模型,该模型包含3个弛豫过程和1个分子间伸缩振动模式.其中,慢速弛豫模式的强度(弛豫时间20~160ps)贡献了主要的介电强度,中间弛豫模式与其密切相关.而快速弛豫模式(弛豫时间约为1ps)只占了大约5%的介电常数.     

Examination of the in vitro degradation of [14C] pentaerythritol tetranitrate in rat and human blood with an improved thin-layer radiochromatographic procedure

Improvements were made on a reported thin-layer radiochromatographic assay for the determination of [14C]pentaerythritol tetranitrate (PETN) and its metabolites in whole blood, using methanol instead of dioxane as the extracting solvent. Recovery of total radioactivity for the entire work-up procedure was greater than 90%, and the distribution of PETN and its metabolites in degraded blood samples was found to be reproducible. This modified method appeared simpler and yielded better recovery of radioactivity than the literature method. In vitro metabolism of [14C]PETN in rat and human blood was examined by incubation of the drug with fresh blood at 37 degrees C for 60 min. In rat blood, the half-life of PETN degradation was about 15 min producing the trinitrate, dinitrate and mononitrate metabolites. Human blood was also capable of degrading PETN in vitro, but at a lower rate than rat blood, yielding only the trinitrate metabolite in quantifiable amounts during the incubation period. Equilibrium of PETN between plasma and red blood cells was observed within 1 min after PETN addition to both rat and human blood. The apparent plasma/red blood cells partition ratios of PETN were 1.1 and 1.7 for rat and human blood, respectively. PETN degradation was approximately ten times slower in rat plasma than in rat blood, suggesting that enzymes in erythrocytes are important for PETN metabolism in rat whole blood.     

Dipole relaxation in normal aliphatic alcohols

A new mechanism of dipole relaxation and molecular motion in normal aliphatic alcohols is suggested, based on the assumption of the existence of consequent reorientations of the molecules in chains It is shown that the experimental data obtained by different methods of dielectric spectroscopy, including TDS. can be described accurately within the framework of the model.     

Vibrational relaxation of the C-X stretching mode in halobenzenes

The Raman linewidths of the carbon-halogen stretching mode in several halobenzenes were measured as a function of temperature and in dilute cyclohexane solution. It was found that, as in an earlier investigation on haloalkanes, the vibrational relaxation efficiencies, τ^?1_iso, vary in the order Cl > Br > I. The temperature dependence of τ_iso, together with the results in dilute solution, provide evidence that dipolar interactions play a significant role in vibrational relaxation in these systems.     

Molecular vibration mode assignment of nematic liquid crystal 5CB on Terahertz spectra

Using DFT/B3LYP/6-311++G** method, the molecular structure and absorption spectra in terahertz (THz) range of liquid crystal 5CB are investigated. In a frequency range 0–15 THz, an assignment of the vibrational modes corresponding to absorption frequencies is performed using potential energy distribution for the first time. It is found that the cyano group radical (–CN) do actively take part in the strongest THz absorption of 1.743, 3.942, 5.169 and 14.769 THz in different vibration modes. The results suggest that the strong polar group should be avoided in designing liquid crystal molecule and mixtures in order to reduce the absorption intensity in THz range.