1
|
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
|
2
|
Majland M, Berg Jensen R, Greisen Højlund M, Thomas Zinner N, Christiansen O. Optimizing the number of measurements for vibrational structure on quantum computers: coordinates and measurement schemes. Chem Sci 2023; 14:7733-7742. [PMID: 37476724 PMCID: PMC10355095 DOI: 10.1039/d3sc01984e] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2023] [Accepted: 06/07/2023] [Indexed: 07/22/2023] Open
Abstract
One of the primary challenges prohibiting demonstrations of practical quantum advantages for near-term devices amounts to excessive measurement overheads for estimating relevant physical quantities such as ground state energies. However, with major differences between the electronic and vibrational structures of molecules, the question of how the resource requirements of computing anharmonic, vibrational states can be reduced remains relatively unexplored compared to its electronic counterpart. Importantly, bosonic commutation relations, distinguishable Hilbert spaces and vibrational coordinates allow manipulations of the vibrational system that can be exploited to minimize resource requirements. In this work, we investigate the impact of different coordinate systems and measurement schemes on the number of measurements needed to estimate anharmonic, vibrational states for a variety of three-mode (six-mode) molecules. We demonstrate an average of 3-fold (1.5-fold), with up to 7-fold (2.5-fold), reduction in the number of measurements required by employing appropriate coordinate transformations, based on an automized construction of qubit Hamiltonians from a conventional vibrational structure program.
Collapse
Affiliation(s)
- Marco Majland
- Kvantify Aps DK-2300 Copenhagen S Denmark
- Department of Physics and Astronomy, Aarhus University DK-8000 Aarhus C Denmark
- Department of Chemistry, Aarhus University DK-8000 Aarhus C Denmark
| | - Rasmus Berg Jensen
- Department of Physics and Astronomy, Aarhus University DK-8000 Aarhus C Denmark
- Department of Chemistry, Aarhus University DK-8000 Aarhus C Denmark
| | | | - Nikolaj Thomas Zinner
- Kvantify Aps DK-2300 Copenhagen S Denmark
- Department of Physics and Astronomy, Aarhus University DK-8000 Aarhus C Denmark
| | - Ove Christiansen
- Kvantify Aps DK-2300 Copenhagen S Denmark
- Department of Chemistry, Aarhus University DK-8000 Aarhus C Denmark
| |
Collapse
|
3
|
Hanson-Heine MWD. Static Electron Correlation in Anharmonic Molecular Vibrations: A Hybrid TAO-DFT Study. J Phys Chem A 2022; 126:7273-7282. [PMID: 36164938 PMCID: PMC9574917 DOI: 10.1021/acs.jpca.2c05881] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
Hybrid thermally-assisted-occupation density functional theory is used to examine the effects of static electron correlation on the prediction of a benchmark set of experimentally observed molecular vibrational frequencies. The B3LYP and B97-1 thermally-assisted-occupation measure of static electron correlation is important for describing the vibrations of many of the molecules that make up several popular test sets of experimental data. Shifts are seen for known multireference systems and for many molecules containing atoms from the second row of the periodic table of elements. Several molecules only show significant shifts in select vibrational modes, and significant improvements are seen for the prediction of hydrogen stretching frequencies throughout the test set.
Collapse
|
4
|
Hino K, Kurashige Y. Matrix Product State Formulation of the MCTDH Theory in Local Mode Representations for Anharmonic Potentials. J Chem Theory Comput 2022; 18:3347-3356. [PMID: 35606892 DOI: 10.1021/acs.jctc.2c00243] [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
The matrix product state formulation of the multiconfiguration time-dependent Hartree theory, MPS-MCTDH, reported previously [Kurashige, J. Chem. Phys. 2018, 19, 194114] is extended to realistic anharmonic potentials with n-mode representations beyond the linear vibronic coupling model. For realistic vibrational potentials, the local mode representation should give a more compact representation of the potentials, i.e., lowering the dimensionality of the entanglements, than the normal coordinates, and the MPS-MCTDH formulation should work more efficiently and maintain the accuracy with a small bond dimension of the MPS ansatz. In fact, it was confirmed that the use of the local coordinates made the interaction matrices diagonal dominant and the number of terms in the n-body expansion of the potentials was significantly reduced. The method was applied to the IR spectrum of the CH2O molecule, the zero-point energies, and the vibrational energy redistribution dynamics of polyenes C2nH2n+2. The results showed that the efficiency of the MPS-MCTDH method is significantly accelerated by the use of local coordinates even if the long-range interactions are included in the potential.
Collapse
Affiliation(s)
- Kentaro Hino
- Department of Chemistry, Graduate School of Science, Kyoto University, Kitashirakawa Oiwake-cho, Sakyo-ku, Kyoto 606-8502, Japan
| | - Yuki Kurashige
- Department of Chemistry, Graduate School of Science, Kyoto University, Kitashirakawa Oiwake-cho, Sakyo-ku, Kyoto 606-8502, Japan
| |
Collapse
|
5
|
Changala PB. Vibronic mean-field and perturbation theory for Jahn-Teller and pseudo-Jahn-Teller molecules. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1962556] [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]
|
6
|
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: 4] [Impact Index Per Article: 1.3] [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
|
7
|
Hanson-Heine MWD. Reduced Two-Electron Interactions in Anharmonic Molecular Vibrational Calculations Involving Localized Normal Coordinates. J Chem Theory Comput 2021; 17:4383-4391. [PMID: 34087068 DOI: 10.1021/acs.jctc.1c00314] [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
Spatially localized vibrational normal mode coordinates are shown to reduce the importance of calculating the full set of two-electron terms in the molecular electronic Schrödinger equation. Electron correlation and dispersion interactions become less significant in (E,E)-1,3,5,7-octatetraene vibrational self-consistent field calculations when displacing remote atoms along multiple coordinates. Electron correlation interactions between spatially remote modes are also found to be less important compared to their corresponding uncorrelated interaction terms. Attenuation of the Coulomb operator indicates that the two-electron terms between remote electrons become less important for accurately describing the strongly contributing mode-coupling terms between sets of localized vibrational modes.
Collapse
|
8
|
Mathea T, Petrenko T, Rauhut G. VCI Calculations Based on Canonical and Localized Normal Coordinates for Non-Abelian Molecules: Accurate Assignment of the Vibrational Overtones of Allene. J Phys Chem A 2021; 125:990-998. [DOI: 10.1021/acs.jpca.0c10429] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Tina Mathea
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Taras Petrenko
- 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
|
9
|
Ollitrault PJ, Baiardi A, Reiher M, Tavernelli I. Hardware efficient quantum algorithms for vibrational structure calculations. Chem Sci 2020; 11:6842-6855. [PMID: 32874524 PMCID: PMC7448527 DOI: 10.1039/d0sc01908a] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2020] [Accepted: 06/02/2020] [Indexed: 11/21/2022] Open
Abstract
We introduce a framework for the calculation of ground and excited state energies of bosonic systems suitable for near-term quantum devices and apply it to molecular vibrational anharmonic Hamiltonians. Our method supports generic reference modal bases and Hamiltonian representations, including the ones that are routinely used in classical vibrational structure calculations. We test different parametrizations of the vibrational wavefunction, which can be encoded in quantum hardware, based either on heuristic circuits or on the bosonic Unitary Coupled Cluster Ansatz. In particular, we define a novel compact heuristic circuit and demonstrate that it provides a good compromise in terms of circuit depth, optimization costs, and accuracy. We evaluate the requirements, number of qubits and circuit depth, for the calculation of vibrational energies on quantum hardware and compare them with state-of-the-art classical vibrational structure algorithms for molecules with up to seven atoms.
Collapse
Affiliation(s)
- Pauline J Ollitrault
- IBM Quantum , IBM Research - Zurich , Säumerstrasse 4 , 8803 Rüschlikon , Switzerland .
- Laboratory of Physical Chemistry , ETH Zürich , Vladimir-Prelog-Weg 2 , 8093 Zürich , Switzerland .
| | - Alberto Baiardi
- Laboratory of Physical Chemistry , ETH Zürich , Vladimir-Prelog-Weg 2 , 8093 Zürich , Switzerland .
| | - Markus Reiher
- Laboratory of Physical Chemistry , ETH Zürich , Vladimir-Prelog-Weg 2 , 8093 Zürich , Switzerland .
| | - Ivano Tavernelli
- IBM Quantum , IBM Research - Zurich , Säumerstrasse 4 , 8803 Rüschlikon , Switzerland .
| |
Collapse
|
10
|
Mathea T, Rauhut G. Assignment of vibrational states within configuration interaction calculations. J Chem Phys 2020; 152:194112. [DOI: 10.1063/5.0009732] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Tina Mathea
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Guntram Rauhut
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| |
Collapse
|
11
|
Klinting EL, Lauvergnat D, Christiansen O. Vibrational Coupled Cluster Computations in Polyspherical Coordinates with the Exact Analytical Kinetic Energy Operator. J Chem Theory Comput 2020; 16:4505-4520. [DOI: 10.1021/acs.jctc.0c00261] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
| | - David Lauvergnat
- Université Paris-Saclay, CNRS, Institut de Chimie Physique, UMR8000, 91405 Orsay, France
| | - Ove Christiansen
- Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark
| |
Collapse
|
12
|
Klinting EL, Christiansen O, König C. Toward Accurate Theoretical Vibrational Spectra: A Case Study for Maleimide. J Phys Chem A 2020; 124:2616-2627. [DOI: 10.1021/acs.jpca.9b11915] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
| | - Ove Christiansen
- Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark
| | - Carolin König
- Institute of Physical Chemistry, Kiel University, Max-Eyth-Straße 1, D-24118 Kiel, Germany
| |
Collapse
|
13
|
Puzzarini C, Bloino J, Tasinato N, Barone V. Accuracy and Interpretability: The Devil and the Holy Grail. New Routes across Old Boundaries in Computational Spectroscopy. Chem Rev 2019; 119:8131-8191. [DOI: 10.1021/acs.chemrev.9b00007] [Citation(s) in RCA: 114] [Impact Index Per Article: 22.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Affiliation(s)
- Cristina Puzzarini
- Dipartimento di Chimica “Giacomo Ciamician”, Università di Bologna, Via F. Selmi 2, I-40126 Bologna, Italy
| | - Julien Bloino
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Nicola Tasinato
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| |
Collapse
|
14
|
Ziegler B, Rauhut G. Localized Normal Coordinates in Accurate Vibrational Structure Calculations: Benchmarks for Small Molecules. J Chem Theory Comput 2019; 15:4187-4196. [DOI: 10.1021/acs.jctc.9b00381] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Benjamin Ziegler
- 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
|
15
|
Sibert EL. Modeling vibrational anharmonicity in infrared spectra of high frequency vibrations of polyatomic molecules. J Chem Phys 2019; 150:090901. [DOI: 10.1063/1.5079626] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Edwin L. Sibert
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| |
Collapse
|
16
|
Panek PT, Hoeske AA, Jacob CR. On the choice of coordinates in anharmonic theoretical vibrational spectroscopy: Harmonic vs. anharmonic coupling in vibrational configuration interaction. J Chem Phys 2019; 150:054107. [PMID: 30736699 DOI: 10.1063/1.5083186] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023] Open
Abstract
By a suitable choice of coordinates, the computational effort required for calculations of anharmonic vibrational spectra can be reduced significantly. By using suitable localized-mode coordinates obtained from an orthogonal transformation of the conventionally used normal-mode coordinates, anharmonic couplings between modes can be significantly reduced. However, such a transformation introduces harmonic couplings between the localized modes. To elucidate the role of these harmonic couplings, we consider the vibrational self-consistent field (VSCF)/vibrational configuration interaction (VCI) calculations for both few-mode model systems and for ethene as a molecular test case. We show that large harmonic couplings can result in significant errors in localized-mode L-VSCF/L-VCI calculations and study the convergence with respect to the size of the VCI excitation space. To further elucidate the errors introduced by harmonic couplings, we discuss the connection between L-VSCF/L-VCI and vibrational exciton models. With the help of our results, we propose an algorithm for the localization of normal modes in suitable subsets that are chosen to strictly limit the errors introduced by the harmonic couplings while still leading to maximally localized modes.
Collapse
Affiliation(s)
- Paweł T Panek
- Technische Universität Braunschweig, Institute of Physical and Theoretical Chemistry, Gaußstraße 17, 38106 Braunschweig, Germany
| | - Adrian A Hoeske
- Technische Universität Braunschweig, Institute of Physical and Theoretical Chemistry, Gaußstraße 17, 38106 Braunschweig, Germany
| | - Christoph R Jacob
- Technische Universität Braunschweig, Institute of Physical and Theoretical Chemistry, Gaußstraße 17, 38106 Braunschweig, Germany
| |
Collapse
|
17
|
Ziegler B, Rauhut G. Rigorous use of symmetry within the construction of multidimensional potential energy surfaces. J Chem Phys 2018; 149:164110. [DOI: 10.1063/1.5047912] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023] Open
Affiliation(s)
- Benjamin Ziegler
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Guntram Rauhut
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| |
Collapse
|
18
|
Micciarelli M, Conte R, Suarez J, Ceotto M. Anharmonic vibrational eigenfunctions and infrared spectra from semiclassical molecular dynamics. J Chem Phys 2018; 149:064115. [DOI: 10.1063/1.5041911] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Marco Micciarelli
- Dipartimento di Chimica, Università degli Studi di Milano, via C. Golgi 19, 20133 Milano, Italy
| | - Riccardo Conte
- Dipartimento di Chimica, Università degli Studi di Milano, via C. Golgi 19, 20133 Milano, Italy
| | - Jaime Suarez
- Dipartimento di Chimica, Università degli Studi di Milano, via C. Golgi 19, 20133 Milano, Italy
| | - Michele Ceotto
- Dipartimento di Chimica, Università degli Studi di Milano, via C. Golgi 19, 20133 Milano, Italy
| |
Collapse
|
19
|
Morzan UN, Alonso de Armiño DJ, Foglia NO, Ramírez F, González Lebrero MC, Scherlis DA, Estrin DA. Spectroscopy in Complex Environments from QM–MM Simulations. Chem Rev 2018; 118:4071-4113. [DOI: 10.1021/acs.chemrev.8b00026] [Citation(s) in RCA: 97] [Impact Index Per Article: 16.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Affiliation(s)
- Uriel N. Morzan
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Diego J. Alonso de Armiño
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Nicolás O. Foglia
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Francisco Ramírez
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Mariano C. González Lebrero
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Damián A. Scherlis
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Darío A. Estrin
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| |
Collapse
|
20
|
Klinting EL, Thomsen B, Godtliebsen IH, Christiansen O. Employing general fit-bases for construction of potential energy surfaces with an adaptive density-guided approach. J Chem Phys 2018; 148:064113. [DOI: 10.1063/1.5016259] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
| | - Bo Thomsen
- Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark
| | | | - Ove Christiansen
- Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark
| |
Collapse
|
21
|
Hanson-Heine MWD. Reduced Basis Set Dependence in Anharmonic Frequency Calculations Involving Localized Coordinates. J Chem Theory Comput 2018; 14:1277-1285. [PMID: 29385338 DOI: 10.1021/acs.jctc.7b01075] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Localized normal coordinates are known to be effective in speeding up anharmonic frequency calculations by reducing the complexity of the nuclear Hamiltonian and wave function. Displacing atoms in localized coordinates can also cause relatively small changes in the electronic structure, which can be exploited for further computational efficiency improvements during ab initio electronic structure calculations of the potential energy surface by reducing the electronic basis set dependence. Three different schemes for reducing the basis set dependence have been investigated in this work. These include combining localized coordinate schemes with general mixed basis sets, distance based force-field reductions, and using coordinate specific basis sets. The importance of accurately describing electronic interactions is found to diminish both for multicoordinate terms involving the displacement of remote atoms and when describing the interactions between more remote atoms within specific coordinates.
Collapse
|
22
|
Bulik IW, Frisch MJ, Vaccaro PH. Vibrational self-consistent field theory using optimized curvilinear coordinates. J Chem Phys 2017; 147:044110. [DOI: 10.1063/1.4995440] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Ireneusz W. Bulik
- Department of Chemistry, Yale University, New Haven, Connecticut 06511, USA
| | - Michael J. Frisch
- Gaussian, Inc., 340 Quinnipiac St. Bldg. 40, Wallingford, Connecticut 06492, USA
| | - Patrick H. Vaccaro
- Department of Chemistry, Yale University, New Haven, Connecticut 06511, USA
| |
Collapse
|
23
|
Thomsen B, Kawakami T, Shigemoto I, Sugita Y, Yagi K. Weight-Averaged Anharmonic Vibrational Analysis of Hydration Structures of Polyamide 6. J Phys Chem B 2017; 121:6050-6063. [DOI: 10.1021/acs.jpcb.7b00372] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Bo Thomsen
- Theoretical
Molecular Science Laboratory, RIKEN, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | - Tomonori Kawakami
- Advanced
Materials Research Laboratories, Toray Industries, Inc., 2-1 Sonoyama 3-chome, Otsu, Shiga 520-0842, Japan
| | - Isamu Shigemoto
- Advanced
Materials Research Laboratories, Toray Industries, Inc., 2-1 Sonoyama 3-chome, Otsu, Shiga 520-0842, Japan
| | - Yuji Sugita
- Theoretical
Molecular Science Laboratory, RIKEN, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
- RIKEN iTHES, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
- RIKEN Advanced Institute for Computational Science, 7-1-26 Minatojima-Minamimachi,
Chuo-ku, Kobe, Hyogo 650-0047, Japan
- RIKEN Quantitative Biology Center, 6-7-1 Minatojima-Minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Kiyoshi Yagi
- Theoretical
Molecular Science Laboratory, RIKEN, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
- RIKEN iTHES, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| |
Collapse
|
24
|
Madsen NK, Godtliebsen IH, Christiansen O. Efficient algorithms for solving the non-linear vibrational coupled-cluster equations using full and decomposed tensors. J Chem Phys 2017; 146:134110. [DOI: 10.1063/1.4979498] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Niels K. Madsen
- Department of Chemistry, Aarhus University, 8000 Aarhus C, Denmark
| | | | - Ove Christiansen
- Department of Chemistry, Aarhus University, 8000 Aarhus C, Denmark
| |
Collapse
|
25
|
Yagi K, Thomsen B. Infrared Spectra of Protonated Water Clusters, H+(H2O)4, in Eigen and Zundel Forms Studied by Vibrational Quasi-Degenerate Perturbation Theory. J Phys Chem A 2017; 121:2386-2398. [DOI: 10.1021/acs.jpca.6b11189] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Kiyoshi Yagi
- Theoretical
Molecular Science Laboratory and ‡iTHES, RIKEN, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | - Bo Thomsen
- Theoretical
Molecular Science Laboratory and ‡iTHES, RIKEN, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| |
Collapse
|
26
|
Rai P, Sargsyan K, Najm H, Hermes MR, Hirata S. Low-rank canonical-tensor decomposition of potential energy surfaces: application to grid-based diagrammatic vibrational Green's function theory. Mol Phys 2017. [DOI: 10.1080/00268976.2017.1288937] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- Prashant Rai
- Sandia National Laboratories, Livermore, CA, USA
| | | | - Habib Najm
- Sandia National Laboratories, Livermore, CA, USA
| | - Matthew R. Hermes
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, IL, USA
| | - So Hirata
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, IL, USA
| |
Collapse
|
27
|
Hanson-Heine MWD. Intermediate vibrational coordinate localization with harmonic coupling constraints. J Chem Phys 2017; 144:204116. [PMID: 27250288 DOI: 10.1063/1.4951011] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Optimized normal coordinates can significantly improve the speed and accuracy of vibrational frequency calculations. However, over-localization can occur when using unconstrained spatial localization techniques. The unintuitive mixtures of stretching and bending coordinates that result can make interpreting spectra more difficult and also cause artificial increases in mode-coupling during anharmonic calculations. Combining spatial localization with a constraint on the coupling between modes can be used to generate coordinates with properties in-between the normal and fully localized schemes. These modes preserve the diagonal nature of the mass-weighted Hessian matrix to within a specified tolerance and are found to prevent contamination between the stretching and bending vibrations of the molecules studied without a priori classification of the different types of vibration present. Relaxing the constraint can also be used to identify which normal modes form specific groups of localized modes. The new coordinates are found to center on more spatially delocalized functional groups than their fully localized counterparts and can be used to tune the degree of vibrational correlation energy during anharmonic calculations.
Collapse
Affiliation(s)
- Magnus W D Hanson-Heine
- School of Chemistry, University of Nottingham, University Park, Nottingham NG7 2RD, United Kingdom
| |
Collapse
|
28
|
Panek PT, Jacob CR. On the benefits of localized modes in anharmonic vibrational calculations for small molecules. J Chem Phys 2017; 144:164111. [PMID: 27131535 DOI: 10.1063/1.4947213] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Anharmonic vibrational calculations can already be computationally demanding for relatively small molecules. The main bottlenecks lie in the construction of the potential energy surface and in the size of the excitation space in the vibrational configuration interaction (VCI) calculations. To address these challenges, we use localized-mode coordinates to construct potential energy surfaces and perform vibrational self-consistent field and L-VCI calculations [P. T. Panek and C. R. Jacob, ChemPhysChem 15, 3365 (2014)] for all vibrational modes of two prototypical test cases, the ethene and furan molecules. We find that the mutual coupling between modes is reduced when switching from normal-mode coordinates to localized-mode coordinates. When using such localized-mode coordinates, we observe a faster convergence of the n-mode expansion of the potential energy surface. This makes it possible to neglect higher-order contributions in the n-mode expansion of the potential energy surface or to approximate higher-order contributions in hybrid potential energy surfaces, which reduced the computational effort for the construction of the anharmonic potential energy surface significantly. Moreover, we find that when using localized-mode coordinates, the convergence with respect to the VCI excitation space proceeds more smoothly and that the error at low orders is reduced significantly. This makes it possible to devise low-cost models for obtaining a first approximation of anharmonic corrections. This demonstrates that the use of localized-mode coordinates can be beneficial already in anharmonic vibrational calculations of small molecules and provides a possible avenue for enabling such accurate calculations also for larger molecules.
Collapse
Affiliation(s)
- Paweł T Panek
- TU Braunschweig, Institute of Physical and Theoretical Chemistry, Hans-Sommer-Str. 10, 38106 Braunschweig, Germany
| | - Christoph R Jacob
- TU Braunschweig, Institute of Physical and Theoretical Chemistry, Hans-Sommer-Str. 10, 38106 Braunschweig, Germany
| |
Collapse
|
29
|
Battocchio G, Madsen NK, Christiansen O. Density matrices and iterative natural modals in vibrational structure theory. Mol Phys 2016. [DOI: 10.1080/00268976.2016.1243263] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
30
|
Monteiro JG, Barbosa AG. VSCF calculations for the intra- and intermolecular vibrational modes of the water dimer and its isotopologs. Chem Phys 2016. [DOI: 10.1016/j.chemphys.2016.09.027] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
31
|
Cheng X, Talbot JJ, Steele RP. Tuning vibrational mode localization with frequency windowing. J Chem Phys 2016; 145:124112. [DOI: 10.1063/1.4963109] [Citation(s) in RCA: 37] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023] Open
Affiliation(s)
- Xiaolu Cheng
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, Salt Lake City, Utah 84112, USA
| | - Justin J. Talbot
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, Salt Lake City, Utah 84112, USA
| | - Ryan P. Steele
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, Salt Lake City, Utah 84112, USA
| |
Collapse
|
32
|
Beć KB, Futami Y, Wójcik MJ, Nakajima T, Ozaki Y. Spectroscopic and Computational Study of Acetic Acid and Its Cyclic Dimer in the Near-Infrared Region. J Phys Chem A 2016; 120:6170-83. [PMID: 27482762 DOI: 10.1021/acs.jpca.6b04470] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Anharmonic vibrational analysis of near-infrared (NIR) spectra of acetic acid was carried out by anharmonic quantum chemical calculation in a wide concentration range of its CCl4 solution. By predicting vibrational spectra of acetic acid for the first time over a wide NIR region, it was possible to elucidate the influence of the formation of acetic acid cyclic dimer on its NIR spectrum. Quantum chemical simulations were based on coupled cluster and density functional theory quantum methods. Additionally, Møller-Plesset perturbation theory was employed for the additional calculation of hydrogen bonding stabilization energies. An anharmonic vibrational analysis was performed with the use of generalized second-order vibrational perturbation theory (GVPT2). A hybrid approach was assumed, in which monomeric species was treated by CCSD(T)/aug-cc-pVDZ (harmonic approximation) and B3LYP/SNSD (anharmonic approximation) methods. For the cyclic dimer, B3LYP and B2PLYP single and double hybrid functionals, paired with an SNSD basis set, were employed. DFT calculations were augmented with additional empirical dispersion correction. It was found that quantum chemically calculated vibrational modes in the NIR region are in a good agreement with experimental data. The results of anharmonic vibrational analysis were supported by a harmonic shift analysis, for elucidating the very strong anharmonic coupling observed between stretching modes of hydrogen bonded bridge in the cyclic dimer. However, the calculated wavenumbers for combination modes of double hydrogen bonded bridge in the cyclic dimer, which are very sensitive to the formation of hydrogen bonding, were found to be underestimated by quantum chemical methods. Therefore, by band fitting, the wavenumbers and shape parameters for these bands were found, and the modeled spectra were adjusted accordingly. A high accuracy of simulated spectra was achieved, and a detailed analysis of the experimental NIR spectra of acetic acid was possible, with successful identification of numerous experimental bands, including those which originate from concentration effects. It was also found that the main spectral features observed in the NIR spectra of carboxylic acid upon the formation of hydrogen bond should be accounted for combination modes of the stretching and bending vibrations of double hydrogen-bonded bridge in the cyclic dimers of acetic acid.
Collapse
Affiliation(s)
- Krzysztof B Beć
- Department of Chemistry, School of Science and Technology, Kwansei Gakuin University , Sanda, Hyogo 669-1337, Japan.,RIKEN , 519-1399 Aramaki-Aoba, Aoba-ku, Sendai, Miyagi 980-0845, Japan
| | - Yoshisuke Futami
- Department of Biological and Chemical Systems Engineering, National Institute of Technology, Kumamoto College , Yatsushiro, Kumamoto 866-8501, Japan
| | - Marek J Wójcik
- Faculty of Chemistry, Jagiellonian University , Ingardena 3, 30-060 Kraków, Poland
| | - Takahito Nakajima
- RIKEN Advanced Institute for Computational Science , 7-1-26, Minatojima-minami-machi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Yukihiro Ozaki
- Department of Chemistry, School of Science and Technology, Kwansei Gakuin University , Sanda, Hyogo 669-1337, Japan
| |
Collapse
|
33
|
Hanson-Heine MWD. Examining the impact of harmonic correlation on vibrational frequencies calculated in localized coordinates. J Chem Phys 2016; 143:164104. [PMID: 26520495 DOI: 10.1063/1.4934234] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Abstract
Carefully choosing a set of optimized coordinates for performing vibrational frequency calculations can significantly reduce the anharmonic correlation energy from the self-consistent field treatment of molecular vibrations. However, moving away from normal coordinates also introduces an additional source of correlation energy arising from mode-coupling at the harmonic level. The impact of this new component of the vibrational energy is examined for a range of molecules, and a method is proposed for correcting the resulting self-consistent field frequencies by adding the full coupling energy from connected pairs of harmonic and pseudoharmonic modes, termed vibrational self-consistent field (harmonic correlation). This approach is found to lift the vibrational degeneracies arising from coordinate optimization and provides better agreement with experimental and benchmark frequencies than uncorrected vibrational self-consistent field theory without relying on traditional correlated methods.
Collapse
Affiliation(s)
- Magnus W D Hanson-Heine
- School of Chemistry, University of Nottingham, University Park, Nottingham NG7 2RD, United Kingdom
| |
Collapse
|
34
|
Zimmerman PM, Smereka P. Optimizing Vibrational Coordinates To Modulate Intermode Coupling. J Chem Theory Comput 2016; 12:1883-91. [PMID: 26914536 DOI: 10.1021/acs.jctc.5b01168] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The choice of coordinate system strongly affects the convergence properties of vibrational structure computations. Two methods for efficient generation of improved vibrational coordinates are presented and justified by analysis of a model anharmonic two-mode Hessian and numerical computations on polyatomic molecules. To produce optimal coordinates, metrics which quantify off-diagonal couplings over a grid of Hessian matrices are minimized through unitary rotations of the vibrational basis. The first proposed metric minimizes the total squared off-diagonal coupling, and the second minimizes the total squared change in off-diagonal coupling. In this procedure certain anharmonic modes tend to localize, for example X-H stretches. The proposed methods do not rely on prior fitting of the potential energy, vibrational structure computations, or localization metrics, so they are unique from previous vibrational coordinate generation algorithms and are generally applicable to polyatomic molecules. Fitting the potential to the approximate n-mode representation in the optimized bases for all-trans polyenes shows that off-diagonal anharmonic couplings are substantially reduced by the new choices of coordinate system. Convergence of vibrational energies is examined in detail for ethylene, and it is shown that coupling-optimized modes converge in vibrational configuration interaction computations to within 1 cm(-1) using only 3-mode couplings, where normal modes require 4-mode couplings for convergence. Comparison of the vibrational configuration interaction convergence with respect to excitation level for the two proposed metrics shows that minimization of the total off-diagonal coupling is most effective for low-cost vibrational structure computations.
Collapse
Affiliation(s)
- Paul M Zimmerman
- Department of Chemistry and ‡Department of Mathematics, University of Michigan , Ann Arbor, Michigan 48109, United States
| | - Peter Smereka
- Department of Chemistry and ‡Department of Mathematics, University of Michigan , Ann Arbor, Michigan 48109, United States
| |
Collapse
|
35
|
König C, Hansen MB, Godtliebsen IH, Christiansen O. FALCON: A method for flexible adaptation of local coordinates of nuclei. J Chem Phys 2016; 144:074108. [DOI: 10.1063/1.4941846] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Carolin König
- Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark
| | | | | | - Ove Christiansen
- Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark
| |
Collapse
|
36
|
Yagi K. Development of Molecular Vibrational Structure Theory with an Explicit Account of Anharmonicity. ACTA ACUST UNITED AC 2016. [DOI: 10.3175/molsci.10.a0085] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
37
|
Klinting EL, König C, Christiansen O. Hybrid Optimized and Localized Vibrational Coordinates. J Phys Chem A 2015; 119:11007-21. [DOI: 10.1021/acs.jpca.5b08496] [Citation(s) in RCA: 32] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
| | - Carolin König
- Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark
| | - Ove Christiansen
- Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark
| |
Collapse
|
38
|
Faucheaux JA, Hirata S. Higher-order diagrammatic vibrational coupled-cluster theory. J Chem Phys 2015; 143:134105. [DOI: 10.1063/1.4931472] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Jacob A. Faucheaux
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, USA
| | - So Hirata
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, USA
- CREST, Japan Science and Technology Agency, 4-1-8 Honcho, Kawaguchi, Saitama 332-0012, Japan
| |
Collapse
|
39
|
Sibaev M, Crittenden DL. The PyPES library of high quality semi-global potential energy surfaces. J Comput Chem 2015; 36:2200-7. [DOI: 10.1002/jcc.24192] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2015] [Revised: 07/19/2015] [Accepted: 08/10/2015] [Indexed: 11/07/2022]
Affiliation(s)
- Marat Sibaev
- Department of Chemistry; University of Canterbury; Christchurch New Zealand
| | | |
Collapse
|
40
|
Czarnecki MA, Morisawa Y, Futami Y, Ozaki Y. Advances in Molecular Structure and Interaction Studies Using Near-Infrared Spectroscopy. Chem Rev 2015; 115:9707-44. [DOI: 10.1021/cr500013u] [Citation(s) in RCA: 146] [Impact Index Per Article: 16.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
| | - Yusuke Morisawa
- Department
of Chemistry, School of Science and Engineering, Kinki University, Higashi-Osaka, Osaka 577-8502, Japan
| | - Yoshisuke Futami
- Department
of Biological and Chemical Systems Engineering, National Institute of Technology, Kumamoto College, Yatsushiro, Kumamoto 866-8501, Japan
| | - Yukihiro Ozaki
- Department
of Chemistry, School of Science and Technology, Kwansei Gakuin University, Sanda, Hyogo 669-1337, Japan
| |
Collapse
|
41
|
Cheng X, Steele RP. Efficient anharmonic vibrational spectroscopy for large molecules using local-mode coordinates. J Chem Phys 2015; 141:104105. [PMID: 25217902 DOI: 10.1063/1.4894507] [Citation(s) in RCA: 84] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023] Open
Abstract
This article presents a general computational approach for efficient simulations of anharmonic vibrational spectra in chemical systems. An automated local-mode vibrational approach is presented, which borrows techniques from localized molecular orbitals in electronic structure theory. This approach generates spatially localized vibrational modes, in contrast to the delocalization exhibited by canonical normal modes. The method is rigorously tested across a series of chemical systems, ranging from small molecules to large water clusters and a protonated dipeptide. It is interfaced with exact, grid-based approaches, as well as vibrational self-consistent field methods. Most significantly, this new set of reference coordinates exhibits a well-behaved spatial decay of mode couplings, which allows for a systematic, a priori truncation of mode couplings and increased computational efficiency. Convergence can typically be reached by including modes within only about 4 Å. The local nature of this truncation suggests particular promise for the ab initio simulation of anharmonic vibrational motion in large systems, where connection to experimental spectra is currently most challenging.
Collapse
Affiliation(s)
- Xiaolu Cheng
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, Salt Lake City, Utah 84112, USA
| | - Ryan P Steele
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, Salt Lake City, Utah 84112, USA
| |
Collapse
|
42
|
Ramakrishnan R, Rauhut G. Semi-quartic force fields retrieved from multi-mode expansions: Accuracy, scaling behavior, and approximations. J Chem Phys 2015; 142:154118. [DOI: 10.1063/1.4918587] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Raghunathan Ramakrishnan
- Institute of Physical Chemistry and National Center for Computational Design and Discovery of Novel Materials (MARVEL), Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Guntram Rauhut
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| |
Collapse
|
43
|
König C, Christiansen O. Automatic determination of important mode–mode correlations in many-mode vibrational wave functions. J Chem Phys 2015; 142:144115. [DOI: 10.1063/1.4916518] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Carolin König
- Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark
| | - Ove Christiansen
- Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark
| |
Collapse
|
44
|
The performance of RI-MP2 based potential energy surfaces in a vibrational self-consistent field treatment. Chem Phys Lett 2015. [DOI: 10.1016/j.cplett.2014.11.050] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
45
|
Largely reduced grid densities in a vibrational self-consistent field treatment do not significantly impact the resultingwavenumbers. Molecules 2014; 19:21253-75. [PMID: 25525825 PMCID: PMC6270979 DOI: 10.3390/molecules191221253] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2014] [Revised: 11/25/2014] [Accepted: 12/09/2014] [Indexed: 11/16/2022] Open
Abstract
Especially for larger molecules relevant to life sciences, vibrational self-consistent field (VSCF) calculations can become unmanageably demanding even when only first and second order potential coupling terms are considered. This paper investigates to what extent the grid density of the VSCF’s underlying potential energy surface can be reduced without sacrificing accuracy of the resulting wavenumbers. Including single-mode and pair contributions, a reduction to eight points per mode did not introduce a significant deviation but improved the computational efficiency by a factor of four. A mean unsigned deviation of 1.3% from the experiment could be maintained for the fifteen molecules under investigation and the approach was found to be applicable to rigid, semi-rigid and soft vibrational problems likewise. Deprotonated phosphoserine, stabilized by two intramolecular hydrogen bonds, was investigated as an exemplary application.
Collapse
|
46
|
Thomsen B, Yagi K, Christiansen O. A simple state-average procedure determining optimal coordinates for anharmonic vibrational calculations. Chem Phys Lett 2014. [DOI: 10.1016/j.cplett.2014.07.043] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
47
|
Panek PT, Jacob CR. Efficient calculation of anharmonic vibrational spectra of large molecules with localized modes. Chemphyschem 2014; 15:3365-77. [PMID: 25080397 DOI: 10.1002/cphc.201402251] [Citation(s) in RCA: 54] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2014] [Indexed: 11/11/2022]
Abstract
The analysis and interpretation of the vibrational spectra of complex (bio)molecular systems, such as polypeptides and proteins, requires support from quantum-chemical calculations. Such calculations are currently restricted to the harmonic approximation. Here, we show how one of the main bottlenecks in such calculations, the evaluation of the potential energy surface, can be overcome by using localized modes instead of the commonly employed normal modes. We apply such local vibrational self-consistent field (L-VSCF) and vibrational configuration interaction (L-VCI) calculations to a cyclic water tetramer and a helical hexa-alanine peptide. The results show that the use of localized modes is equivalent to the commonly used normal modes, but offers several advantages. First, a faster convergence with respect to the excitation level is observed in L-VCI calculations. Second, the localized modes provide a reduced representation of the couplings between modes that show a regular coupling pattern. This can be used to disregard a significant number of small two-mode potentials a priori. Several such reduced coupling approximations are explored, and we show that the number of single-point calculations required to evaluate the potential energy surface can be significantly reduced without introducing noticeable errors in the resulting vibrational spectra.
Collapse
Affiliation(s)
- Paweł T Panek
- Karlsruhe Institute of Technology (KIT), Center for Functional Nanostructures and Institute of Physical Chemistry, Wolfgang-Gaede-Str. 1a, 76131 Karlsruhe (Germany) http://www.christophjacob.eu
| | | |
Collapse
|