The role of electron-nuclear coupling on multi-state photoelectron spectra, scattering processes and phase transitions

Ortho-Para Conversion for H+ + H2 Collision in Low Temperature: A Fully Close-Coupled Time-Dependent Wave Packet Study

The collision-induced rate coefficients of ortho-para conversion for the H+ + H2 reaction provide accurate information to probe the lifetime of cold environments in interstellar media. Rotationally resolved reaction probabilities are calculated at the low collision energy regime (0 < Ecol ≤ 0.3 eV) by employing the coupled three-dimensional (3D) time-dependent wave packet (TDWP) formalism in hyperspherical coordinates on a recently constructed ab initio ground adiabatic potential energy surface of H3+ [J. Chem. Phys. 2014, 141, 204306] for the process H+ + H2 (v = 0, j = 0-5) → H+ + H2 (v' = 0, j'). Cross-sections are then computed from the converged reaction probabilities as a function of total angular momentum (J) over the same energy regime and subsequently employed to obtain the rate constants for the ortho-to-para (O-P) and para-to-ortho (P-O) conversions and their ratio. The ratio of ortho-para conversion shows (a) appropriate convergence with the inclusion of a higher number of initial rotational states as well as a reasonable agreement with the results from a quantum statistical method and (b) a peak at lower temperature that could be due to the available collision energy for transitions involving lower rotational states (j = 0 → j' = 1).

Adiabat-to-Diabat Angle in Seam Space: Renner-Teller-Type and Pseudo-Jahn-Teller-Type Problems

In this study, we examine the adiabat-to-diabat (ATD) angles for trajectories in 2-dimensional vibrational subspace of the seam space of two degenerate states. In circulating around the tangential touching degeneracy center, the ATD angle is changed by2 π ${2\pi }$ or 0, similar to the Renner-Teller problem and the pseudo-Jahn-Teller problem, respectively. These ATD angle profiles may be indistinguishable from those of circulating multiple conical intersections or a pseudo-Jahn-Teller center. Methods to discern those seemingly indistinguishable cases are proposed. A sharp zigzag variation of the ATD angle is seen as a feature for trajectories that graze a pseudo-Jahn-Teller-type tangential touching center, in contrast to the monotonic steep variation for grazing a conical intersection or a Renner-Teller-type tangential touching center.

Coupled 3D (J ≥ 0) Time-Dependent Wave Packet Calculation for the F + H2 Reaction on Accurate Ab Initio Multi-State Diabatic Potential Energy Surfaces

We had calculated adiabatic potential energy surfaces (PESs), nonadiabatic, and spin-orbit (SO) coupling terms among the lowest three electronic states (12A', 22A', and 12A″) of the F + H2 system using the multireference configuration interaction (MRCI) level of theory, and the adiabatic-to-diabatic transformation equations were solved to formulate the diabatic Hamiltonian matrix [J. Chem. Phys. 2020, 153, 174301] for the entire region of the nuclear configuration space. The accuracy of such diabatic PESs is explored by performing scattering calculations to evaluate integral cross sections (ICSs) and rate constants. The nonadiabatic and SO effects are studied by utilizing coupled 3D time-dependent wave packet formalism with zero and nonzero total angular momentum on multiple adiabatic/diabatic surfaces calculation. We depict the convergence profiles of reaction probabilities for the reactive as well as nonreactive processes on various electronic states at different collision energies with respect to total angular momentum including all helicity quantum numbers. Finally, total ICSs are calculated as functions of collision energies for the initial rovibrational state (v = 0, j = 0) of the H2 molecule along with the temperature-dependent rate coefficient, where those quantities are compared with previous theoretical and experimental results.

Photoelectron spectra of benzene: Can path dependent diabatic surfaces provide unique observables?

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̃2E1g, B̃2E2g, and C̃2A2u) of benzene radical cation (C6H6+) 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.

Beyond Born-Oppenheimer Constructed Diabatic Potential Energy Surfaces for HeH2. J Phys Chem A 2023;127:3832-3847. [PMID: 37098130 DOI: 10.1021/acs.jpca.3c01047]

First-principles based beyond Born-Oppenheimer theory has been employed to construct multistate global Potential-Energy Surfaces (PESs) for the HeH₂⁺ system by explicitly incorporating the Nonadiabatic Coupling Terms (NACTs). Adiabatic PESs and NACTs for the lowest four electronic states (1²A', 2²A', 3²A' and 4²A') are evaluated as functions of hyperangles for a grid of fixed values of the hyperradius in hyperspherical coordinates. Conical intersection between different states are validated by integrating the NACTs along appropriately chosen contours. Subsequently, adiabatic-to-diabatic (ADT) transformation angles are determined by solving the ADT equations to construct the diabatic potential matrix for the HeH₂⁺ system which are smooth, single-valued, continuous, and symmetric and are suitable for performing accurate scattering calculations for the titled system.

Jahn-Teller Effect in Orthorhombic Manganites: Ab Initio Hamiltonian and Roto-vibrational Spectrum

For the first time, using three different electronic structure methodologies, namely, CASSCF, RS2c, and MRCI(SD), we construct ab initio adiabatic potential energy surfaces (APESs) and nonadiabatic coupling term (NACT) of two electronic states (⁵E_g) of MnO₆^9- unit, where eight such units share one La atom in LaMnO₃ crystal. While fitting those APESs with analytic functions of normal modes (Q_x, Q_y), an empirical scaling factor is introduced considering the mass ratio of eight MnO₆^9- units with and without one La atom to explore the environmental (mass) effect on MnO₆^9- unit. When the roto-vibrational levels of MnO₆^9- Hamiltonian are calculated, peak positions computed from ab initio constructed excited APESs are found to be enough close with the experimental satellite transitions [ J. Exp. Theor. Phys. 2016, 122, 890-901] endorsing our earlier model results [ J. Chem. Phys. 2019, 150, 064703]. In order to explore the electron-nuclear coupling in an alternate way, theoretically "exact" and numerically "accurate" beyond Born-Oppenheimer (BBO) theory based diabatic potential energy surfaces (PESs) of MnO₆^9- are constructed to generate the photoelectron (PE) spectra. The PE spectral band also exhibits good peak by peak correspondence with the higher satellite transitions in the dielectric function spectra of the LaMnO₃ complex.

Beyond Born-Oppenheimer based diabatic surfaces of 1,3,5-C6H3F3+ to generate the photoelectron spectra using time-dependent discrete variable representation approach

In this article, Beyond Born-Oppenheimer (BBO) treatment is implemented to construct diabatic potential energy surfaces (PESs) of 1,3,5-C₆H₃F₃⁺ over a series [eighteen (18)] of two-dimensional (2D) nuclear planes constituted with eleven normal modes (Q₂, Q_9x, Q_9y, Q_13x, Q_13y, Q_18x, Q_18y, Q_10x, Q_10y, Q_12x and Q_12y) to include all possible nonadiabatic interactions among six coupled electronic states (X̃²E'', , B̃²E' 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-C₆H₃F₃. Our theoretically calculated spectra for X̃²E'', 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.

Conical intersection seams in spin-orbit coupled systems with an even number of electrons: A numerical study based on neural network fit surfaces

In this work, we consider the existence and topography of seams of conical intersections (CIs) for two key singlet-triplet systems, including a uniformly scaled spin-orbit interaction. The basic one triplet and one singlet state system denoted as (S₀,T₁) and the two singlets and one triplet system denoted as (S₀,S₁,T₁) are treated. Essential to this analysis are realistic electronic structure data taken from a recently reported neural network fit for the 1,2¹A and 1³A states of NH₃, including H_sf (spin-free) and H_so (spin-orbit) surfaces derived from high quality ab initio wavefunctions. Three types of seams for the (S₀,S₁,T₁) system are reported, which depend on the choice of the electronic Hamiltonian, H^e. The nonrelativistic CI seam [H^e = H_sf, (S₀,S₁)], the energy minimized nonrelativistic singlet-triplet intersection seam [H^e = H_sf, (S₀,T₁)], and the fully relativistic seam in the spin-diabatic representation (H^e = H_tot = H_sf + H_so) are reported as functions of R(N-H). The derivative couplings are computed using H^e = H_tot and H_sf from the fit data. The line integral of the derivative coupling is employed to juxtapose the geometric phase in the relativistic, H^e = H_tot, and nonrelativistic, H^e = H_sf, cases. It is found for the (S₀,T₁) system that there is no CI in the spin-adiabatic representation, while for the (S₀,S₁,T₁) system, CI can only be formed for two pairs of spin-adiabatic electronic states. The geometric phase effect thus needs to be handled with care when it comes to spin-nonconserving dynamics simulations.

Dynamical calculations of O(3P) + OH(2Π) reaction on the CHIPR potential energy surface using the fully coupled time-dependent wave-packet approach in hyperspherical coordinates

We have carried out quantum dynamics calculations for the O + OH → H + O₂ reaction on the CHIPR [A. J. C. Varandas, J. Chem. Phys., 2013, 138, 134117] potential energy surface (PES) for ground state HO₂ 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 + O₂ → O + OH) and backward (O + OH → H + O₂) reactions of this reactive system, the equilibrium constant of the reversible process [H + O₂ ⇌ O + OH] is calculated as a function of temperature and compared with previous experimental measurements.