**Êtes-vous un étudiant de l'EPFL à la recherche d'un projet de semestre?**

Travaillez avec nous sur des projets en science des données et en visualisation, et déployez votre projet sous forme d'application sur GraphSearch.

Publication# High-order geometric integrators for the variational Gaussian wavepacket dynamics and application to vibronic spectra at finite temperature

Résumé

Molecular quantum dynamics simulations are essential for understanding many fundamental phenomena in physics and chemistry. They often require solving the time-dependent Schrödinger equation for molecular nuclei, which is challenging even for medium-sized systems. Therefore, for the majority of molecules, it is necessary to use approximate methods to evaluate the nuclear quantum dynamics.In this thesis, we comprehensively investigate the semiclassical variational thawed Gaussian wavepacket dynamics. This method is the most accurate of the single-trajectory Gaussian-based techniques. Moreover, it is symplectic, time-reversible, and norm- and energy-conserving. In one-dimensional double-well potentials, we numerically demonstrate that the variational method may include quantum tunneling approximately. To model anharmonic multidimensional systems, we develop a non-separable coupled Morse potential for which the expectation value in a Gaussian wavepacket can be evaluated analytically. To show that the variational method, unlike the exact quantum grid-based techniques, can simulate high-dimensional systems, we do most of our analysis on a twenty-dimensional coupled Morse system.The equations of motion in the variational method require expectation values of the potential energy, its gradient, and its Hessian, which can be prohibitively expensive to be evaluated for real molecular systems. One way to overcome this problem is to Taylor expand the potential energy around the center of the Gaussian wavepacket. Using the local harmonic approximation of the potential in the variational method yields Heller's thawed Gaussian approximation, which conserves neither the symplectic structure nor the effective energy. Therefore, we employ the local cubic approximation of the potential. The local cubic variational thawed Gaussian wavepacket dynamics is symplectic, conserves the effective energy, and improves the accuracy over the local harmonic method.Although they are more accurate than Heller's thawed Gaussian approximation, both the fully variational and local cubic variational methods are also much more expensive. In order to reduce the computational cost, the high-order geometric integrators are extensively investigated in this thesis. We implement the second-order geometric integrator developed by Faou and Lubich and apply various composition methods to obtain high-order geometric integrators. The high-order integrators can be significantly more efficient than the second-order algorithm, while at the same time being much more accurate. In addition, the geometric integrators preserve the geometric properties of the variational methods, such as norm conservation, symplecticity, time reversibility, and, for sufficiently small time steps, energy or effective energy conservation.To include non-zero temperature effects in the simulation of vibronic spectra, we use the concept of thermo-field dynamics. By computing the vibronic spectra of a harmonic potential at different temperatures, we demonstrate the validity of the thermo-field thawed Gaussian wavepacket dynamics. Furthermore, for the calculation of vibronic spectra of systems with low anharmonicity at finite temperatures, we show the superior accuracy of the fully variational and local cubic variational thermo-field thawed Gaussian wavepacket dynamics over the local harmonic thermo-field method.

Source officielle

Cette page est générée automatiquement et peut contenir des informations qui ne sont pas correctes, complètes, à jour ou pertinentes par rapport à votre recherche. Il en va de même pour toutes les autres pages de ce site. Veillez à vérifier les informations auprès des sources officielles de l'EPFL.

MOOCs associés (32)

Concepts associés (46)

Publications associées (52)

Plasma Physics: Introduction

Learn the basics of plasma, one of the fundamental states of matter, and the different types of models used to describe it, including fluid and kinetic.

Plasma Physics: Introduction

Learn the basics of plasma, one of the fundamental states of matter, and the different types of models used to describe it, including fluid and kinetic.

Numerical methods for ordinary differential equations

Numerical methods for ordinary differential equations are methods used to find numerical approximations to the solutions of ordinary differential equations (ODEs). Their use is also known as "numerical integration", although this term can also refer to the computation of integrals. Many differential equations cannot be solved exactly. For practical purposes, however – such as in engineering – a numeric approximation to the solution is often sufficient. The algorithms studied here can be used to compute such an approximation.

Calcul numérique d'une intégrale

En analyse numérique, il existe une vaste famille d’algorithmes dont le but principal est d’estimer la valeur numérique de l’intégrale définie sur un domaine particulier pour une fonction donnée (par exemple l’intégrale d’une fonction d’une variable sur un intervalle). Ces techniques procèdent en trois phases distinctes : Décomposition du domaine en morceaux (un intervalle en sous-intervalles contigus) ; Intégration approchée de la fonction sur chaque morceau ; Sommation des résultats numériques ainsi obtenus.

Time reversibility

A mathematical or physical process is time-reversible if the dynamics of the process remain well-defined when the sequence of time-states is reversed. A deterministic process is time-reversible if the time-reversed process satisfies the same dynamic equations as the original process; in other words, the equations are invariant or symmetrical under a change in the sign of time. A stochastic process is reversible if the statistical properties of the process are the same as the statistical properties for time-reversed data from the same process.

Jiri Vanicek, Roya Moghaddasi Fereidani

Among the single-trajectory Gaussian-based methods for solving the time-dependent Schrödinger equation, the variational Gaussian approximation is the most accurate one. In contrast to Heller’s original thawed Gaussian approximation, it is symplectic, conse ...

2023Accurate simulations of molecular quantum dynamics are crucial for understanding numerous natural processes and experimental results. Yet, such high-accuracy simulations are challenging even for relatively simple systems where the Born-Oppenheimer approxim ...

The explicit split-operator algorithm has been extensively used for solving not only linear but also nonlinear time-dependent Schrödinger equations. When applied to the nonlinear Gross–Pitaevskii equation, the method remains time-reversible, norm-conservin ...

2021