Publication

Which form of the molecular Hamiltonian is the most suitable for simulating the nonadiabatic quantum dynamics at a conical intersection?

Related publications (28)

From low-rank retractions to dynamical low-rank approximation and back

Daniel Kressner, Axel Elie Joseph Séguin, Gianluca Ceruti

In algorithms for solving optimization problems constrained to a smooth manifold, retractions are a well-established tool to ensure that the iterates stay on the manifold. More recently, it has been demonstrated that retractions are a useful concept for ot ...

Springer2024

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

Roya Moghaddasi Fereidani

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 ...

EPFL2024

Applications of the thawed Gaussian approximation to electronic spectroscopy

Eriks Kletnieks

The exploration of electronically excited states and the study of diverse photochemical and photophysical processes are the main goals of molecular electronic spectroscopy. Exact quantum-mechanical simulation of such experiments is, however, beyond current ...

EPFL2024

High-order geometric integrators for the variational Gaussian approximation

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 ...

2023

TimeEvolver: A program for time evolution with improved error bound

Sebastian Zell

We present TimeEvolver, a program for computing time evolution in a generic quantum system. It relies on well-known Krylov subspace techniques to tackle the problem of multiplying the exponential of a large sparse matrix iH, where His the Hamiltonian, with ...

ELSEVIER2022

High-order geometric integrators for nonadiabatic molecular quantum dynamics and their applications to explore conical intersections

Seonghoon Choi

Accurate 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 ...

EPFL2022

Structure-preserving approaches and data-driven closure modeling for model order reduction

Nicolò Ripamonti

In this thesis, we propose model order reduction techniques for high-dimensional PDEs that preserve structures of the original problems and develop a closure modeling framework leveraging the Mori-Zwanzig formalism and recurrent neural networks. Since high ...

EPFL2022

Rank-adaptive structure-preserving model order reduction of Hamiltonian systems

Cecilia Pagliantini, Jan Sickmann Hesthaven, Nicolò Ripamonti

This work proposes an adaptive structure-preserving model order reduction method for finite-dimensional parametrized Hamiltonian systems modeling non-dissipative phenomena. To overcome the slowly decaying Kolmogorov width typical of transport problems, the ...

EDP SCIENCES S A2022

Geometrical Treatise on the Modelling of 3D Particulate Inclusion-Matrix Microstructures with an Application to Historical Stone Masonry Walls

Mahmoud S. M. Shaqfa

As historical stone masonry structures are vulnerable and prone to damage in earthquakes, investigating their structural integrity is important to reduce injuries and casualties while preserving their historical value. Stone masonry is a composite material ...

EPFL2022

High-order geometric integrators for representation-free Ehrenfest dynamics

Jiri Vanicek, Seonghoon Choi

Ehrenfest dynamics is a useful approximation for ab initio mixed quantum-classical molecular dynamics that can treat electronically nonadiabatic effects. Although a severe approximation to the exact solution of the molecular time-dependent Schrödinger equa ...

2021

How important are the residual nonadiabatic couplings for an accurate simulation of nonadiabatic quantum dynamics in a quasidiabatic representation?

Jiri Vanicek, Seonghoon Choi

Diabatization of the molecular Hamiltonian is a standard approach to remove the singularities of nonadiabatic couplings at conical intersections of adiabatic potential energy surfaces. In general, it is impossible to eliminate the nonadiabatic couplings en ...

2021

Probabilistic and Bayesian methods for uncertainty quantification of deterministic and stochastic differential equations

Giacomo Garegnani

In this thesis we explore uncertainty quantification of forward and inverse problems involving differential equations. Differential equations are widely employed for modeling natural and social phenomena, with applications in engineering, chemistry, meteor ...

EPFL2021

Time-reversible and norm-conserving high-order integrators for the nonlinear time-dependent Schrödinger equation: Application to local control theory

Jiri Vanicek, Julien Roulet

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

Wavelet-Fourier CORSING techniques for multidimensional advection-diffusion-reaction equations

Fabio Nobile, Simone Brugiapaglia

We present and analyze a novel wavelet-Fourier technique for the numerical treatment of multidimensional advection–diffusion–reaction equations based on the COmpRessed SolvING (CORSING) paradigm. Combining the Petrov–Galerkin technique with the compressed ...

2021

Symplectic dynamical low rank approximation of wave equations with random parameters

Fabio Nobile, Eleonora Musharbash, Eva Vidlicková

In this paper we propose a dynamical low-rank strategy for the approximation of second order wave equations with random parameters. The governing equation is rewritten in Hamiltonian form and the approximate solution is expanded over a set of 2S dynamical ...

2020

Rapid 3-D Magnetic Integral Field Computation of Current-Carrying Finite Arc Segments With Rectangular Cross Section

Basile Kawkabani, Frédéric Michael Maurer

The computation of 3-D magnetic fields is a demanding task in the analysis of electrical machines and other electromagnetic devices. In this context, integral field calculation provides a smooth solution, high precision and resolution, "on-demand"-calculat ...

Institute of Electrical and Electronics Engineers2020

Functional Inverse Problems on Spheres: Theory, Algorithms and Applications

Matthieu Martin Jean-André Simeoni

Many scientific inquiries in natural sciences involve approximating a spherical field -namely a scalar quantity defined over a continuum of directions- from generalised samples of the latter (e.g. directional samples, local averages, etc). Such an approxim ...

EPFL2020

Dynamics of optical forces and torques in plasmonic systems: a surface integral equation

Olivier Martin, Karim Achouri, Andrei Kiselev

A novel approach is introduced to determine the time evolution of optical forces and torques on arbitrary shape nanostructures by combining Maxwell's stress tensor with the surface integral equation method (SIE). Conventional time averaging of Maxwell’s st ...

SPIEDigitalLibrary2020

Coordinate parameterisation and spectral method optimisation for Beltrami field solver in stellarator geometry

David Pfefferlé, Aravind Kumar Chandiran, Antoine Baillod, Ajay Kumar

The numerical solution of the stepped pressure equilibrium (Hudson et al 2012 Phys. Plasmas 19 112502) requires a fast and robust solver to obtain the Beltrami field in three-dimensional geometry such as stellarators. The spectral method implemented in the ...

2020

Efficient geometric integrators for nonadiabatic quantum dynamics. I. The adiabatic representation

Jiri Vanicek, Seonghoon Choi

Geometric integrators of the Schrödinger equation conserve exactly many invariants of the exact solution. Among these integrators, the split-operator algorithm is explicit and easy to implement but, unfortunately, is restricted to systems whose Hamiltonian ...

2019