1
|
Alaal N, Brorsen KR. Multicomponent heat-bath configuration interaction with the perturbative correction for the calculation of protonic excited states. J Chem Phys 2021; 155:234107. [PMID: 34937361 DOI: 10.1063/5.0076006] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this study, we extend the multicomponent heat-bath configuration interaction (HCI) method to excited states. Previous multicomponent HCI studies have been performed using only the variational stage of the HCI algorithm as they have largely focused on the calculation of protonic densities. Because this study focuses on energetic quantities, a second-order perturbative correction after the variational stage is essential. Therefore, this study implements the second-order Epstein-Nesbet correction to the variational stage of multicomponent HCI for the first time. Additionally, this study introduces a new procedure for calculating reference excitation energies for multicomponent methods using the Fourier-grid Hamiltonian (FGH) method, which should allow the one-particle electronic basis set errors to be better isolated from errors arising from an incomplete description of electron-proton correlation. The excited-state multicomponent HCI method is benchmarked by computing protonic excitations of the HCN and FHF- molecules and is shown to be of similar accuracy to previous excited-state multicomponent methods such as the multicomponent time-dependent density-functional theory and equation-of-motion coupled-cluster theory relative to the new FGH reference values.
Collapse
Affiliation(s)
- Naresh Alaal
- Department of Chemistry, University of Missouri, Columbia, Missouri 65203, USA
| | - Kurt R Brorsen
- Department of Chemistry, University of Missouri, Columbia, Missouri 65203, USA
| |
Collapse
|
2
|
Witzorky C, Paramonov G, Bouakline F, Jaquet R, Saalfrank P, Klamroth T. Gaussian-Type Orbital Calculations for High Harmonic Generation in Vibrating Molecules: Benchmarks for H 2. J Chem Theory Comput 2021; 17:7353-7365. [PMID: 34747605 DOI: 10.1021/acs.jctc.1c00837] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The response of the hydrogen molecular ion, H2+, to few-cycle laser pulses of different intensities is simulated. To treat the coupled electron-nuclear motion, we use adiabatic potentials computed with Gaussian-type basis sets together with a heuristic ionization model for the electron and a grid representation for the nuclei. Using this mixed-basis approach, the time-dependent Schrödinger equation is solved, either within the Born-Oppenheimer approximation or with nonadiabatic couplings included. The dipole response spectra are compared to all-grid-based solutions for the three-body problem, which we take as a reference to benchmark the Gaussian-type basis set approaches. Also, calculations employing the fixed-nuclei approximation are performed, to quantify effects due to nuclear motion. For low intensities and small ionization probabilities, we get excellent agreement of the dynamics using Gaussian-type basis sets with the all-grid solutions. Our investigations suggest that high harmonic generation (HHG) and high-frequency response, in general, can be reliably modeled using Gaussian-type basis sets for the electrons for not too high harmonics. Further, nuclear motion destroys electronic coherences in the response spectra even on the time scale of about 30 fs and affects HHG intensities, which reflect the electron dynamics occurring on the attosecond time scale. For the present system, non-Born-Oppenheimer effects are small. The Gaussian-based, nonadiabatically coupled, time-dependent multisurface approach to treat quantum electron-nuclear motion beyond the non-Born-Oppenheimer approximation can be easily extended to approximate wavefunction methods, such as time-dependent configuration interaction singles (TD-CIS), for systems where no benchmarks are available.
Collapse
Affiliation(s)
- Christoph Witzorky
- Institut für Chemie, Universität Potsdam, Karl-Liebknecht-Str. 24-25, D-14476 Potsdam-Golm, Germany
| | - Guennaddi Paramonov
- Institut für Chemie, Universität Potsdam, Karl-Liebknecht-Str. 24-25, D-14476 Potsdam-Golm, Germany
| | - Foudhil Bouakline
- Institut für Chemie, Universität Potsdam, Karl-Liebknecht-Str. 24-25, D-14476 Potsdam-Golm, Germany
| | - Ralph Jaquet
- Theoretische Chemie, Universität Siegen, D-57068 Siegen, Germany
| | - Peter Saalfrank
- Institut für Chemie, Universität Potsdam, Karl-Liebknecht-Str. 24-25, D-14476 Potsdam-Golm, Germany
| | - Tillmann Klamroth
- Institut für Chemie, Universität Potsdam, Karl-Liebknecht-Str. 24-25, D-14476 Potsdam-Golm, Germany
| |
Collapse
|
3
|
McLaughlin PF, Chakraborty A. Compact Real-Space Representation of Excited States Using Frequency-Dependent Explicitly Correlated Electron–Hole Interaction Kernel. J Chem Theory Comput 2020; 16:5762-5770. [DOI: 10.1021/acs.jctc.9b01238] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Peter F. McLaughlin
- Department of Chemistry, Syracuse University, Syracuse, New York 13244 United States
| | - Arindam Chakraborty
- Department of Chemistry, Syracuse University, Syracuse, New York 13244 United States
| |
Collapse
|
4
|
Xu X, Yang Y. Full-quantum descriptions of molecular systems from constrained nuclear–electronic orbital density functional theory. J Chem Phys 2020; 153:074106. [DOI: 10.1063/5.0014001] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023] Open
Affiliation(s)
- Xi Xu
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, 1101 University Avenue, Madison, Wisconsin 53706, USA
| | - Yang Yang
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, 1101 University Avenue, Madison, Wisconsin 53706, USA
| |
Collapse
|
5
|
Pavošević F, Culpitt T, Hammes-Schiffer S. Multicomponent Quantum Chemistry: Integrating Electronic and Nuclear Quantum Effects via the Nuclear–Electronic Orbital Method. Chem Rev 2020; 120:4222-4253. [DOI: 10.1021/acs.chemrev.9b00798] [Citation(s) in RCA: 53] [Impact Index Per Article: 13.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Fabijan Pavošević
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Tanner Culpitt
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Sharon Hammes-Schiffer
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| |
Collapse
|
6
|
Xu X, Yang Y. Constrained nuclear-electronic orbital density functional theory: Energy surfaces with nuclear quantum effects. J Chem Phys 2020; 152:084107. [DOI: 10.1063/1.5143371] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Affiliation(s)
- Xi Xu
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin-Madison, 1101 University Avenue, Madison, Wisconsin 53706, USA
| | - Yang Yang
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin-Madison, 1101 University Avenue, Madison, Wisconsin 53706, USA
| |
Collapse
|
7
|
Brorsen KR. Quantifying Multireference Character in Multicomponent Systems with Heat-Bath Configuration Interaction. J Chem Theory Comput 2020; 16:2379-2388. [DOI: 10.1021/acs.jctc.9b01273] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Affiliation(s)
- Kurt R. Brorsen
- Department of Chemistry, University of Missouri, Columbia, Missouri 65203, United States
| |
Collapse
|
8
|
Patel P, Wilson AK. Domain-based local pair natural orbital methods within the correlation consistent composite approach. J Comput Chem 2019; 41:800-813. [PMID: 31891196 DOI: 10.1002/jcc.26129] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2019] [Revised: 10/29/2019] [Accepted: 12/01/2019] [Indexed: 01/15/2023]
Abstract
Ab initio composite approaches have been utilized to model and predict main group thermochemistry within 1 kcal mol-1 , on average, from well-established reliable experiments, primarily for molecules with less than 30 atoms. For molecules of increasing size and complexity, such as biomolecular complexes, composite methodologies have been limited in their application. Therefore, the domain-based local pair natural orbital (DLPNO) methods have been implemented within the correlation consistent composite approach (ccCA) framework, namely DLPNO-ccCA, to reduce the computational cost (disk space, CPU (central processing unit) time, memory) and predict energetic properties such as enthalpies of formation, noncovalent interactions, and conformation energies for organic biomolecular complexes including one of the largest molecules examined via composite strategies, within 1 kcal mol-1 , after calibration with 119 molecules and a set of linear alkanes. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Prajay Patel
- Department of Chemistry, Michigan State University, East Lansing, Michigan, 48824
| | - Angela K Wilson
- Department of Chemistry, Michigan State University, East Lansing, Michigan, 48824
| |
Collapse
|
9
|
Pavošević F, Hammes-Schiffer S. Multicomponent coupled cluster singles and doubles and Brueckner doubles methods: Proton densities and energies. J Chem Phys 2019; 151:074104. [PMID: 31438716 DOI: 10.1063/1.5116113] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The nuclear-electronic orbital (NEO) framework enables computationally practical coupled cluster calculations of multicomponent molecular systems, in which all electrons and specified nuclei, typically protons, are treated quantum mechanically. In addition to energies, computing accurate proton densities is essential for the calculation of reliable molecular properties, including vibrationally averaged geometries and vibrational frequencies. Herein, the Lagrangian formalism for the multicomponent coupled cluster with single and double excitations (NEO-CCSD) method is derived and implemented. The multicomponent coupled cluster with double excitations method using optimized Brueckner orbitals, denoted as NEO-BCCD, is also developed. Both of these methods are used to compute the proton densities for two molecular systems. The results illustrate that orbital relaxation effects, which can be included either indirectly with the NEO-CCSD method or directly with the NEO-BCCD method, are critical for computing even qualitatively accurate proton densities. Both methods are also able to provide accurate proton affinities and vibrationally averaged optimized geometries. This Lagrangian formalism will enable the calculation of other properties such as analytical nuclear gradients and Hessians with NEO coupled cluster methods. Moreover, the accuracy of these methods may be improved systematically by the inclusion of higher-order excitations. Thus, this work provides the foundation for a wide range of future methodological developments and applications within the NEO framework.
Collapse
Affiliation(s)
- Fabijan Pavošević
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, USA
| | - Sharon Hammes-Schiffer
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, USA
| |
Collapse
|
10
|
Mejía-Rodríguez D, de la Lande A. Multicomponent density functional theory with density fitting. J Chem Phys 2019; 150:174115. [DOI: 10.1063/1.5078596] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Daniel Mejía-Rodríguez
- Laboratoire de Chimie Physique, Université Paris Sud/CNRS, Université Paris Saclay, 15 Avenue Jean Perrin, 91405 Orsay, France
| | - Aurélien de la Lande
- Laboratoire de Chimie Physique, Université Paris Sud/CNRS, Université Paris Saclay, 15 Avenue Jean Perrin, 91405 Orsay, France
| |
Collapse
|
11
|
Pavošević F, Culpitt T, Hammes-Schiffer S. Multicomponent Coupled Cluster Singles and Doubles Theory within the Nuclear-Electronic Orbital Framework. J Chem Theory Comput 2018; 15:338-347. [DOI: 10.1021/acs.jctc.8b01120] [Citation(s) in RCA: 45] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Fabijan Pavošević
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Tanner Culpitt
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Sharon Hammes-Schiffer
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| |
Collapse
|
12
|
Affiliation(s)
- Edit Mátyus
- Institute of Chemistry, ELTE Eötvös Loránd University, Budapest, Hungary
| |
Collapse
|
13
|
Brorsen KR, Schneider PE, Hammes-Schiffer S. Alternative forms and transferability of electron-proton correlation functionals in nuclear-electronic orbital density functional theory. J Chem Phys 2018; 149:044110. [DOI: 10.1063/1.5037945] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Kurt R. Brorsen
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Ave., Urbana, Illinois 61801, USA
| | - Patrick E. Schneider
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, USA
| | - Sharon Hammes-Schiffer
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, USA
| |
Collapse
|
14
|
Bayne MG, Scher JA, Ellis BH, Chakraborty A. Linked-Cluster Formulation of Electron–Hole Interaction Kernel in Real-Space Representation without Using Unoccupied States. J Chem Theory Comput 2018; 14:3656-3666. [DOI: 10.1021/acs.jctc.8b00123] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Michael G. Bayne
- Department of Chemistry, Syracuse University, Syracuse, New York 13244 United States
| | - Jeremy A. Scher
- Department of Chemistry, Syracuse University, Syracuse, New York 13244 United States
| | - Benjamin H. Ellis
- Department of Chemistry, Syracuse University, Syracuse, New York 13244 United States
| | - Arindam Chakraborty
- Department of Chemistry, Syracuse University, Syracuse, New York 13244 United States
| |
Collapse
|
15
|
Goli M, Shahbazian S. Developing effective electronic-only coupled-cluster and Møller-Plesset perturbation theories for the muonic molecules. Phys Chem Chem Phys 2018; 20:16749-16760. [PMID: 29881845 DOI: 10.1039/c8cp02489h] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Recently we have proposed an effective Hartree-Fock (EHF) theory for the electrons of the muonic molecules that is formally equivalent to the HF theory within the context of the nuclear-electronic orbital theory [Phys. Chem. Chem. Phys., 2018, 20, 4466]. In the present report we extend the muon-specific effective electronic structure theory beyond the EHF level by introducing the effective second order Møller-Plesset perturbation theory (EMP2) and the effective coupled-cluster theory at single and double excitation levels (ECCSD) as well as an improved version including perturbative triple excitations (ECCSD(T)). These theories incorporate electron-electron correlation into the effective paradigm and through their computational implementation, a diverse set of small muonic species is considered as a benchmark at these post-EHF levels. A comparative computational study on this set demonstrates that the muonic bond length is in general non-negligibly longer than corresponding hydrogenic analogs. Next, the developed post-EHF theories are applied for the muoniated N-heterocyclic carbene/silylene/germylene and the muoniated triazolium cation revealing the relative stability of the sticking sites of the muon in each species. The computational results, in line with previously reported experimental data demonstrate that the muon generally prefers to attach to the divalent atom with carbeneic nature. A detailed comparison of these muonic adducts with the corresponding hydrogenic adducts reveals subtle differences that have already been overlooked.
Collapse
Affiliation(s)
- Mohammad Goli
- School of Nano Science, Institute for Research in Fundamental Sciences (IPM), 19395-5531, Tehran, Iran.
| | | |
Collapse
|
16
|
Yang Y, Culpitt T, Hammes-Schiffer S. Multicomponent Time-Dependent Density Functional Theory: Proton and Electron Excitation Energies. J Phys Chem Lett 2018; 9:1765-1770. [PMID: 29553738 DOI: 10.1021/acs.jpclett.8b00547] [Citation(s) in RCA: 43] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
The quantum mechanical treatment of both electrons and protons in the calculation of excited state properties is critical for describing nonadiabatic processes such as photoinduced proton-coupled electron transfer. Multicomponent density functional theory enables the consistent quantum mechanical treatment of more than one type of particle and has been implemented previously for studying ground state molecular properties within the nuclear-electronic orbital (NEO) framework, where all electrons and specified protons are treated quantum mechanically. To enable the study of excited state molecular properties, herein the linear response multicomponent time-dependent density functional theory (TDDFT) is derived and implemented within the NEO framework. Initial applications to FHF- and HCN illustrate that NEO-TDDFT provides accurate proton and electron excitation energies within a single calculation. As its computational cost is similar to that of conventional electronic TDDFT, the NEO-TDDFT approach is promising for diverse applications, particularly nonadiabatic proton transfer reactions, which may exhibit mixed electron-proton vibronic excitations.
Collapse
Affiliation(s)
- Yang Yang
- Department of Chemistry , Yale University , 225 Prospect Street , New Haven , Connecticut 06520 , United States
| | - Tanner Culpitt
- Department of Chemistry , Yale University , 225 Prospect Street , New Haven , Connecticut 06520 , United States
| | - Sharon Hammes-Schiffer
- Department of Chemistry , Yale University , 225 Prospect Street , New Haven , Connecticut 06520 , United States
| |
Collapse
|
17
|
Rayka M, Goli M, Shahbazian S. Toward a muon-specific electronic structure theory: effective electronic Hartree-Fock equations for muonic molecules. Phys Chem Chem Phys 2018; 20:4466-4477. [PMID: 29372727 DOI: 10.1039/c7cp07599e] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
An effective set of Hartree-Fock (HF) equations are derived for electrons of muonic systems, i.e., molecules containing a positively charged muon, conceiving the muon as a quantum oscillator, which are completely equivalent to the usual two-component HF equations used to derive stationary states of the muonic molecules. In these effective equations, a non-Coulombic potential is added to the orthodox coulomb and exchange potential energy terms, which describes the interaction of the muon and the electrons effectively and is optimized during the self-consistent field cycles. While in the two-component HF equations a muon is treated as a quantum particle, in the effective HF equations it is absorbed into the effective potential and practically transformed into an effective potential field experienced by electrons. The explicit form of the effective potential depends on the nature of muon's vibrations and is derivable from the basis set used to expand the muonic spatial orbital. The resulting effective Hartree-Fock equations are implemented computationally and used successfully, as a proof of concept, in a series of muonic molecules containing all atoms from the second and third rows of the Periodic Table. To solve the algebraic version of the equations muon-specific Gaussian basis sets are designed for both muon and surrounding electrons and it is demonstrated that the optimized exponents are quite distinct from those derived for the hydrogen isotopes. The developed effective HF theory is quite general and in principle can be used for any muonic system while it is the starting point for a general effective electronic structure theory that incorporates various types of quantum correlations into the muonic systems beyond the HF equations.
Collapse
Affiliation(s)
- Milad Rayka
- Department of Physics and Department of Physical and Computational Chemistry, Shahid Beheshti University, G. C., Evin, P.O. Box 19395-4716, Tehran, 19839, Iran.
| | | | | |
Collapse
|
18
|
Yang Y, Brorsen KR, Culpitt T, Pak MV, Hammes-Schiffer S. Development of a practical multicomponent density functional for electron-proton correlation to produce accurate proton densities. J Chem Phys 2017; 147:114113. [DOI: 10.1063/1.4996038] [Citation(s) in RCA: 55] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yang Yang
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Ave., Urbana, Illinois 61801, USA
| | - Kurt R. Brorsen
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Ave., Urbana, Illinois 61801, USA
| | - Tanner Culpitt
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Ave., Urbana, Illinois 61801, USA
| | - Michael V. Pak
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Ave., Urbana, Illinois 61801, USA
| | - Sharon Hammes-Schiffer
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Ave., Urbana, Illinois 61801, USA
| |
Collapse
|
19
|
Brorsen KR, Yang Y, Hammes-Schiffer S. Multicomponent Density Functional Theory: Impact of Nuclear Quantum Effects on Proton Affinities and Geometries. J Phys Chem Lett 2017; 8:3488-3493. [PMID: 28686449 DOI: 10.1021/acs.jpclett.7b01442] [Citation(s) in RCA: 78] [Impact Index Per Article: 11.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Nuclear quantum effects such as zero point energy play a critical role in computational chemistry and often are included as energetic corrections following geometry optimizations. The nuclear-electronic orbital (NEO) multicomponent density functional theory (DFT) method treats select nuclei, typically protons, quantum mechanically on the same level as the electrons. Electron-proton correlation is highly significant, and inadequate treatments lead to highly overlocalized nuclear densities. A recently developed electron-proton correlation functional, epc17, has been shown to provide accurate nuclear densities for molecular systems. Herein, the NEO-DFT/epc17 method is used to compute the proton affinities for a set of molecules and to examine the role of nuclear quantum effects on the equilibrium geometry of FHF-. The agreement of the computed results with experimental and benchmark values demonstrates the promise of this approach for including nuclear quantum effects in calculations of proton affinities, pKa's, optimized geometries, and reaction paths.
Collapse
Affiliation(s)
- Kurt R Brorsen
- Department of Chemistry, University of Illinois at Urbana-Champaign , 600 South Mathews Avenue, Urbana, Illinois 61801, United States
| | - Yang Yang
- Department of Chemistry, University of Illinois at Urbana-Champaign , 600 South Mathews Avenue, Urbana, Illinois 61801, United States
| | - Sharon Hammes-Schiffer
- Department of Chemistry, University of Illinois at Urbana-Champaign , 600 South Mathews Avenue, Urbana, Illinois 61801, United States
| |
Collapse
|
20
|
Gharabaghi M, Shahbazian S. Incorporating nuclear vibrational energies into the "atom in molecules" analysis: An analytical study. J Chem Phys 2017; 146:154106. [PMID: 28433028 DOI: 10.1063/1.4979994] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
Abstract
The quantum theory of atoms in molecules (QTAIM) is based on the clamped nucleus paradigm and solely working with the electronic wavefunctions, so does not include nuclear vibrations in the AIM analysis. On the other hand, the recently extended version of the QTAIM, called the multi-component QTAIM (MC-QTAIM), incorporates both electrons and quantum nuclei, i.e., those nuclei treated as quantum waves instead of clamped point charges, into the AIM analysis using non-adiabatic wavefunctions. Thus, the MC-QTAIM is the natural framework to incorporate the role of nuclear vibrations into the AIM analysis. In this study, within the context of the MC-QTAIM, the formalism of including nuclear vibrational energy in the atomic basin energy is developed in detail and its contribution is derived analytically using the recently proposed non-adiabatic Hartree product nuclear wavefunction. It is demonstrated that within the context of this wavefunction, the quantum nuclei may be conceived pseudo-adiabatically as quantum oscillators and both isotropic harmonic and anisotropic anharmonic oscillator models are used to compute the zero-point nuclear vibrational energy contribution to the basin energies explicitly. Inspired by the results gained within the context of the MC-QTAIM analysis, a heuristic approach is proposed within the context of the QTAIM to include nuclear vibrational energy in the basin energy from the vibrational wavefunction derived adiabatically. The explicit calculation of the basin contribution of the zero-point vibrational energy using the uncoupled harmonic oscillator model leads to results consistent with those derived from the MC-QTAIM.
Collapse
Affiliation(s)
- Masumeh Gharabaghi
- Faculty of Chemical and Petroleum Sciences, Shahid Beheshti University, G. C., Evin, P.O. Box 19395-4716, 19839 Tehran, Iran
| | - Shant Shahbazian
- Department of Physics, Shahid Beheshti University, G. C., Evin, P.O. Box 19395-4716, 19839 Tehran, Iran
| |
Collapse
|
21
|
Brorsen KR, Pak MV, Hammes-Schiffer S. Calculation of Positron Binding Energies and Electron–Positron Annihilation Rates for Atomic Systems with the Reduced Explicitly Correlated Hartree–Fock Method in the Nuclear–Electronic Orbital Framework. J Phys Chem A 2017; 121:515-522. [DOI: 10.1021/acs.jpca.6b10124] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Kurt R. Brorsen
- Department of Chemistry, 600 South Mathews Avenue, University of Illinois at Urbana−Champaign, Urbana, Illinois 61801, United States
| | - Michael V. Pak
- Department of Chemistry, 600 South Mathews Avenue, University of Illinois at Urbana−Champaign, Urbana, Illinois 61801, United States
| | - Sharon Hammes-Schiffer
- Department of Chemistry, 600 South Mathews Avenue, University of Illinois at Urbana−Champaign, Urbana, Illinois 61801, United States
| |
Collapse
|
22
|
Udagawa T, Sugiura K, Suzuki K, Tachikawa M. Unusual H/D isotope effect in isomerization and keto–enol tautomerism reactions of pyruvic acid: nuclear quantum effect restricts some rotational isomerization reactions. RSC Adv 2017. [DOI: 10.1039/c6ra28271g] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022] Open
Abstract
H/D isotope effects on isomerization and keto–enol tautomerism reactions of the pyruvic acid molecule have been investigated using the multicomponent B3LYP methods, which can take account of the nuclear quantum effect of protons and deuterons.
Collapse
Affiliation(s)
- Taro Udagawa
- Department of Chemistry and Biomolecular Science
- Faculty of Engineering
- Gifu University
- Gifu 501-1193
- Japan
| | - Keita Sugiura
- Department of Chemistry and Biomolecular Science
- Faculty of Engineering
- Gifu University
- Gifu 501-1193
- Japan
| | - Kimichi Suzuki
- Quantum Chemistry Division
- Graduate School of Science
- Yokohama City University
- Yokohama
- Japan
| | - Masanori Tachikawa
- Quantum Chemistry Division
- Graduate School of Science
- Yokohama City University
- Yokohama
- Japan
| |
Collapse
|
23
|
|
24
|
Affiliation(s)
- Harry F. King
- Department of Chemistry University at Buffalo State University of New York Buffalo, New York 14260, United States
| |
Collapse
|
25
|
Udagawa T, Tachikawa M. Nuclear quantum effect and H/D isotope effect on F + (H2O)n → FH + (H2O)n−1OH (n = 1-3) reactions. J Chem Phys 2016; 145:164310. [DOI: 10.1063/1.4966162] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
26
|
Culpitt T, Brorsen KR, Pak MV, Hammes-Schiffer S. Multicomponent density functional theory embedding formulation. J Chem Phys 2016; 145:044106. [DOI: 10.1063/1.4958952] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Affiliation(s)
- Tanner Culpitt
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Ave, Urbana, Illinois 61801, USA
| | - Kurt R. Brorsen
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Ave, Urbana, Illinois 61801, USA
| | - Michael V. Pak
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Ave, Urbana, Illinois 61801, USA
| | - Sharon Hammes-Schiffer
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Ave, Urbana, Illinois 61801, USA
| |
Collapse
|
27
|
Ellis BH, Aggarwal S, Chakraborty A. Development of the Multicomponent Coupled-Cluster Theory for Investigation of Multiexcitonic Interactions. J Chem Theory Comput 2015; 12:188-200. [DOI: 10.1021/acs.jctc.5b00879] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Benjamin H. Ellis
- Department
of Chemistry, Syracuse University, Syracuse, New York 13244, United States
| | - Somil Aggarwal
- Jamesville-DeWitt
High School, DeWitt, New York 13214, United States
| | - Arindam Chakraborty
- Department
of Chemistry, Syracuse University, Syracuse, New York 13244, United States
| |
Collapse
|
28
|
Yang Y, Kylänpää I, Tubman NM, Krogel JT, Hammes-Schiffer S, Ceperley DM. How large are nonadiabatic effects in atomic and diatomic systems? J Chem Phys 2015; 143:124308. [DOI: 10.1063/1.4931667] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yubo Yang
- Department of Physics, University of Illinois, Urbana, Illinois 61801, USA
| | - Ilkka Kylänpää
- Department of Physics, University of Illinois, Urbana, Illinois 61801, USA
- Department of Physics, Tampere University of Technology, P.O. Box 692, FI-33101 Tampere, Finland
| | - Norm M. Tubman
- Department of Physics, University of Illinois, Urbana, Illinois 61801, USA
| | - Jaron T. Krogel
- Oak Ridge National Laboratory, Materials Sciences & Technology Division, Oak Ridge, Tennessee 37831, USA
| | | | - David M. Ceperley
- Department of Physics, University of Illinois, Urbana, Illinois 61801, USA
| |
Collapse
|
29
|
Brorsen KR, Sirjoosingh A, Pak MV, Hammes-Schiffer S. Nuclear-electronic orbital reduced explicitly correlated Hartree-Fock approach: Restricted basis sets and open-shell systems. J Chem Phys 2015; 142:214108. [DOI: 10.1063/1.4921304] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Kurt R. Brorsen
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Ave., Urbana, Illinois 61801, USA
| | - Andrew Sirjoosingh
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Ave., Urbana, Illinois 61801, USA
| | - Michael V. Pak
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Ave., Urbana, Illinois 61801, USA
| | - Sharon Hammes-Schiffer
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Ave., Urbana, Illinois 61801, USA
| |
Collapse
|