TSH calculations

Trajectory surface-hopping (TSH) calculations have been accomplished for H3+ on DIM hypersurfaces.2,498 A surface-hopping mechanism has also been suggested as applicable to several other systems, including various He2+ charge-transfer reactions,495 the H -H2 reaction,499 and the C+-H2 radiative association process.500... 

Once we know how to calculate and for these five types of processes, we can calculate these quantities for any process by taking advantage of the fact that 0 and H are state properties. Tlie procedure is to construct a hypothetical process path from the initial state to the final state consisting of a series of steps of the given five types. Having done this, we calculate A/ for each of the steps, and then add the A/ s for the steps to calculate A// for the total process. Since H is a state property, tsH calculated for the hypothetical process path—which we constructed for convenience—is the same as for the path actually followed by the process. The same procedure can be followed to calculate Af for any process. 

With respect to sharp oscillations in reaction probability that occur at higher energy due to reflection from the excited singlet, we expect that this is a real effect that should be captured by TSH calculations. Indeed, somewhat analogous... 

A related question is whether there would be differences in reactivity between different tine structure states of 0( P) due to the accidental cancellation that causes state 1 to be largely dec oupled from state 4. TIk presenc.e or absence of this effect in 3D TSH calculations would be crucially dej)endence on the underlying basis set (adiabatic or diabatic) used. For exanijjle, this would not show up if the diabatie basis set that we used to set up our problem were used, but it should in an adiabatic basis. However, the signatures of this effect are subtle, as it mostly leads to a sujjjjression of reactivity in state 2 comjjared to state 1, but there is already suppression close to threshold due to the thicker barrier. 

Total cross sections, translational energy distributions and velocity contour diagrams for many of the products for the H+ + D2 and D+ + HD reactions have been obtained from the crossed ion beam-molecular beam measurements and the TSH calculations making these reactions probably the best understood in the field. This combined experimental and theoretical investigation provided a unique opportunity for a detailed comparison of an a priori theory of chemical reactions with experiments. These comparisons are very encouraging since the TSH theory involves no empirical information or adjustable parameters. 

Recently, Hoffmann and Schatz(7ij have developed a new level of treatment of spin-orbit effects in bimolecular reactions which enables a more sophisticated treatment of intersystem crossing d3mamics than in the past. In this treatment high quality electronic structure n thods are used to determine global surfaces for the reaction and spin-orbit matrix elements, and then trajectory surface hopping (TSH) methods are used to determine properties of the bimolecular collisions such as reactive cross sections and state distribution information. In an application of this theory to the O + H2 reaction, the spin-orbit matrix elements were determined as a function of position, and then TSH calculations were done within a diabatic representation to determine cross sections. Intersystem crossing effects were found to be small for O + H2 due to... 

TSH calculations were done using the fewest switches method 13). )etails are similar to work that was described earlier (77), except that here we have considered calculations in both the diabatic and adiabatic representations. A time step of 1 atomic unit was used for all calculations, and 2000 trajectories were used to determine the reaction probabilities at each energy. Variation of the results with respect to the numerical parameters was within statistical uncertainty. We found the diabatic results to be significantly less sensitive to time step dian the adiabatic results. 

MeONa both attack the FOSS core), the hydrolysis being calculated to be endothermic with a calculated energy significantly dependent on the method of calculation. The reactivity of TsHs appears to be lower than expected when compared with simple silanes such as (MeOlsSi-H which may be due to the steric bulk and relative rigidity of the FOSS cage. ... 

Zeno paradox. On the other hand, recovering these interferences from a single path leads to excessive correlation, as evidenced by the highly oscillatory results obtained with TSH for Tully s third, extended coupling with reflection, model. This is remedied effortlessly in FMS, and one may speculate that FMS will tend to the opposite behavior Interferences that are truly present will tend to be damped if insufficient basis functions are available. This is probably preferable to the behavior seen in TSH, where there is a tendency to accentuate phase interferences and it is often unclear whether the interference effects are treated correctly. This last point can be seen in the results of the second, dual avoided crossing, model, where the TSH results exhibit oscillation, but with the wrong structure at low energies. The correct behavior can be reproduced by the FMS calculations with only ten basis functions [38]. 

Tsh,MD = 0 °C- Pc = 008 mbar. Tsh,SD = 40 °C Below dW data calculated from DR measurements during secondary drying as a function of drying time. ( ) Mean dW of all three runs bars, minimum and maximum of all calculated dW... 

The data and the discussion above should be taken as guide for a process development, not as absolute. The data are from one specific project and rounded for simplicity, e.g. at a given p., Tice will vary with Tsh, but the result will not be round figures as shown the tMD has been calculated with a constant K o = 60 kj/ °C h m2 as shown in Figurel.58 and Table 1.10.1, Ktot decreases with decreasing pressure. The dependence of Tjce on pr can be very different as used in the example (see Section 2.6). 

DR can be calculated from automatic pressure rise measurements similarly to Tice. During SD no time limit for the pressure rise exists, the product has to warm up close to Tsh anyway. The accuracy of the DR measurement depends on the reproducibility of the vacuum gauge and the amount of solids per chamber volume (see Eq. (16a)). As a rule of thumb we can use 1 g solid/L of chamber volume results in DR 0.1%/h 10%, if dt = 60 °C. For the pressure rise time, 60-120 s is recommended depending on the product temperature and how well the rule 1 g/L is met (longer measuring time for ratios <1). 

Table 3.6 Shelf temperatures (Tsh), measured sublimation temperatures (TMD), measured time of sublimation (tMDm), calculated time of sublimation (tMDc), calculated temperature during desorption drying (Ts )) and estimated total drying time (t,ot) /Tables 1 and 2 from [3.80])...

Thus, molecular vibrations play a crucial role in ET processes, and MD calculations involving quantum transitions, such as our MD-TSH method, offer considerable promise as a powerful tool for investigating their role in ET reactions. 

---