1
|
Wang R, Zhao H, Sun Z. Reactant-Product Decoupling Technique Using the Intermediate Coordinate Method. J Phys Chem A 2024; 128:3726-3741. [PMID: 38666315 DOI: 10.1021/acs.jpca.4c01148] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/12/2024]
Abstract
Although the reactant-product decoupling (RPD) technique was proposed over two decades ago, it remains an efficient approach for calculating product state-resolved information on some simple direct reactions using the quantum wave packet method. In the past, usually the RPD technique employed the collocation method to transform the wave function between reactant and product arrangements, which requires quite large computational efforts. In this work, the intermediate coordinate (IC) method is employed to realize the RPD technique. Numerical examples demonstrate that this new IC RPD (IRPD) technique has superior computational efficiency compared with the original method employing the collocation method. Especially, the new IRPD technique significantly saves disk space and computer memory. To illustrate the features of our new method, the total reaction probabilities of the H + H2, H + Br2, and F + H2 reactions with J = 0 and the differential cross sections of the H + H2 and F + H2 reactions at a series of collision energy are calculated and presented. With this efficient and effective new RPD technique, the Li + HF reaction, which involves sharp resonances with long-range wave functions in the van der Waals wells in both the reactant and product arrangements, is also calculated with several J at the product state-resolved level to reveal the ability of the RPD technique for describing resonance wave functions. With these numerical examples, it is found that, for the reaction with resonances, the RPD approach should be applied carefully. Otherwise, it is very possible that the resonances could disappear with the application of the RPD technique.
Collapse
Affiliation(s)
- Ransheng Wang
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Hailin Zhao
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
| | - Zhigang Sun
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
| |
Collapse
|
2
|
Mao Y, Yang Z, Buren B, Chen M. Unveiling Quantum Interference in the D + + H 2 Nonadiabatic Reaction Dynamics at Low Collision Energies. J Phys Chem A 2024; 128:420-430. [PMID: 38174889 DOI: 10.1021/acs.jpca.3c07097] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2024]
Abstract
Fully converged nonadiabatic dynamics calculations of the D+ + H2 → H+ + HD reaction are performed at low temperatures using the time-dependent wave packet approach based on a set of precise 3 × 3 diabatic potential energy surfaces (PESs) ( Phys. Chem. Chem. Phys., 2021, 23, 7735-7747, DOI: 10.1039/D0CP04100A). The D+ + H2 reaction is mediated by a dense manifold of resonances associated with the deep potential well on the ground-state PES. The calculated results show that the nonadiabatic coupling can affect the resonance positions, deviating from the expectation based solely on adiabatic considerations. Furthermore, significant forward-backward asymmetry in total differential cross sections (DCSs) is revealed, which is markedly influenced by nonadiabatic effects. The nonadiabatic effects not only affect the contribution of partial waves in the reaction but also make the interference patterns in the DCSs change significantly.
Collapse
Affiliation(s)
- Ye Mao
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics, Dalian University of Technology, Dalian 116024, PR China
| | - Zijiang Yang
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics, Dalian University of Technology, Dalian 116024, PR China
| | - Bayaer Buren
- School of Science, Shenyang University of Technology, Shenyang 110870, PR China
| | - Maodu Chen
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics, Dalian University of Technology, Dalian 116024, PR China
| |
Collapse
|
3
|
Menéndez M, Garcia E, Lara M, Jambrina PG, Aoiz FJ. Li + HF and Li + HCl Reactions Revisited I: QCT Calculations and Simulation of Experimental Results. J Phys Chem A 2023; 127:6924-6944. [PMID: 37579497 PMCID: PMC10461305 DOI: 10.1021/acs.jpca.3c03763] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2023] [Revised: 07/20/2023] [Indexed: 08/16/2023]
Abstract
The Li + HF and Li + HCl reactions share some common features. They have the same kinematics, relatively small barrier heights, bent transition states, and are both exothermic when the zero point energy is considered. Nevertheless, the pioneering crossed beam experiments by Lee and co-workers in the 80s (Becker et al., J. Chem. Phys. 1980, 73, 2833) revealed that the dynamics of the two reactions differ significantly, especially at low collision energies. In this work, we present theoretical simulations of their results in the laboratory frame (LAB), based on quasiclassical trajectories and obtained using accurate potential energy surfaces. The calculated LAB angular distributions and time-of-flight spectra agree well with the raw experimental data, although our simulations do not reproduce the experimentally derived center-of-mass (CM) differential cross section and velocity distributions. The latter were derived by forward convolution fitting under the questionable assumption that the CM recoil velocity and scattering angle distribution were uncoupled, while our results show that the coupling between them is relevant. Some important insights into the reaction mechanism discussed in the article by Becker et al. had not been contrasted with those that can be extracted from the theoretical results. Among them, the correlation between the angular momenta involved in the reactions has also been examined. Given the kinematics of both systems, the reagent orbital angular momentum, l , is almost completely transformed into the rotation of the product diatom, j'. However, contrary to the coplanar mechanism proposed in the original paper, we find that the initial and final relative orbital angular momenta are not necessarily parallel. Both reactions are found to be essentially direct, although about 15% of the LiFH complexes live longer than 200 fs.
Collapse
Affiliation(s)
- Marta Menéndez
- Departamento
de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Ernesto Garcia
- Departamento
de Química Física, Universidad
del País Vasco (UPV/EHU), 01006 Vitoria, Spain
| | - Manuel Lara
- Departamento
de Química Física Aplicada, Facultad de Ciencias, Universidad Autónoma de Madrid, 28039 Madrid, Spain
| | - Pablo G. Jambrina
- Departamento
de Química Física, Facultad de Ciencias Químicas, Universidad de Salamanca, 37008 Salamanca, Spain
| | - F. Javier Aoiz
- Departamento
de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
4
|
Karabulut E, Celik FA, Korkmaz ET. The long-lived reactive nitrogen species in the troposphere: DFTB model for atmospheric applications. Phys Chem Chem Phys 2023; 25:5569-5581. [PMID: 36727207 DOI: 10.1039/d2cp05344f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
The longest lived reactive NO2 molecule formation in a dry and clean air environment under a high-temperature shock wave was investigated under three basic reactions (R2 for the O + NO system, R6 for the NO + NO3 system, and R7 for the NO + O3 system) in the atmospheric environment. With certain approaches, a DFTB3 model was used, which gave results close to the density functional theory. In the calculations, the related reactions up to 250 ps were examined at individual specific temperatures, and the temperature ranges that contributed to the formation of the NO2 molecule were determined. Moreover, a shock wave with both heating and cooling channels was applied only on R2 to see whether molecular concentrations were in good agreement with atmospheric information. The reaction products were examined under a shock wave of about 20 ps. At the end of the study, the applicability of the DFTB model to atmospheric systems was demonstrated by comparing it with experimental data and information. QCT approach was also used for the calculation of reaction rate constants of only O2-formation on the O + NO system. Here, all systems are focused on nitrogen species containing oxygen. In particular, the highest-population NO molecule that emerged in the lightning flash event was used as the reactant, while systems existing with the longest lived NO2 in the atmosphere after the lightning flash were focused in the product channel. As a result of the study, the hypothesis of geophysicists that almost all NO2 formed in the lightning flash event originates from the NO + O system was disproved. It has been proven that the presence of NO3 molecules that can withstand high temperatures in such systems should be evaluated.
Collapse
Affiliation(s)
- Ezman Karabulut
- Vocational School of Health Services, Bitlis Eren University, 13000 Bitlis, Turkey
| | - Fatih Ahmet Celik
- Faculty of Arts&Sciences, Physics Department, Bitlis Eren University, 13000 Bitlis, Turkey
| | - Ebru Tanboğa Korkmaz
- Faculty of Arts&Sciences, Physics Department, Bitlis Eren University, 13000 Bitlis, Turkey
| |
Collapse
|
5
|
Karabulut E. Oxygen Molecule Formation and the Puzzle of Nitrogen Dioxide and Nitrogen Oxide during Lightning Flash. J Phys Chem A 2022; 126:5363-5374. [PMID: 35920809 DOI: 10.1021/acs.jpca.2c02378] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Unlike the compounds of the natural air atmosphere, the lightning systems are primarily focused on NO(X2Π), NO2(12A'), and O(3P) concentrations that occurred newly and highly in the ground electronic structure. While the NO/NO2 concentrations ratio is about 2000 during the lightning flash, this ratio becomes about 0.8 right after the lightning flash. The reason for this decrease in the ratio is the disappearance of the high temperature that prevents the formation of NO2 (with the combination of NO and O) and of the photon energy that causes its dissociation (NO2 + hv → NO + O) right after the lightning flash. However, this study will focus on the reactions that contribute to the NO concentration, except for the combination of N and O atoms during lightning flash. To do this, it was focused on the reactive scattering states (especially the NO-exchange) of the NO + O collision and the photo-dissociation of NO2, which provide the formation of the NO molecule in the ground electronic state. This case raises important questions. To what extent do the NO-exchange reaction and the photo-dissociation of NO2 contribute to the atmospherically observed NO molecules? or how can the vibrational quantum states of the NO molecules formed by the photo-dissociation be effected on the NO + O1 collision to produce a NO1 molecule? These conditions may contribute to the concentrations of NO high during lightning flashes. Under low collision energy (between 0.1 and 0.3 eV), the NO (v = 0) population dissociated by a photon can act as reactants in the NO-exchange reactive scattering on the doublet electronic state. Since it is assumed that all of the NO2 molecules are due to NO in the lightning flash system, this is one of the reasons that makes the NO population so high during lightning flash. Therefore, in the light of considering that the lightning system supports the formation of highly vibrating molecular groups, it might also support the formation of O2 molecules. In particular, it was shown that the v = 4 quantum state of the NO molecule over the doublet state between collision energies of 0.9-1.5 eV and the v = 5 quantum state of the NO molecule over the quartet state between collision energies of 1.0-1.5 eV contribute to O2 formation.
Collapse
Affiliation(s)
- Ezman Karabulut
- Vocational School of Health Service, Bitlis Eren University, 13000 Bitlis, Turkey
| |
Collapse
|
6
|
Gomez-Carrasco S, Felix-Gonzalez D, Aguado A, Roncero O. Spin-orbit transitions in the N +( 3P JA) + H 2 → NH +(X 2 Π, 4Σ -)+ H( 2S) reaction, using adiabatic and mixed quantum-adiabatic statistical approaches. J Chem Phys 2022; 157:084301. [DOI: 10.1063/5.0102376] [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/14/2022] Open
Abstract
The cross section and rate constants for the title reaction are calculated for all the spin-orbit states of N+(3PJA) using two statistical approaches, one purely adiabatic and the other one mixing quantum capture for the entrance channel and adiabatic treatment for the products channel. This is made by using a symmetry adapted basis set combining electronic (spin and orbital) and nuclear angular momenta in the reactants channel. To this aim, accurate ab initio calculations are performed separately for reactants and products. In the reactants channel, the three lowest electronic states (without spin-orbit couplings) have been diabatized, and the spin-orbit couplings have been introduced through a model localizing the spin-orbit interactions in the N+ atom, which yields accurate results as compared to ab initio calculations including spin-orbit couplings. For the products, eleven purely adiabatic spin-orbit states have been determined with ab initio calculations. The reactive rate constants thus obtained are in very good agreement with the available experimental data for several ortho-H2 fractions, assuming a thermal initial distribution of spin-orbit states. The rate constants for selected spin-orbit JA states are obtained, to provide a proper validation of the spin-orbit effects to obtain the experimental rate constants.
Collapse
Affiliation(s)
| | | | - Alfredo Aguado
- Departamento de Química Física Aplicada, Universidad Autónoma de Madrid, Spain
| | | |
Collapse
|
7
|
Lü Y, Zhang C, Wang H, Guo Q, Li Y. An accurate many-body expansion potential energy surface for AlH 2 (2 2A') and quantum dynamics in Al( 3P) + H 2 ( v0 = 0-3, j0 = 0, 2, 4, 6) collisions. Phys Chem Chem Phys 2022; 24:16637-16646. [PMID: 35766326 DOI: 10.1039/d2cp01802k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
An accurate potential energy surface is constructed for the excited state of AlH2 by fitting extensive ab initio points calculated at the multi-reference configuration interaction level based on aug-cc-pV(Q+d)Z and aug-cc-pV(5+d)Z basis sets. All the calculated energies are corrected via the many-body expansion method and extrapolated to the complete basis set limit. The various topographic features of the new potential energy surface are investigated to demonstrate the correct behavior of Al(3P) + H2(X1Σg+) and AlH(a3Π) + H(2S) dissociation limits. By employing the time-dependent wave packet approach, the integral scattering cross-sections obtained from the Coriolis coupling calculation and the centrifugal sudden approximation, respectively, are compared in detail and show that the former has a higher effect on the reaction. Moreover, the thermal rate constants for Al(3P) + H2 (v0 = 0-3, j0 = 0, 2, 4, 6) in the temperature range of 0-5000 K are calculated, thereby providing insights into the influence of ro-vibrational quantum numbers on the thermal rate constants.
Collapse
Affiliation(s)
- Yanling Lü
- Department of Physics, Liaoning University, Shenyang 110036, China.
| | - ChengYuan Zhang
- Department of Physics, Liaoning University, Shenyang 110036, China.
| | - Hainan Wang
- Department of Physics, Liaoning University, Shenyang 110036, China.
| | - Qiang Guo
- Department of Physics, Liaoning University, Shenyang 110036, China. .,Lvyuan Institute of Energy & Environmental Science and Technology, Liaoning University, Shenyang, 110036, China
| | - Yongqing Li
- Department of Physics, Liaoning University, Shenyang 110036, China. .,Lvyuan Institute of Energy & Environmental Science and Technology, Liaoning University, Shenyang, 110036, China
| |
Collapse
|
8
|
The time-dependent quantum wave packet calculation for reaction S ‒( 2P)+H 2( 1Σ g+)→SH ‒( 1Σ)+H( 2S)
on a new ab initio potential energy surface. CHINESE J CHEM PHYS 2022. [DOI: 10.1063/1674-0068/cjcp2111239] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
|
9
|
Buren B, Chen M, Sun Z, Guo H. Quantum Wave Packet Treatment of Cold Nonadiabatic Reactive Scattering at the State-To-State Level. J Phys Chem A 2021; 125:10111-10120. [PMID: 34767377 DOI: 10.1021/acs.jpca.1c08105] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Cold and ultracold collisions are dominated by quantum effects, such as resonances, tunneling, and nonadiabatic transitions between different electronic states. Due to the extremely long de Broglie wavelength in such processes, quantum reactive scattering is most conveniently characterized using the time-independent close-coupling (TICC) methods. However, the TICC approach is difficult for systems with a large number of channels because of its steep numerical scaling laws. Here, a recently proposed quantum wave packet (WP) approach for solving adiabatic reactive scattering problems at low collision energies is extended to include nonadiabatic transitions. To impose the outgoing boundary conditions, the total scattering wavefunction is split into three parts, the interaction, the asymptotic, and the long-range regions. Each region is associated with a different set of basis functions, which could be optimized separately. In this way, an extremely long grid can be used to accommodate the characteristic long de Broglie wavelengths in the scattering coordinate. The better numerical scaling laws of the WP approach have the potential for handling larger nonadiabatic reactive systems at low temperatures in the future.
Collapse
Affiliation(s)
- Bayaer Buren
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics, Dalian University of Technology, Dalian 116024, China.,Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Maodu Chen
- Key Laboratory of Materials Modification by Laser, Electron, and Ion Beams (Ministry of Education), School of Physics, Dalian University of Technology, Dalian 116024, China
| | - Zhigang Sun
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
10
|
Roncero O, Andrianarijaona V, Aguado A, Sanz-Sanz C. Vibrational effects in the quantum dynamics of the H + D 2+ charge transfer reaction. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1948125] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- O. Roncero
- IFF-CSIC, Instituto de Física Fundamental, Madrid, Spain
| | - V. Andrianarijaona
- Department of Physics, Pacific Union College, Angwin, CA, USA
- Department of Physics and Engineering, Southern Adventist University, Collegedale, TN, USA
| | - A. Aguado
- Unidad Asociada UAM-CSIC, Departamento de Química Física Aplicada, Universidad Autónoma de Madrid, Madrid, Spain
| | - C. Sanz-Sanz
- Unidad Asociada UAM-CSIC, Departamento de Química Física Aplicada, Universidad Autónoma de Madrid, Madrid, Spain
| |
Collapse
|
11
|
Ma H, Zhang C, Song Y, Ma F, Li Y. Accurate High-Level Ab Initio-Based Global Potential Energy Surface and Quantum Dynamics Calculation for the First Excited State of CH 2. J Phys Chem A 2021; 125:5490-5498. [PMID: 34137628 DOI: 10.1021/acs.jpca.1c02413] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A full three-dimensional global potential energy surface (PES), covering the whole configuration space, is reported first for the title system by fitting high-level ab initio energies at the multireference configuration interaction level with the aug-cc-pV6Z basis set. In this work, the many-body expansion method is invoked to fit the innate character of the CH2+(12A″) PES. The topographical features are examined in detail based on the new global PES and in accordance with the other calculations from the ab initio energies, which show the correct behavior at the C+(2P) + H2(X1Σg+) and CH+(a3Π) + H(2S) dissociation limits. Using a time-dependent wave packet method, we provide insights into the dynamics behavior for reaction of C+(2P) + H2(X1Σg+) → CH+(a3Π) + H(2S). The integral cross sections and reaction probabilities increase monotonically in terms of the collision energy.
Collapse
Affiliation(s)
- Hongyu Ma
- Department of Physics, Liaoning University, Shenyang 110036, China
| | - Chengyuan Zhang
- Department of Physics, Liaoning University, Shenyang 110036, China
| | - Yuzhi Song
- Collaborative Innovation Center of Light Manipulations and Applications, Shandong Normal University, Jinan 250358, China
| | - Fengcai Ma
- Department of Physics, Liaoning University, Shenyang 110036, China
| | - Yongqing Li
- Department of Physics, Liaoning University, Shenyang 110036, China
| |
Collapse
|
12
|
Sanz-Sanz C, Aguado A, Roncero O. Near-resonant effects in the quantum dynamics of the H + H 2 + → H 2 + H + charge transfer reaction and isotopic variants. J Chem Phys 2021; 154:104104. [PMID: 33722048 DOI: 10.1063/5.0044320] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The non-adiabatic quantum dynamics of the H + H2 + → H2 + H+ charge transfer reactions, and some isotopic variants, is studied with an accurate wave packet method. A recently developed 3 × 3 diabatic potential model is used, which is based on very accurate ab initio calculations and includes the long-range interactions for ground and excited states. It is found that for initial H2 +(v = 0), the quasi-degenerate H2(v' = 4) non-reactive charge transfer product is enhanced, producing an increase in the reaction probability and cross section. It becomes the dominant channel from collision energies above 0.2 eV, producing a ratio between v' = 4 and the rest of v's, which that increase up to 1 eV. The H + H2 + → H2 + + H exchange reaction channel is nearly negligible, while the reactive and non-reactive charge transfer reaction channels are of the same order, except that corresponding to H2(v' = 4), and the two charge transfer processes compete below 0.2 eV. This enhancement is expected to play an important vibrational and isotopic effect that needs to be evaluated. For the three proton case, the problem of the permutation symmetry is discussed when using reactant Jacobi coordinates.
Collapse
Affiliation(s)
- Cristina Sanz-Sanz
- Unidad Asociada UAM-CSIC, Departamento de Química Física Aplicada, Facultad de Ciencias M-14, Universidad Autónoma de Madrid, 28049 Madrid, Spain
| | - Alfredo Aguado
- Unidad Asociada UAM-CSIC, Departamento de Química Física Aplicada, Facultad de Ciencias M-14, Universidad Autónoma de Madrid, 28049 Madrid, Spain
| | - Octavio Roncero
- Instituto de Física Fundamental (IFF-CSIC), C.S.I.C., Serrano 123, 28006 Madrid, Spain
| |
Collapse
|
13
|
Zhang W, Meng Q, Gao S, Song Y. Theoretical insight into the vibrational excitation effect of the S+(4S) + H2(X1∑+g) reaction. Chem Phys Lett 2021. [DOI: 10.1016/j.cplett.2020.138257] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
14
|
Sathyamurthy N, Mahapatra S. Time-dependent quantum mechanical wave packet dynamics. Phys Chem Chem Phys 2020; 23:7586-7614. [PMID: 33306771 DOI: 10.1039/d0cp03929b] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Starting from a model study of the collinear (H, H2) exchange reaction in 1959, the time-dependent quantum mechanical wave packet (TDQMWP) method has come a long way in dealing with systems as large as Cl + CH4. The fast Fourier transform method for evaluating the second order spatial derivative of the wave function and split-operator method or Chebyshev polynomial expansion for determining the time evolution of the wave function for the system have made the approach highly accurate from a practical point of view. The TDQMWP methodology has been able to predict state-to-state differential and integral reaction cross sections accurately, in agreement with available experimental results for three dimensional (H, H2) collisions, and identify reactive scattering resonances too. It has become a practical computational tool in predicting the observables for many A + BC exchange reactions in three dimensions and a number of larger systems. It is equally amenable to determining the bound and quasi-bound states for a variety of molecular systems. Just as it is able to deal with dissociative processes (without involving basis set expansion), it is able to deal with multi-mode nonadiabatic dynamics in multiple electronic states with equal ease. We present an overview of the method and its strength and limitations, citing examples largely from our own research groups.
Collapse
|
15
|
Goswami S, Sahoo J, Paul SK, Rao TR, Mahapatra S. Effect of Reagent Vibration and Rotation on the State-to-State Dynamics of the Hydrogen Exchange Reaction, H + H 2 → H 2 + H. J Phys Chem A 2020; 124:9343-9359. [PMID: 33124827 DOI: 10.1021/acs.jpca.0c06707] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
State-to-state dynamics of the benchmark hydrogen exchange reaction H + H2 (v = 0-4, j = 0-3) → H2 (v', j') + H is investigated with the aid of the real wave packet approach of Gray and Balint-Kurti (J. Chem. Phys. 1998, 108, 950-962) and electronic ground BKMP2 potential energy surface of Boothroyd et al. (J. Chem. Phys. 1996, 104, 7139-7152). Initial state-selected and product state-resolved reaction probabilities, integral cross section, and product diatom vibrational and rotational level populations at a few collision energies are reported to elucidate the energy disposal mechanism. State-specific thermal rate constants are also calculated and compared with the available literature results. Coriolis coupling terms of the nuclear Hamiltonian are included, and calculations are parallelized over the helicity quantum number, Ω'. Attempts are made, in particular, to study the effect of reagent vibrational and rotational excitations on the dynamical attributes. It is found that the calculations become computationally expensive with reagent vibrational and rotational excitation. Reagent vibrational excitation is found to enhance the reactivity and has significant impact on the energy disposal to the vibrational and rotational degrees of freedom of the product. The interplay of reagent translational and vibrational energy on the product vibrational distribution unfolds an important aspect of the energy disposal mechanism. The effect of reagent rotation on the state-to-state dynamics is found not to be very significant, and the weak effect turns out to be specific to v'.
Collapse
Affiliation(s)
- Sugata Goswami
- School of Chemistry, University of Hyderabad, Hyderabad 500046, India
| | - Jayakrushna Sahoo
- School of Chemistry, University of Hyderabad, Hyderabad 500046, India
| | - Suranjan K Paul
- School of Chemistry, University of Hyderabad, Hyderabad 500046, India
| | - T Rajagopala Rao
- School of Chemistry, University of Hyderabad, Hyderabad 500046, India
| | - S Mahapatra
- School of Chemistry, University of Hyderabad, Hyderabad 500046, India
| |
Collapse
|
16
|
Bulut N, Roncero O, Lique F. Possible Formation and Destruction of the OD + Ions in the Interstellar Medium. J Phys Chem A 2020; 124:6552-6561. [DOI: 10.1021/acs.jpca.0c05021] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Niyazi Bulut
- Department of Physics, Firat University, 23169 Elazig̃, Turkey
| | - Octavio Roncero
- Instituto de Fı́sica Fundamental, CSIC, C/Serrano, 123, 28006 Madrid, Spain
| | - François Lique
- LOMC - UMR 6294, CNRS-Université du Havre, 25 rue Philippe Lebon, BP 1123, F-76063 Le Havre, France
| |
Collapse
|
17
|
Araújo JP, Alves MD, da Silva RS, Ballester MY. A comparative study of analytic representations of potential energy curves for O2, N2, and SO in their ground electronic states. J Mol Model 2019; 25:198. [DOI: 10.1007/s00894-019-4079-3] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/16/2019] [Accepted: 05/28/2019] [Indexed: 11/24/2022]
|
18
|
Gans B, Garcia GA, Boyé-Péronne S, Pratt ST, Guillemin JC, Aguado A, Roncero O, Loison JC. Origin band of the first photoionizing transition of hydrogen isocyanide. Phys Chem Chem Phys 2019; 21:2337-2344. [PMID: 30656348 PMCID: PMC6469576 DOI: 10.1039/c8cp07737a] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The photoelectron spectrum of the X1Σ+ → X+2Σ+ ionizing transition of hydrogen isocyanide (HNC) is measured for the first time at a fixed photon energy (13 eV). The assignment of the spectrum is supported by wave-packet calculations simulating the photoionization transition spectrum and using ab initio calculations of the potential energy surfaces for the three lowest electronic states of the cation. The photoelectron spectrum allows the retrieval of the fundamental of the CN stretching mode of the cationic ground state ([small nu, Greek, tilde]3 = 2260 ± 80 cm-1) and the adiabatic ionization energy of hydrogen isocyanide: IE(HNC) = 12.011 ± 0.010 eV, which is far below that of HCN (IE(HCN) = 13.607 eV). In light of this latter result, the thermodynamics of the HCN+/HNC+ isomers is discussed and a short summary of the values available in the literature is given.
Collapse
Affiliation(s)
- Bérenger Gans
- Institut des Sciences Moléculaires d'Orsay (ISMO), UMR 8214 CNRS, Univ. Paris-Sud, Université Paris-Saclay, F-91405 Orsay, France.
| | | | | | | | | | | | | | | |
Collapse
|
19
|
Xie C, Liu X, Guo H. State-to-state quantum dynamics of the H + LiF → Li + HF reaction on an accurate ab initio potential energy surface. Chem Phys 2018. [DOI: 10.1016/j.chemphys.2018.04.012] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
|
20
|
Jambrina PG, Menéndez M, Aoiz FJ. Angular momentum-scattering angle quantum correlation: a generalized deflection function. Chem Sci 2018; 9:4837-4850. [PMID: 29910936 PMCID: PMC5982215 DOI: 10.1039/c7sc05489k] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/28/2017] [Accepted: 04/25/2018] [Indexed: 11/23/2022] Open
Abstract
A natural generalization of the classical deflection function, the functional dependence of the deflection angle on the angular momentum (or the impact parameter), is the joint probability density function of these two quantities, revealing the correlation between them. It provides, at a glance, detailed information about the reaction mechanisms and how changes in the impact parameter affect the product angular distribution. It is also useful to predict the presence of quantum phenomena such as interference. However, the classical angular momentum-scattering angle correlation function has a limited use whenever quantum effects become important. Rigorously speaking, there is not a quantum equivalent of the classical joint distribution, as the differential cross section depends on the coherences between the different values of J caused by the cross terms in the expansion of partial waves. In this article, we present a simple method to calculate a quantum analog of this correlation, a generalized deflection function that can shed light onto the reaction mechanism using just quantum mechanical results. Our results show that there is a very good agreement between the quantum and classical correlation functions as long as quantum effects are not all relevant. When this is not the case, it will also be shown that the quantum correlation function is most useful to observe the extent of quantum effects such as interference among different reaction mechanisms.
Collapse
Affiliation(s)
- P G Jambrina
- Departamento de Química Física Aplicada , Universidad Autonoma de Madrid , 28049 , Madrid , Spain .
| | - M Menéndez
- Departamento de Química Física I , Facultad de Ciencias Químicas , Universidad Complutense de Madrid , 28040 Madrid , Spain .
| | - F J Aoiz
- Departamento de Química Física I , Facultad de Ciencias Químicas , Universidad Complutense de Madrid , 28040 Madrid , Spain .
| |
Collapse
|
21
|
Liu X, Xie C, Guo H. A new potential energy surface and state-to-state quantum dynamics of the Li + HF → H + LiF reaction. Chem Phys 2018. [DOI: 10.1016/j.chemphys.2018.01.023] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
22
|
He H, Zhu W, Su W, Dong L, Li B. Accurate Quantum Wave Packet Study of the Deep Well D + + HD Reaction: Product Ro-vibrational State-Resolved Integral and Differential Cross Sections. J Phys Chem A 2018; 122:2319-2328. [DOI: 10.1021/acs.jpca.7b08941] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Haixiang He
- School of Chemistry & Chemical Engineering, Guangxi University, Nanning 530004, People’s Republic of China
| | - Weimin Zhu
- School of Chemistry & Chemical Engineering, Guangxi University, Nanning 530004, People’s Republic of China
| | - Wenli Su
- School of Chemistry & Chemical Engineering, Guangxi University, Nanning 530004, People’s Republic of China
| | - Lihui Dong
- School of Chemistry & Chemical Engineering, Guangxi University, Nanning 530004, People’s Republic of China
| | - Bin Li
- School of Chemistry & Chemical Engineering, Guangxi University, Nanning 530004, People’s Republic of China
| |
Collapse
|
23
|
Tan RS, Zhai HC, Yan W, Gao F, Lin SY. A new ab initio potential energy surface of LiClH (1A') system and quantum dynamics calculation for Li + HCl (v = 0, j = 0-2) → LiCl + H reaction. J Chem Phys 2017; 146:164305. [PMID: 28456188 DOI: 10.1063/1.4982066] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
A new ab initio potential energy surface (PES) for the ground state of Li + HCl reactive system has been constructed by three-dimensional cubic spline interpolation of 36 654 ab initio points computed at the MRCI+Q/aug-cc-pV5Z level of theory. The title reaction is found to be exothermic by 5.63 kcal/mol (9 kcal/mol with zero point energy corrections), which is very close to the experimental data. The barrier height, which is 2.99 kcal/mol (0.93 kcal/mol for the vibrationally adiabatic barrier height), and the depth of van der Waals minimum located near the entrance channel are also in excellent agreement with the experimental findings. This study also identified two more van der Waals minima. The integral cross sections, rate constants, and their dependence on initial rotational states are calculated using an exact quantum wave packet method on the new PES. They are also in excellent agreement with the experimental measurements.
Collapse
Affiliation(s)
- Rui Shan Tan
- School of Physics, Shandong University, Jinan 250100, China
| | - Huan Chen Zhai
- School of Physics, Shandong University, Jinan 250100, China
| | - Wei Yan
- School of Physics, Shandong University, Jinan 250100, China
| | - Feng Gao
- School of Physics, Shandong University, Jinan 250100, China
| | - Shi Ying Lin
- School of Physics, Shandong University, Jinan 250100, China
| |
Collapse
|
24
|
Aguado A, Roncero O, Zanchet A, Agúndez M, Cernicharo J. The Photodissociation of HCN and HNC: Effects on the HNC/HCN Abundance Ratio in the Interstellar Medium. THE ASTROPHYSICAL JOURNAL 2017; 838:33. [PMID: 28522878 PMCID: PMC5433558 DOI: 10.3847/1538-4357/aa63ee] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/15/2023]
Abstract
The impact of the photodissociation of HCN and HNC isomers is analyzed in different astrophysical environments. For this purpose, the individual photodissociation cross section of HCN and HNC isomers have been calculated in the 7-13.6 eV photon energy range for a temperature of 10 K. These calculations are based on the ab initio calculation of three-dimensional adiabatic potential energy surfaces of the 21 lower electronic states. The cross sections are then obtained using a quantum wave packet calculation of the rotational transitions needed to simulate a rotational temperature of 10 K. The cross section calculated for HCN shows significant differences with respect to the experimental one, and this is attributed to the need of considering non-adiabatic transitions. Ratios between the photodissociation rates of HCN and HNC under different ultraviolet radiation fields have been computed by renormalizing the rates to the experimental one. It is found that HNC is photodissociated faster than HCN by a factor of 2.2, for the local interstellar radiation field, and 9.2, for the solar radiation field at 1 au. We conclude that to properly describe the HNC/HCN abundance ratio in astronomical environments illuminated by an intense ultraviolet radiation field it is necessary to use different photodissociation rates for each of the two isomers, obtained by integrating the product of the photodissociation cross sections and ultraviolet radiation field over the relevant wavelength range.
Collapse
Affiliation(s)
- Alfredo Aguado
- Departamento de Química Física Aplicada (UAM), Unidad Asociada a IFF-CSIC, Facultad de Ciencias Módulo 14, Universidad Autónoma de Madrid, 28049, Madrid, Spain
| | - Octavio Roncero
- Instituto de Física Fundamental (IFF-CSIC), C.S.I.C., Serrano 123, 28006 Madrid, Spain
| | - Alexandre Zanchet
- Instituto de Física Fundamental (IFF-CSIC), C.S.I.C., Serrano 123, 28006 Madrid, Spain
| | - Marcelino Agúndez
- Instituto de Ciencia de Materiales de Madrid, CSIC, C/ Sor Juana Inés de la Cruz 3, Cantoblanco 28049, Spain
| | - José Cernicharo
- Instituto de Ciencia de Materiales de Madrid, CSIC, C/ Sor Juana Inés de la Cruz 3, Cantoblanco 28049, Spain
| |
Collapse
|
25
|
Wu H, Duan ZX, Yin SH, Zhao GJ. State-resolved dynamics study of the H + HS reaction on the 3A′ and 3A″ states with time-dependent quantum wave packet method. J Chem Phys 2016; 145:124305. [DOI: 10.1063/1.4962543] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Hui Wu
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, University of Chinese Academy of Sciences, Dalian 116023, China
| | - Zhi-Xin Duan
- School of Science, Dalian Jiao Tong University, Dalian 116023, China
| | - Shu-Hui Yin
- Department of Physics, Dalian Maritime University, Dalian 116023, China
| | - Guang-Jiu Zhao
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, University of Chinese Academy of Sciences, Dalian 116023, China
| |
Collapse
|
26
|
Zanchet A, Roncero O, Bulut N. Quantum and quasi-classical calculations for the S⁺ + H₂(v,j) → SH⁺(v',j') + H reactive collisions. Phys Chem Chem Phys 2016; 18:11391-400. [PMID: 27055725 PMCID: PMC4894552 DOI: 10.1039/c6cp00604c] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
State-to-state cross-sections for the S(+) + H2(v,j) → SH(+)(v',j') + H endothermic reaction are obtained using quantum wave packet (WP) and quasi-classical (QCT) methods for different initial ro-vibrational H2(v,j) over a wide range of translation energies. The final state distribution as a function of the initial quantum number is obtained and discussed. Additionally, the effect of the internal excitation of H2 on the reactivity is carefully studied. It appears that energy transfer among modes is very inefficient that vibrational energy is the most favorable for the reaction, and rotational excitation significantly enhances the reactivity when vibrational energy is sufficient to reach the product. Special attention is also paid to an unusual discrepancy between classical and quantum dynamics for low rotational levels while agreement improves with rotational excitation of H2. An interesting resonant behaviour found in WP calculations is also discussed and associated with the existence of roaming classical trajectories that enhance the reactivity of the title reaction. Finally, a comparison with the experimental results of Stowe et al. for S(+) + HD and S(+) + D2 reactions exhibits a reasonably good agreement with those results.
Collapse
Affiliation(s)
- Alexandre Zanchet
- Instituto de Física Fundamental, CSIC, C/Serrano, 123, 28006 Madrid, Spain.
| | - Octavio Roncero
- Instituto de Física Fundamental, CSIC, C/Serrano, 123, 28006 Madrid, Spain.
| | - Niyazi Bulut
- Department of Physics, Firat University, Elazig, Turkey
| |
Collapse
|
27
|
Chenel A, Roncero O, Aguado A, Agúndez M, Cernicharo J. Photodissociation of HCN and HNC isomers in the 7-10 eV energy range. J Chem Phys 2016; 144:144306. [PMID: 27083720 PMCID: PMC4894478 DOI: 10.1063/1.4945389] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The ultraviolet photoabsorption spectra of the HCN and HNC isomers have been simulated in the 7-10 eV photon energy range. For this purpose, the three-dimensional adiabatic potential energy surfaces of the 7 lowest electronic states, and the corresponding transition dipole moments, have been calculated, at multireference configuration interaction level. The spectra are calculated with a quantum wave packet method on these adiabatic potential energy surfaces. The spectra for the 3 lower excited states, the dissociative electronic states, correspond essentially to predissociation peaks, most of them through tunneling on the same adiabatic state. The 3 higher electronic states are bound, hereafter electronic bound states, and their spectra consist of delta lines, in the adiabatic approximation. The radiative lifetime towards the ground electronic states of these bound states has been calculated, being longer than 10 ns in all cases, much longer that the characteristic predissociation lifetimes. The spectra of HCN is compared with the available experimental and previous theoretical simulations, while in the case of HNC there are no previous studies to our knowledge. The spectrum for HNC is considerably more intense than that of HCN in the 7-10 eV photon energy range, which points to a higher photodissociation rate for HNC, compared to HCN, in astrophysical environments illuminated by ultraviolet radiation.
Collapse
Affiliation(s)
- Aurelie Chenel
- Instituto de Física Fundamental (IFF-CSIC), C.S.I.C., Serrano 123, 28006 Madrid, Spain
| | - Octavio Roncero
- Instituto de Física Fundamental (IFF-CSIC), C.S.I.C., Serrano 123, 28006 Madrid, Spain
| | - Alfredo Aguado
- Departamento de Química Física Aplicada (UAM), Unidad Asociada a IFF-CSIC, Facultad de Ciencias Módulo 14, Universidad Autónoma de Madrid, 28049, Madrid, Spain
| | - Marcelino Agúndez
- Instituto de Ciencias de Materiales (CCMM-CSIC), C.S.I.C., Cantoblanco, 28049 Madrid, Spain
| | - José Cernicharo
- Instituto de Ciencias de Materiales (CCMM-CSIC), C.S.I.C., Cantoblanco, 28049 Madrid, Spain
| |
Collapse
|
28
|
Bulut N, Lique F, Roncero O. Exchange and Inelastic OH(+) + H Collisions on the Doublet and Quartet Electronic States. J Phys Chem A 2015. [PMID: 26203890 DOI: 10.1021/acs.jpca.5b05246] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The exchange and inelastic state-to-state cross sections for the OH(+) + H collisions are computed from wave packet calculations using the doublet and quartet ground electronic potential energy surface (PES) correlating to the open shell reactants, for collision energies in the range of 1 meV to 0.7 eV. The doublet PES presents a deep insertion well, of ≈6 eV, but the exchange reaction has a rather low probability, showing that the mechanism is not statistical. This well is also responsible of a rather high rotational energy transfer, which makes the rigid-rotor approach overestimate the cross section for low Δj transitions and for high collisonal energies. The quartet PES, with a much shallower well, also presents a low exchange reaction cross section, but the inelastic state-to-state cross sections are very well reproduced by rigid-rotor calculations. When the electronic partition is used to obtain the total state-to-state cross section, the contribution of the doublet state becomes small, and the resulting total cross sections become close to those obtained for the quartet state. Thus, the total (quartet and doublet) cross sections for this open shell system can be reproduced rather satisfactorily by those obtained with the rigid-rotor approximation on the quartet state. Finally, we compare the new OH(+)-H cross sections with OH(+)-He ones recently computed. We found significant differences, especially for transitions with large Δj showing that specific OH(+)-H calculations had to be performed to accurately analyze the OH(+) emission from interstellar molecular clouds.
Collapse
Affiliation(s)
- Niyazi Bulut
- Firat University , Department of Physics, 23169 Elazig̃, Turkey
| | - François Lique
- LOMC - UMR 6294, CNRS-Université du Havre , 25 rue Philippe Lebon, BP 1123-76 063, Le Havre, France
| | - Octavio Roncero
- Instituto de Física Fundamental, CSIC , C/Serrano, 123, 28006 Madrid, Spain
| |
Collapse
|
29
|
Lara M, Jambrina PG, Aoiz FJ, Launay JM. Cold and ultracold dynamics of the barrierless D+ + H2 reaction: Quantum reactive calculations for ∼R−4 long range interaction potentials. J Chem Phys 2015; 143:204305. [DOI: 10.1063/1.4936144] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Manuel Lara
- Departamento de Química Física Aplicada, Facultad de Ciencias, Universidad Autónoma de Madrid, 28049 Madrid, Spain
| | - P. G. Jambrina
- Departamento de Química Física, Facultad de Química, Universidad Complutense, 28040 Madrid, Spain
| | - F. J. Aoiz
- Departamento de Química Física, Facultad de Química, Universidad Complutense, 28040 Madrid, Spain
| | - J.-M. Launay
- Institut de Physique de Rennes, UMR CNRS 6251, Université de Rennes I, F-35042 Rennes, France
| |
Collapse
|
30
|
Li WT, Chen MD, Sun ZG. Quantum Dynamics of Li+HF/DF Reaction Investigated by a State-to-State Time-dependent Wave Packet Approach. CHINESE J CHEM PHYS 2015. [DOI: 10.1063/1674-0068/28/cjcp1507151] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
|
31
|
Bulut N, Kłos J, Roncero O. Quantum mechanical calculations of state-to-state cross sections and rate constants for the F + DCl → Cl + DF reaction. J Chem Phys 2015; 142:214310. [PMID: 26049499 DOI: 10.1063/1.4922110] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present accurate state-to-state quantum wave packet calculations of integral cross sections and rate constants for the title reaction. Calculations are carried out on the best available ground 1(2)A' global adiabatic potential energy surface of Deskevich et al. [J. Chem. Phys. 124, 224303 (2006)]. Converged state-to-state reaction cross sections have been calculated for collision energies up to 0.5 eV and different initial rotational and vibrational excitations, DCl(v = 0, j = 0 - 1; v = 1, j = 0). Also, initial-state resolved rate constants of the title reaction have been calculated in a temperature range of 100-400 K. It is found that the initial rotational excitation of the DCl molecule does not enhance reactivity, in contract to the reaction with the isotopologue HCl in which initial rotational excitation produces an important enhancement. These differences between the isotopologue reactions are analyzed in detail and attributed to the presence of resonances for HCl(v = 0, j), absent in the case of DCl(v = 0, j). For vibrational excited DCl(v = 1, j), however, the reaction cross section increases noticeably, what is also explained by another resonance.
Collapse
Affiliation(s)
- Niyazi Bulut
- Department of Physics, Firat University, 23169 Elazig˜, Turkey
| | - Jacek Kłos
- Department of Chemistry and Biochemistry, University of Maryland, College Park, Maryland 20742-2021, USA
| | - Octavio Roncero
- Instituto de Física Fundamental (IFF-CSIC), C.S.I.C., Serrano 123, 28006 Madrid, Spain
| |
Collapse
|
32
|
Koner D, Barrios L, González-Lezana T, Panda AN. Wave packet and statistical quantum calculations for the He + NeH⁺ → HeH⁺ + Ne reaction on the ground electronic state. J Chem Phys 2015; 141:114302. [PMID: 25240353 DOI: 10.1063/1.4895567] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A real wave packet based time-dependent method and a statistical quantum method have been used to study the He + NeH(+) (v, j) reaction with the reactant in various ro-vibrational states, on a recently calculated ab initio ground state potential energy surface. Both the wave packet and statistical quantum calculations were carried out within the centrifugal sudden approximation as well as using the exact Hamiltonian. Quantum reaction probabilities exhibit dense oscillatory pattern for smaller total angular momentum values, which is a signature of resonances in a complex forming mechanism for the title reaction. Significant differences, found between exact and approximate quantum reaction cross sections, highlight the importance of inclusion of Coriolis coupling in the calculations. Statistical results are in fairly good agreement with the exact quantum results, for ground ro-vibrational states of the reactant. Vibrational excitation greatly enhances the reaction cross sections, whereas rotational excitation has relatively small effect on the reaction. The nature of the reaction cross section curves is dependent on the initial vibrational state of the reactant and is typical of a late barrier type potential energy profile.
Collapse
Affiliation(s)
- Debasish Koner
- Department of Chemistry, Indian Institute of Technology Guwahati, Guwahati 781039, India
| | - Lizandra Barrios
- Instituto de Física Fundamental, C.S.I.C., Serrano 123, Madrid 28006, Spain
| | | | - Aditya N Panda
- Department of Chemistry, Indian Institute of Technology Guwahati, Guwahati 781039, India
| |
Collapse
|
33
|
Bulut N, Castillo J, Jambrina PG, Kłos J, Roncero O, Aoiz FJ, Bañares L. Accurate Time-Dependent Wave Packet Calculations for the O+ + H2 → OH+ + H Ion–Molecule Reaction. J Phys Chem A 2015; 119:11951-62. [DOI: 10.1021/acs.jpca.5b00815] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- N. Bulut
- Department of Physics, Firat University, 23169 Elazig̃, Turkey
| | - J.F. Castillo
- Departamento
de Química Física I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid (Unidad Asociada I+D+i CSIC), 28040 Madrid, Spain
| | - P. G. Jambrina
- Departamento
de Química Física I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid (Unidad Asociada I+D+i CSIC), 28040 Madrid, Spain
| | - J. Kłos
- Department of Chemistry
and Biochemistry, University of Maryland, College Park, Maryland 20742-2021, United States
| | - O. Roncero
- Departamento de Física Atómica, Molecular y de Agregados,
Instituto de Física Fundamental, CSIC, C/Serrano, 123, 28006 Madrid, Spain
| | - F. J. Aoiz
- Departamento
de Química Física I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid (Unidad Asociada I+D+i CSIC), 28040 Madrid, Spain
| | - L. Bañares
- Departamento
de Química Física I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid (Unidad Asociada I+D+i CSIC), 28040 Madrid, Spain
| |
Collapse
|
34
|
Sahoo T, Ghosh S, Adhikari S, Sharma R, Varandas AJC. Low-temperature D+ + H2 reaction: A time-dependent coupled wave-packet study in hyperspherical coordinates. J Chem Phys 2015; 142:024304. [DOI: 10.1063/1.4905379] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Tapas Sahoo
- Department of Physical Chemistry, Indian Association for the Cultivation of Science, 2A & 2B Raja S. C. Mullick Road, Jadavpur, Kolkata 700032, West Bengal, India
| | - Sandip Ghosh
- Department of Physical Chemistry, Indian Association for the Cultivation of Science, 2A & 2B Raja S. C. Mullick Road, Jadavpur, Kolkata 700032, West Bengal, India
| | - Satrajit Adhikari
- Department of Physical Chemistry, Indian Association for the Cultivation of Science, 2A & 2B Raja S. C. Mullick Road, Jadavpur, Kolkata 700032, West Bengal, India
| | - Rahul Sharma
- Departamento de Química, and Centro de Química, Universidade de Coimbra, 3004-535 Coimbra, Portugal
| | - António J. C. Varandas
- Departamento de Química, and Centro de Química, Universidade de Coimbra, 3004-535 Coimbra, Portugal
| |
Collapse
|
35
|
|
36
|
|
37
|
Yuan J, Cheng D, Chen M. Time-dependent wave packet and quasiclassical trajectory studies of the Au + HD reaction: competition between the reactive channels. RSC Adv 2014. [DOI: 10.1039/c4ra06297c] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
|
38
|
Goswami S, Rao TR, Mahapatra S, Bussery-Honvault B, Honvault P. Time-Dependent Quantum Wave Packet Dynamics of S + OH Reaction on Its Electronic Ground State. J Phys Chem A 2014; 118:5915-26. [DOI: 10.1021/jp504757g] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Sugata Goswami
- School
of Chemistry, University of Hyderabad, Hyderabad 500046, India
| | - T. Rajagopala Rao
- Laboratoire
Interdisciplinaire Carnot de Bourgogne, UMR CNRS 6303, Université Bourgogne, 21078 Dijon Cedex, France
| | - S. Mahapatra
- School
of Chemistry, University of Hyderabad, Hyderabad 500046, India
| | - B. Bussery-Honvault
- Laboratoire
Interdisciplinaire Carnot de Bourgogne, UMR CNRS 6303, Université Bourgogne, 21078 Dijon Cedex, France
| | - P. Honvault
- Laboratoire
Interdisciplinaire Carnot de Bourgogne, UMR CNRS 6303, Université Bourgogne, 21078 Dijon Cedex, France
- UFR
ST, Université de Franche-Comté, 25030 Besançon
Cedex, France
| |
Collapse
|
39
|
Yuan J, Cheng D, Sun Z, Chen M. Time-dependent quantum wave packet and quasiclassical trajectory studies of the Au(2S) + H2(X1∑+g) → AuH(X1∑+g) + H(2S) reaction. Mol Phys 2014. [DOI: 10.1080/00268976.2014.920113] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
40
|
González-Lezana T, Honvault P, Scribano Y. Dynamics of the D(+) + H2 → HD + H(+) reaction at the low energy regime by means of a statistical quantum method. J Chem Phys 2014; 139:054301. [PMID: 23927256 DOI: 10.1063/1.4816638] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The D(+) +H2(v = 0, j = 0, 1) → HD+H(+) reaction has been investigated at the low energy regime by means of a statistical quantum mechanical (SQM) method. Reaction probabilities and integral cross sections (ICSs) between a collisional energy of 10(-4) eV and 0.1 eV have been calculated and compared with previously reported results of a time independent quantum mechanical (TIQM) approach. The TIQM results exhibit a dense profile with numerous narrow resonances down to Ec ~ 10(-2) eV and for the case of H2(v = 0, j = 0) a prominent peak is found at ~2.5 × 10(-4) eV. The analysis at the state-to-state level reveals that this feature is originated in those processes which yield the formation of rotationally excited HD(v' = 0, j' > 0). The statistical predictions reproduce reasonably well the overall behaviour of the TIQM ICSs at the larger energy range (Ec ≥ 10(-3) eV). Thermal rate constants are in qualitative agreement for the whole range of temperatures investigated in this work, 10-100 K, although the SQM values remain above the TIQM results for both initial H2 rotational states, j = 0 and 1. The enlargement of the asymptotic region for the statistical approach is crucial for a proper description at low energies. In particular, we find that the SQM method leads to rate coefficients in terms of the energy in perfect agreement with previously reported measurements if the maximum distance at which the calculation is performed increases noticeably with respect to the value employed to reproduce the TIQM results.
Collapse
|
41
|
Xie TX, Zhang YY, Shi Y, Jin MX. Theoretical Study of Reagent Rotational Excitation Effect on the Stereodynamics of H + LiF→HF+Li Reaction. CHINESE J CHEM PHYS 2014. [DOI: 10.1063/1674-0068/27/01/39-44] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
|
42
|
|
43
|
Ge M, Yang H, Zheng Y. The dynamical study of O(1D) + HCl(v = 0, j = 0) reaction at hyperthermal collision energies. Chem Cent J 2013; 7:177. [PMID: 24237765 PMCID: PMC4176982 DOI: 10.1186/1752-153x-7-177] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2013] [Accepted: 11/07/2013] [Indexed: 11/21/2022] Open
Abstract
Backgrounds The quasi-classical trajectory calculations for O(1D) + HCl → OH + Cl (R1) and O(1D) + HCl → ClO + H (R2) reactions have been performed at hyperthermal collision energies (60.0, 90.0, and 120.0 kal/mol) on the 1A' state. Reaction probabilities and integral cross sections are calculated. The product rotational distributions for the two channels, and the product rotational alignment parameters are investigated. Also, the alignment and the orientation of the products have been predicted through the angular distribution functions (concerning the initial/final velocity vector, and the product rotational angular momentum vector). To have a deeper understanding of the natures of the vector correlation between reagent and product relative velocities, a natural generalization of the differential cross section __PDDCS00, is calculated. Results The OH + Cl channel is the main product channel and is observed to have essentially isotropic rotational distributions. The ClO + H channel is found to be clearly rotationally polarized. Conclusions The dynamical, especially the stereodynamical characters are quite different for the two channels of the title reaction. Most reactions occur directly, except for R2 reaction at the collision energies of 60.0 and 120.0 kcal/mol. The alignment and orientation effects are weak/strong for R1/R2 reaction. The well structure on the potential energy surface and hyperthermal collision energies might result in the dynamical effects.
Collapse
Affiliation(s)
- Meihua Ge
- School of Physics, Shandong University, Jinan 250100, China.
| | | | | |
Collapse
|
44
|
Krasilnikov MB, Popov RS, Roncero O, De Fazio D, Cavalli S, Aquilanti V, Vasyutinskii OS. Polarization of molecular angular momentum in the chemical reactions Li + HF and F + HD. J Chem Phys 2013; 138:244302. [DOI: 10.1063/1.4809992] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
|
45
|
CHENG JIE, YUE XIANFANG. QUASICLASSICAL TRAJECTORY STUDY OF STEREODYNAMICS FOR THE REACTIONS Li+ HF/DF/TF. JOURNAL OF THEORETICAL & COMPUTATIONAL CHEMISTRY 2013. [DOI: 10.1142/s0219633613500089] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Stereodynamics of the reaction Li + HF (v = 0,j = 0) → LiF + H and its isotopic variants on the ground electronic state (12A′) potential energy surface (PES) are studied by employing the quasiclassical trajectory (QCT) method. At a collision energy of 2.2 kcal/mol, product rotational angular momentum distributions, P(θr) and P(ϕr), are calculated in the center-of-mass (CM) frame. The results demonstrate that the product rotational angular momentum j′ is not only aligned along the direction perpendicular to the reagent relative velocity vector k, but also oriented along the negative y-axis. The four generalized polarization-dependent differential cross sections (PDDCSs) are also computed. The PDDCS00 distribution shows a sideways scattering for the reaction Li + HF and a strongly backward scattering for the reaction Li + DF . However, it displays both the forward and backward scatterings for the reaction Li + TF . These features demonstrate that the Li + HF and Li + DF reactions proceed predominantly through the direct reaction mechanism. However, the Li + TF reaction undergoes both the direct and indirect reaction mechanisms. The PDDCS21- distribution indicates that the product angular distributions are anisotropic.
Collapse
Affiliation(s)
- JIE CHENG
- Department of Physics and Information Engineering, Jining University, Qufu 273155, P. R. China
| | - XIAN-FANG YUE
- Department of Physics and Information Engineering, Jining University, Qufu 273155, P. R. China
| |
Collapse
|
46
|
Krasil’nikov MB, Vasyutinskii OS, Roncero O. Angular momentum polarization of the molecules in the Li + HF reaction. RUSSIAN JOURNAL OF PHYSICAL CHEMISTRY B 2013. [DOI: 10.1134/s199079311301003x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
47
|
Jorfi M, González-Lezana T, Zanchet A, Honvault P, Bussery-Honvault B. Quasiclassical Trajectory and Statistical Quantum Calculations for the C + OH → CO + H Reaction on the First Excited 12A″ Potential Energy Surface. J Phys Chem A 2013; 117:1872-9. [DOI: 10.1021/jp309764g] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- M. Jorfi
- Institut de Chimie des Milieux et des Matériaux de Poitiers,
UMR CNRS 6503, Université de Poitiers, 86022 Poitiers Cedex, France
| | - T. González-Lezana
- Instituto de Física Fundamental, C.S.I.C., Serrano 123, Madrid 28006, Spain
| | - A. Zanchet
- Instituto de Física Fundamental, C.S.I.C., Serrano 123, Madrid 28006, Spain
| | - P. Honvault
- Laboratoire Interdisciplinaire Carnot de Bourgogne, UMR CNRS 6303, Université de Bourgogne, 21078 Dijon Cedex,
France
- UFR Sciences et
Techniques, Université de Franche-Comté, 25030 Besançon Cedex, France
| | - B. Bussery-Honvault
- Laboratoire Interdisciplinaire Carnot de Bourgogne, UMR CNRS 6303, Université de Bourgogne, 21078 Dijon Cedex,
France
| |
Collapse
|
48
|
Xu W, Zhang P. Accurate study on the quantum dynamics of the He + HeH(+) (X1Σ+) reaction on a new ab initio potential energy surface for the lowest 1(1)A' electronic singlet state. J Phys Chem A 2013; 117:1406-12. [PMID: 23347266 DOI: 10.1021/jp312084r] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A time-dependent quantum wave packet method is used to investigate the dynamics of the He + HeH(+)(X(1)Σ(+)) reaction based on a new potential energy surface [Liang et al., J. Chem. Phys.2012, 136, 094307]. The coupled channel (CC) and centrifugal-sudden (CS) reaction probabilities as well as the total integral cross sections are calculated. A comparison of the results with and without Coriolis coupling revealed that the number of K states N(K) (K is the projection of the total angular momentum J on the body-fixed z axis) significantly influences the reaction threshold. The effective potential energy profiles of each N(K) for the He + HeH(+) reaction in a collinear geometry indicate that the barrier height gradually decreased with increased N(K). The calculated time evolution of CC and CS probability density distribution over the collision energy of 0.27-0.36 eV at total angular momentum J = 50 clearly suggests a lower reaction threshold of CC probabilities. The CC cross sections are larger than the CS results within the entire energy range, demonstrating that the Coriolis coupling effect can effectively promote the He + HeH(+) reaction.
Collapse
Affiliation(s)
- Wenwu Xu
- State key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Science, Dalian 116023, PR China
| | | |
Collapse
|
49
|
Gómez-Carrasco S, Bulut N, Bañares L, Roncero O. Wave packet calculations on nonadiabatic effects for the O(3P)+HF(1Σ+) reaction under hyperthermal conditions. J Chem Phys 2012; 137:114309. [PMID: 22998264 DOI: 10.1063/1.4753811] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We present wave packet calculations of total and state-to-state reaction probabilities and integral cross sections for the nonadiabatic dynamics of the O((3)P)+HF → F((2)P)+OH((2)Π) reaction at hyperthermal collision energies ranging from 1.2 to 2.4 eV. The validity of the centrifugal sudden approximation is discussed for the title reaction and a comprehensive investigation of the influence of nonadiabatic effects on the dynamics of this reactive system at high (hyperthermal) collision energies is presented. In general, nonadiabatic effects are negligible for averaged observables, such as total reaction probabilities and integral cross sections, but they are clearly observed in detailed observables such as rotationally state-resolved reaction probabilities. A critical discussion of nonadiabatic effects on the dynamics of the title reaction is carried out by comparing with the reverse reaction and the characteristics of the adiabatic and diabatic potential energy surfaces involved.
Collapse
Affiliation(s)
- S Gómez-Carrasco
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad de Salamanca, Plaza de los Caídos s/n, 37008 Salamanca, Spain
| | | | | | | |
Collapse
|
50
|
Yue XF, Wang MS. Isotope effects on product polarization and reaction mechanism in the Li+HF(v=0,j=0)→LiF+H reaction. Chem Phys 2012. [DOI: 10.1016/j.chemphys.2012.07.005] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|