1
|
Schneider M, Rauhut G. Comparison of curvilinear coordinates within vibrational structure calculations based on automatically generated potential energy surfaces. J Chem Phys 2024; 161:094102. [PMID: 39225517 DOI: 10.1063/5.0225991] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2024] [Accepted: 08/12/2024] [Indexed: 09/04/2024] Open
Abstract
For floppy molecules showing internal rotations and/or large amplitude motions, curvilinear internal coordinates are known to be superior to rectilinear normal coordinates within vibrational structure calculations. Due to the myriad definitions of internal coordinates, automated and efficient potential energy surface generators necessitate a high degree of flexibility, supporting the properties arising from these coordinates. Within this work, an approach to deal with these challenges is presented, including key elements, such as the selection of appropriate fit functions, the exploitation of symmetry, the positioning of grid points, or elongation limits for different coordinates. These elements are tested for five definitions of curvilinear coordinates, with three of them being generated in an automated manner. Calculations for semi-rigid molecules, namely H2O, H2CO, CH2F2, and H2CNH, demonstrate the general functionality of the implemented algorithms. Additional calculations for the HOPO molecule highlight the benefits of these curvilinear coordinates for systems with large amplitude motions. This new implementation allowed us to compare the performance of these different coordinate systems with respect to the convergence of the underlying expansion of the potential energy surface and subsequent vibrational configuration interaction calculations.
Collapse
Affiliation(s)
- Moritz Schneider
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Guntram Rauhut
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| |
Collapse
|
2
|
Rey M, Carrington T. Using nested tensor train contracted basis functions with group theoretical techniques to compute (ro)-vibrational spectra of molecules with non-Abelian groups. J Chem Phys 2024; 161:044102. [PMID: 39037133 DOI: 10.1063/5.0219434] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2024] [Accepted: 06/30/2024] [Indexed: 07/23/2024] Open
Abstract
In this paper, we use nested tensor-train contractions to compute vibrational and ro-vibrational energy levels of molecules with five and six atoms. At each step, we fully exploit symmetry by using symmetry adapted basis functions obtained from an irreducible tensor method. Contracted basis functions are determined by diagonalizing reduced dimensional Hamiltonian matrices. The size of matrices of eigenvectors, used to account for coupling between groups of coordinates, is reduced by discarding rows and columns. The size of the matrices that must be diagonalized is thus substantially reduced, making it possible to use direct eigensolvers, even for molecules with five and six atoms. The symmetry-adapted contracted vibrational basis functions have been used to compute J = 0 energy levels of the CH3CN (C3v) and J > 0 levels of CH4.
Collapse
Affiliation(s)
- Michaël Rey
- Groupe de Spectrométrie Moléculaire et Atmosphérique, UMR CNRS 7331, BP 1039, F-51687 Reims Cedex 2, France
| | - Tucker Carrington
- Chemistry Department, Queen's University, Kingston, Ontario K7L 3N6, Canada
| |
Collapse
|
3
|
Sharma M, Banik S, Roy TK. Performance of Effective Harmonic Oscillator Approach for the Calculations of Vibrational Transition Energies of Large Molecules. J Phys Chem A 2024; 128:5762-5776. [PMID: 38979981 DOI: 10.1021/acs.jpca.4c01583] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/10/2024]
Abstract
The accuracy and performance of the effective harmonic oscillator approximation for the description of anharmonic vibrational structure calculations are tested for large molecular systems and compared with experimental values along with vibrational self-consistent field and second-order perturbation theories. The effective harmonic oscillator approach is an effective single-particle approximation where the variational parameters are the centroids and widths of the multidimensional Gaussian product functions posited as the vibrational wave functions. A comprehensive calculation for 849 transitions that include the fundamentals, two and three quanta overtone transitions, and several combination bands of three polyaromatic hydrocarbons and one DNA nucleobase with a total of 231 normal modes are assessed. A comparison of EHO results with the experimental values is done for the polyaromatic hydrocarbons, and a close agreement is found between the two results. It also offers anharmonic eigenstates and eigenfunctions that are nearly identical with vibrational self-consistent field theory. An extensive analysis on the resultant wave functions of the excited states is performed. The overall root-mean-square deviation (RMSD) between these two methods for 849 transitions understudy is only about 8.3 cm-1, suggesting the effective harmonic oscillator as a viable alternative for the reliable calculations of transition energies of large molecular systems.
Collapse
Affiliation(s)
- Mokshi Sharma
- Department of Chemistry and Chemical Sciences, Central University of Jammu, Rahya-Suchani (Bagla), Jammu, Jammu and Kashmir 181143, India
| | - Subrata Banik
- Department of Chemistry, School of Chemical and Biotechnology, SASTRA Deemed University, Thanjavur, Tamil Nadu 613401, India
| | - Tapta Kanchan Roy
- Department of Chemistry and Chemical Sciences, Central University of Jammu, Rahya-Suchani (Bagla), Jammu, Jammu and Kashmir 181143, India
| |
Collapse
|
4
|
Zhang Q, Wang RS, Wang L. Neural canonical transformations for vibrational spectra of molecules. J Chem Phys 2024; 161:024103. [PMID: 38979703 DOI: 10.1063/5.0209255] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2024] [Accepted: 06/19/2024] [Indexed: 07/10/2024] Open
Abstract
The behavior of polyatomic molecules around their equilibrium positions can be regarded as that of quantum-coupled anharmonic oscillators. Solving the corresponding Schrödinger equations enables the interpretation or prediction of the experimental spectra of molecules. In this study, we developed a novel approach for solving the excited states of anharmonic vibrational systems. The normal coordinates of the molecules are transformed into new coordinates through a normalizing flow parameterized by a neural network. This facilitates the construction of a set of orthogonal many-body variational wavefunctions. This methodology has been validated on an exactly solvable 64-dimensional coupled harmonic oscillator, yielding numerical results with a relative error of 10-6. The neural canonical transformations are also applied to calculate the energy levels of two specific molecules, acetonitrile (CH3CN) and ethylene oxide (C2H4O). These molecules involve 12 and 15 vibrational modes, respectively. A key advantage of this approach is its flexibility concerning the potential energy surface, as it requires no specific form. Furthermore, this method can be readily implemented on large-scale distributed computing platforms, making it easy to extend to investigating complex vibrational structures.
Collapse
Affiliation(s)
- Qi Zhang
- Beijing National Laboratory for Condensed Matter Physics and Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
| | - Rui-Si Wang
- Beijing National Laboratory for Condensed Matter Physics and Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Lei Wang
- Beijing National Laboratory for Condensed Matter Physics and Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
- Songshan Lake Materials Laboratory, Dongguan, Guangdong 523808, China
| |
Collapse
|
5
|
Suda K, Yokogawa D. Vibrational Self-Consistent Field (VSCF) and Post-VSCF Method Calculations Combined with the Reference Interaction Site Model Self-Consistent Field Method Coupled with the Constrained Spatial Electron Density Distribution: Applications to NaHCOO in Aqueous Phase. J Chem Theory Comput 2024; 20:4885-4892. [PMID: 38815984 DOI: 10.1021/acs.jctc.4c00360] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/01/2024]
Abstract
Investigating vibrational behavior in solution is crucial for understanding molecular dynamics within a solvent environment. Notably, the analysis of Raman spectra for molecules in solution is important owing to its ability to unveil intricate solute-solvent interactions. Previous studies have effectively employed frequency calculations utilizing the reference interaction site model self-consistent field method in conjunction with constrained spatial electron density distribution (RISM-SCF-cSED) to understand molecular vibrations in solution, primarily focusing on fundamental vibrational modes. However, the oversight of overtones and combination tones in these studies prompted us to combine the vibrational self-consistent field (VSCF) and vibrational second-order Mo̷ller-Plesset perturbation (VMP2) methods with RISM-SCF-cSED to address these aspects theoretically. Illustrating the efficacy of this integrated approach, we computed the Raman spectra of sodium formate (NaHCOO) in water, revealing the necessity of accounting for molecular anharmonicity in solution vibrational analysis. Our findings underscore the potency of VSCF and VMP2 in conjunction with RISM-SCF-cSED as a robust theoretical framework for such calculations.
Collapse
Affiliation(s)
- Kayo Suda
- Graduate School of Arts and Sciences, The University of Tokyo, 3-8-1 Komaba, Meguro-ku, Tokyo 153-8902, Japan
| | - Daisuke Yokogawa
- Graduate School of Arts and Sciences, The University of Tokyo, 3-8-1 Komaba, Meguro-ku, Tokyo 153-8902, Japan
| |
Collapse
|
6
|
Schneider M, Rauhut G. VSCF/VCI theory based on the Podolsky Hamiltonian. J Chem Phys 2024; 160:214118. [PMID: 38832735 DOI: 10.1063/5.0213401] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2024] [Accepted: 05/14/2024] [Indexed: 06/05/2024] Open
Abstract
While the vibrational spectra of semi-rigid molecules can be computed on approaches relying on the Watson Hamiltonian, floppy molecules or molecular clusters are better described by Hamiltonians, which are capable of dealing with any curvilinear coordinates. It is the kinetic energy operator (KEO) of these Hamiltonians, which render the correlated calculations relying on them rather costly. Novel implementation of vibrational self-consistent field theory and vibrational configuration interaction theory on the basis of the Podolsky Hamiltonian are reported, in which the inverse of the metric tensor, i.e., the G matrix, is represented by an n-mode expansion expressed in terms of polynomials. An analysis of the importance of the individual terms of the KEO with respect to the truncation orders of the n-mode expansion is provided. Benchmark calculations have been performed for the cis-HOPO and methanimine, H2CNH, molecules and are compared to experimental data and to calculations based on the Watson Hamiltonian and the internal coordinate path Hamiltonian.
Collapse
Affiliation(s)
- Moritz Schneider
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Guntram Rauhut
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| |
Collapse
|
7
|
Wang XG, Carrington T. A two-step quadrature-based variational calculation of ro-vibrational levels and wavefunctions of CO 2 using a bisector- x molecule-fixed frame. Phys Chem Chem Phys 2024; 26:15181-15191. [PMID: 38752328 DOI: 10.1039/d4cp00655k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/30/2024]
Abstract
In this paper, we propose a new two-step strategy for computing ro-vibrational energy levels and wavefunctions of a triatomic molecule and apply it to CO2. A two-step method [J. Tennyson and B. T. Sutcliffe, Mol. Phys., 1986, 58, 1067] uses a basis whose functions are products of K-dependent "vibrational" functions and symmetric top functions. K is the quantum number for the molecule-fixed z component of the angular momentum. For a linear molecule, a two-step method is efficient because the Hamiltonian used to compute the basis functions includes the largest coupling term. The most important distinguishing feature of the two-step method we propose is that it uses an associated Legendre basis and quadrature rather than a K-dependent discrete variable representation. This reduces the cost of the calculation and simplifies the method. We have computed ro-vibrational energy levels with J up to 100 for CO2, on an accurate available potential energy surface which is known as the AMES-2 PES and present a subset of those levels. We have converged most levels up to 20 000 cm-1 to 0.0001 cm-1.
Collapse
Affiliation(s)
- Xiao-Gang Wang
- Chemistry Department, Queen's University, Kingston, Ontario K7L 3N6, Canada.
| | - Tucker Carrington
- Chemistry Department, Queen's University, Kingston, Ontario K7L 3N6, Canada.
| |
Collapse
|
8
|
Bader F, Lauvergnat D, Christiansen O. Efficient vibrationally correlated calculations using n-mode expansion-based kinetic energy operators. Phys Chem Chem Phys 2024; 26:11469-11481. [PMID: 38546727 DOI: 10.1039/d4cp00423j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/18/2024]
Abstract
Due to its efficiency and flexibility, the n-mode expansion is a frequently used tool for representing molecular potential energy surfaces in quantum chemical simulations. In this work, we investigate the performance of n-mode expansion-based models of kinetic energy operators in general polyspherical coordinate systems. In particular, we assess the operators with respect to accuracy in vibrationally correlated calculations and their effect on potential energy surface construction with the adaptive density guided approach. Our results show that the n-mode expansion-based operator variants are reliable and systematically improvable approximations of the full kinetic energy operator. Moreover, we introduce a workflow to generate the n-mode expanded kinetic energy operators on-the-fly within the adaptive density guided approach. This scheme can be applied in studies of species and coordinate systems, for which an analytical form of the kinetic energy operator is not available.
Collapse
Affiliation(s)
- Frederik Bader
- Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark.
| | - David Lauvergnat
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, Orsay 91405, France.
| | - Ove Christiansen
- Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark.
| |
Collapse
|
9
|
Schröder B, Rauhut G. From the Automated Calculation of Potential Energy Surfaces to Accurate Infrared Spectra. J Phys Chem Lett 2024; 15:3159-3169. [PMID: 38478898 PMCID: PMC10961845 DOI: 10.1021/acs.jpclett.4c00186] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2024] [Revised: 02/20/2024] [Accepted: 02/28/2024] [Indexed: 03/22/2024]
Abstract
Advances in the development of quantum chemical methods and progress in multicore architectures in computer science made the simulation of infrared spectra of isolated molecules competitive with respect to established experimental methods. Although it is mainly the multidimensional potential energy surface that controls the accuracy of these calculations, the subsequent vibrational structure calculations need to be carefully converged in order to yield accurate results. As both aspects need to be considered in a balanced way, we focus on approaches for molecules of up to 12-15 atoms with respect to both parts, which have been automated to some extent so that they can be employed in routine applications. Alternatives to machine learning will be discussed, which appear to be attractive, as long as local regions of the potential energy surface are sufficient. The automatization of these methods is still in its infancy, and the generalization to molecules with large amplitude motions or molecular clusters is far from trivial, but many systems relevant for astrophysical studies are already in reach.
Collapse
Affiliation(s)
- Benjamin Schröder
- Institute
of Physical Chemistry, University of Goettingen, Tammannstrasse 6, Göttingen 37077, Germany
| | - Guntram Rauhut
- Institute
for Theoretical Chemistry, University of
Stuttgart, Pfaffenwaldring 55, Stuttgart 70569, Germany
| |
Collapse
|
10
|
Attal L, Falvo C, Calvo F, Parneix P. Modeling the dynamics of quantum systems coupled to large-dimensional baths using effective energy states. J Chem Phys 2024; 160:044107. [PMID: 38270236 DOI: 10.1063/5.0184299] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2023] [Accepted: 12/14/2023] [Indexed: 01/26/2024] Open
Abstract
The quantum dynamics of a low-dimensional system in contact with a large but finite harmonic bath is theoretically investigated by coarse-graining the bath into a reduced set of effective energy states. In this model, the couplings between the system and the bath are obtained from statistically averaging over the discrete, degenerate effective states. Our model is aimed at intermediate bath sizes in which non-Markovian processes and energy transfer between the bath and the main system are important. The method is applied to a model system of a Morse oscillator coupled to 40 harmonic modes. The results are found to be in excellent agreement with the direct quantum dynamics simulations presented in the work of Bouakline et al. [J. Phys. Chem. A 116, 11118-11127 (2012)], but at a much lower computational cost. Extension to larger baths is discussed in comparison to the time-convolutionless method. We also extend this study to the case of a microcanonical bath with finite initial internal energies. The computational efficiency and convergence properties of the effective bath states model with respect to relevant parameters are also discussed.
Collapse
Affiliation(s)
- Loïse Attal
- Institut des Sciences Moléculaires d'Orsay, CNRS, Université Paris-Saclay, 91405 Orsay, France
| | - Cyril Falvo
- Institut des Sciences Moléculaires d'Orsay, CNRS, Université Paris-Saclay, 91405 Orsay, France
- Université Grenoble-Alpes, CNRS, LiPhy, 38000 Grenoble, France
| | - Florent Calvo
- Université Grenoble-Alpes, CNRS, LiPhy, 38000 Grenoble, France
| | - Pascal Parneix
- Institut des Sciences Moléculaires d'Orsay, CNRS, Université Paris-Saclay, 91405 Orsay, France
| |
Collapse
|
11
|
Houston PL, Qu C, Yu Q, Pandey P, Conte R, Nandi A, Bowman JM. A New A Priori Method to Avoid Calculation of Negligible Hamiltonian Matrix Elements in CI Calculation. J Phys Chem A 2024; 128:479-487. [PMID: 38180902 DOI: 10.1021/acs.jpca.3c06985] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2024]
Abstract
Hamiltonian matrices typically contain many elements that are negligibly small compared to the diagonal elements, even with methods to prune the underlying basis. Because for general potentials the calculation of H-matrix elements is a major part of the computational effort to obtain eigenvalues and eigenfunctions of the Hamiltonian, there is strong motivation to investigate locating these negligible elements without calculating them or at least avoid calculating them. We recently demonstrated an effective means to "learn" negligible elements using machine learning classification (J. Chem. Phys. 2023, 159, 071101). Here we present a simple, new method to avoid calculating them by using a cut-off value for the absolute difference in the quantum numbers for the bra and ket. This method is demonstrated for many of the same case studies as were used in the paper above, namely for realistic H-matrices of H2O, the vinyl radical, C2H3, and glycine, C2H5NO2. The new method is compared to the recently reported machine learning approach. In addition, we point out an important synergy between the two methods.
Collapse
Affiliation(s)
- Paul L Houston
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
- Department of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | - Chen Qu
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Qi Yu
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Priyanka Pandey
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Riccardo Conte
- Dipartimento di Chimica, Università degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| | - Apurba Nandi
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - Joel M Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
12
|
Manzhos S, Ihara M. A controlled study of the effect of deviations from symmetry of the potential energy surface (PES) on the accuracy of the vibrational spectrum computed with collocation. J Chem Phys 2023; 159:211103. [PMID: 38038200 DOI: 10.1063/5.0182373] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2023] [Accepted: 11/12/2023] [Indexed: 12/02/2023] Open
Abstract
Symmetry, in particular permutational symmetry, of a potential energy surface (PES) is a useful property in quantum chemical calculations. It facilitates, in particular, state labelling and identification of degenerate states. In many practically important applications, however, these issues are unimportant. The imposition of exact symmetry and the perception that it is necessary create additional methodological requirements narrowing or complicating algorithmic choices that are thereby biased against methods and codes that by default do not incorporate symmetry, including most off-the-shelf machine learning methods that cannot be directly used if exact symmetry is demanded. By introducing symmetric and unsymmetric errors into the PES of H2CO in a controlled way and computing the vibrational spectrum with collocation using symmetric and nonsymmetric collocation point sets, we show that when the deviations from an ideal PES are random, imposition of exact symmetry does not bring any practical advantages. Moreover, a calculation ignoring symmetry may be more accurate. We also compare machine-learned PESs with and without symmetrization and demonstrate that there is no advantage of imposing exact symmetry for the accuracy of the vibrational spectrum.
Collapse
Affiliation(s)
- Sergei Manzhos
- School of Materials and Chemical Technology, Tokyo Institute of Technology, Ookayama 2-12-1, Meguro-ku, Tokyo 152-8552, Japan
| | - Manabu Ihara
- School of Materials and Chemical Technology, Tokyo Institute of Technology, Ookayama 2-12-1, Meguro-ku, Tokyo 152-8552, Japan
| |
Collapse
|
13
|
De Vos J, Schröder B, Rauhut G. Comprehensive quantum chemical analysis of the (ro)vibrational spectrum of thiirane and its deuterated isotopologue. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2023; 302:123083. [PMID: 37423098 DOI: 10.1016/j.saa.2023.123083] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/31/2023] [Revised: 06/23/2023] [Accepted: 06/26/2023] [Indexed: 07/11/2023]
Abstract
The (ro)vibrational spectra of thiirane, c-C2H4S, and its fully deuterated isotopologue, c-C2D4S, have been studied by means of vibrational configuration interaction theory, VCI, its incremental variant, iVCI, and subsequent variational rovibrational calculations, RVCI, which rely on multidimensional potential energy surfaces of coupled-cluster quality including up to four-mode coupling terms. Accurate geometrical parameters, fundamental vibrational transitions and first overtones, rovibrational spectra and rotational spectroscopic constants have been determined from these calculations and were compared with experimental results whenever available. A number of tentative misassignments in the vibrational spectra could be resolved and most results for the deuterated thiirane are high-level predictions, which may guide experiments to come. Besides this, a new implementation of infrared intensities within the iVCI framework has been tested for the transitions of the title compounds and are compared with results obtained from standard VCI calculations.
Collapse
Affiliation(s)
- John De Vos
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569, Stuttgart, Germany; Department of Chemistry, Ghent University, Krijgslaan 281, 9000, Ghent, Belgium.
| | - Benjamin Schröder
- Institute for Physical Chemistry, University of Göttingen, Tammannstr. 6, 37077, Göttingen, Germany.
| | - Guntram Rauhut
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569, Stuttgart, Germany.
| |
Collapse
|
14
|
Tran HK, Berkelbach TC. Vibrational heat-bath configuration interaction with semistochastic perturbation theory using harmonic oscillator or VSCF modals. J Chem Phys 2023; 159:194101. [PMID: 37965997 PMCID: PMC10653875 DOI: 10.1063/5.0172702] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2023] [Accepted: 10/23/2023] [Indexed: 11/16/2023] Open
Abstract
Vibrational heat-bath configuration interaction (VHCI)-a selected configuration interaction technique for vibrational structure theory-has recently been developed in two independent works [J. H. Fetherolf and T. C. Berkelbach, J. Chem. Phys. 154, 074104 (2021); A. U. Bhatty and K. R. Brorsen, Mol. Phys. 119, e1936250 (2021)], where it was shown to provide accuracy on par with the most accurate vibrational structure methods with a low computational cost. Here, we eliminate the memory bottleneck of the second-order perturbation theory correction using the same (semi)stochastic approach developed previously for electronic structure theory. This allows us to treat, in an unbiased manner, much larger perturbative spaces, which are necessary for high accuracy in large systems. Stochastic errors are easily controlled to be less than 1 cm-1. We also report two other developments: (i) we propose a new heat-bath criterion and an associated exact implicit sorting algorithm for potential energy surfaces expressible as a sum of products of one-dimensional potentials; (ii) we formulate VHCI to use a vibrational self-consistent field (VSCF) reference, as opposed to the harmonic oscillator reference configuration used in previous reports. Our tests are done with quartic and sextic force fields, for which we find that with VSCF, the minor improvements to accuracy are outweighed by the higher computational cost associated the matrix element evaluations. We expect VSCF-based VHCI to be important for more general potential representations, for which the harmonic oscillator basis function integrals are no longer analytic.
Collapse
Affiliation(s)
- Henry K. Tran
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | | |
Collapse
|
15
|
Vogt E, Simkó I, Császár AG, Kjaergaard HG. Quantum Chemical Investigation of the Cold Water Dimer Spectrum in the First OH-Stretching Overtone Region Provides a New Interpretation. J Phys Chem A 2023; 127:9409-9418. [PMID: 37930939 DOI: 10.1021/acs.jpca.3c03705] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2023]
Abstract
Intramolecular vibrational transition wavenumbers and intensities were calculated in the fundamental HOH-bending, fundamental OH-stretching, first OH-stretching-HOH-bending combination, and first OH-stretching overtone (ΔvOH = 2) regions of the water dimer's spectrum. Furthermore, the rotational-vibrational spectrum was calculated in the ΔvOH = 2 region at 10 K, corresponding to the temperature of the existing jet-expansion experiments. The calculated spectrum was obtained by combining results from a full-dimensional (12D) vibrational and a reduced-dimensional vibrational-rotational-tunneling model. The ΔvOH = 2 spectral region is rich in features due to contributions from multiple vibrational-rotational-tunneling sub-bands. Origins of the experimental vibrational bands depend on the assignment of the observed sub-bands. Based on our calculations, we assign the observed sub-bands, and our reassignment leads to new values for the vibrational band origins of the free donor and antisymmetric acceptor OH-stretching first overtones of ∼7227 and ∼7238 cm-1, respectively. The observed bands with origins at 7192.34 and ∼7366 cm-1 are assigned to the symmetric acceptor OH-stretching first overtone and the OH-stretching combination of the donor, respectively.
Collapse
Affiliation(s)
- Emil Vogt
- Department of Chemistry, University of Copenhagen, Universitetsparken 5, Copenhagen Ø DK-2100, Denmark
| | - Irén Simkó
- Laboratory of Molecular Structure and Dynamics, Institute of Chemistry ELTE Eötvös Loránd University, Pázmány Péter Sétány 1/A, Budapest H-1117, Hungary
- HUN-REN-ELTE Complex Chemical Systems Research Group, P.O. Box 32, Budapest 112 H-1518, Hungary
| | - Attila G Császár
- Laboratory of Molecular Structure and Dynamics, Institute of Chemistry ELTE Eötvös Loránd University, Pázmány Péter Sétány 1/A, Budapest H-1117, Hungary
- HUN-REN-ELTE Complex Chemical Systems Research Group, P.O. Box 32, Budapest 112 H-1518, Hungary
| | - Henrik G Kjaergaard
- Department of Chemistry, University of Copenhagen, Universitetsparken 5, Copenhagen Ø DK-2100, Denmark
| |
Collapse
|
16
|
Manzhos S, Ihara M. Neural Network with Optimal Neuron Activation Functions Based on Additive Gaussian Process Regression. J Phys Chem A 2023; 127:7823-7835. [PMID: 37698519 DOI: 10.1021/acs.jpca.3c02949] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/13/2023]
Abstract
Feed-forward neural networks (NNs) are a staple machine learning method widely used in many areas of science and technology, including physical chemistry, computational chemistry, and materials informatics. While even a single-hidden-layer NN is a universal approximator, its expressive power is limited by the use of simple neuron activation functions (such as sigmoid functions) that are typically the same for all neurons. More flexible neuron activation functions would allow the use of fewer neurons and layers and thereby save computational cost and improve expressive power. We show that additive Gaussian process regression (GPR) can be used to construct optimal neuron activation functions that are individual to each neuron. An approach is also introduced that avoids nonlinear fitting of neural network parameters by defining them with rules. The resulting method combines the advantage of robustness of a linear regression with the higher expressive power of an NN. We demonstrate the approach by fitting the potential energy surfaces of the water molecule and formaldehyde. Without requiring any nonlinear optimization, the additive-GPR-based approach outperforms a conventional NN in the high-accuracy regime, where a conventional NN suffers more from overfitting.
Collapse
Affiliation(s)
- Sergei Manzhos
- School of Materials and Chemical Technology, Tokyo Institute of Technology, Ookayama 2-12-1, Meguro-ku, Tokyo 152-8552, Japan
| | - Manabu Ihara
- School of Materials and Chemical Technology, Tokyo Institute of Technology, Ookayama 2-12-1, Meguro-ku, Tokyo 152-8552, Japan
| |
Collapse
|
17
|
Sharma D, Roy TK. Accuracy of Different Electronic Basis Set Families for Anharmonic Molecular Vibrations: A Comprehensive Benchmark Study. J Phys Chem A 2023; 127:7132-7147. [PMID: 37603414 DOI: 10.1021/acs.jpca.3c02874] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/23/2023]
Abstract
In this work, the accuracy and convergence of different electronic basis set families for the computation of anharmonic molecular vibrational spectroscopic calculations are benchmarked. A series of 39 different basis sets from different families following their hierarchy are assessed on VSCF and VSCF-PT2 algorithms with commonly used MP2 and DFT based B3LYP-D potentials for a set of molecular systems. Such an effort has been validated in a previous work ( J. Phys. Chem. A 2020, 124, 9203-9221) with split-valence basis sets for fundamentals and intensities. Here, fundamental transitions, vibrationally excited states, and intensities are compared with the experimental data to estimate the accuracy for a series of Jensen, Dunning, Calendar, Karlsruhe, and Sapporo basis set families. The convergence of basis sets are also compared with the large ANO basis set. Comprehensive statistical error analysis in terms of accuracy and precision was carried out to assess the performance of each basis set. It is observed that the improvement for the calculated harmonic and anharmonic values from the smaller basis sets to the medium (i.e., triple-ξ) is considerable. Beyond this, from medium to large basis sets, the convergence is slow and mostly posits nearly converged values. Basis sets with and without diffuse functions offer characteristically different accuracies and convergence patterns. Finally, recommendations are given on the choice of basis set chosen as black-box which can balance between accuracy and computational time, estimation of the errors, and their selections especially for large molecules.
Collapse
Affiliation(s)
- Dhiksha Sharma
- Department of Chemistry and Chemical Sciences, Central University of Jammu, Rahya-Suchani (Bagla), Jammu, J&K 181143 India
| | - Tapta Kanchan Roy
- Department of Chemistry and Chemical Sciences, Central University of Jammu, Rahya-Suchani (Bagla), Jammu, J&K 181143 India
| |
Collapse
|
18
|
Schröder B. Ab Initio Rovibrational Spectroscopy of the Acetylide Anion. Molecules 2023; 28:5700. [PMID: 37570670 PMCID: PMC10420331 DOI: 10.3390/molecules28155700] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2023] [Revised: 07/22/2023] [Accepted: 07/22/2023] [Indexed: 08/13/2023] Open
Abstract
In this work the rovibrational spectrum of the acetylide anion HCC- is investigated using high-level electronic structure methods and variational rovibrational calculations. Using a composite approach the potential energy surface and dipole surface is constructed from explicitly correlated coupled-cluster accounting for corrections due to core-valence correlation, scalar relativistic effects and higher-order excitation effects. Previous approaches for approximating the latter are critically evaluated. Employing the composite potential, accurate spectroscopic parameters determined from variational calculations are presented. In comparison to the few available reference data the present results show excellent agreement with ground state rotational constants within 0.005% of the experimental value. Intensities determined from the variational calculations suggest the bending fundamental transition ν2 around 510 cm-1 to be the best target for detection. The rather weak CD stretching fundamental ν1 in deuterated isotopologues show a second-order resonance with the (0,20,1) state and the consequences are discussed in some detail. The spectroscopic parameters and band intensities provided for a number of vibrational bands in isotopologues of the acetylide anion should facilitate future spectroscopic investigations.
Collapse
Affiliation(s)
- Benjamin Schröder
- Institute of Physical Chemistry, University of Goettingen, Tammannstr. 6, 37077 Göttingen, Germany
| |
Collapse
|
19
|
Guan J, Lu Y, Sen K, Abdul Nasir J, Desmoutier AW, Hou Q, Zhang X, Logsdail AJ, Dutta G, Beale AM, Strange RW, Yong C, Sherwood P, Senn HM, Catlow CRA, Keal TW, Sokol AA. Computational infrared and Raman spectra by hybrid QM/MM techniques: a study on molecular and catalytic material systems. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2023; 381:20220234. [PMID: 37211033 PMCID: PMC10200352 DOI: 10.1098/rsta.2022.0234] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/23/2022] [Accepted: 04/04/2023] [Indexed: 05/23/2023]
Abstract
Vibrational spectroscopy is one of the most well-established and important techniques for characterizing chemical systems. To aid the interpretation of experimental infrared and Raman spectra, we report on recent theoretical developments in the ChemShell computational chemistry environment for modelling vibrational signatures. The hybrid quantum mechanical and molecular mechanical approach is employed, using density functional theory for the electronic structure calculations and classical forcefields for the environment. Computational vibrational intensities at chemical active sites are reported using electrostatic and fully polarizable embedding environments to achieve more realistic vibrational signatures for materials and molecular systems, including solvated molecules, proteins, zeolites and metal oxide surfaces, providing useful insight into the effect of the chemical environment on the signatures obtained from experiment. This work has been enabled by the efficient task-farming parallelism implemented in ChemShell for high-performance computing platforms. This article is part of a discussion meeting issue 'Supercomputing simulations of advanced materials'.
Collapse
Affiliation(s)
- Jingcheng Guan
- Department of Chemistry, University College London, London WC1H 0AJ, UK
| | - You Lu
- STFC Scientific Computing, Daresbury Laboratory, Keckwick Lane, Daresbury, Warrington WA4 4AD, UK
| | - Kakali Sen
- STFC Scientific Computing, Daresbury Laboratory, Keckwick Lane, Daresbury, Warrington WA4 4AD, UK
| | - Jamal Abdul Nasir
- Department of Chemistry, University College London, London WC1H 0AJ, UK
| | | | - Qing Hou
- Department of Chemistry, University College London, London WC1H 0AJ, UK
- Institute of Photonic Chips, University of Shanghai for Science of Technology, Shanghai 201512, People’s Republic of China
| | - Xingfan Zhang
- Department of Chemistry, University College London, London WC1H 0AJ, UK
| | - Andrew J. Logsdail
- Cardiff Catalysis Institute, School of Chemistry, Cardiff University, Park Place, Cardiff CF10 3AT, UK
| | - Gargi Dutta
- Department of Chemistry, University College London, London WC1H 0AJ, UK
- Department of Physics, Balurghat College, Balurghat 733101, West Bengal, India
| | - Andrew M. Beale
- Department of Chemistry, University College London, London WC1H 0AJ, UK
- Research Complex at Harwell, Rutherford Appleton Laboratory, Harwell Oxford, Didcot OX11 0FA, UK
| | - Richard W. Strange
- School of Life Sciences, University of Essex, Wivenhoe Park, Colchester, Essex CO4 3SQ, UK
| | - Chin Yong
- STFC Scientific Computing, Daresbury Laboratory, Keckwick Lane, Daresbury, Warrington WA4 4AD, UK
| | - Paul Sherwood
- Department of Chemistry, Lancaster University, Lancaster LA1 4YB, UK
| | - Hans M. Senn
- School of Chemistry, University of Glasgow, Joseph Black Building, Glasgow G12 8QQ, UK
| | - C. Richard A. Catlow
- Department of Chemistry, University College London, London WC1H 0AJ, UK
- Cardiff Catalysis Institute, School of Chemistry, Cardiff University, Park Place, Cardiff CF10 3AT, UK
- Research Complex at Harwell, Rutherford Appleton Laboratory, Harwell Oxford, Didcot OX11 0FA, UK
| | - Thomas W. Keal
- STFC Scientific Computing, Daresbury Laboratory, Keckwick Lane, Daresbury, Warrington WA4 4AD, UK
| | - Alexey A. Sokol
- Department of Chemistry, University College London, London WC1H 0AJ, UK
| |
Collapse
|
20
|
Wang Y, Chen Z, Yang Y. Calculating Vibrational Excited State Absorptions with Excited State Constrained Minimized Energy Surfaces. J Phys Chem A 2023. [PMID: 37335973 DOI: 10.1021/acs.jpca.3c01420] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/21/2023]
Abstract
The modeling and interpretation of vibrational spectra are crucial for studying reaction dynamics using vibrational spectroscopy. Most prior theoretical developments focused on describing fundamental vibrational transitions while fewer developments focused on vibrational excited state absorptions. In this study, we present a new method that uses excited state constrained minimized energy surfaces (CMESs) to describe vibrational excited state absorptions. The excited state CMESs are obtained similarly to the previous ground state CMES development in our group but with additional wave function orthogonality constraints. Using a series of model systems, including the harmonic oscillator, Morse potential, double-well potential, quartic potential, and two-dimensional anharmonic potential, we demonstrate that this new procedure provides good estimations of the transition frequencies for vibrational excited state absorptions. These results are significantly better than those obtained from harmonic approximations using conventional potential energy surfaces, demonstrating the promise of excited state CMES-based methods for calculating vibrational excited state absorptions in real systems.
Collapse
Affiliation(s)
- Yiwen Wang
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin─Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Zehua Chen
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin─Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Yang Yang
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin─Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| |
Collapse
|
21
|
Wang XG, Carrington T. Computing excited OH stretch states of water dimer in 12D using contracted intermolecular and intramolecular basis functions. J Chem Phys 2023; 158:084107. [PMID: 36859104 DOI: 10.1063/5.0139586] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023] Open
Abstract
Due to the ubiquity and importance of water, water dimer has been intensively studied. Computing the (ro-)vibrational spectrum of water dimer is challenging. The potential has eight wells separated by low barriers, which makes harmonic approximations of limited utility. A variational approach is imperative, but difficult because there are 12 coupled vibrational coordinates. In this paper, we use a product contracted basis whose functions are products of intramolecular and intermolecular functions computed using an iterative eigensolver. An intermediate matrix F facilitates calculating matrix elements. Using F, it is possible to do calculations on a general potential without storing the potential on the full quadrature grid. We find that surprisingly many intermolecular functions are required. This is due to the importance of coupling between inter- and intra-molecular coordinates. The full G16 symmetry of water dimer is exploited. We calculate, for the first time, monomer excited stretch states and compare P(1) transition frequencies with their experimental counterparts. We also compare with experimental vibrational shifts and tunneling splittings. Surprisingly, we find that the largest tunneling splitting, which does not involve the interchange of the two monomers, is smaller in the asymmetric stretch excited state than in the ground state. Differences between levels we compute and those obtained with a [6+6]D adiabatic approximation [Leforestier et al. J. Chem. Phys. 137 014305 (2012)] are ∼0.6 cm-1 for states without monomer excitation, ∼4 cm-1 for monomer excited bend states, and as large as ∼10 cm-1 for monomer excited stretch states.
Collapse
Affiliation(s)
- Xiao-Gang Wang
- Chemistry Department, Queen's University, Kingston, Ontario K7L 3N6, Canada
| | - Tucker Carrington
- Chemistry Department, Queen's University, Kingston, Ontario K7L 3N6, Canada
| |
Collapse
|
22
|
Zapata Trujillo JC, McKemmish LK. Model Chemistry Recommendations for Scaled Harmonic Frequency Calculations: A Benchmark Study. J Phys Chem A 2023; 127:1715-1735. [PMID: 36753303 DOI: 10.1021/acs.jpca.2c06908] [Citation(s) in RCA: 12] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/09/2023]
Abstract
Despite the widespread popularity of scaled harmonic frequency calculations to predict experimental fundamental frequencies in chemistry, sparse benchmarking is available to guide users on the appropriate level of theory and basis set choices (model chemistry) or deep understanding of expected errors. An updated assessment of the best approach for scaling to minimize errors is also overdue. Here, we assess the performance of over 600 popular, contemporary, and robust model chemistries in the calculation of scaled harmonic frequencies, evaluating different scaling factor types and their implications in the scaled harmonic frequencies and model chemistry performance. We can summarize our results into three main findings: (1) Using model-chemistry-specific scaling factors optimized for three different frequency regions (low (<1,000 cm-1), mid (1,000-2,000 cm-1), and high (>2,000 cm-1)) results in substantial improvements in the agreement between the scaled harmonic and experimental frequencies compared to other choices. (2) Larger basis sets and more robust levels of theory generally lead to superior performance; however, the particular model chemistry choice matters and poor choices lead to significantly reduced accuracies. (3) Outliers are expected in routine calculations regardless of the model chemistry choice. Our benchmarking results here do not consider the intensity of vibrational transitions; however, we draw upon previous benchmarking results for dipole moments that highlight the importance of diffuse functions (i.e., augmented basis sets) in high-quality intensity predictions. In terms of specific recommendations, overall, the highest accuracy model chemistries are double-hybrid density functional approximations with a non-Pople augmented triple-ζ basis set, which can produce median frequency errors down to 7.6 cm-1 (DSD-PBEP86/def2-TZVPD), which is very close to the error in the harmonic approximation, i.e., the anharmonicity error. Double-ζ basis sets should not be used with double-hybrid functionals as there is no improvement compared to hybrid functional results (unlike for double-hybrid triple-ζ model chemistries). Note that 6-311G* and 6-311+G* basis sets perform like a double-ζ basis set for vibrational frequencies. After scaling, all studied hybrid functionals with non-Pople triple-ζ basis sets will produce median errors of less than 15 cm-1, with the best result of 9.9 cm-1 with B97-1/def2-TZVPD. Appropriate matching of double-ζ basis sets with hybrid functionals can produce high-quality results, but the precise choice of functional and basis set is more important. The B97-1, TPSS0-D3(BJ), or ωB97X-D hybrid density functionals with 6-31G*, pc-1, or pcseg-1 are recommended for fast routine calculations, all delivering median errors of 11-12 cm-1. Note that dispersion corrections are not easily available for B97-1; given its strong performance here, we recommend these be added to major programs in coming updates.
Collapse
Affiliation(s)
| | - Laura K McKemmish
- School of Chemistry, University of New South Wales, 2052 Sydney, NSW, Australia
| |
Collapse
|
23
|
Fayaz A, Banik S, Kanchan Roy T. The importance of electron correlations on vibrational anharmonicities and potential energy surfaces. COMPUT THEOR CHEM 2023. [DOI: 10.1016/j.comptc.2023.114059] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
|
24
|
Schneider M, Rauhut G. Quantum chemical rovibrational analysis of aminoborane and its isotopologues. J Comput Chem 2023; 44:298-306. [PMID: 35582830 DOI: 10.1002/jcc.26893] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2022] [Revised: 04/22/2022] [Accepted: 04/28/2022] [Indexed: 12/31/2022]
Abstract
Aminoborane, H2 NBH2 and its isotopologues, H2 N10 BH2 , D2 NBD2 , and D2 N10 BD2 , have been studied by high-level ab initio methods. All calculations rely on multidimensional potential energy surfaces and dipole moment surfaces including high-order mode coupling terms, which have been obtained from electronic structure calculations at the level of explicitly correlated coupled-cluster theory, CCSD(T)-F12, or the distinguishable cluster approximation, DCSD. Subsequent vibrational structure calculations based on second-order vibrational perturbation theory, VPT2, and vibrational configuration interaction theory, VCI, were used to determine rotational constants, centrifugal distortion constants, vibrationally averaged geometrical parameters and (an)harmonic vibrational frequencies. The impact of core-correlation effects is discussed in detail. Rovibrational VCI calculations were used to simulate the gas phase spectra of these species and an in-depth analysis of the ν7 band of aminoborane is provided. Color-coding is used to reveal the identity of the individual progressions of the rovibrational transitions for this particular mode.
Collapse
Affiliation(s)
- Moritz Schneider
- Institute for Theoretical Chemistry, University of Stuttgart, Stuttgart, Germany
| | - Guntram Rauhut
- Institute for Theoretical Chemistry, University of Stuttgart, Stuttgart, Germany
| |
Collapse
|
25
|
Manzhos S, Ihara M. The loss of the property of locality of the kernel in high-dimensional Gaussian process regression on the example of the fitting of molecular potential energy surfaces. J Chem Phys 2023; 158:044111. [PMID: 36725493 DOI: 10.1063/5.0136156] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023] Open
Abstract
Kernel-based methods, including Gaussian process regression (GPR) and generally kernel ridge regression, have been finding increasing use in computational chemistry, including the fitting of potential energy surfaces and density functionals in high-dimensional feature spaces. Kernels of the Matern family, such as Gaussian-like kernels (basis functions), are often used which allow imparting to them the meaning of covariance functions and formulating GPR as an estimator of the mean of a Gaussian distribution. The notion of locality of the kernel is critical for this interpretation. It is also critical to the formulation of multi-zeta type basis functions widely used in computational chemistry. We show, on the example of fitting of molecular potential energy surfaces of increasing dimensionality, the practical disappearance of the property of locality of a Gaussian-like kernel in high dimensionality. We also formulate a multi-zeta approach to the kernel and show that it significantly improves the quality of regression in low dimensionality but loses any advantage in high dimensionality, which is attributed to the loss of the property of locality.
Collapse
Affiliation(s)
- Sergei Manzhos
- School of Materials and Chemical Technology, Tokyo Institute of Technology, Ookayama 2-12-1, Meguro-ku, Tokyo 152-8552, Japan
| | - Manabu Ihara
- School of Materials and Chemical Technology, Tokyo Institute of Technology, Ookayama 2-12-1, Meguro-ku, Tokyo 152-8552, Japan
| |
Collapse
|
26
|
Ulenikov ON, Bekhtereva ES, Gromova OV, Fomchenko AL, Nikolaeva NI, Glushkov PA, Bauerecker S. On the method and results of calculation of vibrational tetrahedral sub-level energies and resonance interactions in vibrational spectra of high symmetry molecules: CH 4 and GeH 4 as an application of the XY 4 (T d) molecule. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2023; 284:121796. [PMID: 36095857 DOI: 10.1016/j.saa.2022.121796] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/09/2022] [Revised: 08/22/2022] [Accepted: 08/25/2022] [Indexed: 06/15/2023]
Abstract
We report here the analytical description of one of the important problems in the study of XY4 (Td) molecules, namely, description of vibrational tetrahedral sub-level structures and resonance interactions caused by the high symmetry of a molecule. The results obtained are applied to description of the vibrational energy spectrum of the CH4 and GeH4 molecules.
Collapse
Affiliation(s)
- O N Ulenikov
- Research School of High-Energy Physics, National Research Tomsk Polytechnic University, Tomsk, 634050, Russia.
| | - E S Bekhtereva
- Research School of High-Energy Physics, National Research Tomsk Polytechnic University, Tomsk, 634050, Russia
| | - O V Gromova
- Research School of High-Energy Physics, National Research Tomsk Polytechnic University, Tomsk, 634050, Russia
| | - A L Fomchenko
- Research School of High-Energy Physics, National Research Tomsk Polytechnic University, Tomsk, 634050, Russia
| | - N I Nikolaeva
- Research School of High-Energy Physics, National Research Tomsk Polytechnic University, Tomsk, 634050, Russia
| | - P A Glushkov
- Research School of High-Energy Physics, National Research Tomsk Polytechnic University, Tomsk, 634050, Russia
| | - S Bauerecker
- Institut fur Physikalische und Theoretische Chemie, Technische Universitat Braunschweig, D - 38106, Braunschweig, Germany
| |
Collapse
|
27
|
Tschöpe M, Rauhut G. Convergence of series expansions in rovibrational configuration interaction (RVCI) calculations. J Chem Phys 2022; 157:234105. [PMID: 36550038 DOI: 10.1063/5.0129828] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022] Open
Abstract
Rotational and rovibrational spectra are a key in astrophysical studies, atmospheric science, pollution monitoring, and other fields of active research. The ab initio calculation of such spectra is fairly sensitive with respect to a multitude of parameters and all of them must be carefully monitored in order to yield reliable results. Besides the most obvious ones, i.e., the quality of the multidimensional potential energy surface and the vibrational wavefunctions, it is the representation of the μ-tensor within the Watson Hamiltonian, which has a significant impact on the desired line lists or simulated spectra. Within this work, we studied the dependence of high-resolution rovibrational spectra with respect to the truncation order of the μ-tensor within the rotational contribution and the Coriolis coupling operator of the Watson operator. Moreover, the dependence of the infrared intensities of the rovibrational transitions on an n-mode expansion of the dipole moment surface has been investigated as well. Benchmark calculations are provided for thioformaldehyde, which has already served as a test molecule in other studies and whose rovibrational spectrum was found to be fairly sensitive. All calculations rely on rovibrational configuration interaction theory and the discussed high-order terms of the μ-tensor are a newly implemented feature, whose theoretical basics are briefly discussed.
Collapse
Affiliation(s)
- Martin Tschöpe
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Guntram Rauhut
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| |
Collapse
|
28
|
Yu Q, Hammes-Schiffer S. Multidimensional Quantum Dynamical Simulation of Infrared Spectra under Polaritonic Vibrational Strong Coupling. J Phys Chem Lett 2022; 13:11253-11261. [PMID: 36448842 DOI: 10.1021/acs.jpclett.2c03245] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/17/2023]
Abstract
Recent experimental and theoretical studies demonstrate that the chemical reactivity of molecules can be modified inside an optical cavity. Here, we provide a theoretical framework for conducting multidimensional quantum simulations of the infrared (IR) spectra for molecules interacting with cavity modes. A single water molecule under polaritonic vibrational strong coupling serves as an illustrative example. Combined with accurate potential energy and dipole moment surfaces, our cavity vibrational self-consistent field/virtual state configuration interaction (cav-VSCF/VCI) approach can predict the IR spectra when the molecule is inside or outside the cavity. The spectral signatures of Rabi splittings and shifts of certain bands are found to be strongly dependent on the frequency and polarization direction of the cavity modes. Analyses of the simulated spectra show that polaritonic vibrational strong coupling can induce unconventional couplings among the molecule's vibrational modes, suggesting that intramolecular vibrational energy transfer can be significantly accelerated by the cavity.
Collapse
Affiliation(s)
- Qi Yu
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | | |
Collapse
|
29
|
Kelemen AK, Luber S. On the vibrations of formic acid predicted from first principles. Phys Chem Chem Phys 2022; 24:28109-28120. [PMID: 36385362 PMCID: PMC9710498 DOI: 10.1039/d2cp04417j] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2022] [Accepted: 11/09/2022] [Indexed: 08/26/2023]
Abstract
In this article, we review recent first principles, anharmonic studies on the molecular vibrations of gaseous formic acid in its monomer form. Transitions identified as fundamentals for both cis- and trans form reported in these studies are collected and supported by results from high-resolution experiments. Attention is given to the effect of coordinate coupling on the convergence of the computed vibrational states.
Collapse
Affiliation(s)
- Anna Klára Kelemen
- Department of Chemistry, University of Zurich, CH-8057 Zurich, Switzerland
| | - Sandra Luber
- Department of Chemistry, University of Zurich, CH-8057 Zurich, Switzerland.
| |
Collapse
|
30
|
Conte R, Nandi A, Qu C, Yu Q, Houston PL, Bowman JM. Semiclassical and VSCF/VCI Calculations of the Vibrational Energies of trans- and gauche-Ethanol Using a CCSD(T) Potential Energy Surface. J Phys Chem A 2022; 126:7709-7718. [PMID: 36240438 DOI: 10.1021/acs.jpca.2c06322] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A recent full-dimensional Δ-Machine learning potential energy surface (PES) for ethanol is employed in semiclassical and vibrational self-consistent field (VSCF) and virtual-state configuration interaction (VCI) calculations, using MULTIMODE, to determine the anharmonic vibrational frequencies of vibration for both the trans and gauche conformers of ethanol. Both semiclassical and VSCF/VCI energies agree well with the experimental data. We find significant mixing between the VSCF basis states due to Fermi resonances between bending and stretching modes. The same effects are also accurately described by the full-dimensional semiclassical calculations. These are the first high-level anharmonic calculations using a PES, in particular a "gold-standard" CCSD(T) one.
Collapse
Affiliation(s)
- Riccardo Conte
- Dipartimento di Chimica, Università degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| | - Apurba Nandi
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Chen Qu
- Independent Researcher, Toronto, Ontario M9B0E3, Canada
| | - Qi Yu
- Department of Chemistry Yale University, New Haven, Connecticut 06520, United States
| | - Paul L Houston
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States.,Department of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | - Joel M Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
31
|
Tyuterev V, Tashkun S, Rey M, Nikitin A. High-order contact transformations of molecular Hamiltonians: general approach, fast computational algorithm and convergence of ro-vibrational polyad models. Mol Phys 2022. [DOI: 10.1080/00268976.2022.2096140] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/16/2022]
Affiliation(s)
- Vladimir Tyuterev
- Groupe de Spectrométrie Moléculaire et Atmosphérique, Université de Reims, Reims, France
- Laboratory of Molecular Quantum Mechanics and Radiative Transfer, Tomsk State University, Tomsk, Russia
| | - Sergey Tashkun
- Laboratory of Theoretical Spectroscopy, V.E. Zuev Institute of Atmospheric Optics, Siberian Branch, Russian Academy of Sciences, Tomsk, Russia
| | - Michael Rey
- Groupe de Spectrométrie Moléculaire et Atmosphérique, Université de Reims, Reims, France
| | - Andrei Nikitin
- Laboratory of Theoretical Spectroscopy, V.E. Zuev Institute of Atmospheric Optics, Siberian Branch, Russian Academy of Sciences, Tomsk, Russia
| |
Collapse
|
32
|
Carvalho JR, Vidal LN. Calculation of absolute Raman scattering cross-sections using vibrational self-consistent field/vibrational configuration interaction wave functions. J Comput Chem 2022; 43:1484-1494. [PMID: 35731622 DOI: 10.1002/jcc.26951] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2022] [Revised: 05/21/2022] [Accepted: 06/02/2022] [Indexed: 11/12/2022]
Abstract
In the present study, the differential scattering cross-sections, depolarization ratios and Raman shifts of small molecular systems are obtained from configuration iteration wave functions of vibrational self-consistent field (VSCF) states. The transition polarizabilities were modeled using the Placzek approximation, neglecting those contributions not arising from the electric dipole mechanism. This theoretical approach is considered a good approximation for samples that absorb in the UV range if the excitation radiation falls in the visible region, as is the case of the molecules selected for the present study, namely: water, methane, and acetylene. Potential energy and electronic polarizability surfaces are calculated by the CCSD(T) and CC3 methods with aug-cc-p(C)V(T,Q,5)Z basis sets. The vibrational Hamiltonian includes the vibrational angular momentum contribution of the Watson kinetic energy operator. As expected, due to the variational nature of the VSCF and vibrational configuration interaction (VCI) methods, the Raman transition wavenumbers are substantially improved over the harmonic predictions. Surprisingly, the scattering cross-sections obtained using the harmonic approximation or the VSCF method better agrees with the experimental values than those cross-sections predicted using VCI wave functions. The more significant deviations of the VCI results from the experimental reference may be related to the significant uncertainties of the measured cross-sections. Still, it may also indicate that the VCI Raman transition moments may require a more accurate description of the electronic polarizability surface. Finally, the depolarization ratios calculated for H2 O and C2 D2 using harmonic and VCI wave functions have similar accuracy, whereas, for C2 H2 and C2 HD, the VCI results are more accurate.
Collapse
Affiliation(s)
- Jhonatas R Carvalho
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas, USA
| | - Luciano N Vidal
- Departamento Acadêmico de Química e Biologia, Universidade Tecnológica Federal do Paraná, Curitiba, Paraná, Brazil
| |
Collapse
|
33
|
Töpfer K, Upadhyay M, Meuwly M. Quantitative molecular simulations. Phys Chem Chem Phys 2022; 24:12767-12786. [PMID: 35593769 PMCID: PMC9158373 DOI: 10.1039/d2cp01211a] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2022] [Accepted: 04/30/2022] [Indexed: 11/21/2022]
Abstract
All-atom simulations can provide molecular-level insights into the dynamics of gas-phase, condensed-phase and surface processes. One important requirement is a sufficiently realistic and detailed description of the underlying intermolecular interactions. The present perspective provides an overview of the present status of quantitative atomistic simulations from colleagues' and our own efforts for gas- and solution-phase processes and for the dynamics on surfaces. Particular attention is paid to direct comparison with experiment. An outlook discusses present challenges and future extensions to bring such dynamics simulations even closer to reality.
Collapse
Affiliation(s)
- Kai Töpfer
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland.
| | - Meenu Upadhyay
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland.
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland.
| |
Collapse
|
34
|
Kędziera D, Rauhut G, Császár AG. Structure, energetics, and spectroscopy of the chromophores of HHe+n, H 2He+n, and He+n clusters and their deuterated isotopologues. Phys Chem Chem Phys 2022; 24:12176-12195. [PMID: 35543594 DOI: 10.1039/d1cp05535f] [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
The linear molecular ions H2He+, HHe+2, and He+3 are the central units (chromophores) of certain He-solvated complexes of the H2He+n, HHe+n, and He+n families, respectively. These are complexes which do exist, according to mass-spectrometry studies, up to very high n values. Apparently, for some of the H2He+n and He+n complexes, the linear symmetric tetratomic H2He+2 and the diatomic He+2 cations, respectively, may also be the central units. In this study, definitive structures, relative energies, zero-point vibrational energies, and (an)harmonic vibrational fundamentals, and, in some cases, overtones and combination bands, are established mostly for the triatomic chromophores. The study is also extended to the deuterated isotopologues D2He+, DHe+2, and D2He+2. To facilitate and improve the electronic-structure computations performed, new atom-centered, fixed-exponent, Gaussian-type basis sets called MAX, with X = T(3), Q(4), P(5), and H(6), are designed for the H and He atoms. The focal-point-analysis (FPA) technique is employed to determine definitive relative energies with tight uncertainties for reactions involving the molecular ions. The FPA results determined include the 0 K proton and deuteron affinities of the 4He atom, 14 875(9) cm-1 [177.95(11) kJ mol-1] and 15 229(8) cm-1 [182.18(10) kJ mol-1], respectively, the dissociation energies of the He+2 → He+ + He, HHe+2 → HHe+ + He, and He+3 → He+2 + He reactions, 19 099(13) cm-1 [228.48(16) kJ mol-1], 3948(7) cm-1 [47.23(8) kJ mol-1], and 1401(12) cm-1 [16.76(14) kJ mol-1], respectively, the dissociation energy of the DHe+2 → DHe+ + He reaction, 4033(6) cm-1 [48.25(7) kJ mol-1], the isomerization energy between the two linear isomers of the [H, He, He]+ system, 3828(40) cm-1 [45.79(48) kJ mol-1], and the dissociation energies of the H2He+ → H+2 + He and the H2He+2 → H2He+ + He reactions, 1789(4) cm-1 [21.40(5) kJ mol-1] and 435(6) cm-1 [5.20(7) kJ mol-1], respectively. The FPA estimates of the first dissociation energy of D2He+ and D2He+2 are 1986(4) cm-1 [23.76(5) kJ mol-1] and 474(5) cm-1 [5.67(6) kJ mol-1], respectively. Determining the vibrational fundamentals of the triatomic chromophores with second-order vibrational perturbation theory (VPT2) and vibrational configuration interaction (VCI) techniques, both built around the Eckart-Watson Hamiltonian, proved unusually challenging. For the species studied, VPT2 has difficulties yielding dependable results, in some cases even for the fundamentals of the H-containing molecular cations, while carefully executed VCI computations yield considerably improved spectroscopic results. In a few cases unusually large anharmonic corrections to the fundamentals, on the order of 15% of the harmonic value, have been observed.
Collapse
Affiliation(s)
- Dariusz Kędziera
- Faculty of Chemistry, Nicolaus Copernicus University in Torun, 87-100 Torun, Poland.
| | - Guntram Rauhut
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany.
| | - Attila G Császár
- Laboratory of Molecular Structure and Dynamics, Institute of Chemistry, ELTE Eötvös Loránd University and MTA-ELTE Complex Chemical Systems Research Group, Pázmány Péter sétány 1/A, H-1117 Budapest, Hungary.
| |
Collapse
|
35
|
Begušić T, Tapavicza E, Vaníček J. Applicability of the Thawed Gaussian Wavepacket Dynamics to the Calculation of Vibronic Spectra of Molecules with Double-Well Potential Energy Surfaces. J Chem Theory Comput 2022; 18:3065-3074. [PMID: 35420803 DOI: 10.1021/acs.jctc.2c00030] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Simulating vibrationally resolved electronic spectra of anharmonic systems, especially those involving double-well potential energy surfaces, often requires expensive quantum dynamics methods. Here, we explore the applicability and limitations of the recently proposed single-Hessian thawed Gaussian approximation for the simulation of spectra of systems with double-well potentials, including 1,2,4,5-tetrafluorobenzene, ammonia, phosphine, and arsine. This semiclassical wavepacket approach is shown to be more robust and to provide more accurate spectra than the conventional harmonic approximation. Specifically, we identify two cases in which the Gaussian wavepacket method is especially useful due to the breakdown of the harmonic approximation: (i) when the nuclear wavepacket is initially at the top of the potential barrier but delocalized over both wells, e.g., along a low-frequency mode, and (ii) when the wavepacket has enough energy to classically go over the low potential energy barrier connecting the two wells. The method is efficient and requires only a single classical ab initio molecular dynamics trajectory, in addition to the data required to compute the harmonic spectra. We also present an improved algorithm for computing the wavepacket autocorrelation function, which guarantees that the evaluated correlation function is continuous for arbitrary size of the time step.
Collapse
Affiliation(s)
- Tomislav Begušić
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Enrico Tapavicza
- Department of Chemistry and Biochemistry, California State University, Long Beach, 1250 Bellflower Boulevard, Long Beach, California 90840-9507, United States
| | - Jiří Vaníček
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
36
|
Schröder B, Rauhut G. Comparison of body definitions for incremental vibrational configuration interaction theory (iVCI). J Chem Phys 2022; 156:174103. [DOI: 10.1063/5.0085082] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Within incremental vibrational configuration interaction theory (iVCI), the vibrational state energy is determined by means of a many-body expansion, i.e., it is a sum of terms of increasing order, which allow for an embarrassingly parallel evaluation. The convergence of this expansion depends strongly on the definition of the underlying bodies, which essentially decompose the correlation space into fragments. The different definitions considered here comprise mode-based bodies, excitation level-based bodies, and energy-based bodies. An analysis of the convergence behavior revealed that accounting for resonances within these definitions is mandatory and leads to a substantial improvement of the convergence, that is, the expansions can be truncated at lower orders. Benchmark calculations and systematic comparisons of the different body definitions for a small set of molecules, i.e., ketene, ethene, and diborane, have been conducted to study the overall performance of these iVCI implementations with respect to accuracy and central processing unit time.
Collapse
Affiliation(s)
- Benjamin Schröder
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Guntram Rauhut
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| |
Collapse
|
37
|
Motta M, Rice JE. Emerging quantum computing algorithms for quantum chemistry. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2022. [DOI: 10.1002/wcms.1580] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Mario Motta
- IBM Quantum, IBM Research‐Almaden San Jose California USA
| | - Julia E. Rice
- IBM Quantum, IBM Research‐Almaden San Jose California USA
| |
Collapse
|
38
|
Vogt E, Simkó I, Császár AG, Kjaergaard HG. Reduced-dimensional vibrational models of the water dimer. J Chem Phys 2022; 156:164304. [PMID: 35490001 DOI: 10.1063/5.0090013] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A model based on the finite-basis representation of a vibrational Hamiltonian expressed in internal coordinates is developed. The model relies on a many-mode, low-order expansion of both the kinetic energy operator and the potential energy surface (PES). Polyad truncations and energy ceilings are used to control the size of the vibrational basis to facilitate accurate computations of the OH stretch and HOH bend intramolecular transitions of the water dimer (H2 16O)2. Advantages and potential pitfalls of the applied approximations are highlighted. The importance of choices related to the treatment of the kinetic energy operator in reduced-dimensional calculations and the accuracy of different water dimer PESs are discussed. A range of different reduced-dimensional computations are performed to investigate the wavenumber shifts in the intramolecular transitions caused by the coupling between the intra- and intermolecular modes. With the use of symmetry, full 12-dimensional vibrational energy levels of the water dimer are calculated, predicting accurately the experimentally observed intramolecular fundamentals. It is found that one can also predict accurate intramolecular transition wavenumbers for the water dimer by combining a set of computationally inexpensive reduced-dimensional calculations, thereby guiding future effective-Hamiltonian treatments.
Collapse
Affiliation(s)
- Emil Vogt
- Department of Chemistry, University of Copenhagen, Universitetsparken 5, DK-2100 Copenhagen Ø, Denmark
| | - Irén Simkó
- Laboratory of Molecular Structure and Dynamics, Institute of Chemistry, ELTE Eötvös Loránd University, Pázmány Péter Sétány 1/A, H-1117 Budapest, Hungary
| | - Attila G Császár
- Laboratory of Molecular Structure and Dynamics, Institute of Chemistry, ELTE Eötvös Loránd University, Pázmány Péter Sétány 1/A, H-1117 Budapest, Hungary
| | - Henrik G Kjaergaard
- Department of Chemistry, University of Copenhagen, Universitetsparken 5, DK-2100 Copenhagen Ø, Denmark
| |
Collapse
|
39
|
Erfort S, Tschoepe M, Rauhut G. Efficient and Automated Quantum Chemical Calculation of Rovibrational Nonresonant Raman Spectra. J Chem Phys 2022; 156:124102. [DOI: 10.1063/5.0087359] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
An outline of a newly developed program for the simulation of rovibrational nonresonant Raman spectra is presented. This program is an extension of our recently developed code for rovibrational infrared spectra [J. Chem Phys. 152 (2020) 244104] and relies on vibrational wavefunctions from variational configuration interaction theory to allow for an almost fully automated calculation of such spectra in pure ab initio fashion. Due to efficient contraction schemes this program requires modest computational resources and it can be controlled by only a few lines of input. As the required polarizability surfaces are also computed in an automated fashion, this implementation enables the routine application to small molecules. For demonstrating its capabilities, benchmark calculations for water H216O are compared to reference data and spectra for the beryllium dihydride dimer, Be2H4 (D2h), are predicted. The inversion symmetry of the D2h systems lead to complementary infrared and Raman spectra, which are needed both for a comprehensive investigation of this system.
Collapse
Affiliation(s)
- Sebastian Erfort
- Institute for Theoretical Chemistry, University of Stuttgart Faculty of Chemistry, Germany
| | | | - Guntram Rauhut
- Institut fuer Theoretische Chemie, University of Stuttgart Faculty of Chemistry, Germany
| |
Collapse
|
40
|
Roy TK. Performance of Vibrational Self-Consistent Field Theory for Accurate Potential Energy Surfaces: Fundamentals, Excited States, and Intensities. J Phys Chem A 2022; 126:608-622. [PMID: 35050620 DOI: 10.1021/acs.jpca.1c09989] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The performance of vibrational structure calculations beyond harmonic approximation in the framework of the vibrational self-consistent field method with second-order perturbation corrections (VSCF-PT2) is investigated in conjunction with very accurate potential energy surfaces (PESs) given by various coupled-cluster electronic structure theories. The quality of anharmonic calculations depends on the accuracy of the underlying multidimensional PES obtained from its functional form, which is given by the level of electronic structure theory. Two such highest levels of typical coupled-cluster electronic structure methods, CCSD and the ″gold standard″ CCSD(T), along with their variants such as CCD, CR-CCL (completely renormalized CR-CC(2,3) approach), and CCSD(TQ) are tested for the construction of accurate anharmonic potentials without any fitting or ad hoc scaling and using cc-pVTZ basis sets. The accuracy of VSCF-PT2 theory in comparison to experimental values is tested for a series of 16 molecules with 135 fundamental bands, 64 overtones, and combination bands and also for 39 intensities. It is found that CCD and CCSD bind the potential tighter than CCSD(T) and the computed VSCF-PT2 transitions are more blue-shifted showing higher deviation from the experiment. In general, VSCF-PT2 results computed at the CCSD(T) potential offer a good cost/accuracy ratio, with the mean absolute deviation and the mean absolute percentage error with the experiment being ∼16 cm-1 and 1.38, respectively, for fundamentals. Additionally, while the CR-CCL and CCSD(TQ) methods offer similar levels of accuracies as compared to CCSD(T), the former offers a better accuracy/cost ratio than the latter and is a suitable alternative to CCSD(T).
Collapse
Affiliation(s)
- Tapta Kanchan Roy
- Department of Chemistry and Chemical Sciences, Central University of Jammu, Rahya-Suchani (Bagla), Samba, Jammu and Kashmir 181143, India
| |
Collapse
|
41
|
Manzhos S, Sasaki E, Ihara M. Easy representation of multivariate functions with low-dimensional terms via Gaussian process regression kernel design: applications to machine learning of potential energy surfaces and kinetic energy densities from sparse data. MACHINE LEARNING: SCIENCE AND TECHNOLOGY 2022. [DOI: 10.1088/2632-2153/ac4949] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
Abstract
We show that Gaussian process regression (GPR) allows representing multivariate functions with low-dimensional terms via kernel design. When using a kernel built with high-dimensional model representation (HDMR), one obtains a similar type of representation as the previously proposed HDMR-GPR scheme while being faster and simpler to use. We tested the approach on cases where highly accurate machine learning is required from sparse data by fitting potential energy surfaces and kinetic energy densities.
Collapse
|
42
|
Mathea T, Petrenko T, Rauhut G. Advances in vibrational configuration interaction theory - part 2: Fast screening of the correlation space. J Comput Chem 2022; 43:6-18. [PMID: 34651704 DOI: 10.1002/jcc.26764] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2021] [Revised: 09/20/2021] [Accepted: 09/20/2021] [Indexed: 11/11/2022]
Abstract
For larger molecules, the computational demands of configuration selective vibrational configuration interaction theory (cs-VCI) are usually dominated by the configuration selection process, which commonly is based on second order vibrational Møller-Plesset perturbation (VMP2) theory. Here we present two techniques, which lead to substantial accelerations of such calculations while retaining the desired high accuracy of the final results. The first one introduces the concept of configuration classes, which allows for a highly efficient exploitation of the analogs of the Slater-Condon rules in vibrational structure calculations with large correlation spaces. The second approach uses a VMP2 like vector for augmenting the targeted vibrational wavefunction within the selection of configurations and thus avoids any intermediate diagonalization steps. The underlying theory is outlined and benchmark calculations are provided for highly correlated vibrational states of several molecules.
Collapse
Affiliation(s)
- Tina Mathea
- Institute for Theoretical Chemistry, University of Stuttgart, Stuttgart, Germany
| | - Taras Petrenko
- Institute for Theoretical Chemistry, University of Stuttgart, Stuttgart, Germany
| | - Guntram Rauhut
- Institute for Theoretical Chemistry, University of Stuttgart, Stuttgart, Germany
| |
Collapse
|
43
|
Huang Q, Matsuda Y, Eguchi R, Fujii A, Kuo J. Understanding Fermi resonances behind the complex vibrational spectra of the methyl groups in simple alcohol, thiol, and their ethers. J CHIN CHEM SOC-TAIP 2022. [DOI: 10.1002/jccs.202100281] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Affiliation(s)
- Qian‐Rui Huang
- Institute of Atomic and Molecular Sciences Academia Sinica Taipei Taiwan
| | - Yoshiyuki Matsuda
- Department of Chemistry, Graduate School of Science Tohoku University Sendai Japan
| | - Riku Eguchi
- Department of Chemistry, Graduate School of Science Tohoku University Sendai Japan
| | - Asuka Fujii
- Department of Chemistry, Graduate School of Science Tohoku University Sendai Japan
| | - Jer‐Lai Kuo
- Institute of Atomic and Molecular Sciences Academia Sinica Taipei Taiwan
| |
Collapse
|
44
|
Manzhos S, Ihara M. Computational vibrational spectroscopy of molecule-surface interactions: what is still difficult and what can be done about it. Phys Chem Chem Phys 2022; 24:15158-15172. [DOI: 10.1039/d2cp01389d] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Interactions of molecules with solid surfaces are responsible for key functionalities for a range of currently actively pursued technologies, including heterogeneous catalysis for synthesis or decomposition of molecules, sensitization, surface...
Collapse
|
45
|
Omodemi O, Sprouse S, Herbert D, Kaledin M, Kaledin AL. On the Cartesian Representation of the Molecular Polarizability Tensor Surface by Polynomial Fitting to Ab Initio Data. J Chem Theory Comput 2021; 18:37-45. [PMID: 34958587 DOI: 10.1021/acs.jctc.1c01015] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We describe an approach to constructing an analytic Cartesian representation of the molecular dipole polarizability tensor surface in terms of polynomials in interatomic distances with a training set of ab initio data points obtained from a molecular dynamics (MD) simulation or by any other available means. The proposed formulation is based on a perturbation treatment of the unmodified point dipole polarizability model of Applequist [ J. Am. Chem. Soc. 1972, 94, 2952] and is shown here to be, by construction (i) free of short-range or other singularities or discontinuities, (ii) symmetric and translationally invariant, and (iii) nonreliant on a body-fixed coordinate system. Permutational invariance of like nuclei is demonstrated to be readily applicable, making this approach useful for highly fluxional and reactive systems. Derivation of the method is described in detail, adding brief didactic numerical examples of H2 and H2O and concluding with an MD simulation of the Raman spectrum of H5O2+ at 300 K with the polarizability tensor fitted to CCSD(T)/aug-cc-pVTZ data obtained using the HBB-4B potential [ J. Chem. Phys. 2005, 122, 044308].
Collapse
Affiliation(s)
- Oluwaseun Omodemi
- Department of Chemistry & Biochemistry, Kennesaw State University, 370 Paulding Avenue NW, Box # 1203, Kennesaw, Georgia 30144, United States
| | - Sarah Sprouse
- Department of Chemistry & Biochemistry, Kennesaw State University, 370 Paulding Avenue NW, Box # 1203, Kennesaw, Georgia 30144, United States
| | - Destyni Herbert
- Department of Chemistry & Biochemistry, Kennesaw State University, 370 Paulding Avenue NW, Box # 1203, Kennesaw, Georgia 30144, United States
| | - Martina Kaledin
- Department of Chemistry & Biochemistry, Kennesaw State University, 370 Paulding Avenue NW, Box # 1203, Kennesaw, Georgia 30144, United States
| | - Alexey L Kaledin
- Cherry L. Emerson Center for Scientific Computation, Emory University, 1515 Dickey Drive, Atlanta, Georgia 30322, United States
| |
Collapse
|
46
|
Simmons J, Carrington Jr. T. Using collocation and solutions for a sum-of-product potential to compute vibrational energy levels for general potentials. Chem Phys Lett 2021. [DOI: 10.1016/j.cplett.2021.138967] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
47
|
Han J, Grofe A, Gao J. Variational Energy Decomposition Analysis of Charge-Transfer Interactions between Metals and Ligands in Carbonyl Complexes. Inorg Chem 2021; 60:14060-14071. [PMID: 34460236 DOI: 10.1021/acs.inorgchem.1c01367] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Variational energy decomposition analyses have been presented to quantify the σ-dative, ligand-to-metal forward charge transfer (CT) and the π-conjugative, metal-to-ligand backward charge delocalization on a series of isolelectronic transition-metal carbonyl complexes M(CO)6, including M = Ti2-, V-, Cr, Mn+, and Fe2+. Although the qualitative features of these energy terms are understood, well-defined quantitative studies have been scarce. Consistent with early findings, electrostatic and Pauli exchange effects play a key role in σ-donation, resulting in blue shifts in ligand vibrational frequency in the complex geometries. Excluding chemical bonding interactions between the CO ligand and the metal fragments in the energy decomposition analysis, we found that loosely bound electrostatic complexes can be formed at a longer metal-to-ligand distance due to the exponential decay of Pauli exchange. In all complexes, the overall binding stabilization can be attributed to CT effects, with opposing trends between σ-donation and π-back bonding that follows an order of Ti2- (4.4) > V1- (2.6) > Cr (1.5) > Mn1+ (1.1) > Fe2+ (0.5) in π-to-σ CT ratio. These electronic and energetic features are mirrored in the vibrational frequency shifts induced by different factors. The present investigation may help stimulate the use of energy decomposition techniques to understand the structure and activity of metallocatalysts using density functional theory.
Collapse
Affiliation(s)
- Jingting Han
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, Jilin University, Changchun, Jilin Province 130023, China
| | - Adam Grofe
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, Jilin University, Changchun, Jilin Province 130023, China.,Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, China
| | - Jiali Gao
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, China.,Beijing University Shenzhen Graduate School, Shenzhen 518055, China.,Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| |
Collapse
|
48
|
Yagi K, Sugita Y. Anharmonic Vibrational Calculations Based on Group-Localized Coordinates: Applications to Internal Water Molecules in Bacteriorhodopsin. J Chem Theory Comput 2021; 17:5007-5020. [PMID: 34296615 PMCID: PMC10986902 DOI: 10.1021/acs.jctc.1c00060] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
An efficient anharmonic vibrational method is developed exploiting the locality of molecular vibration. Vibrational coordinates localized to a group of atoms are employed to divide the potential energy surface (PES) of a system into intra- and inter-group contributions. Then, the vibrational Schrödinger equation is solved based on a PES, in which the inter-group coupling is truncated at the harmonic level while accounting for the intra-group anharmonicity. The method is applied to a pentagonal hydrogen bond network (HBN) composed of internal water molecules and charged residues in a membrane protein, bacteriorhodopsin. The PES is calculated by the quantum mechanics/molecular mechanics (QM/MM) calculation at the level of B3LYP-D3/aug-cc-pVDZ. The infrared (IR) spectrum is computed using a set of coordinates localized to each water molecule and amino acid residue by second-order vibrational quasi-degenerate perturbation theory (VQDPT2). Benchmark calculations show that the proposed method yields the N-D/O-D stretching frequencies with an error of 7 cm-1 at the cost reduced by more than five times. In contrast, the harmonic approximation results in a severe error of 150 cm-1. Furthermore, the size of QM regions is carefully assessed to find that the QM regions should include not only the pentagonal HBN itself but also its HB partners. VQDPT2 calculations starting from transient structures obtained by molecular dynamics simulations have shown that the structural sampling has a significant impact on the calculated IR spectrum. The incorporation of anharmonicity, sufficiently large QM regions, and structural samplings are of essential importance to reproduce the experimental IR spectrum. The computational spectrum paves the way for decoding the IR signal of strong HBNs and helps elucidate their functional roles in biomolecules.
Collapse
Affiliation(s)
- Kiyoshi Yagi
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | - Yuji Sugita
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
- Computational
Biophysics Research Team, RIKEN Center for
Computational Science, 7-1-26 Minatojima-Minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
- Laboratory
for Biomolecular Function Simulation, RIKEN
Center for Biosystems Dynamics Research, 1-6-5 Minatojima-Minamimachi,
Chuo-ku, Kobe, Hyogo 650-0047, Japan
| |
Collapse
|
49
|
Qu C, Houston PL, Conte R, Nandi A, Bowman JM. MULTIMODE Calculations of Vibrational Spectroscopy and 1d Interconformer Tunneling Dynamics in Glycine Using a Full-Dimensional Potential Energy Surface. J Phys Chem A 2021; 125:5346-5354. [PMID: 34110169 DOI: 10.1021/acs.jpca.1c03738] [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/30/2022]
Abstract
A full-dimensional, permutationally invariant polynomial potential energy surface for glycine recently reported (R. Conte et al., J. Chem. Phys. 2020, 153, 244301) is used with the code MULTIMODE to determine the IR absorption spectra for Conformers I and II using a new separable dipole moment function. The calculated spectra agree well with the experimental ones. The full-dimensional nature of the potential allows us also to examine dynamical results, such as tunneling rates. Remarkably, using a one-dimensional path based on the potential energy surface to estimate the tunneling rate from Conformer VI to Conformer I, good agreement is found with the recent experimental measurement. Finally a brief comparison of our potential energy surface with a recently reported sGDML one is made.
Collapse
Affiliation(s)
- Chen Qu
- Department of Chemistry & Biochemistry, University of Maryland, College Park, Maryland 20742, United States
| | - Paul L Houston
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States.,Department of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | - Riccardo Conte
- Dipartimento di Chimica, Università degli Studi di Milano, via Golgi 19, 20133 Milano, Italy
| | - Apurba Nandi
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Joel M Bowman
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
50
|
Kuo JL. Disentangling the Complex Vibrational Spectra of Hydrogen-Bonded Clusters of 2-Pyridone with Ab Initio Structural Search and Anharmonic Analysis. J Phys Chem A 2021; 125:4306-4312. [PMID: 33998811 DOI: 10.1021/acs.jpca.1c00587] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Vibrational spectra of 1:1 clusters of 2-pyridone (2PY) with water, ammonia, and other hydrogen bond-forming molecules have been measured by several experimental groups over the past two decades. Complex vibrational signatures associated with the N-H stretching fundamental at 3 μm are often observed. Several anharmonic coupling schemes have been proposed; however, the origin of these commonly seen complex features remains unclear. In this work, we present our theoretical analysis on the structure and vibrational spectra of these clusters using ab initio random search and ab initio anharmonic algorithms, respectively. Low-energy conformers were found to be hydrogen-bonded clusters and their vibrational spectra at 3 μm were simulated with ab initio anharmonic algorithms. We demonstrate that simple anharmonic mechanisms of Fermi resonance (FR), coupling between NH stretching modes, and overtone/combinations of skeleton modes of 2PY can lead to the complex vibrational signatures observed experimentally. Since this vibrational coupling scheme is inherent to the cis-amides with adjacent N-H and C═O groups when a hydrogen bond is formed with 2PY as the donor and acceptor, we believe that such a phenomenon is general to other hydrogen-bonding systems with the same functional group.
Collapse
Affiliation(s)
- Jer-Lai Kuo
- Institute of Atomic and Molecular Sciences, Academia Sinica, Taipei 10617, Taiwan
| |
Collapse
|