Atomistic-scale simulations of the initial chemical events in the thermal initiation of triacetonetriperoxide

To study the initial chemical events related to the detonation of triacetonetriperoxide (TATP), we have performed a series of molecular dynamics (MD) simulations. In these simulations we used the ReaxFF reactive force field, which we have extended to reproduce the quantum mechanics (QM)-derived relative energies of the reactants, products, intermediates, and transition states related to the TATP unimolecular decomposition. We find excellent agreement between the QM-predicted reaction products and those observed from 100 independent ReaxFF unimolecular MD cookoff simulations. Furthermore, the primary reaction products and average initiation temperature observed in these 100 independent unimolecular cookoff simulations match closely with those observed from a TATP condensed-phase cookoff simulation, indicating that unimolecular decomposition dominates the thermal initiation of the TATP condensed phase. Our simulations demonstrate that thermal initiation of condensed-phase TATP is entropy-driven (rather than enthalpy-driven), since the initial reaction (which mainly leads to the formation of acetone, O(2), and several unstable C(3)H(6)O(2) isomers) is almost energy-neutral. The O(2) generated in the initiation steps is subsequently utilized in exothermic secondary reactions, leading finally to formation of water and a wide range of small hydrocarbons, acids, aldehydes, ketones, ethers, and alcohols.     

Dissociative recombination of the weakly bound NO-dimer cation: cross sections and three-body dynamics

Dissociative recombination (DR) of the dimer ion (NO)(2) (+) has been studied at the heavy-ion storage ring CRYRING at the Manne Siegbahn Laboratory, Stockholm. The experiments were aimed at determining details on the strongly enhanced thermal rate coefficient for the dimer, interpreting the dissociation dynamics of the dimer ion, and studying the degree of similarity to the behavior in the monomer. The DR rate reveals that the very large efficiency of the dimer rate with respect to the monomer is limited to electron energies below 0.2 eV. The fragmentation products reveal that the breakup into the three-body channel NO+O+N dominates with a probability of 0.69+/-0.02. The second most important channel yields NO+NO fragments with a probability of 0.23+/-0.03. Furthermore, the dominant three-body breakup yields electronic and vibrational ground-state products, NO(upsilon=0)+N((4)S)+O((3)P), in about 45% of the cases. The internal product-state distribution of the NO fragment shows a similarity with the product-state distribution as predicted by the Franck-Condon overlap between a NO moiety of the dimer ion and a free NO. The dissociation dynamics seem to be independent of the NO internal energy. Finally, the dissociation dynamics reveal a correlation between the kinetic energy of the NO fragment and the degree of conservation of linear momentum between the O and N product atoms. The observations support a mechanism in which the recoil takes place along one of the NO bonds in the dimer.     

Ab initio calculation of the NH(3sigma-)-NH(3sigma-) interaction potentials in the quintet, triplet, and singlet states

We present the ab initio potential-energy surfaces of the NH-NH complex that correlate with two NH molecules in their 3sigma- electronic ground state. Three distinct potential-energy surfaces, split by exchange interactions, correspond to the coupling of the S(A) = 1 and S(B) = 1 electronic spins of the monomers to dimer states with S = 0, 1, and 2. Exploratory calculations on the quintet (S = 2), triplet (S = 1), and singlet (S = 0) states and their exchange splittings were performed with the valence bond self-consistent-field method that explicitly accounts for the nonorthogonality of the orbitals on different monomers. The potential surface of the quintet state, which can be described by a single Slater determinant reference function, was calculated at the coupled cluster level with single and double excitations and noniterative treatment of the triples. The triplet and singlet states require multiconfiguration reference wave functions and the exchange splittings between the three potential surfaces were calculated with the complete active space self-consistent-field method supplemented with perturbative configuration interaction calculations of second and third orders. Full potential-energy surfaces were computed as a function of the four intermolecular Jacobi coordinates, with an aug-cc-pVTZ basis on the N and H atoms and bond functions at the midpoint of the intermolecular vector R. An analytical representation of these potentials was given by expanding their dependence on the molecular orientations in coupled spherical harmonics, and representing the dependence of the expansion coefficients on the intermolecular distance R by the reproducing kernel Hilbert space method. The quintet surface has a van der Waals minimum of depth D(e) = 675 cm(-1) at R(e) = 6.6a0 for a linear geometry with the two NH electric dipoles aligned. The singlet and triplet surfaces show similar, slightly deeper, van der Waals wells, but when R is decreased the weakly bound NH dimer with S = 0 and S = 1 converts into the chemically bound N2H2 diimide (also called diazene) molecule with only a small energy barrier to overcome.     

End-evaporation dynamics revisited

We present analytical results on the so-called end-evaporation kinetics in equilibrium polymeric systems following a temperature jump (T jump). A T jump prepares the system with a nonequilibrium length distribution, after which it relaxes back to its equilibrium state. Starting from a master equation, we develop a mean-field analytical theory based on a generating function approach, which allows explicit approximate expressions for the monomer and dimer concentrations to be derived in a discrete setting; the concentrations of the other chains as well as the average chain length were shown to be entirely expressible in terms of the monomer and dimer concentrations. We find that the calculated monomer and dimer concentrations as well as the average chain length are in good agreement with numerical simulation results and do not suffer from some of the defects of earlier continuum theories. Furthermore, the relaxation was shown to take place in three different stages. The first stage comprises the very fast relaxation of the monomers to almost their equilibrium concentration; the other polymer chains have hardly relaxed. During the second stage, which is highly nonlinear, a redistribution of material at practically constant monomer density takes place. Only in the final stage of the relaxation process the chain concentrations approach their true equilibrium values. In this stage there are only very small shifts in the concentrations of chains, which are governed by extremely slow "indirect" monomer-mediated processes.     

Phase behavior and interfacial properties of nonadditive mixtures of Onsager rods

Within a second virial theory, we study bulk phase diagrams as well as the free planar isotropic-nematic interface of binary mixtures of nonadditive thin and thick hard rods. For species of the same type, the excluded volume is determined only by the dimensions of the particles, whereas for dissimilar ones it is taken to be larger or smaller than that, giving rise to a nonadditivity that can be positive or negative. We argue that such a nonadditivity can result from modeling of soft interactions as effective hard-core interactions. The nonadditivity enhances or reduces the fractionation at isotropic-nematic (IN) coexistence and may induce or suppress a demixing of the high-density nematic phase into two nematic phases of different composition (N(1) and N(2)), depending on whether the nonadditivity is positive or negative. The interfacial tension between coexisting isotropic and nematic phases shows an increase with increasing fractionation at the IN interface, and complete wetting of the IN(2) interface by the N(1) phase upon approach of the triple-point coexistence. In all explored cases bulk and interfacial properties of the nonadditive mixtures exhibit a striking and quite unexpected similarity with the properties of additive mixtures of different diameter ratio.     

Furanosyl Oxocarbenium Ion Conformational Energy Landscape Maps as a Tool to Study the Glycosylation Stereoselectivity of 2-Azidofuranoses, 2-Fluorofuranoses and Methyl Furanosyl Uronates

The 3D shape of glycosyl oxocarbenium ions determines their stability and reactivity and the stereochemical course of S_N1 reactions taking place on these reactive intermediates is dictated by the conformation of these species. The nature and configuration of functional groups on the carbohydrate ring affect the stability of glycosyl oxocarbenium ions and control the overall shape of the cations. We herein map the stereoelectronic substituent effects of the C2-azide, C2-fluoride and C4-carboxylic acid ester on the stability and reactivity of the complete suite of diastereoisomeric furanoses by using a combined computational and experimental approach. Surprisingly, all furanosyl donors studied react in a highly stereoselective manner to provide the 1,2-cis products, except for the reactions in the xylose series. The 1,2-cis selectivity for the ribo-, arabino- and lyxo-configured furanosides can be traced back to the lowest-energy ³E or E₃ conformers of the intermediate oxocarbenium ions. The lack of selectivity for the xylosyl donors is related to the occurrence of oxocarbenium ions adopting other conformations.     

Destructive adsorption of CCl4 over lanthanum-based solids: linking activity to acid-base properties

The relative activities of a low-surface crystalline and high-surface amorphous LaOCl, further denoted as S1 and S2, have been compared for the destructive adsorption of CCl4. It was found that the intrinsic activity of S2 is higher than that of S1. Both samples were characterized with X-ray diffraction (XRD), X-ray photoelectron spectroscopy (XPS), N2-physisorption, and Raman and infrared (IR) spectroscopy. IR was used in combination with CO2, CO, and methanol as probe molecules. The CO2 experiments showed that different carbonate species are formed on both materials. For S1, a high surface concentration of bidentate carbonate species and a lower concentration of monodentate carbonate were observed. In the case of S2, bulk carbonates were present together with bridged carbonates. CO adsorption shows that S2 and S1 have very similar Lewis acid sites. However, methanol adsorption experiments showed that S2 had a higher number of stronger Lewis acid sites than S1 and that twofold coordinated methoxy species were more strongly bound than threefold coordinated methoxy species. Because of the analogy between methanol dissociation and the removal of the first chlorine atom in the destructive adsorption of CCl4, the sites enabling twofold coordination were likely to be the same Lewis acid sites actively involved in the destructive adsorption of CCl4. La2O3 was less active than the two LaOCl materials, and therefore, the intrinsic activity of the catalyst increases as the strength of the Lewis acid sites increases. S2 contains more chlorine at the surface than S1, which is expressed by the higher number of sites enabling twofold coordination. Moreover, this explains the difference in destructive adsorption capacity for CCl4 that was observed for the samples S1 and S2. Since LaCl3, being the most acidic phase, is not active for the destructive adsorption of CCl4, basic oxygen atoms, however, remain needed to stabilize the reaction intermediate CCl3 as La-O-CCl3.     

Excitation decay pathways of Lhca proteins: a time-resolved fluorescence study

Light-harvesting complex I (LHCI), which serves as a peripheral antenna for photosystem I (PSI) in green plants, consists mainly of four polypeptides, Lhca1-4. We report room temperature emission properties of individual reconstituted monomeric Lhca proteins (Lhca1, -2, -3, and -4) and dimeric Lhca1/4, performed by steady-state and time-resolved fluorescence techniques. The emission quantum yields of the samples are approximately 0.12, 0.085, 0.081, 0.041, and 0.063 for Lhca1, -2, -3, -4, and the -1/4 dimer, respectively, which is considerably lower than the value of 0.22 found for light-harvesting complex II (LHCII), the main peripheral antenna complex of photosystem II in green plants. The decay components of LHCI proteins can be divided in two categories: Lhca1 and Lhca3 have decay times of 1.1-1.6 ns and 3.3-3.6 ns, and Lhca2 and Lhca4 have decay times of 0.7-0.9 ns and 3.1-3.2 ns. These categories seem to correlate with the pigment composition of the samples. All decay times are faster than that observed previously for LHCII. When the absolute emission yields and the lifetimes of the Lhca samples are combined, the overall emission properties of the individual Lhca proteins are expressed in terms of their emitting dipole moment strength. In the samples without extreme red states, that is, Lhca1 and Lhca2, the emitting dipole moment has a value close to unity (relative to monomeric chlorophyll in acetone), which is similar to that for LHCII, whereas, in the samples with the red-most state (F-730), that is, Lhca3, -4, and the -1/4 dimer, the emitting dipole moment has a value less than unity (0.6-0.8), which can be explained by mixing the red-most (exciton) state with a dark charge-transfer state, as suggested in previous PSI red pigment studies. In addition, we find a lifetime component of approximately 50-150 ps in all red-pigment-containing samples, which cannot be due to "slow" energy transfer, but is instead assigned to an unrelaxed state of the pigment-protein, which, on this time-scale, is converted into the final emitting state.     

Green's-function reaction dynamics: a particle-based approach for simulating biochemical networks in time and space

We have developed a new numerical technique, called Green's-function reaction dynamics (GFRD), that makes it possible to simulate biochemical networks at the particle level and in both time and space. In this scheme, a maximum time step is chosen such that only single particles or pairs of particles have to be considered. For these particles, the Smoluchowski equation can be solved analytically using Green's functions. The main idea of GFRD is to exploit the exact solution of the Smoluchoswki equation to set up an event-driven algorithm, which combines in one step the propagation of the particles in space with the reactions between them. The event-driven nature allows GFRD to make large jumps in time and space when the particles are far apart from each other. Here, we apply the technique to a simple model of gene expression. The simulations reveal that spatial fluctuations can be a major source of noise in biochemical networks. The calculations also show that GFRD is highly efficient. Under biologically relevant conditions, GFRD is up to five orders of magnitude faster than conventional particle-based techniques for simulating biochemical networks in time and space. GFRD is not limited to biochemical networks. It can also be applied to a large number of other reaction-diffusion problems.