1
|
Yu Q, Yang J, Zhang HR, Liang PY, Gao G, Yuan Y, Dou W, Zhou PP. Investigations of the reaction mechanism of sodium with hydrogen fluoride to form sodium fluoride and the adsorption of hydrogen fluoride on sodium fluoride monomer and tetramer. J Mol Model 2024; 30:26. [PMID: 38191945 DOI: 10.1007/s00894-023-05821-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2023] [Accepted: 12/27/2023] [Indexed: 01/10/2024]
Abstract
CONTEXT The reaction between Na and HF is a typical harpooning reaction which is of great interest due to its significance in understanding the elementary chemical reaction kinetics. This work aims to investigate the detailed reaction mechanisms of sodium with hydrogen fluoride and the adsorption of HF on the resultant NaF as well as the (NaF)4 tetramer. The results suggest that the reaction between Na and HF leads to the formation of sodium fluoride salt NaF and hydrogen gas. Na interacts with HF to form a complex HF···Na, and then the approaching of F atom of HF to Na results in a transition state H···F···Na. Accompanied by the broken of H-F bond, the bond forms between F and Na atoms as NaF, then the product NaF is yielded due to the removal of H atom. The resultant NaF can further form (NaF)4 tetramer. The interaction of NaF with HF leads to the complex NaF···HF; the form I as well as II of (NaF)4 can interact with HF to produce two complexes (i.e., (NaF)4(I-1)···HF, (NaF)4(I-2)···HF, (NaF)4(II-1)···HF and (NaF)4(II-2)···HF), but the form III of (NaF)4 can interact with HF to produce only one complex (NaF)4(III)···HF. These complexes were explored in terms of noncovalent interaction (NCI) and quantum theory of atoms in molecules (QTAIM) analyses. NCI analyses confirm the existences of attractive interactions in the complexes HF···Na, NaF···HF, (NaF)4(I-1)···HF, (NaF)4(I-2)···HF, (NaF)4(II-1)···HF and (NaF)4(II-2)···HF, and (NaF)4(III)···HF. QTAIM analyses suggest that the F···Na interaction forms in the HF···Na complex while the F···H hydrogen bonds form in NaF···HF, (NaF)4(I-1)···HF, (NaF)4(I-2)···HF, (NaF)4(II-1)···HF and (NaF)4(II-2)···HF, and (NaF)4(III)···HF complexes. Natural bond orbital (NBO) analyses were also applied to analyze the intermolecular donor-acceptor orbital interactions in these complexes. These results would provide valuable insight into the chemical reaction of Na and HF and the adsorption interaction between sodium fluoride salt and HF. METHODS The calculations were carried out at the M06-L/6-311++G(2d,2p) level of theory which were performed using the Gaussian16 program. Intrinsic reaction coordinate (IRC) calculations were carried out at the same level of theory to confirm that the obtained transition state was true. The molecular surface electrostatic potential (MSEP) was employed to understand how the complex forms. Quantum theory of atoms in molecules (QTAIM) and noncovalent interaction (NCI) analysis was used to know the topology parameters at bond critical points (BCPs) and intermolecular interactions in the complex and intermediate. The topology parameters and the BCP plots were obtained by the Multiwfn software.
Collapse
Affiliation(s)
- Qinwei Yu
- State Key Laboratory of Fluorine & Nitrogen Chemicals, Xi'an Modern Chemistry Research Institute, Xi'an, 710065, People's Republic of China.
| | - Jianming Yang
- State Key Laboratory of Fluorine & Nitrogen Chemicals, Xi'an Modern Chemistry Research Institute, Xi'an, 710065, People's Republic of China.
| | - Hai-Rong Zhang
- Key Laboratory of Advanced Catalysis of Gansu Province, Advanced Catalysis Center, College of Chemistry and Chemical Engineering, Lanzhou University, 222 South Tianshui Road, 730000, Lanzhou, People's Republic of China
| | - Peng-Yu Liang
- Key Laboratory of Advanced Catalysis of Gansu Province, Advanced Catalysis Center, College of Chemistry and Chemical Engineering, Lanzhou University, 222 South Tianshui Road, 730000, Lanzhou, People's Republic of China
| | - Ge Gao
- Key Laboratory of Advanced Catalysis of Gansu Province, Advanced Catalysis Center, College of Chemistry and Chemical Engineering, Lanzhou University, 222 South Tianshui Road, 730000, Lanzhou, People's Republic of China
| | - Yongna Yuan
- School of Information Science & Engineering, Lanzhou University, Lanzhou, Gansu, 730000, People's Republic of China
| | - Wei Dou
- Key Laboratory of Advanced Catalysis of Gansu Province, Advanced Catalysis Center, College of Chemistry and Chemical Engineering, Lanzhou University, 222 South Tianshui Road, 730000, Lanzhou, People's Republic of China
| | - Pan-Pan Zhou
- Key Laboratory of Advanced Catalysis of Gansu Province, Advanced Catalysis Center, College of Chemistry and Chemical Engineering, Lanzhou University, 222 South Tianshui Road, 730000, Lanzhou, People's Republic of China.
| |
Collapse
|
2
|
Yan W, Tan RS, Lin SY. New ab initio potential energy surface of NaFH (1 A') system and quantum dynamics studies for the Na + HF ( v, j) → NaF + H reaction. RSC Adv 2023; 13:15506-15513. [PMID: 37223419 PMCID: PMC10201549 DOI: 10.1039/d3ra01885g] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2023] [Accepted: 05/15/2023] [Indexed: 05/25/2023] Open
Abstract
A global potential energy surface (PES) for the electronic ground state of the Na + HF reactive system is constructed by three-dimensional cubic spline interpolation of 37 000 ab initio points obtained using the multireference configuration interaction method including the Davidson's correction (MRCI + Q) with auc-cc-pV5Z basis set. The endoergicity, well depth and properties of the separated diatomic molecules are in good agreement with experimental estimations. Quantum dynamics calculations have been performed and compared with those of the previous MRCI PES as well as experimental values. The better agreement between theory and experiment indicates the accuracy of the new PES.
Collapse
Affiliation(s)
- Wei Yan
- School of Science, Shandong Jianzhu University Jinan 250101 China
| | - Rui Shan Tan
- School of Science, Shandong Jianzhu University Jinan 250101 China
| | - Shi Ying Lin
- School of Physics, Shandong University Jinan 250100 China
| |
Collapse
|
3
|
Accurate Adiabatic and Diabatic Potential Energy Surfaces for the Reaction of He + H 2. BIOMED RESEARCH INTERNATIONAL 2022; 2022:7552881. [PMID: 35757471 PMCID: PMC9225863 DOI: 10.1155/2022/7552881] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/10/2022] [Accepted: 05/12/2022] [Indexed: 11/17/2022]
Abstract
The accurate adiabatic and diabatic potential energy surfaces, which are for the two lowest states of He + H2, are presented in this study. The Molpro 2012 software package is used, and the large basis sets (aug-cc-pV5Z) are selected. The high-level MCSCF/MRCI method is employed to calculate the adiabatic potential energy points of the title reaction system. The triatomic reaction system is described by Jacobi coordinates, and the adiabatic potential energy surfaces are fitted accurately using the B-spline method. The equilibrium structures and electronic energies for the H2 are provided, and the corresponding different levels of vibrational energies of the ground state are deduced. To better express the diabatic process of the whole reaction, avoid crossing points being calculated and conical intersection also being optimized. Meanwhile, the diabatic potential energy surfaces of the reaction process are constructed. This study will be helpful for the analysis of histopathology and for the study in biological and medical mechanisms.
Collapse
|
4
|
The Reaction Mechanism Study for the F 3 System. BIOMED RESEARCH INTERNATIONAL 2022; 2022:7088063. [PMID: 35528182 PMCID: PMC9071887 DOI: 10.1155/2022/7088063] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/02/2022] [Accepted: 03/25/2022] [Indexed: 11/17/2022]
Abstract
In order to study the F3 system, an accurate global adiabatic potential energy surface is reduced in the present work. The high-level ab initio (MCSCF/MRCI level) methods with big basis set aVQZ are used to calculate 27690 potential energy points in the MOLPRO quantum chemistry package using the Jacobi coordinate. Meanwhile, the B-spline fit method is used to reduce the global potential energy surface in this present work. The shallow well complexes are found in the present work when the angles θ = 30°, 60°, and 90°. Analysing the global potential energy surfaces can get the conclusion that reactants should overcome at least 0.894 eV energy to cross the transition state and reach products. This study will be helpful for the analysis in histopathology and for the study of biological and medical mechanisms.
Collapse
|
5
|
Yin R, Gao N, Cao J, Li Y, Wang D, Huang X. Global accurate diabatic potential surfaces for the reaction H + Li 2. RSC Adv 2020; 10:39226-39240. [PMID: 35518445 PMCID: PMC9057374 DOI: 10.1039/d0ra05777k] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2020] [Accepted: 09/21/2020] [Indexed: 12/25/2022] Open
Abstract
The adiabatic potential energies for the lowest three states of a Li2H system are calculated with a high level ab initio method (MCSCF/MRCI) with a large basis set (aV5Z). The accurate three dimensional B-spline fitting method is used to map the global adiabatic potential energy surfaces, using the existing adiabatic potential energies, for the lowest two adiabatic states of the title reaction system. The different vibrational states and corresponding energies are studied for the diatomic molecule of reactant and products. In order to clearly understand the nonadiabatic process, the avoided crossing area and conical intersection are carefully studied. For further study of the nonadiabatic dynamic reaction, the diabatic potential energy surfaces are deduced in the present work.
Collapse
Affiliation(s)
- Ruilin Yin
- Institute of Theoretical Chemistry, Jilin University Changchun People's Republic of China
| | - Nan Gao
- Department of Thoracic Surgery, China-Japan Union Hospital of Jilin University Changchun People's Republic of China
| | - Jing Cao
- Institute of Theoretical Chemistry, Jilin University Changchun People's Republic of China
| | - Yanchun Li
- Institute of Theoretical Chemistry, Jilin University Changchun People's Republic of China
| | - Dequan Wang
- Institute of Theoretical Chemistry, Jilin University Changchun People's Republic of China
| | - Xuri Huang
- Institute of Theoretical Chemistry, Jilin University Changchun People's Republic of China
| |
Collapse
|