A QM ϩ MM direct chemical dynamics simulation was performed to study collisions of protonated octaglycine, gly 8 -H ϩ , with the diamond {111} surface at an initial collision energy E i of 100 eV and incident angle i of 0°and 45°. The semiempirical model AM1 was used for the gly 8 -H ϩ intramolecular potential, so that its fragmentation could be studied. Shattering dominates gly 8 -H ϩ fragmentation at i ϭ 0°, with 78% of the ions dissociating in this way. At i ϭ 45°shattering is much less important. For i ϭ 0°there are 304 different pathways, many related by their backbone cleavage patterns. For the i ϭ 0°fragmentations, 59% resulted from both a-x and b-y cleavages, while for i ϭ 45°70% of the fragmentations occurred with only a-x cleavage. For i ϭ 0°, the average percentage energy transfers to the internal degrees of freedom of the ion and the surface, and the energy remaining in ion translation are 45%, 26%, and 29%. For 45°these percentages are 26%, 12%, and 62%. The percentage energy-transfer to ⌬E int for i ϭ 0°is larger than that reported in previous experiments for collisions of des-Arg 1 -bradykinin with a diamond surface at the same i . This difference is discussed in terms of differences between the model diamond surface used in the simulations and the diamond surface prepared for the experiments. (J Am Soc Mass Spectrom 2009, 20, 939 -948)
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.
Algorithms are presented for sampling quantum microcanonical ensembles for a potential energy minimum and for the conical intersection at the minimum energy crossing point of two coupled electronic states. These ensembles may be used to initialize trajectories for chemical dynamics simulations. The unimolecular dynamics of a microcanonical ensemble about a potential energy minimum may be compared with the dynamics predicted by quantum Rice-Ramsperger-Kassel-Marcus (RRKM) theory. If the dynamics is non-RRKM, it will be of particular interest to determine which states have particularly long lifetimes. Initializing a microcanonical ensemble for the electronically excited state at a conical intersection is a model for electronic nonadiabatic dynamics. The trajectory surface-hopping approach may be used to study the ensuing chemical dynamics. A strength of the model is that zero-point energy conditions are included for the initial nonadiabatic dynamics at the conical intersection.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.