1
|
Kim SS, Rhee YM. Potential energy interpolation with target-customized weighting coordinates: application to excited-state dynamics of photoactive yellow protein chromophore in water. Phys Chem Chem Phys 2024; 26:9021-9036. [PMID: 38440829 DOI: 10.1039/d3cp05643k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/06/2024]
Abstract
Interpolation of potential energy surfaces (PESs) can provide a practical route to performing molecular dynamics simulations with a reliability matching a high-level quantum chemical calculation. An obstacle to its widespread use is perhaps the lack of general and optimal interpolation settings that can be applied in a black-box manner for any given molecular system. How to set up the weights for interpolation is one such task, and we still need to diversify the approaches in order to treat various systems. Here, we develop a new interpolation weighting scheme, which allows us to choose the weighting coordinates in a system-specific manner, by amplifying the contribution from specific internal coordinates. The new weighting scheme with an appropriate selection of coordinates is proved to be effective in reducing the interpolation error along the reaction pathway. As a demonstration, we consider the photoactive yellow protein chromophore system, as it constitutes itself as an interesting target that bears long-standing questions related to excited-state dynamics inside protein environments. We build its two-state diabatic interpolated PES with the new weighting scheme. We indeed see the utility of our scheme by conducting nonadiabatic molecular dynamics simulations with the required semi-global PES based on a limited number of data points.
Collapse
Affiliation(s)
- Seung Soo Kim
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea.
| | - Young Min Rhee
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea.
| |
Collapse
|
2
|
Ootani Y, Satoh A, Harabuchi Y, Taketsugu T. Trajectory on-the-fly molecular dynamics approach to tunneling splitting in the electronic excited state: A case of tropolone. J Comput Chem 2020; 41:1549-1556. [PMID: 32239685 DOI: 10.1002/jcc.26199] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2020] [Revised: 03/10/2020] [Accepted: 03/16/2020] [Indexed: 11/10/2022]
Abstract
The semiclassical tunneling method is applied to evaluate the tunneling splitting of tropolone due to the intramolecular proton transfer in the electronic excited state, first time, in a framework of the trajectory on-the-fly molecular dynamics (TOF-MD) approach. To prevent unphysical zero-point vibrational energy transfer among the normal modes of vibration, quantum zero-point vibrational energies are assigned only to the vibrational modes related to intramolecular proton transfer, whereas the remaining modes are treated as bath modes. Practical ways to determine the tunnel-initiating points and tunneling path are introduced. It is shown that the tunneling splitting decreases as the bath-mode energy increases. The experimental tunneling splitting value is well reproduced by the present TOF-MD approach based on the Wentzel-Kramers-Brillouin (WKB) approximation.
Collapse
Affiliation(s)
- Yusuke Ootani
- Institute for Materials Research, Tohoku University, Sendai, Japan
| | - Aya Satoh
- Department of Chemistry, Faculty of Science, Hokkaido University, Sapporo, Japan
| | - Yu Harabuchi
- Department of Chemistry, Faculty of Science, Hokkaido University, Sapporo, Japan
| | - Tetsuya Taketsugu
- Department of Chemistry, Faculty of Science, Hokkaido University, Sapporo, Japan.,Institute for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Sapporo, Japan
| |
Collapse
|
3
|
Lin SY, Zhang P, Zhang JZ. Hybrid many-body-expansion/Shepard-interpolation method for constructing ab initio potential energy surfaces for quantum dynamics calculations. Chem Phys Lett 2013. [DOI: 10.1016/j.cplett.2012.11.089] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
|
4
|
Ootani Y, Taketsugu T. Ab initio molecular dynamics approach to tunneling splitting in polyatomic molecules. J Comput Chem 2011; 33:60-5. [DOI: 10.1002/jcc.21943] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2011] [Revised: 08/15/2011] [Accepted: 08/17/2011] [Indexed: 11/08/2022]
|
5
|
|
6
|
Matanović I, Doslić N, Kühn O. Ground and asymmetric CO-stretch excited state tunneling splittings in the formic acid dimer. J Chem Phys 2007; 127:014309. [PMID: 17627348 DOI: 10.1063/1.2748048] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Abstract
There has been some controversy concerning the assignment of measured tunneling splittings for the formic acid dimer in the vibrational ground state and the asymmetric CO-stretching excited state. The discussion is intimately related to the question whether the fundamental excitation of the CO-vibration promotes or hinders tunneling. Here we will address this issue on the basis of a five-dimensional reaction space Hamiltonian which includes three large amplitude coordinates as well as two harmonic modes whose linear superposition reproduces the asymmetric CO-vibrational mode. Within density functional theory using the B3LYP functional together with a 6-311++G(3df,3pd) basis set we obtain a ground state tunneling splitting which is about 2.4 larger than the one for the CO-stretching excited state.
Collapse
Affiliation(s)
- I Matanović
- Department of Physical Chemistry, Rudjer Bosković Institute, 10000 Zagreb, Croatia
| | | | | |
Collapse
|
7
|
Tishchenko O, Truhlar DG. Global Potential Energy Surfaces with Correct Permutation Symmetry by Multiconfiguration Molecular Mechanics. J Chem Theory Comput 2007; 3:938-48. [DOI: 10.1021/ct600315h] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Oksana Tishchenko
- Chemistry Department and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431
| | - Donald G. Truhlar
- Chemistry Department and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431
| |
Collapse
|
8
|
Miani A, Carloni P, Raugei S. A novel parametrization scheme for classical and quantum mechanical simulations of large, floppy molecular systems. Chem Phys Lett 2006. [DOI: 10.1016/j.cplett.2006.06.063] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
9
|
Mil'nikov GV, Yagi K, Taketsugu T, Nakamura H, Hirao K. Simple and accurate method to evaluate tunneling splitting in polyatomic molecules. J Chem Phys 2006; 120:5036-45. [PMID: 15267369 DOI: 10.1063/1.1647052] [Citation(s) in RCA: 56] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Abstract
A practical and accurate semiclassical method for calculating the tunneling splitting of the ground state in polyatomic molecules is presented based on a recent version of the instanton theory [J. Chem. Phys. 115, 6881 (2001)]. The method uses ab initio quantum chemical data for the potential energy surface without any concomitant extrapolation and requires only a small number of ab initio data points to get convergence even for large molecules. This enables one to use an advanced level of electronic structure theory and achieve a high accuracy of the result. The method is applied to the 9-atomic malonaldehyde molecule by making use of the potential energy surface at the level of CCSD(T) with the hybrid basis set of aug-cc-pVTZ (for oxygen atoms and the transferred hydrogen atom) and cc-pVTZ (for other atoms).
Collapse
Affiliation(s)
- Gennady V Mil'nikov
- Department of Theoretical Studies, Institute for Molecular Science, Myodaiji, Okazaki 444-8585, Japan.
| | | | | | | | | |
Collapse
|
10
|
Giese K, Kühn O. The all-Cartesian reaction plane Hamiltonian: Formulation and application to the H-atom transfer in tropolone. J Chem Phys 2005; 123:054315. [PMID: 16108647 DOI: 10.1063/1.1978869] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this work we present an all-Cartesian reaction surface approach, where the large amplitude coordinates span the so-called reaction plane, that is, the unique plane defined by the two minima and the saddle-point structure of an isomerization reaction. Orthogonal modes are treated within harmonic approximation which gives the total Hamiltonian an almost separable form that is suitable for multidimensional quantum dynamics calculations. The reaction plane Hamiltonian is constructed for the H-atom transfer in tropolone as an example for a system with an intramolecular O...H-O hydrogen bond. We find ground-state tunneling splittings of 3.5 and 0.16 cm(-1) for the normal and deuterated species, respectively. We calculated infrared-absorption spectra for a four-dimensional model focusing on the low-frequency region. Here, we identify a reaction mode which is closely connected to the tautomerization that is reflected in the increase of tunneling splitting to 18 cm(-1) upon excitation.
Collapse
Affiliation(s)
- Kai Giese
- Institut für Chemie, Physikalische und Theoretische Chemie, Freie Universität Berlin, Takustrasse 3, D-14195 Berlin, Germany
| | | |
Collapse
|
11
|
Ishida T, Schatz GC. A local interpolation scheme using no derivatives in potential sampling: application to O(1D) + H2 system. J Comput Chem 2003; 24:1077-86. [PMID: 12759907 DOI: 10.1002/jcc.10252] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
We recently proposed a local interpolation scheme, in which interpolant moving least squares (IMLS) and Shepard interpolation are employed to describe potential energy surfaces. This IMLS/Shepard scheme is used to interpolate quantum chemical potential energy surfaces for which analytical derivatives are not available. In this study, we apply the scheme to the highly exothermic O((1)D) + H(2) --> H + OH reaction and compare it with results based on Shepard interpolation using second-order Taylor expansions. An analytical surface is used to define the potential function so that errors in the interpolation function may accurately be determined. We find that the present scheme reproduces the correct reactive cross-sections more accurately than the Shepard scheme, and with rms errors for energy and gradients that are significantly smaller than those from Shepard interpolation. This occurs even though the present scheme does not utilize derivative and Hessian information, whereas the Shepard interpolation does. The Bayesian approach proposed by Bettens and Collins does not improve the IMLS/Shepard results significantly, although it does the Shepard-only approach. The accuracy of the IMLS/Shepard scheme is surprising, but can be explained by the more global nature of the interpolation.
Collapse
Affiliation(s)
- Toshimasa Ishida
- Research Center for Molecular-scale Nanoscience Center, Institute for Molecular Science, Okazaki 444-8585, Japan.
| | | |
Collapse
|
12
|
Yagi K, Taketsugu T, Hirao K. Generation of full-dimensional potential energy surface of intramolecular hydrogen atom transfer in malonaldehyde and tunneling dynamics. J Chem Phys 2001. [DOI: 10.1063/1.1418436] [Citation(s) in RCA: 90] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
13
|
Hoon Choi T, Tae Park S, Soo Kim M. Theoretical and experimental studies of the dissociation dynamics of methaniminium cation, CH2NH2+→CHNH++H2: Reaction path bifurcation. J Chem Phys 2001. [DOI: 10.1063/1.1355309] [Citation(s) in RCA: 16] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
14
|
Yagi K, Taketsugu T, Hirao K, Gordon MS. Direct vibrational self-consistent field method: Applications to H2O and H2CO. J Chem Phys 2000. [DOI: 10.1063/1.481881] [Citation(s) in RCA: 145] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
15
|
Collins MA, Zhang DH. Application of interpolated potential energy surfaces to quantum reactive scattering. J Chem Phys 1999. [DOI: 10.1063/1.480344] [Citation(s) in RCA: 71] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|