Variational computation of anharmonic ground and excited vibrational eigenstates using bound Quartic Force Fields: Application with MCTDH and ElVibRot
Abstract
In this work we introduce the use of Quartic Force fields (QFF) potential expansions in the context of variational calculations. Such potentials are commonly employed in molecular Vibrational Second-Order Perturbation Theory (VPT2) studies, for which equations explicitly dependent on the QFF parameters exist. However, QFF are unbound potentials for more or less large displacements from the reference point and, most of the time, this prevents their use in conjunction with variational wavepacket-based calculations. In this work, we propose a general correction to QFFs and introduce a fully automated numerical approach to avoid their unbound character. Our corrected potentials, bound QFF (bQFF), do not exhibit appreciable modification of the local topography around the region of interest for infrared spectroscopy. As a consequence of this, we can affirm that the vibrational eigenstructure (eigenvalues, eigenstates) remains essentially unaltered by our correction. To illustrate their numerical stability, we have interfaced our bQFF routines in combination with to two well-established quantum simulation software packages MCTDH and \textsc{ElVibRot} which feature variational approaches. More specifically, our bQFFs are separable and hence directly expressible as MCTDH operators. Furthermore, concerning the size of our bQFF expansion, we show that it is possible to tensor-decompose our bQFF in Canonical Polyadic form (CP-bQFF). We use the Monte Carlo Canonical Polyadic decomposition algorithm for this. CP-bQFF results are virtually identical to uncompressed bQFF, but the computational efficiency is largely improved. Our approach paves the way for the automated variational study of anharmonic eigenstates in molecular systems within the reach of QFF-based potentials, using either time-dependent or time-independent schemes.
Related
- Anharmonic Vibrational Structure and Bound States of the He-CO2 van der Waals Complex
- An anharmonic liquid-entropy functional from the Mori-Zwanzig memory kernel
- AnhDis: A Scheme for Anharmonic Nuclear Ensembles and Absorption Spectra
- Effective Harmonic and Anharmonic Models for Vibronic Spectra
- Cartesian and Curvilinear Anharmonic Frequencies and Intensities: Perturb-then-Diagonalize Polyad Transport and Adiabatic Discrete-Variable Representation
BibTeX
@misc{aljaaidi2026variational,
title = {{Variational computation of anharmonic ground and excited vibrational eigenstates using bound Quartic Force Fields: Application with MCTDH and ElVibRot}},
author = {Ali Al-Jaaidi and David Lauvergnat and Daniel Peláez},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.18654},
eprint = {2609.18654},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.18654},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Variational computation of anharmonic ground and excited vibrational eigenstates using bound Quartic Force Fields: Application with MCTDH and ElVibRot
AU - Ali Al-Jaaidi
AU - David Lauvergnat
AU - Daniel Peláez
PY - 2026
DA - 2026/09/29
PB - arXiv
DO - 10.48550/arxiv.2609.18654
UR - https://arxiv.org/abs/2609.18654
AB - In this work we introduce the use of Quartic Force fields (QFF) potential expansions in the context of variational calculations. Such potentials are commonly employed in molecular Vibrational Second-Order Perturbation Theory (VPT2) studies, for which equations explicitly dependent on the QFF parameters exist. However, QFF are unbound potentials for more or less large displacements from the reference point and, most of the time, this prevents their use in conjunction with variational wavepacket-based calculations. In this work, we propose a general correction to QFFs and introduce a fully automated numerical approach to avoid their unbound character. Our corrected potentials, bound QFF (bQFF), do not exhibit appreciable modification of the local topography around the region of interest for infrared spectroscopy. As a consequence of this, we can affirm that the vibrational eigenstructure (eigenvalues, eigenstates) remains essentially unaltered by our correction. To illustrate their numerical stability, we have interfaced our bQFF routines in combination with to two well-established quantum simulation software packages MCTDH and \textsc{ElVibRot} which feature variational approaches. More specifically, our bQFFs are separable and hence directly expressible as MCTDH operators. Furthermore, concerning the size of our bQFF expansion, we show that it is possible to tensor-decompose our bQFF in Canonical Polyadic form (CP-bQFF). We use the Monte Carlo Canonical Polyadic decomposition algorithm for this. CP-bQFF results are virtually identical to uncompressed bQFF, but the computational efficiency is largely improved. Our approach paves the way for the automated variational study of anharmonic eigenstates in molecular systems within the reach of QFF-based potentials, using either time-dependent or time-independent schemes.
ER -