Abstract
Transiently trapped quantum states along the reaction coordinate in the transition-state region of a chemical reaction are normally called Feshbach resonances or dynamical resonances. Feshbach resonances trapped in the HF–OH interaction well have been discovered in an earlier photodetchment study of FH2O−; however, it is not clear whether these resonances are accessible by the F + H2O reaction. Here we report an accurate state-to-state quantum dynamics study of the F + H2O → HF + OH reaction on an accurate newly constructed potential energy surface. Pronounced oscillatory structures are observed in the total reaction probabilities, in particular at collision energies below 0.2 eV. Detailed analysis reveals that these oscillating structures originate from the Feshbach resonance states trapped in the hydrogen bond well on the HF(v′ = 2)-OH vibrationally adiabatic potentials, producing mainly HF(v′ = 1) product. Therefore, the resonances observed in the photodetchment study of FH2O− are accessible to the reaction.
Similar content being viewed by others
Introduction
Feshbach resonances of a chemical reaction are transiently trapped quantum states along the reaction coordinate in the transition-state region. They not only have a profound influence on both the rate and product distribution, but also provide sensitive probes to the potential energy surface (PES) of the reaction. Over the past decades, great efforts have been devoted to detecting resonances in chemical reactions and to studying their structures and dynamics1,2,3. Thanks to the close interaction between theory and experiment, an accurate physical picture of dynamical resonances in the F + H2 (HD) → HF + H(D) reaction has been well established4,5,6,7,8,9,10,11,12. Detailed quantum dynamics studies on increasingly more accurate PESs uncovered that the interaction between H(D) and HF in the post-barrier region softens the HF bond and creates a peculiar well on the HF(v′ = 3)-H(D) vibrationally adiabatic potential (VAP), which supports resonance states producing mainly HF(v′ = 2) product8,9.
Similar to the F + H2 system, the F + H2O → HF + OH reaction is highly exothermic (ΔH0 = −17.6 kcal/mol) and has a low and “early” barrier on the ground electronic state. There is a relatively deep van der Waals (vdw) well (~3 kcal/mol) in the entrance channel and a hydrogen bond (HB) well (~6 kcal/mol) in the exit channel. The reaction is of great importance in atmospheric13 and astrochemistry14, and has attracted a large number of experimental and theoretical studies. The reaction dynamics was first investigated experimentally by Setser and colleagues15,16, who found the HF product is vibrationally excited and rotationally cold. After that, Nesbitt and colleagues17,18 observed a highly nascent vibrational population inversion for the HF(v = 1) product. Theoretically, as the number of the electrons is relatively small, this system represents an attractive prototype for the electronic structure and reaction dynamics calculations. Guo and colleagues19 developed a full-dimensional global PES for the F + H2O system using a multireference configuration interaction (MRCI) method and applied a standard quasi-classical trajectory (QCT) study. In agreement with experiment, the product state distribution is dominated by HF vibrational excited states. The subsequent QCT and quantum mechanical (QM) calculations found that the pre-reaction vdw well enhance reactivity at low energies due to stereodynamics20 and the vibrational excitation of H2O is much more efficient at promoting the reactivity than the translational energy21. On the other hand, Guo and colleagues22,23,24 also found the coupled-cluster method with higher excitations than singles and doubles gives a better description than the MRCI method for the barrier height. They scaled the original MRCI PES with an external correlation scaling method and provided both kinetic and dynamical data for the reaction on the scaled-MRCI PESs22,23,24.
In 2014, Otto et al. carried out a joint experimental and theoretical study on the FH2O− photodetachment dynamics. Both the measured and calculated photoelectron spectra of the anion were found to possess both broad and sharp features25,26,27,28. The former were attributed to short-lived reactive Feshbach resonances and the latter to long-lived nonreactive Feshbach resonances, both supported by product well28. At about the same time, Zhao and Guo29 reported a state-to-state quantum dynamics study on the reaction by using transition-state wave packet method for collision energy up to 2 kcal/mol on the scaled-MRCI PES. Strong oscillatory structures were observed for the J = 0 quantum reaction probabilities as a function of collision energy. However, they attributed these oscillations to tunneling resonances supported by the pre-reaction vdw well29. Therefore, the following questions clearly emerge: Can the reactive Feshbach resonances observed in the photodetachmment study be observed in the F + H2O reaction? And how these resonances affect the reaction?
Here we report an accurate state-to-state quantum scattering study of the title reaction on an accurate newly constructed neural network (NN) PES by using the time-dependent wave packet (TDWP) method30,31. The total reaction probabilities for collision energy below 0.2 eV exhibit many pronounced oscillatory structures especially. Analysis of scattering wave functions and VAPs reveal that these oscillating structures originate from the Feshbach resonance states trapped in the HF(v′ = 2)-OH VAP well, producing mainly the HF(v′ = 1) product. This suggests the Feshbach resonances observed in the photodetachment study can be accessed by the neutral reaction.
Results
Total reaction probabilities and the product-state distributions
The converged total reaction probabilities for the total angular momentum J = 0 in Fig. 1a exhibit many sharp and narrow oscillatory structures especially for collision energy below 0.2 eV. Overall, the reaction probabilities on the present PES resemble those obtained on the scaled-MRCI permutation invariant polynomial NN PES with spin–orbit (SO) corrections24 rather well, except the oscillatory structures and a shift of energy by ~0.015 eV apparently due to different barrier heights for these two PES, but larger than those obtained on the original MRCI PES20 (Supplementary Fig. 3). Figure 1b shows the product HF vibrational state-resolved reaction probabilities as a function of collision energy. The HF(v = 1) state has the largest population in the entire energy region considered here, indicating an inverted HF vibrational state distribution. Before the opening of the HF(v = 2) channel, the HF(v = 0) and HF(v = 1) probabilities have rich oscillating structures. The HF(v = 2) channel opens at the collision energy of 0.191 eV and also has some small oscillating structures in probability. Meanwhile the HF(v = 0) and HF(v = 1) probabilities become rather smooth at high collision energies.
To compare with the experimental and previous QCT results, we calculated the rovibrational state-resolved integral cross-sections based on the reaction probabilities for total angular momentum J = 0 by using the J-shifting approximation, i.e., PJ>0(E) = PJ=0(E − ΔE), where ΔE = B*J(J + 1). The rotational constant is approximately calculated by B* = (2μR)−1 (R*)2 with R* being the length of the scattering coordinate at the transition-state32,33. Figure 2a shows the normalized vibrational populations and Fig. 2b–d show the normalized rotational populations for each v′(v′ = 1–3), of product HF at the collision energy of 0.234 eV (5.4 kcal/mol). The present QM vibrational state distributions agree quite well with the experiment, while the previous QCT results overestimate the population for v′ = 2. The calculated rotational state distributions, which changes from a broad one with a maximum at j = 6 for v′ = 0 to a narrow one with a maximum at j = 0 for v′ = 2, are also consistent with the experiment.
Analysis of scattering wave functions and VAP
To investigate the dynamical original of the oscillatory structures, we extract the scattering wave function at the collision energy of 0.126 eV for J = 0. Inspection of the two-dimensional contour of the wave function in the product region shown in Fig. 3a reveals there exist two nodes along the H–F coordinate (correlating to the HF product) inside the HB well, but only one node for the outgoing wave function with a sudden decrease in amplitude. This sudden change of nodal structure and amplitude in wave function is the characteristic feature of a Feshbach resonance, indicating the resonance is trapped in HF(v′ = 2) vibrational adiabatic well and produces mainly HF(v′ = 1) product.
To further confirm the mechanism of the resonance, we calculated the VAP for HF(v′ = 0–3) in the product region. As shown in Fig. 3b, the calculated VAP for all these vibrational states of HF have HB wells, which may support resonance states for the system. The HF(v′ = 2)-OH VAP has an asymptotic value of 0.193 eV measured from the ground state of the reactants and has a well with a well depth of 0.26 eV at a separation distance between HF and OH equal to 5.0 a0, i.e., the bottom of well is lower than the ground state of the reactants by 0.067 eV. The calculations show that the zero-point energy of the HF–OH in the well is also lower than the ground state of reactant. This means the F + H2O reaction can access to the Feshbach resonances supported by HF(v′ = 2)-OH VAP well even with zero collision energy and the oscillating structures on the HF(v′ = 0) and HF(v′ = 1) probabilities in Fig. 1b originate from these Feshbach resonance states. In addition, the small oscillating structures on the HF(v′ = 2) probabilities at the collision energy above 0.191 eV apparently originate from the bending excited HF(v′ = 2)-OH VAP supported resonance states.
As can be seen from Fig. 3b, the minimum positions for the HF(v′)-OH VAP wells for v′ = 0-2 all locate R~5.0 a0, whereas the well depth increases with v′, from 0.19 eV for v′ = 0, to 0.22 eV for v′ = 1, and 0.26 eV for v′ = 2. Therefore, these VAP wells are HB wells for the HF and OH products mainly through dipole–dipole interaction. With increase of the HF vibrational state, HF diploe moment increases, resulting in a deeper well depth. The Feshbash resonances for the F + H2O reaction are trapped in HF(v′ = 2)-OH HB well.
It is well known that in the F + H2 reaction, the interaction between H and HF in the transition-state region softens the HF bond, manifesting a larger anharmonicity for HF vibration. This anharmonicity considerably lowers the energy level for the HF(v′ = 3) vibrationally excited state, resulting in a peculiar well on the H-HF(v′ = 3) VAP right after the barrier position. This well is substantially deeper and locates much closer to barrier than the H-HF(v′ = 0–2) VAP wells, which originate from the vdw interaction between H and HF. Therefore, the resonance supporting HF(v′ = 2)-OH well in the F + H2O reaction is very different in nature from the H-HF(v′ = 3) VAP well, which support resonance states in the F + H2(v = 0) reaction8,34,35.
Relation between FH2O- photoelectron spectra and neutral reaction probabilities
The photodetachment of the FH2O− anion is also investigated quantum mechanically in full-dimensions on the new accurate NN PES. The anion PES is a newly constructed one by this group based on 11,936 CCSD(T)-F12a/(O and F: cc-pCVTZ-F12, H: cc-pVTZ-F12)36 ab initio energies with fundamental invariants NN (FI-NN) method37, yielding a root-mean-square-error of 0.043 meV. The ground rovibrational eigenstate on the anion PES is calculated and vertically projected onto the neutral PES for a continuous propagation using the diatom–diatom Jacobi coordinates. Figure 4a presents the total photoelectron spectra of the FH2O− system, showing groups of sharp peaks, denoted as S1~S3, accompanied with broad features in the high-energy wings (B1~B3). The present spectra are very similar to that of Ma and Guo28, except for the position of the B3 peak, which comes from the difference of the two neutral PESs. Figure 4a also shows the reactive flux spectra to both the HF + OH and F + H2O products. As can be seen, F + H2O channel becomes dominant after opening. Figure 4b shows the product HF vibrational state-resolved reactive flux spectra to HF + OH. According to the VAP in Fig. 3b, the sharp peaks correspond to Feshbach resonance states trapped in the well, which can only decay into lower HF vibrational states by vibrational predissociation, whereas the broad features correspond to direct dissociation products. To compare with the probabilities of the neutral reaction, we present the spectra only for energy between 1.5 and 1.9 eV in Fig. 4c. As can be seen, the reaction probabilities resemble the spectrum very well, except the oscillatory structures in the reaction probabilities are much richer, indicating the resonance mechanism of the two processes is the same.
Discussion
Therefore, the polyatomic F + H2O reaction also possesses Feshbach resonances trapped in VAP well in product region, similar to the triatomic F + H2 reaction. However, different from the F + H2 reaction in which the HF bond softening creates a peculiar well on the H-HF(v′ = 3) VAP, the VAP well for the F + H2O reaction is mainly HB in nature. Earlier photodetachment study of FH2O− discovered Feshbach resonances trapped in the HF + OH interaction well. In this study, we unambiguously show that the Feshbach resonances observed in the photodetachment study actually can be accessed by the F + H2O reaction. More efforts, in particular more joint efforts between theory and experiment should be devoted to studying these Feshbach resonances in the reaction at a quantitative level of accuracy, and to pushing forward our understanding of Feshbach resonances from triatomic reactions to polyatomic reactions.
Methods
Potential energy surface
The PES used in the calculation is a full-dimensional global one constructed using NN method based on ~26,429 AE-UCCSD(T)-F12a/CVTZ-F12 ab initio energies. The correction terms considering the influence of a larger basis set as well as SO couplings were further implemented. With an overall fitting error of ~1.78 meV measured in terms of root-mean-square-error (Supplementary Table 1), the PES is expected to have a quantitative level of accuracy. The static barrier height for the PES is 1.922 kcal/mol (and raised to 2.262 kcal/mol with SO energy included), lower than the value of 3.838 kcal/mol for the MRCI PES19 and very close to the value of 1.919 suggested by Dawes and collegues22.
Quantum dynamics calculations
Our state-to-state quantum calculation (J = 0) is carried out by using the product coordinates based30,38 scheme-based TDWP method with the coordinate systems shown in Supplementary Fig. 1. Test calculations reveal that the potential-averaged five-dimensional method with the nonreactive OH frozen at the ground vibrational state is of quantitative accuracy for the total reaction probabilities when compared with the full six-dimensional results as shown in Supplementary Fig. 2. The details of the PES and dynamics calculation are provided in Supplementary Information.
Data availability
The data that support the findings of this study are available from the corresponding author upon request.
Code availability
The accurate NN PESs developed in this work are available from the corresponding authors upon request.
References
Schatz, G. C. Detecting resonances. Science 288, 1599–1600 (2000).
Fernández-Alonso, F. & Zare, R. N. Scattering resonances in the simplest chemical reaction. Annu. Rev. Phys. Chem. 53, 67–99 (2002).
Liu, K. Quantum dynamical resonances in chemical reactions: from A + BC to polyatomic systems. Adv. Chem. Phys. 149, 1 (2012).
Schatz, G. C., Bowman, J. M. & Kuppermann, A. Large quantum effects in the collinear F+H2→FH+H reaction. J. Chem. Phys. 58, 4023 (1973).
Wu, S.-F., Johnson, B. R. & Levine, R. D. Quantum mechanical computational studies of chemical reactions: III. Collinear A + BC reaction with some model potential energy surfaces. Mol. Phys. 25, 839–856 (1973).
Neumark, D. M., Wodtke, A. M., Robinson, G. N., Hayden, C. C. & Lee, Y. T. Experimental investigation of resonances in reactive scattering: the F+H2 reaction. Phys. Rev. Lett. 53, 226–229 (1984).
Skodje, R. T. et al. Resonance-mediated chemical reaction: F+HD→HF+D. Phys. Rev. Lett. 85, 1206–1209 (2000).
Qiu, M. et al. Observation of Feshbach resonances in the F+H2→HF+H reaction. Science 311, 1440–1443 (2006).
Ren, Z. et al. Probing the resonance potential in the F atom reaction with hydrogen deuteride with spectroscopic accuracy. Proc. Natl Acad. Sci. USA 105, 12662 (2008).
Wang, T. et al. Dynamical resonances accessible only by reagent vibrational excitation in the F+HD→HF+D reaction. Science 342, 1499–1502 (2013).
Kim, J. B. et al. Spectroscopic observation of resonances in the F+H2 reaction. Science 349, 510–513 (2015).
Yang, T. et al. Enhanced reactivity of fluorine with para-hydrogen in cold interstellar clouds by resonance-induced quantum tunnelling. Nat. Chem. 11, 744–749 (2019).
Ricaud, P. & Lefèvre, F. Fluorine and the Environment: Atmospheric Chemistry, Emissions and Lithosphere (Elsevier, Amsterdam, 2006).
Neufeld, D. A., Zmuidzinas, J., Schilke. P. & Phillips, T. G. Discovery of interstellar hydrogen fluoride. Astrophys. J. 488, L141–L144 (1997).
Duewer, W. H. & Setser, D. W. Infrared chemiluminescence and energy partitioning from reactions of fluorine atoms with hydrides of carbon, silicon, oxygen, sulfur, nitrogen, and phosphorus. J. Chem. Phys. 58, 2310 (1973).
Agrawalla, B. S. & Setser, D. W. Infrared chemiluminescence and laser-induced fluorescence studies of energy disposal by reactions of fluorine and chlorine atoms with hydrogen sulfide, deuterium sulfide, hydrogen selenide, water, water-d2, and methanol. J. Phys. Chem. 90, 2450–2462 (1986).
Zolot, A. M. & Nesbitt, D. J. Crossed jet reactive scattering dynamics of F+H2O→HF(v,J)+OH:HF(v,J) product quantum state distributions under single-collision conditions. J. Chem. Phys. 129, 184305 (2008).
Ziemkiewicz, M. & Nesbitt, D. J. Nonadiabatic reactive scattering in atom+triatom systems: nascent rovibronic distributions in F+H2O→HF+OH. J. Chem. Phys. 131, 2054309 (2009).
Li, J., Dawes, R. & Guo, H. An ab initio based full-dimensional global potential energy surface for FH2O(X2A′) and dynamics for the F+H2O→HF+HO reaction. J. Chem. Phys. 137, 094304 (2012).
Li, J., Jiang, B. & Guo, H. Enhancement of bimolecular reactivity by a pre-reaction van der Waals complex: the case of F+H2O→HF+HO. Chem. Sci. 4, 629–632 (2013).
Li, J., Jiang, B. & Guo, H. Reactant vibrational excitations are more effective than translational energy in promoting an early-barrier reaction F+H2O→HF+HO. J. Am. Chem. Soc. 135, 982–985 (2013).
Nguyen, T. L., Li, J., Dawes, R., Stanton, J. F. & Guo, H. Accurate determination of barrier height and kinetics for the F+H2O→HF+HO reaction. J. Phys. Chem. A 117, 8864–8872 (2013).
Li, J. & Guo, H. Quasi-classical trajectory study of F+H2O→HF+HO reaction: Influence of barrier height, reactant rotational excitation, and isotopic substitution. Chin. J. Chem. Phys. 26, 627–634 (2013).
Li, J., Jiang, B. & Guo, H. Permutation invariant polynomial neural network approach to fitting potential energy surfaces. II. Four-atom systems. J. Chem. Phys. 139, 204103 (2013).
Yang, X., Wang, X.-B. & Wang, L.-S. Photodetachment of F−(H2O)n (n = 1–4): observation of charge-transfer states [F−(H2O)n +] and the transition state of F+H2O hydrogen abstraction reaction. J. Chem. Phys. 115, 2889–2892 (2001).
Ishikawa, Y., Nakajima, T., Yanai, T. & Hirao, K. Ab initio direct molecular dynamics study of the fragmentation of F(H2O) complex generated by photodetachment of F−(H2O) anion complex. Chem. Phys. Lett. 363, 458–464 (2002).
Otto, R. et al. Imaging dynamics on the F+H2O→HF+OH potential energy surfaces from wells to barriers. Science 343, 396–399 (2014).
Ma, J. & Guo, H. Reactive and nonreactive Feshbach resonances accessed by photodetachment of FH2O−. J. Phys. Chem. Lett. 6, 4822–4826 (2015).
Zhao, B. & Guo, H. Modulations of transition-state control of state-to-state dynamics in the F+H2O→HF+HO reaction. J. Phys. Chem. Lett. 6, 676–680 (2015).
Zhang, D. H. State-to-state quantum reactive scattering for four-atom chemical reactions: differential cross section for the H+H2O→H2+OH abstraction reaction. J. Chem. Phys. 125, 133102 (2006).
Liu, S. & Zhang, D. H. A local mode picture for H atom reaction with vibrationally excited H2O: a full dimensional state-to-state quantum dynamics investigation. Chem. Sci. 7, 261–265 (2016).
Bowman, J. M. Reduced dimensionality theory of quantum reactive scattering. J. Phys. Chem. 95, 4960–4968 (1991).
Song, H. & Guo, H. Vibrational and rotational mode specificity in the Cl+H2O→HCl+OH reaction: a quantum dynamical study. J. Phys. Chem. A. 119, 6188–6194 (2015).
Zhang, D. H. & Guo, H. Recent advances in quantum dynamics of bimolecular reactions. Annu. Rev. Phys. Chem. 67, 6.1–6.24 (2016).
Yang, T. et al. Extremely short-lived reaction resonances in Cl+HD(v=1)→DCl+H due to chemical bond softening. Science 347, 60 (2015).
Sarka, J., Lauvergnat, D., Brites, V., Császár, A. G. & Léonard, C. Rovibrational energy levels of the F−(H2O) and F−(D2O) complexes. Phys. Chem. Chem. Phys. 18, 17678–17690 (2016).
Shao, K., Chen, J., Zhao, Z. & Zhang, D. H. Fitting potential energy surfaces with fundamental invariant neural network. J. Chem. Phys. 145, 071101 (2016).
Fu, B. & Zhang, D. H. A time-dependent quantum dynamical study of the H + HBr reaction. J. Phys. Chem. A 111, 9516–9521 (2007).
Acknowledgements
This work was supported by the National Natural Science Foundation of China (grant numbers 21773235, 21590804, and 21688102), the Chinese Academy of Sciences (XDB17010200), and DICP (Grant: DICP ZZBS201611).
Author information
Authors and Affiliations
Contributions
D.H.Z. and S.L. conceived and supervised the research. X.Z. and J.C. developed the PES. S.L. performed the quantum dynamical calculation. L.L. simulated the photoelectron spectra. S.L. and D.H.Z. wrote the manuscript.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Peer review information Nature Communications thanks William Hase, and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Zhang, X., Li, L., Chen, J. et al. Feshbach resonances in the F + H2O → HF + OH reaction. Nat Commun 11, 223 (2020). https://doi.org/10.1038/s41467-019-14097-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-019-14097-y
This article is cited by
-
Observation of resonances in the transition state region of the F + NH3 reaction using anion photoelectron spectroscopy
Nature Chemistry (2023)
-
Atomistic dynamics of elimination and nucleophilic substitution disentangled for the F− + CH3CH2Cl reaction
Nature Chemistry (2021)
-
Mode specificity of a multi-channel reaction prototype: F + CH3OH → HF + CH3O/CH2OH
Theoretical Chemistry Accounts (2020)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.