Thermal decomposition of monomethylhydrazine: Shock tube experiments and kinetic modeling

The thermal decomposition of gaseous monomethylhydrazine (MMH) was studied by recording MMH absorption at 220 nm of the reacting gas behind a reflected shock wave at temperatures of 900–1370 K, pressures of 140–450 kPa, and in mixtures containing 97.5–99 mol% argon. Based on previous work (Sun and Law; J Phys Chem A 2007, 111(19), 3748–3760), a kinetic mechanism was developed over extended temperature and pressure ranges to model these experimental data. Specifically, the temperature and pressure dependence of the unimolecular rate coefficients on the dissociation of MMH and the associated radicals were calculated by the QRRK/Master equation analysis at temperatures of 300–2000 K and pressures of 1–100 atm based on published thermochemical and kinetic parameters. They were then fitted using the Troe formalism and incorporated in the kinetic model. This unadjusted model was then used to predict the MMH decomposition profiles at different temperatures and pressures for seven groups of MMH/Ar mixtures and the half‐life decomposition times from shock tube experiments. Good agreement was observed below 940 K and above 1150 K for the diluted MMH/Ar mixtures. The model predictions further show that the overall MMH decomposition rate follows first‐order kinetics, and that the N–N bond scission is the most sensitive reaction path for the modeling of the homogeneous decomposition of MMH at elevated pressures. However, the model predictions deviate from the experimental data with the incubation period of ca. 100 μs observed in the 1030–1090 K temperature range, and it also predicts longer ignition delays for highly concentrated MMH/Ar mixtures. The discrepancy between the model predictions and experimental data at these special conditions of MMH decomposition was analyzed. © 2008 Wiley Periodicals, Inc. Int J Chem Kinet 41: 176–186, 2009     

Weak interactions and molecular recognition in systems involving electron transfer proteins

Electrostatic interactions and other weak interactions between amino acid side chains on protein surfaces play important roles in molecular recognition, and the mechanism of their intermolecular interactions has gained much interest. We established that charged peptides are useful for investigating the molecular recognition character of proteins and their molecular interaction induced structural changes. Positively charged lysine peptides competitively inhibited electron transfer from reduced cytochrome f (cyt f or cytochrome c (cyt c) to oxidized plastocyanin (PC), due to neutralization of the negatively charged site of PC by formation of PC-lysine peptide complexes. Lysine peptides also inhibited electron transfer from cyt c to cytochrome c peroxidase. Likewise, negatively charged aspartic acid peptides interacted with the positively charged sites of cytfand cyt c, and competitively inhibited electron transfer from reduced cytfor cyt c to oxidized PC and from [Fe(CN)6]4- to oxidized cyt c. Changes in the geometry and a shift to a higher redox potential of the active site Cu of PC on oligolysine binding were detected by spectroscopic and electrochemical measurements, owing to the absence of absorption in the visible region for lysine peptides. Structural and redox potential changes were also observed for cyt f and cyt c by interaction with aspartic acid peptides.     

Classical and interdisciplinary approaches to the design of organic and hybrid molecular systems

This review summarizes the recent advances in the construction of organic and hybrid systems having considerable potential for practical applications in pharmacology, medicine, materials science, petrochemistry, and chemistry of high-energy compounds. The methodological classification of the approaches into classical and interdisciplinary is proposed. Examples of the efficient application of these approaches for investigating complex molecular systems, as well as for developing new methods of synthesis are given.     

Perturbation theory for degenerate states of atomic and molecular systems

The first few terms of the perturbation expansions for the function and the energy shift of a degenerate state of an arbitrary quantum mechanical system are obtained using the adiabatic formula. It is shown how the expansion for the secular operator may be obtained from the expansion for the function. The results are used to calculate energies of the ground and some excited states and multiplet splittings of some beryllium-like ions.     

Thermodynamic and kinetic theory of ionic micellar systems: 2. Statistical-thermodynamic relations

The role of electrochemical potentials in the grand canonical ensemble of ionic micellar systems is characterized. The notion of relative electrochemical potentials is introduced with allowance for the electroneutrality condition. Fundamental relations and primary statistical-thermodynamic relations are derived for ideal and real ionic micellar systems with participation of electrochemical potentials, in which inaccuracies observed in published literature, are eliminated. A differential equation for the osmotic pressure of ionic aggregated system is obtained. Relations that link the work of the aggregation of ionic micelle with chemical and electrochemical potentials and aggregation numbers are established. Separate contributions to the work of aggregation are commented on.     

Operator methods in full configuration interaction theory for molecular systems

We discuss modern trends in the theory and practice of full configuration interaction calculations. We pay the most attention to the wave operator method, in which the wave function is considered as the kernel of a many-particle operator. The corresponding operator equation, equivalent to the Schrödinger equation, automatically leads to a convenient matrix algorithm. We also discuss an alternative approach based on the pairing operator, generalizing the construction of the wave function in the method of one-particle spin-pairing amplitudes.Kharkov University. Translated from Teoreticheskaya i Éksperimental'naya Khimiya, Vol. 27, No. 4, pp. 413–426, July–August, 1991. Original article submitted January 14, 1991.     

Koopmans' theorem for large molecular systems within density functional theory

It is shown that in density functional theory (DFT), Koopmans' theorem for a large molecular system can be stated as follows: The ionization energy of the system equals the negative of the highest occupied molecular orbital (HOMO) energy plus the Coulomb electrostatic energy of removing an electron from the system, or equivalently, the ionization energy of an N-electron system is the negative of the arithmetic average of the HOMO energy of this system and the lowest unoccupied molecular orbital (LUMO) energy of the (N - 1)-electron system. Relations between this DFT Koopmans' theorem and its existing counterparts in the literature are discussed. Some of the previous results are generalized and some are simplified. DFT calculation results of a fullerene molecule, a finite single-walled carbon nanotube and a finite boron nitride nanotube are presented, indicating that this Koopmans' theorem approximately holds, even if the orbital relaxation is taken into consideration.     

Thermodynamic and kinetic theory of ionic micellar systems: 1. Work of aggregation

The previously developed approach to the thermodynamics of molecular aggregates based on the methods of the nucleation theory is extended to including ionic aggregative systems with the choice of the real concentration of monomeric ions of surfactants as a standard concentration of aggregates. Results are also generalized to multicomponent systems and a new expression for the work of aggregation is derived. Three models of ionic aggregates are analyzed, i.e., bare aggregate (containing no counterions), fully dressed aggregate (with maximal number of counterions), and dressed aggregate with arbitrary numbers of counterions. In the consideration of the ionic aggregate as a complex ion that yields simpler ions upon dissociation, the notion of the average activity coefficient of the ionic aggregate is introduced. Characteristics of real dressed ionic aggregates are studied. Restrictions in the use of the model of the dressed ionic aggregate are revealed.     

Experimental evidence of the role of viscosity in the molecular kinetic theory of dynamic wetting

We report an experimental study of the dynamics of spontaneous spreading of aqueous glycerol drops on glass. For a range of glycerol concentrations, we follow the evolution of the radius and contact angle over several decades of time and investigate the influence of solution viscosity. The application of the molecular kinetic theory to the resulting data allows us to extract the coefficient of contact-line friction ζ, the molecular jump frequency κ(0), and the jump length λ for each solution. Our results show that the modified theory, which explicitly accounts for the effect of viscosity, can successfully be applied to droplet spreading. The viscosity affects the jump frequency but not the jump length. In combining these data, we confirm that the contact-line friction of the solution/air interface against the glass is proportional to the viscosity and exponentially dependent on the work of adhesion.     

Classical theory of collisional entropy production

A classical dynamical theory of elementary collision processes is formulated in analogy to the quantum theory of the dynamical scattering matrix, which can be defined for a pure quantum stationary scattering state. The elements of this matrix are probability amplitudes for transitions between internal states defined for given values of a reaction coordinate. The squared magnitudes of these amplitudes, modeled in the proposed classical theory, define normalized internal state population distributions suitable for information theoretical analysis. Statistical entropy and surprisal are defined as dynamical functions of a reaction coordinate. This formalism differs fundamentally from concepts based on the classical Liouville equation.     

Classical description of activated conformational processes in molecular systems coupled to solvent degrees of freedom

Extended Fokker-Planck (FP) equations are generalized stochastic equations which describe the evolution of a set of coordinates, loosely referred to as solute, coupled with a relevant set of solvent variables. They are useful for the analysis of molecular dynamics in liquids, when a time-scale separation between probe motions and relaxation times of interactions with surroundings particles cannot be performed, because of persistence of slowly fluctuating components. In this article, we focus attention on a model system, made up of an angular coordinate and its conjugate momentum, submitted to a bistable potential, and coupled to a dissipative harmonic mode, to investigate the influence of polar solvents on reactive dynamics. The results are appropriate to describe dielectric effects, solvent-controlled conformational changes involving charge transfer which occur in photophysical processes, and the dynamic Stokes shifts observed in time-resolved fluorescence experiments. © 1996 John Wiley & Sons, Inc.     

Differential mobility spectrometry of chlorocarbons with a micro-fabricated drift tube

Chlorocarbons were ionized through gas phase chemistry at ambient pressure in air and resultant ions were characterized using a micro-fabricated drift tube with differential mobility spectrometry (DMS). Positive and negative product ions were characterized simultaneously in a single drift tube equipped with a 3 mCi (63)Ni ion source at 50 degrees C and drift gas of air with 1 ppm moisture. Scans of compensation voltage for most chlorocarbons produced differential mobility spectra with Cl(-) as the sole product ion and a few chlorocarbons produced adduct ions, M (.-) Cl(-). Detection limits were approximately 20-80 pg for gas chromatography-DMS measurements. Chlorocarbons also yielded positive ions through chemical ionization in air and differential mobility spectra showed peaks with characteristic compensation voltages for each substance. Field dependence of mobility was determined for positive and negative ions of each substance and confirmed characteristic behavior for each ion. A DMS analyzer with a membrane inlet was used to continuously monitor effluent from columns of bentonite or synthetic silica beads to determine breakthrough volumes of individual chlorocarbons. These findings suggest a potential of DMS for monitoring subsurface environments either on site or perhaps in situ.     

Extraction of kinetic information from single-molecule experiments.

We discuss two possible approaches for extracting kinetic information from single-molecule experiments. The first approach is based on computing correlation functions from measured fluorescence signals, and the second on studying the statistics of on and off times of the same fluorescence signal. We show that in both cases it is possible to extract kinetic information about the nature of intramolecular fluctuations of the single molecule. We show that for single-molecule kinetics the intramolecular fluctuations produce stochastic memory effects which lead to new dynamic features that do not exist in traditional chemical kinetics. In particular, we investigate a new type of chemical oscillations in correlation functions observed experimentally by Edman and Rigler (Proc. Natl. Acad. Sci. USA 2000, 97, 8266).     

Orientation and alignment in charge exchange processes involving sodium atoms studied by semi-classical molecular theory

The semi-classical molecular model is applied to study the charge exchange processes in the H⁺-Na 3p and Li⁺-Na 3p systems in the keV energy range. The dependence of the charge exchange on the orientation and the alignment of the initial or final state is obtained for the transition probabilities and for the differential cross sections. The 14 state present model for the H⁺-Na system is in good agreement with the experimental differential cross sections. The alignment is explained by orientation occurring in the transfer region. The cross sections predicted with a 28 state model for Li⁺-Na exhibit similar behaviour as for H⁺-Na.     

Variational transition state theory as a tool to determine kinetic selectivity in reactions involving a valley-ridge inflection point

Variational transition state theory has been used to calculate the kinetic isotope effects affecting product ratios in the reaction between (1)O(2) and d(6)-tetramethylethylene. The minimum energy path on the potential energy surface for this process reaches a valley-ridge inflection point and then bifurcates leading to the two final products. Using canonical variational transition state theory, two distinct dynamical bottlenecks were located corresponding to the H- and the D-abstraction, respectively. The calculated KIE at 263 K turns out to be 1.126. Analogously, a H/T KIE of 1.17 at the same temperature has been found for the reaction of (1)O(2) with the tritiated derivative of tetramethylethylene.     

Shock tube pyrolysis of pyrrole and kinetic modeling

The kinetics of pyrolysis of pyrrole dilute in argon have been studied in a single pulse shock tube, using capillary column GC, together with GC/MS and FTIR for product identification, over the temperature range 1200–1700 K, total pressures of 7.5–13.5 atm and nominal mixture compositions of pyrrole of 5000 and 700 ppm (nominal concentrations of 5 × 10^?7 and 7 × 10^?8 mol cm^?3). Time-resolved measurements of the rate of disappearance of pyrrole behind reflected shock waves have been made by absorption spectroscopy at 230 nm, corresponding to the lowest ¹π* ← ¹π transition of pyrrole at pressures of 20 atm and mixture compositions between 1000–2000 ppm pyrrole (1.7–3.0 × 10^?7 mol cm^?3) over the temperature range of 1300 to 1700 K. At the lower end of the studied temperature range, the isomers of pyrrole, allyl cyanide and cis- and trans-crotononitrile, were the principal products, together with hydrogen cyanide and propyne/allene. At elevated temperatures, acetylene, acetonitrile, cyanoacetylene, and hydrogen became important products. The rate of overall disappearance of pyrrole, as measured by absorption spectrometry, was found to be first order in pyrrole concentration, with a rate constant k_dis(pyrrole) = 10^14.1±0.7 exp(?74.1 ± 3.0 kcal mol^?1/RT) s^?1 between 1350–1600 K and at a pressure of 20 atm. First order dependence of pyrrole decomposition and major product formation was also observed in the single pulse experiments over the range of mixture compositions studied. A 75-step reaction model is presented and shown to substantially fit the observed temperature profiles of the major product species and the reactant profile. In the model the initiation reaction is postulated to be the reversible formation of pyrrolenine, (2H-pyrrole). Pyrrolenine can undergo ring scission at the C2? N bond forming a biradical which can rearrange to form allyl cyanide and crotononitrile or undergo decomposition to form HCN and C₃H₄ or acetylene and a precursor of acetonitrile. The model predicts an overall rate of disappearance of pyrrole in agreement with the experimental measurements.     

Extended kinetic theory of polymeric fluids

The setting of classical configuration space kinetic theory of polymeric fluids is extended by adopting the field of internal momentum into the set of independent state variables. The internal momentum predicated by the extended theory compares well with the internal momentum seen in molecular simulations. The extended theory is also shown to provide an intrinscially consistent formulation of the time evolution of macromolecules with internal friction.