An ab initio investigation of non-adiabatic couplings and conical intersections among the lowest five electronic states of the NO3 radical

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.

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 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.

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

We present first principle based beyond Born-Oppenheimer (BBO) theory and its applications on various models as well as realistic spectroscopic and scattering processes, where the Jahn-Teller (JT) theory is brought in conjunction with the BBO approach on the phase transition of lanthanide complexes. Over one and half decades, our development of BBO theory is demonstrated with ab initio calculations on representative molecules of spectroscopic interest (NO2 radical, Na3 and K3 clusters, NO3 radical, C6H6+ and 1,3,5-C6H3F3+ radical cations) as well as triatomic reactive scattering processes (H+ + H2 and F + H2). Such an approach exhibits the effect of JT, Renner-Teller (RT) and pseudo Jahn-Teller (PJT) type of interactions. While implementing the BBO theory, we generate highly accurate diabatic potential energy surfaces (PESs) to carry out quantum dynamics calculation and find excellent agreement with experimental photoelectron spectra of spectroscopic systems and cross-sections/rate constants of scattering processes. On the other hand, such electron-nuclear couplings incorporated through JT theory play a crucial role in dictating higher energy satellite transitions in the dielectric function spectra of the LaMnO3 complex. Overall, this article thoroughly sketches the current perspective of the BBO approach and its connection with JT theory with various applications on physical and chemical processes.

Beyond Born-Oppenheimer constructed diabatic potential energy surfaces for F + H2 reaction

First principles based beyond Born-Oppenheimer theory has been implemented on the F + H₂ system for constructing multistate global diabatic Potential Energy Surfaces (PESs) through the incorporation of Nonadiabatic Coupling Terms (NACTs) explicitly. The spin-orbit (SO) coupling effect on the collision process of the F + H₂ reaction has been included as a perturbation to the non-relativistic electronic Hamiltonian. Adiabatic PESs and NACTs for the lowest three electronic states (1²A', 2²A', and 1²A″) are determined in hyperspherical coordinates as functions of hyperangles for a grid of fixed values of the hyperradius. Jahn-Teller (JT) type conical intersections between the two A' states translate along C_2v and linear geometries in F + H₂. In addition, A' and A″ states undergo Renner-Teller (RT) interaction at collinear configurations of this system. Both JT and RT couplings are validated by integrating NACTs along properly chosen contours. Subsequently, we have solved adiabatic-to-diabatic transformation (ADT) equations to evaluate the ADT angles for constructing the diabatic potential matrix of F + H₂, including the SO coupling terms. The newly calculated diabatic PESs are found to be smooth, single-valued, continuous, and symmetric and can be invoked for performing accurate scattering calculations on the F + H₂ system.

Topological Effects in Vibronically Coupled Degenerate Electronic States: A Case Study on Nitrate and Benzene Radical Cation

We carry out detailed investigation for topological effects of two molecular systems, NO₃ radical and C₆H₆ ⁺ (Bz⁺) radical cation, where the dressed adiabatic, dressed diabatic, and adiabatic-via-dressed diabatic potential energy curves (PECs) are generated employing ab initio calculated adiabatic and diabatic potential energy surfaces (PESs). We have implemented beyond Born-Oppenheimer (BBO) theory for constructing smooth, single-valued, and continuous diabatic PESs for five coupled electronic states [J. Phys. Chem. A 2017, 121, 6314-6326]. In the case of NO₃ radical, the nonadiabatic coupling terms (NACTs) among the low-lying five electronic states, namely, X̃ ²A₂ ^' (1²B₂), A~ ²E″ (1²A₂ and 1²B₁), and B~ ²E' (1²A₁ and 2²B₂), bear the signature of Jahn-Teller (JT) interactions, pseudo JT (PJT) interactions, and accidental conical intersections (CIs). Similarly, Bz⁺ radical cation also exhibits JT, PJT, and accidental CIs in the interested domain of nuclear configuration space. In order to generate dressed PECs, two components of degenerate in-plane asymmetric stretching modes are selectively chosen for both the molecular species (Q _3x -Q _3y pair for NO₃ radical and Q _16x -Q _16y pair for Bz⁺ radical cation). The JT coupling between the electronic states is essentially originated through the asymmetric stretching normal mode pair, where the coupling elements exhibit symmetric and nonlinear functional behavior along Q _3x and Q _16x normal modes.

Aqueous Solvation of SmI3: A Born-Oppenheimer Molecular Dynamics Density Functional Theory Cluster Approach

We report the results of Born-Oppenheimer molecular dynamics (BOMD) simulations on the aqueous solvation of the SmI₃ molecule and of the bare Sm³⁺ cation at room temperature using the cluster microsolvation approach including 37 and 29 water molecules, respectively. The electronic structure calculations were done using the M062X hybrid exchange-correlation functional in conjunction with the 6-31G** basis sets for oxygen and hydrogen. For the iodine and samarium atoms, the Stuttgart-Köln relativistic effective-core potentials were utilized with their associated valence basis sets. When SmI₃ is embedded in the microsolvation environment, we find that substitution of the iodine ions by water molecules around Sm(III) cannot be achieved due to an insufficient number of explicit water molecules to fully solvate the four separate metal and halogen ions. Therefore, we studied the solvation dynamics of the bare Sm³⁺ cation with a 29-water molecule model cluster. Through the Sm-O radial distribution function and the evolution of the Sm-O distances, the present study yields a very tightly bound first rigid Sm(III) solvation shell from 2.3 to 2.9 Å whose integration leads to a coordination number of 9 water molecules and a second softer solvation sphere from 3.9 to 5 Å with 12 water molecules. No water exchange processes were found. The theoretical EXAFS spectrum is in excellent agreement with the experimental spectrum for Sm(III) in liquid water. The strong differences between the solvation patterns of Sm(III) vs Sm(II) are discussed in detail.

Five Electronic State Beyond Born-Oppenheimer Equations and Their Applications to Nitrate and Benzene Radical Cation

We present explicit form of Adiabatic to Diabatic Transformation (ADT) equations and expressions of non-adiabatic coupling terms (NACTs) for a coupled five-state electronic manifold in terms of ADT angles between electronic wave functions. ADT matrices eliminate the numerical instability arising from singularity of NACTs and transform the adiabatic Schrödinger equation to its diabatic form. Two real molecular systems NO₃ and C₆H₆⁺ (Bz⁺) are selectively chosen for the demonstration of workability of those equations. We examine the NACTs among the lowest five electronic states of the NO₃ radical [X̃²A₂^' (1²B₂), òE″ (1²A₂ and 1²B₁) and B̃²E' (1²A₁ and 2²B₂)], in which all types of non-adiabatic interactions, that is, Jahn-Teller (JT) interactions, Pseudo Jahn-Teller (PJT) interactions, and accidental conical intersections (CIs) are present. On the other hand, lowest five electronic states of Bz⁺ [X̃²E_1g (1²B_3g and 1²B_2g), B̃²E_2g (1²A_g and 1²B_1g), and C̃²A_2u (1²B_1u)] depict similar kind of complex feature of non-adiabatic effects. For NO₃ radical, the two components of degenerate in-plane asymmetric stretching mode are taken as a plane of nuclear configuration space (CS), whereas in case of Bz⁺, two pairs are chosen: One is the pair of components of degenerate in-plane asymmetric stretching mode, and the other one is constituted with one of the components each from out-of-plane degenerate bend and in-plane degenerate asymmetric stretching modes. We calculate ab initio adiabatic potential energy surfaces (PESs) and NACTs among the lowest five electronic states at the CASSCF level using MOLPRO quantum chemistry package. Subsequently, the ADT is performed using those newly developed equations to validate the positions of the CIs, evaluate the ADT angles and construct smooth, symmetric, and continuous diabatic PESs for both the molecular systems.