Search

We construct the effective Hartree potential for H2 on Cu(111) as introduced in our earlier work [Dutta et al., J. Chem. Phys. 154, 104103 (2021), and Dutta et al., J. Chem. Phys. 157, 194112 (2022)] starting from the same gas–metal interaction potential obtained for 0 K. Unlike in that work, we now explicitly account for surface expansion at 925 K and investigate different models to describe the surface vibrational modes: (i) a cluster model yielding harmonic normal modes at 0 K and (ii) slab models resulting in phonons at 0 and 925 K according to the quasi-harmonic approximation—all consistently calculated at the density functional theory level with the same exchange–correlation potential. While performing dynamical calculations for the H2(v = 0, j = 0)–Cu(111) system employing Hartree potential constructed with 925 K phonons and surface temperature, (i) the calculated chemisorption probabilities are the highest compared to the other approaches over the energy domain and (ii) the threshold for the reaction probability is the lowest, in close agreement with the experiment. Although the survival probabilities (v′ = 0) depict the expected trend (lower in magnitude), the excitation probabilities (v′ = 1) display a higher magnitude since the 925 K phonons and surface temperature are more effective for the excitation process compared to the phonons/normal modes obtained from the other approaches investigated to describe the surface. [ABSTRACT FROM AUTHOR]

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]

The "accuracy" of diabatic Hamiltonian matrix elements solely depends on the level of ab initio calculation and on the included set of coupled electronic manifold forming a "true" sub‐Hilbert space (SHS). If ab initio calculations are performed with enough accuracy, existence of SHS can be validated only when the Curl Condition attains zero within numerical precision. For the first time, we depict a theoretical and algorithmic formulation on the calculation of Curl Condition to demonstrate the existence of SHS on molecular processes or chemical reactions. [ABSTRACT FROM AUTHOR]

In this article, Beyond Born–Oppenheimer (BBO) treatment is implemented to construct diabatic potential energy surfaces (PESs) of 1,3,5-C6H3F3+ over a series [eighteen (18)] of two-dimensional (2D) nuclear planes constituted with eleven normal modes (Q2, Q9x, Q9y, Q13x, Q13y, Q18x, Q18y, Q10x, Q10y, Q12x and Q12y) to include all possible nonadiabatic interactions among six coupled electronic states (X~;2E′′, , B~;2E′ and). We had formulated explicit expressions of adiabatic to diabatic transformation (ADT) equations [S. Mukherjee, J. Dutta, B. Mukherjee, S. Sardar and S. Adhikari, J. Chem. Phys., 2019, 150, 064308] for the same system forming six state sub-Hilbert space and at present, these ADT equations are solved by incorporating MRCI level ab initio adiabatic PESs and CP-MCSCF calculated nonadiabatic coupling terms (NACTs) to derive diabatic PESs and couplings. Such single-valued, smooth, symmetric and continuous diabatic surface matrices are utilized to carry out multi-state multi-mode nuclear dynamics with the help of time-dependent discrete variable representation (TDDVR) methodology to compute the photoelectron (PE) spectra of 1,3,5-C6H3F3. Our theoretically calculated spectra for X~;2E′′, and states using BBO treatment and TDDVR dynamics show peak by peak correspondence with the experimental results as well as better than the findings of the multi-configuration time-dependent Hartree (MCTDH) method. [ABSTRACT FROM AUTHOR]

We have carried out quantum dynamics calculations for the O + OH → H + O2 reaction on the CHIPR [A. J. C. Varandas, J. Chem. Phys., 2013, 138, 134117] potential energy surface (PES) for ground state HO2 using the fully coupled 3D time-dependent wavepacket formalism [S. Adhikari and A. J. C. Varandas, Comput. Phys. Commun., 2013, 184, 270] in hyperspherical coordinates. Reaction probabilities for J > 0 are calculated for different initial rotational states of the OH radical (v = 0; j = 0, 1). State-to-state as well as total integral cross sections and rate-coefficients are evaluated and compared with previous theoretical calculations and available experimental studies. Using the rate constant for the forward (hereinafter considered to be H + O2 → O + OH) and backward (O + OH → H + O2) reactions of this reactive system, the equilibrium constant of the reversible process [H + O2 O + OH] is calculated as a function of temperature and compared with previous experimental measurements. [ABSTRACT FROM AUTHOR]