1
|
Schatz GC, Wodtke AM, Yang X. Spiers Memorial Lecture: New directions in molecular scattering. Faraday Discuss 2024. [PMID: 38764350 DOI: 10.1039/d4fd00015c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/21/2024]
Abstract
The field of molecular scattering is reviewed as it pertains to gas-gas as well as gas-surface chemical reaction dynamics. We emphasize the importance of collaboration of experiment and theory, from which new directions of research are being pursued on increasingly complex problems. We review both experimental and theoretical advances that provide the modern toolbox available to molecular-scattering studies. We distinguish between two classes of work. The first involves simple systems and uses experiment to validate theory so that from the validated theory, one may learn far more than could ever be measured in the laboratory. The second class involves problems of great complexity that would be difficult or impossible to understand without a partnership of experiment and theory. Key topics covered in this review include crossed-beams reactive scattering and scattering at extremely low energies, where quantum effects dominate. They also include scattering from surfaces, reactive scattering and kinetics at surfaces, and scattering work done at liquid surfaces. The review closes with thoughts on future promising directions of research.
Collapse
Affiliation(s)
- George C Schatz
- Dept of Chemistry, Northwestern University, Evanston, Illinois 60208, USA
| | - Alec M Wodtke
- Institute for Physical Chemistry, Georg August University, Goettingen, Germany
- Max Planck Institute for Multidisciplinary Natural Sciences, Goettingen, Germany.
- International Center for the Advanced Studies of Energy Conversion, Georg August University, Goettingen, Germany
| | - Xueming Yang
- Dalian Institute for Chemical Physics, Chinese Academy of Sciences, Dalian, China
- Department of Chemistry, College of Science, Southern University of Science and Technology, Shenzhen, China
| |
Collapse
|
2
|
Powell A, Gerrits N, Tchakoua T, Somers MF, Busnengo HF, Meyer J, Kroes GJ, Doblhoff-Dier K. Best-of-Both-Worlds Predictive Approach to Dissociative Chemisorption on Metals. J Phys Chem Lett 2024; 15:307-315. [PMID: 38169287 PMCID: PMC10788952 DOI: 10.1021/acs.jpclett.3c02972] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2023] [Revised: 12/21/2023] [Accepted: 12/28/2023] [Indexed: 01/05/2024]
Abstract
Predictive capability, accuracy, and affordability are essential features of a theory that is capable of describing dissociative chemisorption on a metal surface. This type of reaction is important for heterogeneous catalysis. Here we present an approach in which we use diffusion Monte Carlo (DMC) to pin the minimum barrier height and construct a density functional that reproduces this value. This predictive approach allows the construction of a potential energy surface at the cost of density functional theory while retaining near DMC accuracy. Scrutinizing effects of energy dissipation and quantum tunneling, dynamics calculations suggest the approach to be of near chemical accuracy, reproducing molecular beam sticking experiments for the showcase H2 + Al(110) system to ∼1.4 kcal/mol.
Collapse
Affiliation(s)
- Andrew
D. Powell
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RA Leiden, The Netherlands
| | - Nick Gerrits
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RA Leiden, The Netherlands
| | - Theophile Tchakoua
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RA Leiden, The Netherlands
| | - Mark F. Somers
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RA Leiden, The Netherlands
| | - Heriberto F. Busnengo
- Instituto
de Física Rosario (IFIR), CONICET-UNR, 2000 Rosario, Argentina
- Facultad
de Ciencias Exatas, Ingeniería y
Agrimensura, UNR, 2000 Rosario, Argentina
| | - Jörg Meyer
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RA Leiden, The Netherlands
| | - Geert-Jan Kroes
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RA Leiden, The Netherlands
| | - Katharina Doblhoff-Dier
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, 2300 RA Leiden, The Netherlands
| |
Collapse
|
3
|
Fereidounpour P, Ramazani S. Collision of C2 with NH and deuterated analogs in the interstellar medium. Scattering calculations. Mol Phys 2023. [DOI: 10.1080/00268976.2023.2194458] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/30/2023]
|
4
|
Díaz C, Gravielle MS. Grazing incidence fast atom and molecule diffraction: theoretical challenges. Phys Chem Chem Phys 2022; 24:15628-15656. [PMID: 35730987 DOI: 10.1039/d2cp01246d] [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
This perspective article reviews the state-of-the-art of grazing incidence fast atom and molecule diffraction (GIFAD and GIFMD) simulations and addresses the main challenges that theorists, aiming to provide useful inputs in this topic, are facing. We first discuss briefly the methods used to build accurate potential energy surfaces describing the interaction between the projectile and the surface. Subsequently, we focus on the dynamics simulation methods for GIFAD, a phenomenon that has received a lot of experimental attention since 2007, when the first measurements were published. Following this experimental effort, theorists have developed and adapted a bunch of methods able to simulate, analyze and extract information from the experimental outputs. We review these methods, from the very simple ones based on classical dynamics to the full quantum ones, paying special attention to more versatile semiclassical approaches, which include quantum ingredients in the dynamics at a computational cost only slightly higher than that required in classical dynamics. Within the semiclassical framework it is possible, for example, to include in the dynamics the surface phonons and the projectile coherence, two factors that may have a relevant influence on the experimental measurements, at a reasonable computational cost. Finally, we address GIFMD, a phenomenon that has received much less attention and for which there is still a lot of room for research. We review the few examples of GIFMD available in the literature, and we discuss new phenomena associated with the molecular internal degrees of freedom, which may have some impact in other closely related fields, such as molecular reactivity on metal surfaces. Finally, we point out opened questions, raised from the comparisons between theoretical and experimental results, which claim for further experimental efforts.
Collapse
Affiliation(s)
- Cristina Díaz
- Departamento de Química Física, Facultad de CC. Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain.
| | - María Silvia Gravielle
- Instituto de Astronomía y Física del Espacio (IAFE, UBA-CONICET), Ciudad Universitaria, C1428EGA, Buenos Aires, Argentina.
| |
Collapse
|
5
|
Kwon HY, Morrow Z, Kelley CT, Jakubikova E. Interpolation Methods for Molecular Potential Energy Surface Construction. J Phys Chem A 2021; 125:9725-9735. [PMID: 34730973 DOI: 10.1021/acs.jpca.1c06812] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The concept of a potential energy surface (PES) is one of the most important concepts in modern chemistry. A PES represents the relationship between the chemical system's energy and its geometry (i.e., atom positions) and can provide useful information about the system's chemical properties and reactivity. Construction of accurate PESs with high-level theoretical methodologies, such as density functional theory, is still challenging due to a steep increase in the computational cost with the increase of the system size. Thus, over the past few decades, many different mathematical approaches have been applied to the problem of the cost-efficient PES construction. This article serves as a short overview of interpolative methods for the PES construction, including global polynomial interpolation, trigonometric interpolation, modified Shepard interpolation, interpolative moving least-squares, and the automated PES construction derived from these.
Collapse
Affiliation(s)
- Hyuk-Yong Kwon
- Department of Chemistry, North Carolina State University, Raleigh, North Carolina 27695, United States
| | - Zachary Morrow
- Department of Mathematics, North Carolina State University, Raleigh, North Carolina 27695, United States
| | - C T Kelley
- Department of Mathematics, North Carolina State University, Raleigh, North Carolina 27695, United States
| | - Elena Jakubikova
- Department of Chemistry, North Carolina State University, Raleigh, North Carolina 27695, United States
| |
Collapse
|
6
|
Zhou X, Zhang Y, Yin R, Hu C, Jiang B. Neural Network Representations for Studying
Gas‐Surface
Reaction Dynamics: Beyond the
Born‐Oppenheimer
Static Surface Approximation
†. CHINESE J CHEM 2021. [DOI: 10.1002/cjoc.202100303] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Affiliation(s)
- Xueyao Zhou
- Hefei National Laboratory for Physical Science at the Microscale, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes, Department of Chemical Physics University of Science and Technology of China Hefei Anhui 230026 China
| | - Yaolong Zhang
- Hefei National Laboratory for Physical Science at the Microscale, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes, Department of Chemical Physics University of Science and Technology of China Hefei Anhui 230026 China
| | - Rongrong Yin
- Hefei National Laboratory for Physical Science at the Microscale, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes, Department of Chemical Physics University of Science and Technology of China Hefei Anhui 230026 China
| | - Ce Hu
- Hefei National Laboratory for Physical Science at the Microscale, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes, Department of Chemical Physics University of Science and Technology of China Hefei Anhui 230026 China
| | - Bin Jiang
- Hefei National Laboratory for Physical Science at the Microscale, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes, Department of Chemical Physics University of Science and Technology of China Hefei Anhui 230026 China
| |
Collapse
|
7
|
Kroes GJ. Computational approaches to dissociative chemisorption on metals: towards chemical accuracy. Phys Chem Chem Phys 2021; 23:8962-9048. [PMID: 33885053 DOI: 10.1039/d1cp00044f] [Citation(s) in RCA: 33] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We review the state-of-the-art in the theory of dissociative chemisorption (DC) of small gas phase molecules on metal surfaces, which is important to modeling heterogeneous catalysis for practical reasons, and for achieving an understanding of the wealth of experimental information that exists for this topic, for fundamental reasons. We first give a quick overview of the experimental state of the field. Turning to the theory, we address the challenge that barrier heights (Eb, which are not observables) for DC on metals cannot yet be calculated with chemical accuracy, although embedded correlated wave function theory and diffusion Monte-Carlo are moving in this direction. For benchmarking, at present chemically accurate Eb can only be derived from dynamics calculations based on a semi-empirically derived density functional (DF), by computing a sticking curve and demonstrating that it is shifted from the curve measured in a supersonic beam experiment by no more than 1 kcal mol-1. The approach capable of delivering this accuracy is called the specific reaction parameter (SRP) approach to density functional theory (DFT). SRP-DFT relies on DFT and on dynamics calculations, which are most efficiently performed if a potential energy surface (PES) is available. We therefore present a brief review of the DFs that now exist, also considering their performance on databases for Eb for gas phase reactions and DC on metals, and for adsorption to metals. We also consider expressions for SRP-DFs and briefly discuss other electronic structure methods that have addressed the interaction of molecules with metal surfaces. An overview is presented of dynamical models, which make a distinction as to whether or not, and which dissipative channels are modeled, the dissipative channels being surface phonons and electronically non-adiabatic channels such as electron-hole pair excitation. We also discuss the dynamical methods that have been used, such as the quasi-classical trajectory method and quantum dynamical methods like the time-dependent wave packet method and the reaction path Hamiltonian method. Limits on the accuracy of these methods are discussed for DC of diatomic and polyatomic molecules on metal surfaces, paying particular attention to reduced dimensionality approximations that still have to be invoked in wave packet calculations on polyatomic molecules like CH4. We also address the accuracy of fitting methods, such as recent machine learning methods (like neural network methods) and the corrugation reducing procedure. In discussing the calculation of observables we emphasize the importance of modeling the properties of the supersonic beams in simulating the sticking probability curves measured in the associated experiments. We show that chemically accurate barrier heights have now been extracted for DC in 11 molecule-metal surface systems, some of which form the most accurate core of the only existing database of Eb for DC reactions on metal surfaces (SBH10). The SRP-DFs (or candidate SRP-DFs) that have been derived show transferability in many cases, i.e., they have been shown also to yield chemically accurate Eb for chemically related systems. This can in principle be exploited in simulating rates of catalyzed reactions on nano-particles containing facets and edges, as SRP-DFs may be transferable among systems in which a molecule dissociates on low index and stepped surfaces of the same metal. In many instances SRP-DFs have allowed important conclusions regarding the mechanisms underlying observed experimental trends. An important recent observation is that SRP-DFT based on semi-local exchange DFs has so far only been successful for systems for which the difference of the metal work function and the molecule's electron affinity exceeds 7 eV. A main challenge to SRP-DFT is to extend its applicability to the other systems, which involve a range of important DC reactions of e.g. O2, H2O, NH3, CO2, and CH3OH. Recent calculations employing a PES based on a screened hybrid exchange functional suggest that the road to success may be based on using exchange functionals of this category.
Collapse
Affiliation(s)
- Geert-Jan Kroes
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands.
| |
Collapse
|
8
|
Masoumpour MS, Daryanavard M. Ab Initio-Based Global Potential Energy Surface and Reaction Dynamics for H 2S + C. J Phys Chem A 2020; 124:7901-7910. [DOI: 10.1021/acs.jpca.0c05856] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
| | - Marzieh Daryanavard
- Department of chemistry, Estahban Higher Education Center, Estahban 74519-44655, Iran
| |
Collapse
|
9
|
Lomaka A, Tamm T. Linearization of moment tensor potentials for multicomponent systems with a preliminary assessment for short-range interaction energy in water dimer and trimer. J Chem Phys 2020; 152:164115. [DOI: 10.1063/5.0007473] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Andre Lomaka
- Department of Chemistry and Biotechnology, Tallinn University of Technology, Tallinn, Estonia
| | - Toomas Tamm
- Department of Chemistry and Biotechnology, Tallinn University of Technology, Tallinn, Estonia
| |
Collapse
|
10
|
Del Cueto M, Muzas AS, Frankcombe TJ, Martín F, Díaz C. Prominent out-of-plane diffraction in helium scattering from a methyl-terminated Si(111) surface. Phys Chem Chem Phys 2019; 21:15879-15887. [PMID: 31286123 DOI: 10.1039/c9cp02141h] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Due to their electrochemical and oxidative stability, organic-terminated semiconductor surfaces are well suited to applications in, for example, photoelectrodes and electrochemical cells, which explains the lively interest in their detailed characterization. Helium atom scattering (HAS) is a useful tool to carry out such characterization. Here, we have simulated HAS in He/CH3-Si(111) based on density functional theory (DFT) potential energy surfaces (PESs) and multi-configuration time-dependent Hartree (MCTDH) dynamics. Our analysis of HAS shows that most diffraction taking place in this system corresponds to high-order out-of-plane peaks. This is a general trend that does not depend on the specific features of the simulations, such as the inclusion or not of the van der Waals long-range effects. This is the first and only He-surface system for which such huge out-of-plane diffraction has been described. This striking theoretical finding should encourage new experimental developments to confirm this previously unreported effect.
Collapse
Affiliation(s)
- Marcos Del Cueto
- Departamento de Química Módulo 13, Universidad Autónoma de Madrid, 28049 Madrid, Spain.
| | - Alberto S Muzas
- University of New South Wales, Canberra, ACT 2600, Australia
| | | | - Fernando Martín
- Departamento de Química Módulo 13, Universidad Autónoma de Madrid, 28049 Madrid, Spain. and Condensed Matter Physics Center (IFIMAC), Universidad Autónoma de Madrid, 28049 Madrid, Spain and Instituto Madrileño de Estudios Avanzados en Nanociencia (IMDEA-Nanociencia), Cantoblanco, 28049 Madrid, Spain
| | - Cristina Díaz
- Departamento de Química Módulo 13, Universidad Autónoma de Madrid, 28049 Madrid, Spain. and Condensed Matter Physics Center (IFIMAC), Universidad Autónoma de Madrid, 28049 Madrid, Spain and Institute for Advanced Research in Chemistry (IAdChem), Universidad Autónoma de Madrid, 28049 Madrid, Spain
| |
Collapse
|
11
|
Jiang B, Guo H. Dynamics in reactions on metal surfaces: A theoretical perspective. J Chem Phys 2019; 150:180901. [PMID: 31091904 DOI: 10.1063/1.5096869] [Citation(s) in RCA: 31] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Recent advances in theoretical characterization of reaction dynamics on metal surfaces are reviewed. It is shown that the widely available density functional theory of metals and their interactions with molecules have enabled first principles theoretical models for treating surface reaction dynamics. The new theoretical tools include methods to construct high-dimensional adiabatic potential energy surfaces, to characterize nonadiabatic processes within the electronic friction models, and to describe dynamics both quantum mechanically and classically. Three prototypical surface reactions, namely, dissociative chemisorption, Eley-Rideal reactions, and recombinative desorption, are surveyed with a focus on some representative examples. While principles governing gas phase reaction dynamics may still be applicable, the presence of the surface introduces a higher level of complexity due to strong interaction between the molecular species and metal substrate. Furthermore, most of these reactive processes are impacted by energy exchange with surface phonons and/or electron-hole pair excitations. These theoretical studies help to interpret and rationalize experimental observations and, in some cases, guide experimental explorations. Knowledge acquired in these fundamental studies is expected to impact many practical problems in a wide range of interfacial processes.
Collapse
Affiliation(s)
- Bin Jiang
- Hefei National Laboratory for Physical Science at the Microscale, Department of Chemical Physics, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes, University of Science and Technology of China, Hefei 230026, China
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA
| |
Collapse
|
12
|
Nour Ghassemi E, Somers M, Kroes GJ. Test of the Transferability of the Specific Reaction Parameter Functional for H 2 + Cu(111) to D 2 + Ag(111). THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2018; 122:22939-22952. [PMID: 30344838 PMCID: PMC6189907 DOI: 10.1021/acs.jpcc.8b05658] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/13/2018] [Revised: 09/04/2018] [Indexed: 06/08/2023]
Abstract
The accurate description of the dissociative chemisorption of a molecule on a metal surface requires a chemically accurate description of the molecule-surface interaction. Previously, it was shown that the specific reaction parameter approach to density functional theory (SRP-DFT) enables accurate descriptions of the reaction of dihydrogen with metal surfaces in, for instance, H2 + Pt(111), H2 + Cu(111), and H2 + Cu(100). SRP-DFT likewise allowed a chemically accurate description of dissociation of methane on Ni(111) and Pt(111), and the SRP functional for CH4 + Ni(111) was transferable to CH4 + Pt(111), where Ni and Pt belong to the same group. Here, we investigate whether the SRP density functional derived for H2 + Cu(111) also gives chemically accurate results for H2 + Ag(111), where Ag belongs to the same group as Cu. To do this, we have performed quasi-classical trajectory calculations using the six-dimensional potential energy surface of H2 + Ag(111) within the Born-Oppenheimer static surface approximation. The computed reaction probabilities are compared with both state-resolved associative desorption and molecular beam sticking experiments. Our results do not yet show transferability, as the computed sticking probabilities and initial-state selected reaction probabilities are shifted relative to experiment to higher energies by about 2-3 kcal/mol. The lack of transferability may be due to the different character of the SRP functionals for H2 + Cu and CH4 + group 10 metals, the latter containing a van der Waals correlation functional and the former not.
Collapse
|
13
|
Klinting EL, Thomsen B, Godtliebsen IH, Christiansen O. Employing general fit-bases for construction of potential energy surfaces with an adaptive density-guided approach. J Chem Phys 2018; 148:064113. [DOI: 10.1063/1.5016259] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
| | - Bo Thomsen
- Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark
| | | | - Ove Christiansen
- Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark
| |
Collapse
|
14
|
Frankcombe TJ. Interpolating DFT Data for 15D Modeling of Methane Dissociation on an fcc Metal. INT J CHEM KINET 2018. [DOI: 10.1002/kin.21157] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Terry J. Frankcombe
- School of Physical, Environmental and Mathematical Sciences; University of New South Wales; PO Box 7916 Canberra BC 2610 Australia
- Research School of Chemistry; Australian National University; Canberra ACT 2601 Australia
| |
Collapse
|
15
|
Nour Ghassemi E, Wijzenbroek M, Somers MF, Kroes GJ. Chemically accurate simulation of dissociative chemisorption of D2 on Pt(1 1 1). Chem Phys Lett 2017. [DOI: 10.1016/j.cplett.2016.12.059] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
16
|
Del Cueto M, Muzas AS, Somers MF, Kroes GJ, Díaz C, Martín F. Exploring surface landscapes with molecules: rotationally induced diffraction of H 2 on LiF(001) under fast grazing incidence conditions. Phys Chem Chem Phys 2017. [PMID: 28621794 DOI: 10.1039/c7cp02904g] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Atomic diffraction by surfaces under fast grazing incidence conditions has been used for almost a decade to characterize surface properties with more accuracy than with more traditional atomic diffraction methods. From six-dimensional solutions of the time-dependent Schrödinger equation, we show that diffraction of H2 molecules under fast grazing incidence conditions could be even more informative for the characterization of ionic surfaces, due to the large anisotropic electrostatic interaction between the quadrupole moment of the molecule and the electric field created by the ionic crystal. Using the LiF(001) surface as a benchmark, we show that fast grazing incidence diffraction of H2 strongly depends on the initial rotational state of the molecule, while rotationally inelastic processes are irrelevant. We demonstrate that, as a result of the anisotropy of the impinging projectile, initial rotational excitation leads to an increase in intensity of high-order diffraction peaks at incidence directions that satisfy precise symmetry constraints, thus providing a more detailed information on the surface characteristics than that obtained from low-order atomic diffraction peaks under fast grazing incidence conditions. As quadrupole-ion surface potentials are expected to accurately represent the interaction between H2 and any surface with a marked ionic character, our analysis should be of general applicability to any of such surfaces. Finally, we show that a density functional theory description of the molecule-ion surface potential catches the main features observed experimentally.
Collapse
Affiliation(s)
- M Del Cueto
- Departamento de Química Módulo 13, Universidad Autónoma de Madrid, 28049 Madrid, Spain.
| | | | | | | | | | | |
Collapse
|
17
|
Kim CW, Rhee YM. Constructing an Interpolated Potential Energy Surface of a Large Molecule: A Case Study with Bacteriochlorophyll a Model in the Fenna-Matthews-Olson Complex. J Chem Theory Comput 2016; 12:5235-5246. [PMID: 27760297 DOI: 10.1021/acs.jctc.6b00647] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Constructing a reliable potential energy surface (PES) is a key step toward computationally studying the chemical dynamics of any molecular system. The interpolation scheme is a useful tool that can closely follow the accuracy of quantum chemical means at a dramatically reduced computational cost. However, applying interpolation to building a PES of a large molecule is not a straightforward black-box approach, as it frequently encounters practical difficulties associated with its large dimensionality. Here, we present detailed courses of applying interpolation toward building a PES of a large chromophore molecule. We take the example of S0 and S1 electronic states of bacteriochlorophyll a (BChla) molecules in the Fenna-Matthews-Olson light harvesting complex. With a reduced model molecule that bears BChla's main π-conjugated ring, various practical approaches are designed for improving the PES quality in a stable manner and for fine-tuning the final surface such that the surface can be adopted for long time molecular dynamics simulations. Combined with parallel implementation, we show that interpolated mechanics/molecular mechanics (IM/MM) simulations of the entire complex in the nanosecond time scale can be conducted readily without any practical issues. With 1500 interpolation data points for each chromophore unit, the PES error relative to the reference quantum chemical calculation is found to be ∼0.15 eV in the thermally accessible region of the conformational space, together with ∼0.01 eV error in S0 - S1 transition energies. The performance issue related to the use of a large interpolation database within the framework of our parallel routines is also discussed.
Collapse
Affiliation(s)
- Chang Woo Kim
- Center for Self-assembly and Complexity, Institute for Basic Science (IBS) , Pohang 37673, Korea.,Department of Chemistry, Pohang University of Science and Technology (POSTECH) , Pohang 37673, Korea
| | - Young Min Rhee
- Center for Self-assembly and Complexity, Institute for Basic Science (IBS) , Pohang 37673, Korea.,Department of Chemistry, Pohang University of Science and Technology (POSTECH) , Pohang 37673, Korea
| |
Collapse
|
18
|
Wijzenbroek M, Helstone D, Meyer J, Kroes GJ. Dynamics of H2 dissociation on the close-packed (111) surface of the noblest metal: H2 + Au(111). J Chem Phys 2016; 145:144701. [DOI: 10.1063/1.4964486] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Mark Wijzenbroek
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Darcey Helstone
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Jörg Meyer
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Geert-Jan Kroes
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
19
|
Nihill KJ, Hund ZM, Muzas A, Díaz C, Del Cueto M, Frankcombe T, Plymale NT, Lewis NS, Martín F, Sibener SJ. Experimental and theoretical study of rotationally inelastic diffraction of H2(D2) from methyl-terminated Si(111). J Chem Phys 2016; 145:084705. [PMID: 27586939 DOI: 10.1063/1.4961257] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Fundamental details concerning the interaction between H2 and CH3-Si(111) have been elucidated by the combination of diffractive scattering experiments and electronic structure and scattering calculations. Rotationally inelastic diffraction (RID) of H2 and D2 from this model hydrocarbon-decorated semiconductor interface has been confirmed for the first time via both time-of-flight and diffraction measurements, with modest j = 0 → 2 RID intensities for H2 compared to the strong RID features observed for D2 over a large range of kinematic scattering conditions along two high-symmetry azimuthal directions. The Debye-Waller model was applied to the thermal attenuation of diffraction peaks, allowing for precise determination of the RID probabilities by accounting for incoherent motion of the CH3-Si(111) surface atoms. The probabilities of rotationally inelastic diffraction of H2 and D2 have been quantitatively evaluated as a function of beam energy and scattering angle, and have been compared with complementary electronic structure and scattering calculations to provide insight into the interaction potential between H2 (D2) and hence the surface charge density distribution. Specifically, a six-dimensional potential energy surface (PES), describing the electronic structure of the H2(D2)/CH3-Si(111) system, has been computed based on interpolation of density functional theory energies. Quantum and classical dynamics simulations have allowed for an assessment of the accuracy of the PES, and subsequently for identification of the features of the PES that serve as classical turning points. A close scrutiny of the PES reveals the highly anisotropic character of the interaction potential at these turning points. This combination of experiment and theory provides new and important details about the interaction of H2 with a hybrid organic-semiconductor interface, which can be used to further investigate energy flow in technologically relevant systems.
Collapse
Affiliation(s)
- Kevin J Nihill
- The James Franck Institute and Department of Chemistry, The University of Chicago, 929 E. 57th Street, Chicago, Illinois 60637, USA
| | - Zachary M Hund
- The James Franck Institute and Department of Chemistry, The University of Chicago, 929 E. 57th Street, Chicago, Illinois 60637, USA
| | - Alberto Muzas
- Departamento de Química, Módulo 13, Universidad Autónoma de Madrid, 28049 Madrid, Spain
| | - Cristina Díaz
- Departamento de Química, Módulo 13, Universidad Autónoma de Madrid, 28049 Madrid, Spain
| | - Marcos Del Cueto
- Departamento de Química, Módulo 13, Universidad Autónoma de Madrid, 28049 Madrid, Spain
| | - Terry Frankcombe
- School of Physical, Environmental and Mathematical Sciences, University of New South Wales, Canberra ACT 2610, Australia
| | - Noah T Plymale
- Division of Chemistry and Chemical Engineering, Beckman Institute and Kavli Nanoscience Institute, California Institute of Technology, 210 Noyes Laboratory, 127-72, Pasadena, California 91125, USA
| | - Nathan S Lewis
- Division of Chemistry and Chemical Engineering, Beckman Institute and Kavli Nanoscience Institute, California Institute of Technology, 210 Noyes Laboratory, 127-72, Pasadena, California 91125, USA
| | - Fernando Martín
- Departamento de Química, Módulo 13, Universidad Autónoma de Madrid, 28049 Madrid, Spain
| | - S J Sibener
- The James Franck Institute and Department of Chemistry, The University of Chicago, 929 E. 57th Street, Chicago, Illinois 60637, USA
| |
Collapse
|
20
|
Liu T, Zhang Z, Fu B, Yang X, Zhang DH. A seven-dimensional quantum dynamics study of the dissociative chemisorption of H 2O on Cu(111): effects of azimuthal angles and azimuthal angle-averaging. Chem Sci 2016; 7:1840-1845. [PMID: 29899905 PMCID: PMC5965058 DOI: 10.1039/c5sc03689e] [Citation(s) in RCA: 53] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2015] [Accepted: 11/25/2015] [Indexed: 01/04/2023] Open
Abstract
We report the first seven-dimensional quantum dynamics study for the dissociative chemisorption of H2O on Cu(111) using the time-dependent wave-packet approach, based on an accurate nine-dimensional potential energy surface (PES), which is newly developed by neural network fitting to ∼80 000 density functional theory points. This seven-dimensional quantum model allows the examination of the influence of azimuthal angles and also the investigation of the quantitative relationship between the seven-dimensional (7D) dissociation probabilities and those results calculated by the six-dimensional (6D) model with the flat surface approximation. The reactivity strongly depends on the azimuthal rotations due to different barrier heights. Very large differences are seen between the 7D dissociation probabilities and the 6D results with fixed azimuthal angles, at different fixed sites of impact, indicating that the 6D model by neglecting the azimuthal rotation can introduce substantial errors in calculating dissociation probabilities and the 7D quantum dynamics is essential to investigate the dissociation process. A new azimuthal angle-averaging approach is proposed that the 7D dissociation probability can be well reproduced by averaging 6D results over 18 azimuthal angles, in particular in low kinetic energy regions.
Collapse
Affiliation(s)
- Tianhui Liu
- Department of Chemical Physics , University of Science and Technology of China , Hefei , China 230026
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry , Dalian Institute of Chemical Physics , Chinese Academy of Sciences , Dalian , China 116023 . ;
| | - Zhaojun Zhang
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry , Dalian Institute of Chemical Physics , Chinese Academy of Sciences , Dalian , China 116023 . ;
| | - Bina Fu
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry , Dalian Institute of Chemical Physics , Chinese Academy of Sciences , Dalian , China 116023 . ;
| | - Xueming Yang
- Department of Chemical Physics , University of Science and Technology of China , Hefei , China 230026
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry , Dalian Institute of Chemical Physics , Chinese Academy of Sciences , Dalian , China 116023 . ;
| | - Dong H Zhang
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry , Dalian Institute of Chemical Physics , Chinese Academy of Sciences , Dalian , China 116023 . ;
| |
Collapse
|
21
|
Kroes GJ, Díaz C. Quantum and classical dynamics of reactive scattering of H2 from metal surfaces. Chem Soc Rev 2016; 45:3658-700. [DOI: 10.1039/c5cs00336a] [Citation(s) in RCA: 120] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
State-of-the-art theoretical models allow nowadays an accurate description of H2/metal surface systems and phenomena relative to heterogeneous catalysis. Here we review the most relevant ones investigated during the last 10 years.
Collapse
Affiliation(s)
- Geert-Jan Kroes
- Leiden Institute of Chemistry
- Gorlaeus Laboratories
- Leiden University
- 2300 RA Leiden
- The Netherlands
| | - Cristina Díaz
- Departamento de Química
- Módulo 13
- Universidad Autónoma de Madrid
- 28049 Madrid
- Spain
| |
Collapse
|
22
|
Jiang B, Yang M, Xie D, Guo H. Quantum dynamics of polyatomic dissociative chemisorption on transition metal surfaces: mode specificity and bond selectivity. Chem Soc Rev 2016; 45:3621-40. [DOI: 10.1039/c5cs00360a] [Citation(s) in RCA: 122] [Impact Index Per Article: 15.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Recent advances in quantum dynamical characterization of polyatomic dissociative chemisorption on accurate global potential energy surfaces are critically reviewed.
Collapse
Affiliation(s)
- Bin Jiang
- Department of Chemistry and Chemical Biology
- University of New Mexico
- Albuquerque
- USA
- Department of Chemical Physics
| | - Minghui Yang
- Key Laboratory of Magnetic Resonance in Biological Systems
- State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics
- Wuhan Centre for Magnetic Resonance
- Wuhan Institute of Physics and Mathematics
- Chinese Academy of Sciences
| | - Daiqian Xie
- Institute of Theoretical and Computational Chemistry
- Key Laboratory of Mesoscopic Chemistry
- School of Chemistry and Chemical Engineering
- Nanjing University
- Nanjing 210093
| | - Hua Guo
- Department of Chemistry and Chemical Biology
- University of New Mexico
- Albuquerque
- USA
| |
Collapse
|
23
|
D'Arcy JH, Jordan MJT, Frankcombe TJ, Collins MA. H2 Adsorption in a Porous Crystal: Accurate First-Principles Quantum Simulation. J Phys Chem A 2015; 119:12166-81. [PMID: 26322374 DOI: 10.1021/acs.jpca.5b06074] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A general method is presented for constructing, from ab initio quantum chemistry calculations, the potential energy surface (PES) for H2 absorbed in a porous crystalline material. The method is illustrated for the metal-organic framework material MOF-5. Rigid body quantum diffusion Monte Carlo simulations are used in the construction of the PES and to evaluate the quantum ground state of H2 in MOF-5, the zero-point energy, and the enthalpy of adsorption at 0 K.
Collapse
Affiliation(s)
- Jordan H D'Arcy
- School of Chemistry, The University of Sydney , Sydney NSW 2006, Australia
| | | | - Terry J Frankcombe
- Research School of Chemistry, Australian National University , Canberra ACT 0200, Australia
| | - Michael A Collins
- Research School of Chemistry, Australian National University , Canberra ACT 0200, Australia
| |
Collapse
|
24
|
Jiang B, Guo H. Permutation invariant polynomial neural network approach to fitting potential energy surfaces. III. Molecule-surface interactions. J Chem Phys 2014; 141:034109. [DOI: 10.1063/1.4887363] [Citation(s) in RCA: 106] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Bin Jiang
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA
| |
Collapse
|
25
|
Frankcombe TJ. Using Hessian update formulae to construct modified Shepard interpolated potential energy surfaces: Application to vibrating surface atoms. J Chem Phys 2014; 140:114108. [DOI: 10.1063/1.4868637] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
26
|
Boereboom JM, Wijzenbroek M, Somers MF, Kroes GJ. Towards a specific reaction parameter density functional for reactive scattering of H2 from Pd(111). J Chem Phys 2013; 139:244707. [DOI: 10.1063/1.4851355] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|