Juxtaposing the fourth order vibrational operator perturbation theory CVPT(4) and the adaptive VCI (A-VCI): Accuracy, vibrational resonances and polyads of C2H4 and C2D4
Vincent Le Bris , Sergey V. Krasnoshchekov , Egor O. Dobrolyubov , Ilya M. Efremov , Igor V. Polyakov , Olivier Coulaud , Didier Bégué
{"title":"Juxtaposing the fourth order vibrational operator perturbation theory CVPT(4) and the adaptive VCI (A-VCI): Accuracy, vibrational resonances and polyads of C2H4 and C2D4","authors":"Vincent Le Bris , Sergey V. Krasnoshchekov , Egor O. Dobrolyubov , Ilya M. Efremov , Igor V. Polyakov , Olivier Coulaud , Didier Bégué","doi":"10.1016/j.jqsrt.2025.109586","DOIUrl":null,"url":null,"abstract":"<div><div><em>Ab initio</em> prediction of anharmonic vibrational spectra produces an increasing computational overhead for larger molecules, requesting a balance between an accuracy and resources. Two complementary fundamental quantum mechanical approaches, the perturbative and variational, have various strong and weak features, depending on a specific target problem. The vibrational perturbation theory (VPT) treats weak couplings and strong resonances separately, relying on somewhat artificial criteria. In contrast, the more precise but computationally intense variational configuration interaction (VCI) method treats all couplings in universal manner. The active ongoing development of approaches to solving vibrational problems requires an update of comparative benchmarks, helping to choose the best theoretical tools for a particular target. In this work, the performance of two particular modern implementations of these methods was juxtaposed: the second and fourth order operator canonical perturbation theory CVPT(2,4) and a recently proposed adaptive vibrational configuration interaction method (A-VCI). Two practically important C<sub>2</sub>H<sub>4</sub> and C<sub>2</sub>D<sub>4</sub> molecules and an accurate CCSD(T)/cc-pVQZ four-body sextic normal mode PES were employed for benchmarking. The comprehensive picture of vibrational resonances and the polyad quantum number was revealed in the form: <span><math><mrow><mi>P</mi><mrow><mo>[</mo><msub><mrow><mi>C</mi></mrow><mrow><mn>2</mn></mrow></msub><msub><mrow><mi>H</mi></mrow><mrow><mn>4</mn></mrow></msub><mo>]</mo></mrow><mo>=</mo><mn>4</mn><mrow><mo>(</mo><msub><mrow><mi>ν</mi></mrow><mrow><mn>1</mn></mrow></msub><mo>+</mo><msub><mrow><mi>ν</mi></mrow><mrow><mn>5</mn></mrow></msub><mo>+</mo><msub><mrow><mi>ν</mi></mrow><mrow><mn>9</mn></mrow></msub><mo>+</mo><msub><mrow><mi>ν</mi></mrow><mrow><mn>11</mn></mrow></msub><mo>)</mo></mrow><mo>+</mo><mn>2</mn><mrow><mo>(</mo><msub><mrow><mi>ν</mi></mrow><mrow><mn>2</mn></mrow></msub><mo>+</mo><msub><mrow><mi>ν</mi></mrow><mrow><mn>3</mn></mrow></msub><mo>+</mo><msub><mrow><mi>ν</mi></mrow><mrow><mn>6</mn></mrow></msub><mo>+</mo><msub><mrow><mi>ν</mi></mrow><mrow><mn>12</mn></mrow></msub><mo>)</mo></mrow><mo>+</mo><msub><mrow><mi>ν</mi></mrow><mrow><mn>4</mn></mrow></msub><mo>+</mo><msub><mrow><mi>ν</mi></mrow><mrow><mn>7</mn></mrow></msub><mo>+</mo><msub><mrow><mi>ν</mi></mrow><mrow><mn>8</mn></mrow></msub><mo>+</mo><msub><mrow><mi>ν</mi></mrow><mrow><mn>10</mn></mrow></msub></mrow></math></span>. A new quadratic resonance criterion <span><math><mrow><msub><mrow><mi>X</mi></mrow><mrow><mi>j</mi></mrow></msub><mo>=</mo><mrow><mo>|</mo></mrow><msubsup><mrow><mi>V</mi></mrow><mrow><mi>j</mi></mrow><mrow><mn>2</mn></mrow></msubsup><mo>/</mo><mrow><mo>(</mo><mrow><msub><mrow><mo>∑</mo></mrow><mrow><mi>k</mi></mrow></msub><mrow><mo>(</mo><msub><mrow><mi>n</mi></mrow><mrow><mi>k</mi></mrow></msub><mo>−</mo><msub><mrow><mi>m</mi></mrow><mrow><mi>k</mi></mrow></msub><mo>)</mo></mrow><msub><mrow><mi>ω</mi></mrow><mrow><mi>k</mi></mrow></msub></mrow><mo>)</mo></mrow><mrow><mo>|</mo></mrow></mrow></math></span> is proposed and its efficiency in elucidating polyad structures is demonstrated. A striking observation was made that CVPT(2) often produces better predictions of fundamental frequencies (e.g., RMSD for C<sub>2</sub>H<sub>4</sub> fundamentals is 2.7 cm<sup>−1</sup>, and for C<sub>2</sub>D<sub>4</sub> is 1.9 cm<sup>−1</sup>), while CVPT(4) demonstrates an excellent level of correlation with A-VCI results for both fundamentals and two-quanta states.</div></div>","PeriodicalId":16935,"journal":{"name":"Journal of Quantitative Spectroscopy & Radiative Transfer","volume":"346 ","pages":"Article 109586"},"PeriodicalIF":1.9000,"publicationDate":"2025-07-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Journal of Quantitative Spectroscopy & Radiative Transfer","FirstCategoryId":"101","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0022407325002481","RegionNum":3,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q2","JCRName":"OPTICS","Score":null,"Total":0}
引用次数: 0
Abstract
Ab initio prediction of anharmonic vibrational spectra produces an increasing computational overhead for larger molecules, requesting a balance between an accuracy and resources. Two complementary fundamental quantum mechanical approaches, the perturbative and variational, have various strong and weak features, depending on a specific target problem. The vibrational perturbation theory (VPT) treats weak couplings and strong resonances separately, relying on somewhat artificial criteria. In contrast, the more precise but computationally intense variational configuration interaction (VCI) method treats all couplings in universal manner. The active ongoing development of approaches to solving vibrational problems requires an update of comparative benchmarks, helping to choose the best theoretical tools for a particular target. In this work, the performance of two particular modern implementations of these methods was juxtaposed: the second and fourth order operator canonical perturbation theory CVPT(2,4) and a recently proposed adaptive vibrational configuration interaction method (A-VCI). Two practically important C2H4 and C2D4 molecules and an accurate CCSD(T)/cc-pVQZ four-body sextic normal mode PES were employed for benchmarking. The comprehensive picture of vibrational resonances and the polyad quantum number was revealed in the form: . A new quadratic resonance criterion is proposed and its efficiency in elucidating polyad structures is demonstrated. A striking observation was made that CVPT(2) often produces better predictions of fundamental frequencies (e.g., RMSD for C2H4 fundamentals is 2.7 cm−1, and for C2D4 is 1.9 cm−1), while CVPT(4) demonstrates an excellent level of correlation with A-VCI results for both fundamentals and two-quanta states.
期刊介绍:
Papers with the following subject areas are suitable for publication in the Journal of Quantitative Spectroscopy and Radiative Transfer:
- Theoretical and experimental aspects of the spectra of atoms, molecules, ions, and plasmas.
- Spectral lineshape studies including models and computational algorithms.
- Atmospheric spectroscopy.
- Theoretical and experimental aspects of light scattering.
- Application of light scattering in particle characterization and remote sensing.
- Application of light scattering in biological sciences and medicine.
- Radiative transfer in absorbing, emitting, and scattering media.
- Radiative transfer in stochastic media.