1
|
Litman Y, Pós ES, Box CL, Martinazzo R, Maurer RJ, Rossi M. Dissipative tunneling rates through the incorporation of first-principles electronic friction in instanton rate theory. II. Benchmarks and applications. J Chem Phys 2022; 156:194107. [PMID: 35597654 DOI: 10.1063/5.0088400] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In Paper I [Litman et al., J. Chem. Phys. (in press) (2022)], we presented the ring-polymer instanton with explicit friction (RPI-EF) method and showed how it can be connected to the ab initio electronic friction formalism. This framework allows for the calculation of tunneling reaction rates that incorporate the quantum nature of the nuclei and certain types of non-adiabatic effects (NAEs) present in metals. In this paper, we analyze the performance of RPI-EF on model potentials and apply it to realistic systems. For a 1D double-well model, we benchmark the method against numerically exact results obtained from multi-layer multi-configuration time-dependent Hartree calculations. We demonstrate that RPI-EF is accurate for medium and high friction strengths and less accurate for extremely low friction values. We also show quantitatively how the inclusion of NAEs lowers the crossover temperature into the deep tunneling regime, reduces the tunneling rates, and, in certain regimes, steers the quantum dynamics by modifying the tunneling pathways. As a showcase of the efficiency of this method, we present a study of hydrogen and deuterium hopping between neighboring interstitial sites in selected bulk metals. The results show that multidimensional vibrational coupling and nuclear quantum effects have a larger impact than NAEs on the tunneling rates of diffusion in metals. Together with Paper I [Litman et al., J. Chem. Phys. (in press) (2022)], these results advance the calculations of dissipative tunneling rates from first principles.
Collapse
Affiliation(s)
- Y Litman
- MPI for the Structure and Dynamics of Matter, Luruper Chaussee 149, 22761 Hamburg, Germany
| | - E S Pós
- MPI for the Structure and Dynamics of Matter, Luruper Chaussee 149, 22761 Hamburg, Germany
| | - C L Box
- Department of Chemistry, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - R Martinazzo
- Department of Chemistry, Università degli Studi di Milano, Via Golgi 19, 20133 Milano, Italy
| | - R J Maurer
- Department of Chemistry, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - M Rossi
- MPI for the Structure and Dynamics of Matter, Luruper Chaussee 149, 22761 Hamburg, Germany
| |
Collapse
|
2
|
Litman Y, Pós ES, Box CL, Martinazzo R, Maurer RJ, Rossi M. Dissipative tunneling rates through the incorporation of first-principles electronic friction in instanton rate theory. I. Theory. J Chem Phys 2022; 156:194106. [PMID: 35597633 DOI: 10.1063/5.0088399] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Reactions involving adsorbates on metallic surfaces and impurities in bulk metals are ubiquitous in a wide range of technological applications. The theoretical modeling of such reactions presents a formidable challenge for theory because nuclear quantum effects (NQEs) can play a prominent role and the coupling of the atomic motion with the electrons in the metal gives rise to important non-adiabatic effects (NAEs) that alter atomic dynamics. In this work, we derive a theoretical framework that captures both NQEs and NAEs and, due to its high efficiency, can be applied to first-principles calculations of reaction rates in high-dimensional realistic systems. More specifically, we develop a method that we coin ring polymer instanton with explicit friction (RPI-EF), starting from the ring polymer instanton formalism applied to a system-bath model. We derive general equations that incorporate the spatial and frequency dependence of the friction tensor and then combine this method with the ab initio electronic friction formalism for the calculation of thermal reaction rates. We show that the connection between RPI-EF and the form of the electronic friction tensor presented in this work does not require any further approximations, and it is expected to be valid as long as the approximations of both underlying theories remain valid.
Collapse
Affiliation(s)
- Y Litman
- MPI for the Structure and Dynamics of Matter, Luruper Chaussee 149, 22761 Hamburg, Germany
| | - E S Pós
- MPI for the Structure and Dynamics of Matter, Luruper Chaussee 149, 22761 Hamburg, Germany
| | - C L Box
- Department of Chemistry, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - R Martinazzo
- Department of Chemistry, Università degli Studi di Milano, Via Golgi 19, 20133 Milano, Italy
| | - R J Maurer
- Department of Chemistry, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - M Rossi
- MPI for the Structure and Dynamics of Matter, Luruper Chaussee 149, 22761 Hamburg, Germany
| |
Collapse
|
3
|
Quartapelle Procopio E, Benincori T, Appoloni G, Mussini PR, Arnaboldi S, Carbonera C, Cirilli R, Cominetti A, Longo L, Martinazzo R, Panigati M, Pò R. A family of solution-processable macrocyclic and open-chain oligothiophenes with atropoisomeric scaffolds: structural and electronic features for potential energy applications. NEW J CHEM 2017. [DOI: 10.1039/c7nj01501a] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
Cyclic and linear inherently chiral oligothiophenes have been employed as donors in bulk heterojunction solar cells with fullerenes as acceptors.
Collapse
|
4
|
Casolo S, Tantardini GF, Martinazzo R. Hydrogen Recombination and Dimer Formation on Graphite from Ab Initio Molecular Dynamics Simulations. J Phys Chem A 2016; 120:5032-40. [DOI: 10.1021/acs.jpca.5b12761] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- S. Casolo
- Dipartimento
di Chimica, Università degli Studi di Milano, via Golgi
19, 20133 Milan, Italy
| | - G. F. Tantardini
- Dipartimento
di Chimica, Università degli Studi di Milano, via Golgi
19, 20133 Milan, Italy
- Istituto di Scienze
e Tecnologie Molecolari, CNR-ISTM, via Golgi 19, 20133 Milan, Italy
| | - R. Martinazzo
- Dipartimento
di Chimica, Università degli Studi di Milano, via Golgi
19, 20133 Milan, Italy
- Istituto di Scienze
e Tecnologie Molecolari, CNR-ISTM, via Golgi 19, 20133 Milan, Italy
| |
Collapse
|
5
|
Bonfanti M, Tantardini GF, Hughes KH, Martinazzo R, Burghardt I. Compact MCTDH Wave Functions for High-Dimensional System-Bath Quantum Dynamics. J Phys Chem A 2012; 116:11406-13. [DOI: 10.1021/jp3064504] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- M. Bonfanti
- Dipartimento di
Chimica, Università degli Studi di Milano, v. Golgi 19, 20133 Milano, Italy
| | - G. F. Tantardini
- Dipartimento di
Chimica, Università degli Studi di Milano, v. Golgi 19, 20133 Milano, Italy
| | - K. H. Hughes
- School of Chemistry, Bangor University, Bangor, Gwynedd LL57 2UW, United
Kingdom
| | - R. Martinazzo
- Dipartimento di
Chimica, Università degli Studi di Milano, v. Golgi 19, 20133 Milano, Italy
| | - I. Burghardt
- Institute of Physical and Theoretical
Chemistry, Goethe University Frankfurt,
Max-von-Laue-Strasse 7, 60438, Frankfurt, Germany
| |
Collapse
|
6
|
Abstract
We investigated binding of hydrogen atoms to small polycyclic aromatic hydrocarbons (PAHs)--i.e., graphene dots with hydrogen-terminated edges--using density functional theory and correlated wavefunction techniques. We considered a number of PAHs with three to seven hexagonal rings and computed binding energies for most of the symmetry unique sites, along with the minimum energy paths for significant cases. The chosen PAHs are small enough to not present radical character at their edges, yet show a clear preference for adsorption at the edge sites which can be attributed to electronic effects. We show how the results, as obtained at different levels of theory, can be rationalized in detail with the help of a few simple concepts derivable from a tight-binding model of the π electrons.
Collapse
Affiliation(s)
- M Bonfanti
- Dipartimento di Chimica Fisica ed Elettrochimica, Università degli Studi di Milano, v. Golgi 19, 20133 Milan, Italy
| | | | | | | | | |
Collapse
|
7
|
|
8
|
Martinazzo R, Hughes KH, Burghardt I. Unraveling a Brownian particle's memory with effective mode chains. Phys Rev E Stat Nonlin Soft Matter Phys 2011; 84:030102. [PMID: 22060317 DOI: 10.1103/physreve.84.030102] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/20/2011] [Indexed: 05/31/2023]
Abstract
Memory effects in quantum dynamical processes involving structured environments are presently difficult, if not impossible, to investigate using standard approaches. Progress can be made by transforming the environmental variables to a suitable chain representation which effectively performs a Markovian embedding of the dynamics. Here, we show that this effective-mode chain representation provides a unique way of unraveling the memory kernel κ(t) as a function of time. Truncated or Markov-closed chains with n effective modes exactly reproduce κ(t) to the 4nth order in time, up to an irrelevant constant of order κ(0)/n. These favorable convergence properties pave the way for efficient quantum simulations of fast (non-Markovian) processes by reduced dynamical models.
Collapse
Affiliation(s)
- R Martinazzo
- Dipartimento di Chimica-Fisica ed Elettrochimica, Università degli Studi di Milano, v. Golgi 19, I-20133 Milano, Italy.
| | | | | |
Collapse
|
9
|
López-López S, Martinazzo R, Nest M. Benchmark calculations for dissipative dynamics of a system coupled to an anharmonic bath with the multiconfiguration time-dependent Hartree method. J Chem Phys 2011; 134:094102. [DOI: 10.1063/1.3556940] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
10
|
Martinazzo R, Vacchini B, Hughes KH, Burghardt I. Communication: Universal Markovian reduction of Brownian particle dynamics. J Chem Phys 2011; 134:011101. [DOI: 10.1063/1.3532408] [Citation(s) in RCA: 88] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
11
|
Martinazzo R, Assoni S, Marinoni G, Tantardini GF. Hot-atom versus Eley–Rideal dynamics in hydrogen recombination on Ni(100). I. The single-adsorbate case. J Chem Phys 2004; 120:8761-71. [PMID: 15267808 DOI: 10.1063/1.1695316] [Citation(s) in RCA: 33] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We compare the efficiency of the Eley-Rideal (ER) reaction with the formation of hot-atom (HA) species in the simplest case, i.e., the scattering of a projectile off a single adsorbate, considering the Hydrogen and Hydrogen-on-Ni(100) system. We use classical mechanics and the accurate embedded diatomics-in-molecules potential to study the collision system over a wide range of collision energies (0.10-1.50 eV), both with a rigid and a nonrigid Ni substrate and for impact on the occupied and neighboring empty cells. In the rigid model metastable and truly bound hot-atoms occur and we find that the cross section for the formation of bound hot-atoms is considerably higher than that for the ER reaction over the whole range of collision energies examined. Metastable hot-atoms form because of the inefficient energy transfer to the adsorbate and have lifetimes of the order 0.1-0.7 ps, depending on the collision energy. When considering the effects of lattice vibrations we find, on average, a consistent energy transfer to the substrate, say 0.1-0.2 eV, which forced us to devise a two-step dynamical model to get rid of the problems associated with the use of periodic boundary conditions. Results for long-lived HA formation due to scattering on the occupied cell at a surface temperature of 120 K agree well with those of the rigid model, suggesting that in the above process the substrate plays only a secondary role and further calculations at surface temperatures of 50 and 300 K are in line with these findings. However, considerably high cross sections for formation of long-lived hot-atoms result also from scattering off the neighboring cells where the energy transfer to the lattice cannot be neglected. Metastable hot-atoms are reduced in number and have usually lifetimes shorter than those of the rigid-model, say less than 0.3 ps. In addition, ER cross sections are only slightly affected by the lattice motion and show a little temperature dependence. Finally, we find also that absorption and reflection strongly depend on the correct consideration of lattice vibrations and the occurrence of trapping.
Collapse
Affiliation(s)
- R Martinazzo
- Department of Physical Chemistry and Electrochemistry, University of Milan, V. Golgi 19, 20133 Milan, Italy
| | | | | | | |
Collapse
|
12
|
Calderoni G, Cargnoni F, Martinazzo R, Raimondi M. Potential energy surface, bound states, and rotational inelastic cross sections of the He-CH[sub 4] system: A theoretical investigation. J Chem Phys 2004; 121:8261-70. [PMID: 15511146 DOI: 10.1063/1.1791111] [Citation(s) in RCA: 17] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We determined two potential energy surfaces (PES) for the He-CH(4) system by means of MP4 and Valence Bond (VB) calculations. The MP4 potential is similar to the one commonly adopted for this system [U. Buck, K. H. Kohl, A. Kolhase, M. Faubel, and U. Staemmler, Mol. Phys. 55, 1255 (1985)], while the VB PES is slightly more attractive. To evaluate the reliability of these potentials, we investigated the scattering properties by performing close coupling calculations, and concluded that: (i) the available experimental data do not permit the ranking among the PES considered; (ii) some theoretical predictions differ considerably from the experimental data, and these discrepancies cannot entirely be ascribed to the inaccuracy of the ab initio calculations; (iii) the scattering properties at low energy might discriminate between the MP4 and VB potentials.
Collapse
Affiliation(s)
- G Calderoni
- Dipartimento di Chimica Fisica ed Elettrochimica, Università degli Studi di Milano, Via Golgi 19, 20133 Milan, Italy
| | | | | | | |
Collapse
|
13
|
Martinazzo R, Tantardini GF, Bodo E, Gianturco FA. Accurate potential energy surfaces for the study of lithium–hydrogen ionic reactions. J Chem Phys 2003. [DOI: 10.1063/1.1621852] [Citation(s) in RCA: 62] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
14
|
Satta M, Bodo E, Martinazzo R, Gianturco FA. Photoexcitation of LiH2+ from selected initial states: A time-dependent model. J Chem Phys 2002. [DOI: 10.1063/1.1482695] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023] Open
|
15
|
Bodo E, Gianturco FA, Martinazzo R, Raimondi M. Reactive Behavior of the [LiH2]+ System I. Evaluation of the Lower-lying Electronic Potentials for the Collinear Geometries. J Phys Chem A 2001. [DOI: 10.1021/jp0123435] [Citation(s) in RCA: 30] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Affiliation(s)
- E. Bodo
- University of Rome “La Sapienza”, Department of Chemistry, P.le A. Moro 5, 00185, Roma, and University of Milan, Department of Physical Chemistry and Electrochemistry and CNR center CSR SRC, V. Golgi 19, 20133, Milan
| | - F. A. Gianturco
- University of Rome “La Sapienza”, Department of Chemistry, P.le A. Moro 5, 00185, Roma, and University of Milan, Department of Physical Chemistry and Electrochemistry and CNR center CSR SRC, V. Golgi 19, 20133, Milan
| | - R. Martinazzo
- University of Rome “La Sapienza”, Department of Chemistry, P.le A. Moro 5, 00185, Roma, and University of Milan, Department of Physical Chemistry and Electrochemistry and CNR center CSR SRC, V. Golgi 19, 20133, Milan
| | - M. Raimondi
- University of Rome “La Sapienza”, Department of Chemistry, P.le A. Moro 5, 00185, Roma, and University of Milan, Department of Physical Chemistry and Electrochemistry and CNR center CSR SRC, V. Golgi 19, 20133, Milan
| |
Collapse
|
16
|
Bodo E, Gianturco FA, Martinazzo R. Reactive Behavior of the [LiH2]+ System II. Collision-Induced Dissociation and Collinear Reaction Dynamics of LiH++H from Quantum Time Dependent Calculations. J Phys Chem A 2001. [DOI: 10.1021/jp012344x] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Affiliation(s)
- E. Bodo
- University of Rome “La Sapienza”, Department of Chemistry, P.le Aldo Moro 5, 00185, Rome, Italy, and University of Milan, Department of Physical Chemistry and Electrochemistry and CNR center CSR SRC, V. Golgi 19, 20133, Milan
| | - F. A. Gianturco
- University of Rome “La Sapienza”, Department of Chemistry, P.le Aldo Moro 5, 00185, Rome, Italy, and University of Milan, Department of Physical Chemistry and Electrochemistry and CNR center CSR SRC, V. Golgi 19, 20133, Milan
| | - R. Martinazzo
- University of Rome “La Sapienza”, Department of Chemistry, P.le Aldo Moro 5, 00185, Rome, Italy, and University of Milan, Department of Physical Chemistry and Electrochemistry and CNR center CSR SRC, V. Golgi 19, 20133, Milan
| |
Collapse
|
17
|
Bodo E, Gianturco F, Martinazzo R, Raimondi M. Possible reaction paths in the LiH+2 chemistry: a computational analysis of the interaction forces. Chem Phys 2001. [DOI: 10.1016/s0301-0104(01)00424-4] [Citation(s) in RCA: 36] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
18
|
Bodo E, Gianturco FA, Martinazzo R, Paesani F, Raimondi M. Testing van der Waals interactions with quantum dynamics: Repulsive anisotropy and well depth in the LiH+He system. J Chem Phys 2000. [DOI: 10.1063/1.1311801] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
19
|
Bodo E, Gianturco FA, Martinazzo R, Forni A, Raimondi M. Spatial Energetics of Protonated LiH: Lower-Lying Potential Energy Surfaces from Valence Bond Calculations. J Phys Chem A 2000. [DOI: 10.1021/jp0022510] [Citation(s) in RCA: 29] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|