Subhankar Sardar

The multistate multimode vibronic dynamics of benzene radical cation with a realistic model Hamiltonian using a parallelized algorithm of the quantumclassical approach

We demonstrate the workability of a parallelized algorithm of the time-dependent discrete variable representation (TDDVR) method to explore the detailed dynamical aspects of vibronic interaction in two three-state model Hamiltonians (X (2)E(1g), B (2)E(2g), C (2)A(2u) and B (2)E(2g), D (2)E(1u), E (2)B(2u)) of benzene radical cation along with a preliminary investigation on its five electronic states (X (2)E(1g), B (2)E(2g), C (2)A(2u), D (2)E(1u), and E(2)B(2u)). Since those electronic states are interconnected through a series of conical intersections, we have used six and nine vibronically important modes for the three- and five-state Hamiltonians, respectively, in order to perform the quantum dynamics on such system. The population profiles calculated by using our TDDVR approach show reasonably good agreement with the results obtained by exact quantum mechanical (multiconfiguration time-dependent Hartree) method, whereas the corresponding (calculated) photoabsorption spectra originating from various electronic states agree well with the experimental ones. It is important to note that the parallelized algorithm of our TDDVR approach reduces the computation cost by more than an order of magnitude compared to its serial analog. The TDDVR approach appears to be a good compromise between accuracy and speed for such large molecular system, where quantum mechanical description is needed in a restricted region.

A quantum-classical approach to the molecular dynamics of butatriene cation with a realistic model Hamiltonian

We are investigating the molecular dynamics of the butatriene cation after excitation from the ground state (X(2)B(2g)) to the first excited electronic state (A(2)B(2u)) by using the time-dependent discrete variable representation (TDDVR) method. The investigation is being carried out with a realistic 18-mode model Hamiltonian consisting of all the vibrational degrees of freedom of the butatriene molecule. First, we perform the simulation on a basic five mode model, and then by including additional thirteen modes as bath on the basic model. This sequential inclusion of bath modes demonstrates the effect of so called weak modes on the subsystem, where the calculations of energy and population transfer from the basic model to the bath quantify the same effect. The spectral profile obtained by using the TDDVR approach shows reasonably good agreement with the results calculated by the quantum mechanical approach/experimental measurement. It appears that the TDDVR approach for those large systems where quantum mechanical description is needed in a restricted region, is a good compromise between accuracy and speed.

The effect of phonon modes on the H2(v, j)/D2(v, j)–Cu(1nn) scattering processes

We include the effect of the phonon modes originating from the three layers of Cu(1nn) surface atoms on the dynamics of incoming molecular [H(2)(v, j)/D(2)(v, j)] degrees of freedom (DOFs) through a mean-field approach, where the surface temperature is incorporated into the effective potential by considering Bose-Einstein probability (BEP) factor for the initial state distribution of the surface modes calculated within harmonic approximation. Such time and temperature dependent effective Hamiltonian is further subdivided assuming a weak coupling between the two sets of molecular DOFs, namely, (x, y, z, Z) and (X, Y), respectively, in particular, to reduce the computational cost and the corresponding coupled quantum dynamical equations of motion have been formulated in terms of Time Dependent Discrete Variable Representation (TDDVR) approach. We demonstrate the workability of TDDVR method to investigate the scattering of H(2)(v, j) on Cu(1nn) surface by calculating the reaction probabilities and scattering cross-sections. Calculated results show that the phonon modes affect (a) the state-to-state transition probabilities of the scattered H(2) molecule substantially but chemisorption and physisorption processes negligibly and (b) the reaction probability of the incoming D(2) molecule noticeably.

A parallelised quantum-classical approach to the molecular dynamics of allene ( ) radical cation

We are investigating the molecular dynamics of the allene system using a parallelised Time Dependent Discrete Variable Representation (TDDVR) methodology by employing a corresponding cation ( ) where its two electronic surfaces (A 2 E and B 2 B 2) are vibronically coupled with each other. In fact the allene radical cation exhibits a three-surface system due to the presence of degeneracy in the A 2 E state. Our initial investigation is carried out on a linear vibronically coupled model Hamiltonian consisting of 11 vibrationally active modes with two potential energy surfaces. We included both the bilinear and quadratic coupled terms in the effective three-surface Hamiltonian of the radical cation. The spectral profiles obtained from the higher order Hamiltonian show better agreement with the experimental spectrum than the result corresponding to the linearly coupled one. Along with the spectral calculation, we also analyse the nuclear cum population dynamics of the same system. The TDDVR calculated spectral profile as well as the population dynamics show reasonably good agreement with the results calculated by an exact quantum mechanical (multiconfiguration time-dependent Hartree, MCTDH) approach as well as experimental measurement. It appears that the TDDVR approach for those large systems where the quantum mechanical description is needed in a restricted region, is a good compromise between accuracy and speed.

Multi-state multi-mode nuclear dynamics on three isomers of C6H4F+2 using parallelized TDDVR approach

We have performed molecular dynamics on the three isomers of the difluorobenzene radical cation (C(6)H(4)F(2)(+)) after excitation from the ground state to a specific higher electronically excited state by using our recently implemented parallelized Time-Dependent Discrete Variable Representation (TDDVR) methodology. A five-state eleven-mode realistic model Hamiltonian for o-C(6)H(4)F(2)(+) and two separate five-state ten-mode Hamiltonians for m- and p-isomer of the same radical cation are considered, where those five electronic states are interconnected through several conical intersections in the vicinity of the Franck-Condon (FC) region and thus the dynamics for each case become complex. The photoelectron, mass analyzed threshold ionization spectra and population profiles obtained by using our TDDVR approach show reasonably good agreement with the results obtained by multiconfiguration time dependent Hartree (MCTDH) method. It is worthwhile to mention that the parallelized TDDVR algorithm reduces the computation time by more than an order of magnitude compared to its serial analog and, therefore, such approach appears to be a good compromise between accuracy and speed for a large molecular system.

The effect of phonon modes on the D2(v=0, j=0)–Cu(111) scattering processes

We include the phonon modes originating from the three layers of Cu(111) surface atoms on the dynamics of incoming molecular [D2(v, j)] degrees of freedom (DOFs) through a mean-field approach, where the surface temperature is incorporated into the effective potential by considering the Bose–Einstein probability factor for the initial state distribution of the surface modes calculated within the harmonic approximation. Such a time- and temperature-dependent effective Hamiltonian is further subdivided assuming a weak coupling between two sets of molecular DOFs, namely (x, y, z, Z) and (X, Y), respectively, in particular, to reduce the computational cost, and the corresponding coupled quantum dynamical equations of motion have been formulated in terms of the time-dependent discrete variable representation (TDDVR) approach. We demonstrate the applicability of the TDDVR method to investigate the collision of H2(v, j) on the Cu(100) surface by calculating the reaction probabilities and scattering cross-sections. Calculated results for the D2(v=0, j=0)–Cu(111) system show that the phonon modes affect the state-to-state transition probabilities of the scattered D2 molecule substantially and chemisorption–physisorption processes noticeably.

Ab initio constructed diabatic surfaces of NO2 and the photodetachment spectra of its anion.

A thorough investigation has been performed for electronic structure, topological effect, and nuclear dynamics of NO2 molecule, where the adiabatic potential energy surfaces (PESs), conical intersections between the ground (X(2)A1) and the first excited state (A(2)B2), and the corresponding non-adiabatic coupling terms between those states are recalculated [Chem. Phys. 416, 11 (2013)] to achieve enough accuracy in dynamics. We employ beyond Born-Oppenheimer theory for these two state sub-Hilbert space to carry out adiabatic to diabatic transformation (ADT) to obtain the ADT angles and thereby, to construct single-valued, smooth, and continuous diabatic PESs. The analytic expressions for the adiabatic PESs and ADT angles are provided to represent a two-state three-mode diabatic Hamiltonian of NO2 for performing nuclear dynamics to calculate the photo-electron spectra of its anion. It appears that not only Jahn-Teller type coupling but also Renner-Teller interaction contributes significantly on the overall spectrum. The coupling between the electronic states (X(2)A1 and A(2)B2) of NO2 is essentially through the asymmetric stretching mode, where the functional form of such interaction is distinctly symmetric and non-linear.

Multisurface Multimode Molecular Dynamical Simulation of Naphthalene and Anthracene Radical Cations by Using Nearly Linear Scalable Time-Dependent Discrete Variable Representation Method

The major portion of the algorithm of the time-dependent discrete variable representation (TDDVR) method is recently parallelized using the shared-memory parallelization scheme with the aim of performing dynamics on relatively large molecular systems. Because of the astronomical importance of naphthalene and anthracene, we have investigated their radical cations as models for theoretical simulation of complex photoelectron spectra and nonradiative decay process using the newly implemented parallel TDDVR code. The strong vibronic coupling among the six lowest doublet electronic states makes these polynuclear hydrocarbons dynamically important. The aim of the present investigation is to show the efficiency of our current TDDVR algorithm to perform dynamics on large dimensional quantum systems in vibronically coupled electronic manifold. Both the sequential and the parallelized TDDVR algorithms are almost linear scalable for an increase in number of processors. Because a significant speed-up is achieved by cycling in the correct way over arrays, all of the simulations are performed within a reasonable wall clock time. Our theoretical spectra well reproduce the features of the corresponding experimental analog. The dynamical outcomes, for example, population, photoelectron spectra, and diffused interstellar bands, etc., of our quantum-classical approach show good agreement with the findings of the well-established quantum dynamical method, that is, multi configuration time-dependent Hartree (MCTDH) approach.

Classical blended quantum formulation of the parallel TDDVR method to the dynamics on furan

A major portion (∼98%) of the quantum-classical molecular dynamics approach, namely, Time-Dependent Discrete Variable Representation (TDDVR) method, is recently parallelized using shared-memory parallelization scheme with the aim of performing dynamics on relatively large molecular systems. Therefore, we have chosen the furan as a model system for dynamical simulation. Four lowest singlet excited electronic states 1A2(3s), 1B2(V), 1A1(V*), and 1B1(3p) of furan are vibronically coupled through several conical intersections in the vicinity of Frank-Condon region. The major focus of the present paper is to explore the efficiency of our newly implemented parallelized TDDVR algorithm to perform dynamics on large dimensional quantum systems in vibronically coupled electronic manifold. The present version of parallelized TDDVR algorithm show closely linear speed up with increasing number of computing processors. As a significant speed up is achieved by cycling in the correct way over arrays, all dynamical simulations are performed within a reasonable wall clock time. The photoelectron spectra calculated by the TDDVR method show peak by peak correspondence with the experimental spectra. The TDDVR calculated quantum-classical dynamical outcomes, viz., population and photoelectron spectra, etc. show good agreement with the findings of the well-established quantum dynamical method, i.e. Multi Configuration Time-Dependent Hartree (MCTDH) approach.

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, NO3 radical and C6H6+ (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. A2017,121, 6314–6326]. In the case of NO3 radical, the nonadiabatic coupling terms (NACTs) among the low-lying five electronic states, namely, X̃2A2′ (12B2), A~2E″ (12A2 and 12B1), and B~2E′ (12A1 and 22B2), 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 (Q3x–Q3y pair for NO3 radical and Q16x–Q16y 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 Q3x and Q16x normal modes.