Direct Dynamics Simulation of the Lifetime of Trimethylene

Post-spin crossing dynamics determine the regioselectivity in open-shell singlet biradical recombination

Comprehensive computational studies reveal unique dynamic effects in a multi-spin-state reaction that determine the regioselectivity of a biradical recombination process.

Gas-phase reactivity of [Ca(formamide)]2+ complex: an example of different dynamical behaviours

In the present contribution, we have summarized our recent work on the comprehension of [Ca(formamide)]²⁺ complex gas-phase unimolecular dissociation. By using different theoretical approaches, we were able to revise the original (and typical for such kind of problems) understanding given in terms of stationary points on the potential energy surface, which did not provide a satisfactory explanation of the experimentally observed reactivity. In particular, we point out how non-statistical and non-intrinsic reaction coordinate mechanisms are of fundamental importance.This article is part of the themed issue 'Theoretical and computational studies of non-equilibrium and non-statistical dynamics in the gas phase, in the condensed phase and at interfaces'.

Quantum chemical investigation of the thermal denitrogenation of 1-pyrazoline

CASPT2//CASSCF(12,12) calculations reveal that a synchronous path is favourable for the denitrogenation of 1-pyrazoline.

Photodissociation of acryloyl chloride at 193 nm: interpretation of the product energy distributions, and new elimination pathways

The use of an automated TS search method leads to the finding of novel HCl elimination pathways.

HCN elimination from vinyl cyanide: product energy partitioning, the role of hydrogen–deuterium exchange reactions and a new pathway

A schematic diagram of HCN elimination channels from vinyl cyanide including a new CCdiss pathway.

Internal energy of HCl upon photolysis of 2-chloropropene at 193 nm investigated with time-resolved Fourier-transform spectroscopy and quasiclassical trajectories

Following photodissociation of 2-chloropropene (H(2)CCClCH(3)) at 193 nm, vibration-rotationally resolved emission spectra of HCl (upsilon < or = 6) in the spectral region of 1900-2900 cm(-1) were recorded with a step-scan time-resolved Fourier-transform spectrometer. All vibrational levels show a small low-J component corresponding to approximately 400 K and a major high-J component corresponding to 7100-18,700 K with average rotational energy of 39+/-(3)(11) kJ mol(-1). The vibrational population of HCl is inverted at upsilon = 2, and the average vibrational energy is 86+/-5 kJ mol(-1). Two possible channels of molecular elimination producing HCl + propyne or HCl + allene cannot be distinguished positively based on the observed internal energy distribution of HCl. The observed rotational distributions fit qualitatively with the distributions of both channels obtained with quasiclassical trajectories (QCTs), but the QCT calculations predict negligible populations for states at small J. The observed vibrational distribution agrees satisfactorily with the total QCT distribution obtained as a weighted sum of contributions from both four-center elimination channels. Internal energy distributions of HCl from 2-chloropropene and vinyl chloride are compared.

Self-consistent polarization neglect of diatomic differential overlap: application to water clusters

Semiempirical self-consistent field (SCF) methods based on the neglect of diatomic differential overlap (NDDO) formalism have the ability to treat the formation and breaking of chemical bonds but have been found to poorly describe hydrogen bonding and weak electrostatic complexes. In contrast, most empirical potentials are not able to describe bond breaking and formation but have the ability to add missing elements of hydrogen bonding by using classical electrostatic interactions. We present a new method which combines aspects of both NDDO-based SCF techniques and classical descriptions of polarization to describe the diffuse nature of the electronic wavefunction in a self-consistent manner. We develop the "self-consistent polarization neglect of diatomic differential overlap" (SCP-NDDO) theory with the additional description of molecular dispersion developed as a second-order perturbation theory expression. The current study seeks to model water-water interactions as a test case. To this end, we have parametrized the method to accurate ab initio complete basis set limit estimates of small water cluster binding energies of Xantheas and co-workers [J. Chem. Phys. 116, 1493 (2002); 120, 823 (2004)]. Overall agreement with the ab initio binding energies (n=2-6, and 8) is achieved with a rms error of 0.19 kcal/mol. We achieve noticeable improvements in the structure, vibrational frequencies, and energetic predictions of water clusters (n< or =21) relative to standard NDDO-based methods.

Comparison of levels of electronic structure theory in direct dynamics simulations of C2H5F --> HF + C2H4 product energy partitioning

Direct dynamics simulations at the MP2/6-311++G** level of theory were performed to study C(2)H(5)F --> HF + C(2)H(4) product energy partitioning. The simulation results are compared with experiment and a previous MP2/6-31G* simulation. The current simulation with the larger basis set releases more energy to HF vibration and less to HF + C(2)H(4) relative translation as compared to the previous simulation with the 6-31G* basis set. The HF rotation and vibration energy distributions determined from the current simulation are in overall very good agreement with previous experimental studies of C(2)H(5)F dissociation by chemical activation and IRMPA. A comparison of the simulations with experiments suggests there may be important mass effects for energy partitioning in HX elimination from haloalkanes. The transition state (TS) structures and energies calculated with MP2 and the 6-31G* and 6-311++G** basis sets are compared with those calculated using CCD, CCSD, CCSD(T), and the 6-311++G** basis set.

Use of a single trajectory to study product energy partitioning in unimolecular dissociation: mass effects for halogenated alkanes

A single trajectory (ST) direct dynamics approach is compared with quasiclassical trajectory (QCT) direct dynamics calculations for determining product energy partitioning in unimolecular dissociation. Three comparisons are made by simulating C(2)H(5)F-->HF + C(2)H(4) product energy partitioning for the MP26-31G(*) and MP26-311 + + G(**) potential energy surfaces (PESs) and using the MP26-31G(*) PES for C(2)H(5)F dissociation as a model to simulate CHCl(2)CCl(3)-->HCl + C(2)Cl(4) dissociation and its product energy partitioning. The trajectories are initiated at the transition state with fixed energy in reaction-coordinate translation E(t) (double dagger). The QCT simulations have zero-point energy (ZPE) in the vibrational modes orthogonal to the reaction coordinate, while there is no ZPE for the STs. A semiquantitative agreement is obtained between the ST and QCT average percent product energy partitionings. The ST approach is used to study mass effects for product energy partitioning in HX(X = F or Cl) elimination from halogenated alkanes by using the MP26-31G(*) PES for C(2)H(5)F dissociation and varying the masses of the C, H, and F atoms. There is, at most, only a small mass effect for partitioning of energy to HX vibration and rotation. In contrast, there are substantial mass effects for partitioning to relative translation and the polyatomic product's vibration and rotation. If the center of mass of the polyatomic product is located away from the C atom from which HX recoils, the polyatomic has substantial rotation energy. Polyatomic products, with heavy atoms such as Cl atoms replacing the H atoms, receive substantial vibration energy that is primarily transferred to the wag-bend motions. For E(t) (double dagger) of 1.0 kcalmol, the ST calculations give average percent partitionings to relative translation, polyatomic vibration, polyatomic rotation, HX vibration, and HX rotation of 74.9%, 6.8%, 1.5%, 14.4%, and 2.4% for C(2)H(5)F dissociation and 39.7%, 38.1%, 0.2%, 16.1%, and 5.9% for a model of CHCl(2)CCl(3) dissociation.

Interpolating moving least-squares methods for fitting potential energy surfaces: Computing high-density potential energy surface data from low-density ab initio data points

A highly accurate and efficient method for molecular global potential energy surface (PES) construction and fitting is demonstrated. An interpolating-moving-least-squares (IMLS)-based method is developed using low-density ab initio Hessian values to compute high-density PES parameters suitable for accurate and efficient PES representation. The method is automated and flexible so that a PES can be optimally generated for classical trajectories, spectroscopy, or other applications. Two important bottlenecks for fitting PESs are addressed. First, high accuracy is obtained using a minimal density of ab initio points, thus overcoming the bottleneck of ab initio point generation faced in applications of modified-Shepard-based methods. Second, high efficiency is also possible (suitable when a huge number of potential energy and gradient evaluations are required during a trajectory calculation). This overcomes the bottleneck in high-order IMLS-based methods, i.e., the high cost/accuracy ratio for potential energy evaluations. The result is a set of hybrid IMLS methods in which high-order IMLS is used with low-density ab initio Hessian data to compute a dense grid of points at which the energy, Hessian, or even high-order IMLS fitting parameters are stored. A series of hybrid methods is then possible as these data can be used for neural network fitting, modified-Shepard interpolation, or approximate IMLS. Results that are indicative of the accuracy, efficiency, and scalability are presented for one-dimensional model potentials as well as for three-dimensional (HCN) and six-dimensional (HOOH) molecular PESs.

Direct dynamics simulations using Hessian-based predictor-corrector integration algorithms

In previous research [J. Chem. Phys. 111, 3800 (1999)] a Hessian-based integration algorithm was derived for performing direct dynamics simulations. In the work presented here, improvements to this algorithm are described. The algorithm has a predictor step based on a local second-order Taylor expansion of the potential in Cartesian coordinates, within a trust radius, and a fifth-order correction to this predicted trajectory. The current algorithm determines the predicted trajectory in Cartesian coordinates, instead of the instantaneous normal mode coordinates used previously, to ensure angular momentum conservation. For the previous algorithm the corrected step was evaluated in rotated Cartesian coordinates. Since the local potential expanded in Cartesian coordinates is not invariant to rotation, the constants of motion are not necessarily conserved during the corrector step. An approximate correction to this shortcoming was made by projecting translation and rotation out of the rotated coordinates. For the current algorithm unrotated Cartesian coordinates are used for the corrected step to assure the constants of motion are conserved. An algorithm is proposed for updating the trust radius to enhance the accuracy and efficiency of the numerical integration. This modified Hessian-based integration algorithm, with its new components, has been implemented into the VENUS/NWChem software package and compared with the velocity-Verlet algorithm for the H(2)CO-->H(2)+CO, O(3)+C(3)H(6), and F(-)+CH(3)OOH chemical reactions.

Detailed Dynamics of the Photodissociation of Cyclobutane

Semiclassical electron-radiation-ion dynamics simulations are reported for the photodissociation of cyclobutane into two molecules of ethylene. The results clearly show the formation of the tetramethylene intermediate diradical, with dissociation completed in approximately 400 fs. In addition, the potential energy surfaces of the electronic ground state and lowest excited-state were calculated at the complete-active-space self-consistent-field/multireference second-order perturbation theory (CASSCF/MRPT2) level with 6-31G* basis sets, along the reaction path determined by the dynamics simulations. There are well-defined energy minima and maxima in the intermediate state region. It is found that both C-C-C bond bending and rotation of the molecule (around the central C-C bond) have important roles in determining the features of the potential energy surfaces for the intermediate species. Finally, the simulations and potential energy surface calculations are applied together in a discussion of the full mechanism for cyclobutane photodissociation.

Rovibrational distributions of HF in the photodissociation of vinyl fluoride at 193 nm: a direct MP2 quasiclassical trajectory study

Quasiclassical trajectory calculations were performed to calculate rovibrational distributions of the nascent HF fragment in the photodissociation of vinyl fluoride at 193 nm. The trajectories were initiated at the transition states of the four-center (4C) and three-center (3C) HF elimination channels, using a microcanonical, quasiclassical normal-mode sampling. In general, the calculated distributions are in reasonably good agreement with experiment. In particular, the trajectory distributions show bimodal character, although not as pronounced as that observed experimentally. The calculations predict that the 3C and 4C distributions are rather similar to each other, which suggests that the low-J and high-J components of the rotational distributions cannot be specifically assigned to each of these channels.

Quasiclassical trajectory calculations on the photodissociation of CF2CHCl at 193 nm: product energy distributions for the HF and HCl eliminations

Quasiclassical trajectory calculations were carried out to determine product energy distributions for the HCl and HF eliminations that take place in the photodissociation of 2-chloro-1,1-difluoroethylene at 193 nm. The trajectories were initiated at the transition states of the HCl and HF elimination channels under microcanonical, quasiclassical conditions, and were propagated with the energies and gradients taken directly from density functional theory calculations. Good agreement with experiment is found, except for the translational energy distribution of the HF elimination channel and the average vibrational energy of the HCl fragment. Possible sources of disagreement are discussed.

Thermal Isomerizations of 2-d-1-(E)-Propenylcyclobutanes to 4-d-3-Methylcyclohexenes

Racemic trans-2-d-1-(E)-propenylcyclobutane at 276 degrees C in the gas phase fragments to give ethylenes and pentadienes, equilibrates with its cis isomer, and rearranges to mixtures of 4-d- and 6-d-3-methylcyclohexenes through [1,3] carbon shifts. The time-dependent distributions of deuterium-labeled isomers of propenylcyclobutanes and 3-methylcyclohexenes reveal a significant secondary deuterium kinetic isotope effect favoring C1-C4 over C1-C2 bond breaking (kH/kD = 1.16 +/- 0.02) and a 72:28 preference for structural isomerizations giving (si + ar) rather than (sr + ai) products through conformationally flexible short-lived diradical intermediates.

A Nested Molecule-Independent Neural Network Approach for High-Quality Potential Fits†

It is shown that neural networks (NNs) are efficient and effective tools for fitting potential energy surfaces. For H2O, a simple NN approach works very well. To fit surfaces for HOOH and H2CO, we develop a nested neural network technique in which we first fit an approximate NN potential and then use another NN to fit the difference of the true potential and the approximate potential. The root-mean-square error (RMSE) of the H2O surface is 1 cm(-1). For the 6-D HOOH and H2CO surfaces, the nested approach does almost as well attaining a RMSE of 2 cm(-1). The quality of the NN surfaces is verified by calculating vibrational spectra. For all three molecules, most of the low-lying levels are within 1 cm(-1) of the exact results. On the basis of these results, we propose that the nested NN approach be considered a method of choice for both simple potentials, for which it is relatively easy to guess a good fitting function, and complicated (e.g., double well) potentials for which it is much harder to deduce an appropriate fitting function. The number of fitting parameters is only moderately larger for the 6-D than for the 3-D potentials, and for all three molecules, decreasing the desired RMSE increases only slightly the number of required fitting parameters (nodes). NN methods, and in particular the nested approach we propose, should be good universal potential fitting tools.

Insight into Selected Reactions in Low-Temperature Dimethyl Ether Combustion from Born−Oppenheimer Molecular Dynamics

Dimethyl ether is under consideration as an alternative diesel fuel. Its combustion chemistry is as yet ill-characterized. Here we use Born-Oppenheimer molecular dynamics (BOMD) based on DFT-B3LYP forces to investigate the short-time dynamics of selected features of the low-temperature dimethyl ether (DME) oxidation potential energy surface. Along the chain propagation pathway, we run BOMD simulations from the transition state involving the decomposition of (*)CH(2)OCH(2)OOH to two CH(2)=O and an (*)OH radical. We predict that formaldehyde C-O stretch overtones are excited, consistent with laser photolysis experiments. We also predict that O-H overtones are excited for the (*)OH formed from (*)CH(2)OCH(2)OOH dissociation. We also investigate short-time dynamics involved in chain branching. First, we examine the isomerization transition state of (*)OOCH(2)OCH(2)OOH --> HOOCH(2)OCHOOH. The latter species is predicted to be a short-lived metastable radical that decomposes within 500 fs to hydroperoxymethyl formate (HPMF; HOOCH(2)OC(=O)H) and the first (*)OH of chain branching. The dissociation of HOOCH(2)OCHOOH exhibits non-RRKM behavior in its lifetime profile, which may be due to conformational constraints or slow intramolecular vibrational energy transfer (IVR) from the nascent H-O bond to the opposite end of the radical, where O-O scission occurs to form HPMF and (*)OH. In a few trajectories, we see HOOCH(2)OCHOOH recross back to (*)OOCH(2)OCH(2)OOH because the isomerization is endothermic, with only an 8 kcal/mol barrier to recrossing. Therefore, some inhibition of chain-branching may be due to recrossing. Second, trajectories run from the transition state leading to the direct decomposition of HPMF (an important source of the second (*)OH radical in chain branching) to HCO, (*)OH, and HC(=O)OH show that these products can recombine to form many other possible products. These products include CH(2)OO + HC(=O)OH, H(2)O + CO + HC(=O)OH, HC(=O)OH + HC(=O)OH, and HC(=O)C(=O)H + H(2)O, which (save CH(2)OO + HC(=O)OH) are all more thermodynamically stable than the original HCO + (*)OH + HC(=O)OH products. Moreover, the multitude of extra products suggest that standard statistical rate theories cannot completely describe the reaction kinetics of significantly oxygenated compounds such as HPMF. These secondary products consume the second (*)OH required for explosive combustion, suggesting an inhibition of DME fuel combustion is likely.

Molecular Dynamics of Complex Gas-Phase Reactive Systems by Time-Dependent Groups

A novel way of assembling the total potential for performing molecular dynamical studies of complex gas-phase reactive chemical systems is introduced. The method breaks the calculation of the total potential and gradients of the potential into time-dependent groups that are governed by spatial cutoffs. These groups evolve during the course of the simulation and their number may increase or diminish as the dynamics of the system determine. In an effort to extend the simulation time of these complex reactive processes and to use high levels of theory when necessary, multiple levels of theory may be used over the groups for the calculation of both the intragroup and intergroup interactions. Representative simulations are performed to illustrate the method and a computationally facile method of obtaining the groups of a simulation are also discussed.

Factors affecting the selection of products from a photochemically generated singlet biradical

The chemistries of a monoradical of the ultrafast "radical-clock" type and a structurally related singlet biradical, generated by Norrish type II photochemistry, are compared. The monoradical is found to undergo the characteristic ring-opening reaction of its class at about 10(10) s(-1) at room temperature. However, the singlet biradical shows no evidence of the analogous ring-opening reaction. The contrasting chemistry is traced not to a fundamental difference in electronic structure of the two intermediates, but rather to a steric interaction that the biradical alone would have to suffer during the ring opening. Although the magnitude of the steric hindrance is small (estimated 15-20 kJ mol(-1)), it is enough to shut down the reaction, because the biradical has other facile product-forming reactions available.

Ab initiodirect dynamics trajectory simulation of C2H5F→C2H4+HF product energy partitioning

Direct dynamics classical trajectory simulations were performed to study product energy partitioning in C(2)H(5)F-->C(2)H(4)+HF dissociation. The intrinsic reaction coordinate potential energy curve, reaction energetics, and transition state (TS) properties were calculated for this reaction at different levels of electronic structure theory, and MP2/6-31G( *) was chosen as a meaningful and practical method for performing the direct dynamics. The trajectories show that the HF bond, uncoupled from the other degrees of freedom, is formed within the first 10 fs as the system moves from the TS towards products. The populations of the HF vibration states, determined from the simulations, decrease monotonically as found from experiments. However, the simulation's populations for the low and high energy vibration states are larger and smaller, respectively, than the experimental results. The HF rotational temperature found from the simulations is in agreement with experiment. Increasing the TS's excess energy gives higher rotational temperatures for both C(2)H(4) and HF. Energy is partitioned to the products from both the excess energy in the TS and the potential energy release in the exit channel. Partitioning from these two energy sources is distinguished by varying the TS's excess energy. An analysis of the simulation's energy disposal shows that the fractions of the excess energy partitioned to relative translation, C(2)H(4) vibration, C(2)H(4) rotation, HF vibration, and HF rotation, are 0.17, 0.64, 0.076, 0.067, and 0.046, respectively, and are in good agreement with previous simulations on empirical potentials and experiments. The partitioning found for the potential energy release is 81%, <0.05%, 5%, 11%, and 3% to relative translation, C(2)H(4) vibration, C(2)H(4) rotation, HF vibration, and HF rotation. This result is substantially different than the deduction from experiments, which summarizes the partitioning as 20%, 45%, 24%, and <12% to relative translation, C(2)H(4) vibration+rotation, HF vibration, and HF rotation. Possible origins of the difference between the simulations and experiments in the release of the potential energy is discussed.

Quasiclassical trajectory study of the collision-induced dissociation of CH[sub 3]SH[sup +]+Ar

Quasiclassical trajectory calculations were carried out to study the dynamics of energy transfer and collision-induced dissociation (CID) of CH(3)SH(+) + Ar at collision energies ranging from 4.34 to 34.7 eV. The relative abundances calculated for the most relevant product ions are found to be in good agreement with experiment, except for the lowest energies investigated. In general, the dissociation to form CH(3)(+) + SH is the dominant channel, even though it is not among the energetically favored reaction pathways. The results corroborate that this selective dissociation observed upon collisional activation arises from a more efficient translational to vibrational energy transfer for the low-frequency C-S stretching mode than for the high-frequency C-H stretching modes, together with weak couplings between the low- and high-frequency modes of vibration. The calculations suggest that CID takes place preferentially by a direct CH(3)(+) + SH detachment, and more efficiently when the Ar atom collides with the methyl group-side of CH(3)SH(+).