1
|
Sardar S. An Exhaustive Quantum-Classical Study of C 6F 6+ Using the Newly Formulated Parallel TDDVR Method. J Phys Chem A 2024; 128:5777-5795. [PMID: 38979982 DOI: 10.1021/acs.jpca.4c02108] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/10/2024]
Abstract
We recently implemented our parallelized quantum-classical dynamical approach, known as the Time-Dependent Discrete Variable Representation (TDDVR) method, which is applied to the spectroscopically important hexafluorobenzene (HFBz) radical cation, where several conical intersections exist in their seven lowest excited electronic states (S11B2u, S21E1g, S31B1u, S41E1u, and S51A2u) considering degeneracy among potential energy surfaces (PESs), to demonstrate their various dynamical aspects. This new parallel version shows almost linear scalability with increasing number of computing processors. To get photoelectron (PE) spectra, Mass-Analyzed Threshold Ionization (MATI) spectra, population dynamics, and many other dynamical observables, the first-principles dynamics is applied at the state-of-the-art level to the corresponding Hamiltonian, where the Jahn-Teller (JT) and pseudo-Jahn-Teller (PJT) type interactions are involved in those coupled seven electronic states. The quantum-classical method is used to thoroughly analyze the effects of these couplings on the nuclear dynamics of the involved electronic states, and the findings are compared with those observables obtained from experiments. Intrinsic dynamical properties are explained using the reduced densities of the wave packet (WP) in a coupled electronic manifold. The PE and MATI spectra of HFBz computed using TDDVR are found to be in good agreement with earlier experimental data and other theoretically simulated spectra.
Collapse
Affiliation(s)
- Subhankar Sardar
- Department of Chemistry, Bhatter College, Dantan, P.O. Dantan, Paschim Medinipur, Pin 721426, India
| |
Collapse
|
2
|
Suzuki K, Kanno M, Koseki S, Kono H. A Structure-Based Gaussian Expansion for Quantum Reaction Dynamics in Molecules: Application to Hydrogen Tunneling in Malonaldehyde. J Phys Chem A 2023; 127:4152-4165. [PMID: 37129441 DOI: 10.1021/acs.jpca.2c09088] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/03/2023]
Abstract
We developed an approximate method for quantum reaction dynamics simulations, namely, a structure-based Gaussian (SBG) expansion approach, where SBG bases for the expansion of the wave function Ψ, expressed by a product of single-atom Cartesian Gaussians centered at the positions of respective nuclei, are mainly placed around critical structures on reaction pathways such as on the intrinsic reaction coordinate (IRC) through a transition state. In the present approach, the "pseudo-lattice points" at which SBGs are deployed are selected in a perturbative manner so as to make moderate the expansion length. We first applied the SBG idea to a two-dimensional quadruple-well model and obtained accurate tunneling splitting values between the lowest four states. We then applied it to hydrogen tunneling in malonaldehyde and achieved a tunneling splitting of 27.1 cm-1 with only 875 SBGs at the MP2/6-31G(d,p) level of theory, in good agreement with 25 cm-1 by the more elaborate multiconfiguration time-dependent Hartree method. Reasonable results were also obtained for singly and doubly deuterated malonaldehyde. We analyzed the tunneling states by utilizing expansion coefficients of individual SBGs and found that 40-45% of the SBGs in Ψ are nonplanar structures and SBGs away from the IRC contribute a little to hydrogen transfer.
Collapse
Affiliation(s)
- Kazuma Suzuki
- Department of Chemistry, Graduate School of Science, Tohoku University, Sendai 980-8578, Japan
| | - Manabu Kanno
- Department of Chemistry, Graduate School of Science, Tohoku University, Sendai 980-8578, Japan
| | - Shiro Koseki
- Department of Chemistry, Graduate School of Science, Osaka Prefecture University, Osaka 599-8531, Japan
| | - Hirohiko Kono
- Department of Chemistry, Graduate School of Science, Tohoku University, Sendai 980-8578, Japan
| |
Collapse
|
3
|
Dutta J, Ravi S, Mukherjee S, Ojha AK, Adhikari S. Jahn-Teller Effect in Orthorhombic Manganites: Ab Initio Hamiltonian and Roto-vibrational Spectrum. J Phys Chem A 2022; 126:691-709. [PMID: 35089047 DOI: 10.1021/acs.jpca.1c08912] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
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 (5Eg) of MnO69- unit, where eight such units share one La atom in LaMnO3 crystal. While fitting those APESs with analytic functions of normal modes (Qx, Qy), an empirical scaling factor is introduced considering the mass ratio of eight MnO69- units with and without one La atom to explore the environmental (mass) effect on MnO69- unit. When the roto-vibrational levels of MnO69- 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 MnO69- 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 LaMnO3 complex.
Collapse
Affiliation(s)
- Joy Dutta
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata-700 032, India
| | - Satyam Ravi
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata-700 032, India.,School of Advance Science and Languages, VIT Bhopal University, Bhopal-466114, India
| | - Soumya Mukherjee
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata-700 032, India
| | - Avik Kumar Ojha
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata-700 032, India
| | - Satrajit Adhikari
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata-700 032, India
| |
Collapse
|
4
|
Mukherjee S, Ravi S, Dutta J, Sardar S, Adhikari S. Beyond Born-Oppenheimer based diabatic surfaces of 1,3,5-C 6H 3F 3+ to generate the photoelectron spectra using time-dependent discrete variable representation approach. Phys Chem Chem Phys 2022; 24:2185-2202. [PMID: 35006221 DOI: 10.1039/d1cp04733g] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
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.
Collapse
Affiliation(s)
- Soumya Mukherjee
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata - 700032, India.
| | - Satyam Ravi
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata - 700032, India. .,School of Advance Science and Languages, VIT Bhopal University, Bhopal - 466114, India
| | - Joy Dutta
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata - 700032, India.
| | - Subhankar Sardar
- Department of Chemistry, Bhatter College, Dantan, Paschim Medinipur - 721426, India
| | - Satrajit Adhikari
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata - 700032, India.
| |
Collapse
|
5
|
Mukherjee S, Ravi S, Naskar K, Sardar S, Adhikari S. A beyond Born–Oppenheimer treatment of C6H6+ radical cation for diabatic surfaces: Photoelectron spectra of its neutral analog using time-dependent discrete variable representation. J Chem Phys 2021; 154:094306. [DOI: 10.1063/5.0040361] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Soumya Mukherjee
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| | - Satyam Ravi
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| | - Koushik Naskar
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| | - Subhankar Sardar
- Department of Chemistry, Bhatter College, Dantan, Paschim Medinipur 721426, India
| | - Satrajit Adhikari
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| |
Collapse
|
6
|
A comparative multi-state multi-dimensional quantum-classical dynamics on compact polycyclic aromatic hydrocarbons (CPAHs) by parallel TDDVR method. COMPUT THEOR CHEM 2020. [DOI: 10.1016/j.comptc.2020.113032] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
7
|
Choi S, Vaníček J. A time-reversible integrator for the time-dependent Schrödinger equation on an adaptive grid. J Chem Phys 2019; 151:234102. [PMID: 31864241 DOI: 10.1063/1.5127856] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
One of the most accurate methods for solving the time-dependent Schrödinger equation uses a combination of the dynamic Fourier method with the split-operator algorithm on a tensor-product grid. To reduce the number of required grid points, we let the grid move together with the wavepacket but find that the naïve algorithm based on an alternate evolution of the wavefunction and grid destroys the time reversibility of the exact evolution. Yet, we show that the time reversibility is recovered if the wavefunction and grid are evolved simultaneously during each kinetic or potential step; this is achieved by using the Ehrenfest theorem together with the splitting method. The proposed algorithm is conditionally stable, symmetric, and time-reversible and conserves the norm of the wavefunction. The preservation of these geometric properties is shown analytically and demonstrated numerically on a three-dimensional harmonic model and collinear model of He-H2 scattering. We also show that the proposed algorithm can be symmetrically composed to obtain time-reversible integrators of an arbitrary even order. We observed 10 000-fold speedup by using the tenth-order instead of the second-order method to obtain a solution with a time discretization error below 10-9. Moreover, using the adaptive grid instead of the fixed grid resulted in a 64-fold reduction in the required number of grid points in the harmonic system and made it possible to simulate the He-H2 scattering for six times longer while maintaining reasonable accuracy. Applicability of the algorithm to high-dimensional quantum dynamics is demonstrated using the strongly anharmonic eight-dimensional Hénon-Heiles model.
Collapse
Affiliation(s)
- Seonghoon Choi
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Jiří Vaníček
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
8
|
Mandal S, Ghosh S, Sardar S, Adhikari S. The TDDVR approach for molecular photoexcitation, molecule–surface and triatomic reactive scattering processes. INT REV PHYS CHEM 2018. [DOI: 10.1080/0144235x.2018.1548103] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
Affiliation(s)
- Souvik Mandal
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, India
| | - Sandip Ghosh
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, India
| | | | - Satrajit Adhikari
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, India
| |
Collapse
|
9
|
Sardar S. Classical blended quantum formulation of the parallel TDDVR method to the dynamics on furan. JOURNAL OF THEORETICAL & COMPUTATIONAL CHEMISTRY 2015. [DOI: 10.1142/s021963361550039x] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
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.
Collapse
Affiliation(s)
- Subhankar Sardar
- Department of Chemistry, Bhatter College, Dantan, P. O. Dantan, Paschim Medinipur, Pin-721426, West Bengal, India
| |
Collapse
|