Quantum Mechanical and Molecular Mechanics (MM3) Studies of Hydrazines

Anharmonic Molecular Mechanics: Ab Initio Based Morse Parametrizations for the Popular MM3 Force Field

Methodologies for creating reactive potential energy surfaces from molecular mechanics force-fields are becoming increasingly popular. To date, molecular mechanics force-fields in biochemistry and small molecule organic chemistry tend to use harmonic expressions to treat bonding stretches, which is a poor approximation in reactive and nonequilibirum molecular dynamics simulations since bonds are often displaced significantly from their equilibrium positions. For such applications there is need for a better treatment of anharmonicity. In this contribution, Morse bonding potentials have been extensively parametrized for the atom types in the MM3 force field of Allinger and co-workers using high level CCSD(T)(F12*) energies. To our knowledge this is among the first instances of a comprehensive parametrization of Morse potentials in a popular organic chemistry force field. In the context of molecular dynamics simulations, these data will: (1) facilitate the fitting of reactive potential energy surfaces using empirical valence bond approaches and (2) enable more accurate treatments of energy transfer.

Gas-phase kinetics study of reaction of OH radical with CH3NHNH2 by second-order multireference perturbation theory

The gas-phase kinetics of H-abstraction reactions of monomethylhydrazine (MMH) by OH radical was investigated by second-order multireference perturbation theory and two-transition-state kinetic model. It was found that the abstractions of the central and terminal amine H atoms by the OH radical proceed through the formation of two hydrogen bonded preactivated complexes with energies of 6.16 and 5.90 kcal mol(-1) lower than that of the reactants, whereas the abstraction of methyl H atom is direct. Due to the multireference characters of the transition states, the geometries and ro-vibrational frequencies of the reactant, transition states, reactant complexes, and product complexes were optimized by the multireference CASPT2/aug-cc-pVTZ method, and the energies of the stationary points of the potential energy surface were refined at the QCISD(T)/CBS level via extrapolation of the QCISD(T)/cc-pVTZ and QCISD(T)/cc-pVQZ energies. It was found that the abstraction reactions of the central and two terminal amine H atoms of MMH have the submerged energy barriers with energies of 2.95, 2.12, and 1.24 kcal mol(-1) lower than that that of the reactants respectively, and the abstraction of methyl H atom has a real energy barrier of 3.09 kcal mol(-1). Furthermore, four MMH radical-H(2)O complexes were found to connect with product channels and the corresponding transition states. Consequently, the rate coefficients of MMH + OH for the H-abstraction of the amine H atoms were determined on the basis of a two-transition-state model, with the total energy E and angular momentum J conserved between the two transition-state regions. In units of cm(3) molecule(-1) s(-1), the rate coefficient was found to be k(1) = 3.37 × 10(-16)T(1.295) exp(1126.17/T) for the abstraction of the central amine H to form the CH(3)N(•)NH(2) radical, k(2) = 2.34 × 10(-17)T(1.907) exp(1052.26/T) for the abstraction of the terminal amine H to form the trans-CH(3)NHN(•)H radical, k(3) = 7.41 × 10(-20)T(2.428) exp(1343.20/T) for the abstraction of the terminal amine H to form the cis-CH(3)NHN(•)H radical, and k(4) = 9.13 × 10(-21)T(2.964) exp(-114.09/T) for the abstraction of the methyl H atom to form the C(•)H(2)NHNH(2) radical, respectively. Assuming that the rate coefficients are additive, the total rate coefficient of these theoretical predictions quantitatively agrees with the measured rate constant at temperatures of 200-650 K, with no adjustable parameters.

Theoretical studies of the in-solution isomeric protonation of non-aromatic six-member rings with two nitrogens

For exploring the preferred site for hydrogen bond formation, theoretical calculations have been performed for a number of six-member, nonaromatic rings allowing for alternative protonation on the ring nitrogens. Gas-phase protonation studies for test molecules indicate that the B3LYP/aug-cc-pvtz and QCISD(T)(CBS) calculations approach the experimental values within about 1 kcal/mol with considerable improvement for relative enthalpies and free energies. Relative free energies calculated at the IEF-PCM/B3LYP/aug-cc-pvtz level predict favorable protonation on the tertiary rather than on the secondary nitrogen both in aqueous solution and in a dichloromethane solvent for saturated rings. Protonation on a nitrogen atom next to a C═C bond is disfavored due to a large increase in internal energy. Monte Carlo simulations considering a counterion and Ewald summation for the long-range electrostatic effects for a 0.1 molar model system predict ΔG(solv)/MC values generally less negative than from the IEF-PCM calculations. These results make the protonation on the tertiary nitrogen even more favored. The solute-solvent pair-energy distribution depends sensitively on the applied model. In conclusion, the freely moving anion has been considered as the most relevant model with overall neutrality for the system and applying the least restrictions.

Theoretical and REMPI spectroscopic study on phenylhydrazine and phenylhydrazine-(Ar)n (n=1, 2) van der Waals complexes

Phenylhydrazine and its van der Waals complexes with one or two argon atoms were investigated with theoretical calculations and resonant two photon ionization (R2PI) spectroscopy. The ab initio and DFT calculations found a conversion of the orbital hybridization of the Nbeta atom from sp3-like in the S0 state to sp2-like in the S1 state, suggesting that the lone pair electrons of the Nbeta atom are involved in a super p-p-pi conjugation over the skeleton of phenylhydrazine in the S1 state. The structural change of the hydrazino group in the S1<--S0 electronic transition was reflected by the vibrational excitations of the hydrazino group observed in the 1C-R2PI spectrum. The band origin of the S1<--S0 transition is determined to be 33610 cm(-1) and the adiabatic ionization energy (IE) of phenylhydrazine, measured by 2C-R2PI spectroscopy, is 62829+/-15 cm(-1). The S1<--S0 electronic transitions of phenylhydrazine-Ar and phenylhydrazine-Ar2 complexes were also observed in the 1C-R2PI spectrum, and their band origins are, respectively, red-shifted by 39 and 80 cm(-1) from that of phenylhydrazine.

Computational Methods in Organic Thermochemistry. 2. Enthalpies and Free Energies of Formation for Functional Derivatives of Organic Hydrocarbons

Standard state enthalpies and free energies of formation for nitrogen-, oxygen-, sulfur-, fluorine-, chlorine-, and silicon-containing compounds can be computed with reasonable accuracy (usually within four and often two kJ/mol) using the G3 and G3MP2 model chemistries. In several of the families, compounds with as many as 10 carbon atoms have been computed. Larger errors are found in the free energies of these longer chain molecules which can be reduced by compensating for the presence of multiple conformers having a significant population at 298.15 K. In some instances, those substances showing large deviations are found to have experimental energies that may be erroneous.

Thermochemical and Kinetic Analysis of the Thermal Decomposition of Monomethylhydrazine: An Elementary Reaction Mechanism

The reaction kinetics for the thermal decomposition of monomethylhydrazine (MMH) was studied with quantum Rice-Ramsperger-Kassel (QRRK) theory and a master equation analysis for pressure falloff. Thermochemical properties were determined by ab initio and density functional calculations. The entropies, S degrees (298.15 K), and heat capacities, Cp degrees (T) (0 < or = T/K < or = 1500), from vibrational, translational, and external rotational contributions were calculated using statistical mechanics based on the vibrational frequencies and structures obtained from the density functional study. Potential barriers for internal rotations were calculated at the B3LYP/6-311G(d,p) level, and hindered rotational contributions to S degrees (298.15 K) and Cp degrees (T) were calculated by solving the Schrödinger equation with free rotor wave functions, and the partition coefficients were treated by direct integration over energy levels of the internal rotation potentials. Enthalpies of formation, DeltafH degrees (298.15 K), for the parent MMH (CH3NHNH2) and its corresponding radicals CH3N*NH2, CH3NHN*H, and C*H2NHNH2 were determined to be 21.6, 48.5, 51.1, and 62.8 kcal mol(-1) by use of isodesmic reaction analysis and various ab initio methods. The kinetic analysis of the thermal decomposition, abstraction, and substitution reactions of MMH was performed at the CBS-QB3 level, with those of N-N and C-N bond scissions determined by high level CCSD(T)/6-311++G(3df,2p)//MPWB1K/6-31+G(d,p) calculations. Rate constants of thermally activated MMH to dissociation products were calculated as functions of pressure and temperature. An elementary reaction mechanism based on the calculated rate constants, thermochemical properties, and literature data was developed to model the experimental data on the overall MMH thermal decomposition rate. The reactions of N-N and C-N bond scission were found to be the major reaction paths for the modeling of MMH homogeneous decomposition at atmospheric conditions.

Origin of Rotational Barriers of the N−N Bond in Hydrazine: NBO Analysis

Hydrazine passes through two transition states, TS1 (phi = 0 degrees ) and TS2 (phi = 180 degrees ), in the course of internal rotation around its N-N bond. The origin of the corresponding rotational barriers in hydrazine has been extensively studied by experimental and theoretical methods. Here, we used natural bond orbital (NBO) analysis and energy decomposition of rotational barrier energy (DeltaE(barrier)) to understand the origin of the torsional potential energy profile of this molecule. DeltaE(barrier) was dissected into structural (DeltaE(struc)), steric exchange (DeltaE(steric)), and hyperconjugative (DeltaE(deloc)) energy contributions. In both transition states, the major barrier-forming contribution is DeltaE(deloc). The TS2 barrier is lowered by pyramidalization of nitrogen atoms through lowering DeltaE(struc), not by N-N bond lengthening through lowering DeltaE(steric). Higher pyramidality of nitrogen atoms of TS2 than that of TS1 explains well why the N-N bond of TS2 is longer than that of TS1. Finally, the steric repulsion between nitrogen lone pairs does not determine the rotational barrier; nuclear-nuclear Coulombic repulsion between outer H/H atoms in TS1 plays an important role in increasing DeltaE(struc). Taken together, we explain the reason for the different TS1 and TS2 barriers. We show that NBO analysis is a useful tool for understanding structures and potential energy surfaces of compounds containing the N-N bond.