Search

While carrying out Beyond Born–Oppenheimer theory based diabatization, the solutions of adiabatic-to-diabatic transformation equations depend on the paths of integration over two-dimensional cross-sections of multi-dimensional space of nuclear degrees of freedom. It is shown that such path-dependent solutions leading to diabatic potential energy surface matrices computed along any two different paths are related through an orthogonal matrix, and thereby, those surface matrices should provide unique observables. While exploring the numerical validity of the theoretical framework, we construct diabatic Hamiltonians for the five low-lying electronic states ( X ̃ 2 E 1 g , B ̃ 2 E 2 g , and C ̃ 2 A 2 u ) of benzene radical cation ( C 6 H 6 + ) along three different approaches of contour integration over two dimensional nuclear planes constituted by seven non-adiabatically active normal modes. Three different diabatic surface matrices are further employed to generate the photoelectron spectra of the benzene molecule (C6H6). It is interesting to note that the spectral peak positions and intensity patterns for all three cases are almost close to each other and also exhibit very good agreement with the experimental results. [ABSTRACT FROM AUTHOR]

We have carried out fully close-coupled three dimensional quantum mechanical wave packet dynamical calculations for the reaction He + H 2 + → HeH + + H on the ground electronic adiabatic potential energy surface and on the lowest two electronic states of newly constructed ab initio calculated diabatic potential energy surfaces for the system [Naskar et al., J. Phys. Chem. A 127, 3832 (2023)]. With the reactant diatom ( H 2 + ) in its roto-vibrational ground state (v = 0, j = 0), the calculations have been carried out in hyperspherical coordinates to obtain the reaction attributes. Convergence profiles of the reaction probability with respect to the total angular momentum quantum number at different collision energies are presented for the title reaction. State-to-state as well as initial state selected integral reaction cross sections are calculated from the fully converged reaction probabilities over a range of collision energies. The integral cross section values computed using the two-state diabatic potential energy surfaces are significantly lower than those obtained using the ground electronic state adiabatic potential energy surface and are in much better agreement with the available experimental results than the latter for total energy greater than 1.1 eV. Therefore, it becomes clear that it is important to include the nonadiabatic coupling terms for a quantitative prediction of the dynamical observables. [ABSTRACT FROM AUTHOR]

The effect of surface mode vibrations on the reactive scattering of D2, initialized in the ground rovibrational state (v = 0, j = 0), from a Cu(111) surface is investigated for different surface temperature situations. We adopt a time and temperature dependent effective Hamiltonian [Dutta et al., J. Chem. Phys. 154, 104103 (2021)] constructed by combining the linearly coupled many oscillator model [Sahoo et al., J. Chem. Phys. 136, 084306 (2012)] and the static corrugation model [M. Wijzenbroek and M. F. Somers, J. Chem. Phys. 137, 054703 (2012)] potential within the mean-field approach. Such an effective Hamiltonian is employed for six-dimensional quantum dynamical calculations to obtain temperature dependent reaction and state-to-state scattering probability profiles as a function of incidence energy of colliding D2 molecules. As reported in the experimental studies, the movements of surface atoms modify the dissociative scattering dynamics at higher surface temperature by exhibiting vibrational quantum and surface atoms' recoil effects in the low and high collision energy domains, respectively. Finally, we compare our present theoretical results with the experimental and other theoretical outcomes, as well as discuss the novelty of our findings. [ABSTRACT FROM AUTHOR]

[Display omitted] We employ trajectory surface hopping (TSH) formalism [Chem. Phys. 349 (2008) 334, J. Chem. Phys. 142 (2015) 144307] to compute cross-sections or/and rate constants of D+ + H 2 (v = 0 , j = 0) and H + H 2 + (v = 0 , j = 0) reactions initiating from ground and first excited electronic state of H 3 + , respectively. While solving Hamilton's equations over adiabatic potential energy surfaces (PESs), switching probabilities are calculated by ab initio nonadiabatic coupling terms for the first time. During hopping, momenta are adjusted to fulfil total energy conservation, where a scaling factor is introduced to explore the conservation of linear momentum. For non-charge transfer process, quasi-classical cross-sections match well with quantum results obtained from diabatic PESs of H 3 + and experimental findings, but deviations are observed at the energy domain of conical intersection driven hopping for charge-transfer processes. Calculated results indicate the origin of limitations of TSH algorithm as well as depict computational advantage with respect to quantum dynamics. [ABSTRACT FROM AUTHOR]