Source author record

Caroline Lasser

Caroline Lasser appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

14works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2022arXiv

An Error Bound for the Time-Sliced Thawed Gaussian Propagation Method

We study the time-sliced thawed Gaussian propagation method, which was recently proposed for solving the time-dependent Schrödinger equation. We introduce a triplet of quadrature-based analysis, synthesis and re-initialization operators to give a rigorous mathematical formulation of the method. Further, we derive combined error bounds for the discretization of the wave packet transform and the time-propagation of the thawed Gaussian basis functions. Numerical experiments in 1D illustrate the theoretical results.

preprint2022arXiv

Laser-induced dynamic alignment of the HD molecule without the Born-Oppenheimer approximation

Laser-induced molecular alignment is well understood within the framework of the Born-Oppenheimer (BO) approximation Without the BO approximation, however, the concept of molecular structure is lost, making alignment hard to define precisely. In this work, we demonstrate the emergence of alignment from the first-ever non-BO quantum dynamics simulations, using the HD molecule exposed to ultrashort laser pulses as a few-body test case We extract the degree of alignment from the non-BO wave function by means of an operator expressed in terms of pseudo-proton coordinates that mimics the BO-based definition of alignment The only essential approximation, in addition to the semiclassical electric-dipole approximation for the matter-field interaction, is the choice of time-independent explicitly correlated Gaussian basis functions. We use a variational, electric-field-dependent basis-set construction procedure, which allows us to keep the basis-set dimension low whilst capturing the main effects of electric polarization on the nuclear and electronic degrees of freedom. The basis-set construction procedure is validated by comparing with virtually exact grid-based simulations for two one-dimensional model systems: laser-driven electron dynamics in a soft attractive Coulomb potential and nuclear rovibrational dynamics in a Morse potential.

preprint2021arXiv

Various variational approximations of quantum dynamics

We investigate variational principles for the approximation of quantum dynamics that apply for approximation manifolds that do not have complex linear tangent spaces. The first one, dating back to McLachlan (1964) minimizes the residuum of the time-dependent Schrödinger equation, while the second one, originating from the lecture notes of Kramer--Saraceno (1981), imposes the stationarity of an action functional. We characterize both principles in terms of metric and a symplectic orthogonality conditions, consider their conservation properties, and derive an elementary a-posteriori error estimate. As an application, we revisit the time-dependent Hartree approximation and frozen Gaussian wave packets.

preprint2016arXiv

Discretising the Herman--Kluk Propagator

The Herman--Kluk propagator is a popular semi-classical approximation of the unitary evolution operator in quantum molecular dynamics. In this paper we formulate the Herman--Kluk propagator as a phase space integral and discretise it by Monte Carlo and quasi-Monte Carlo quadrature. Then, we investigate the accuracy of a symplectic time discretisation by combining backward error analysis with Fourier integral operator calculus. Numerical experiments for two- and six-dimensional model systems support our theoretical results.

preprint2016arXiv

Quantum expectations via spectrograms

We discuss a new phase space method for the computation of quantum expectation values in the high frequency regime. Instead of representing a wavefunction by its Wigner function, which typically attains negative values, we define a new phase space density by adding a first-order Hermite spectrogram term as a correction to the Husimi function. The new phase space density yields accurate approximations of the quantum expectation values as well as allows numerical sampling from non-negative densities. We illustrate the new method by numerical experiments in up to $128$ dimensions.

preprint2015arXiv

A new Phase Space Density for Quantum Expectations

We introduce a new density for the representation of quantum states on phase space. It is constructed as a weighted difference of two smooth probability densities using the Husimi function and first-order Hermite spectrograms. In contrast to the Wigner function, it is accessible by sampling strategies for positive densities. In the semiclassical regime, the new density allows to approximate expectation values to second order with respect to the high frequency parameter and is thus more accurate than the uncorrected Husimi function. As an application, we combine the new phase space density with Egorov's theorem for the numerical simulation of time-evolved quantum expectations by an ensemble of classical trajectories. We present supporting numerical experiments in different settings and dimensions.

preprint2014arXiv

Corrections to Wigner type phase space methods

Over decades, the time evolution of Wigner functions along classical Hamiltonian flows has been used for approximating key signatures of molecular quantum systems. Such approximations are for example the Wigner phase space method, the linearized semiclassical initial value representation, or the statistical quasiclassical method. The mathematical backbone of these approximations is Egorov's theorem. In this paper, we reformulate the well-known second order correction to Egorov's theorem as a system of ordinary differential equations and derive an algorithm with improved asymptotic accuracy for the computation of expectation values. For models with easily evaluated higher order derivatives of the classical Hamiltonian, the new algorithm's corrections are computationally less expensive than the leading order Wigner method. Numerical test calculations for a two-dimensional torsional system confirm the theoretical accuracy and efficiency of the new method.

preprint2014arXiv

Hagedorn wavepackets in time-frequency and phase space

The Hermite functions are an orthonormalbasis of the space of square integrable functions with favourable approximation properties. Allowing for a flexible localization in position and momentum, the Hagedorn wavepackets generalize the Hermite functions also to several dimensions. Using Hagedorn's raising and lowering operators, we derive explicit formulas and recurrence relations for the Wigner and FBI transform of the wavepackets and show their relation to the Laguerre polyomials.

preprint2014arXiv

Landau-Zener type surface hopping algorithms

A class of surface hopping algorithms is studied comparing two recent Landau-Zener (LZ) formulas for the probability of nonadiabatic transitions. One of the formulas requires a diabatic representation of the potential matrix while the other one depends only on the adiabatic potential energy surfaces. For each classical trajectory, the nonadiabatic transitions take place only when the surface gap attains a local minimum. Numerical experiments are performed with deterministically branching trajectories and with probabilistic surface hopping. The deterministic and the probabilistic approach confirm the good agreement of both the LZ probabilities as well the good approximation of the reference solution computed solving the Schroedinger equation via a grid based pseudo-spectral method. Visualizations of position expectations and superimposed surface hopping trajectories with reference position densities illustrate the effective dynamics of the investigated algorithms.

preprint2014arXiv

Nonadiabatic nuclear dynamics of the ammonia cation studied by surface hopping classical trajectory calculations

The Landau--Zener (LZ) type classical-trajectory surface-hopping algorithm is applied to the nonadiabatic nuclear dynamics of the ammonia cation after photoionization of the ground-state neutral molecule to the excited states of the cation. The algorithm employs the recently proposed formula for nonadiabatic LZ transition probabilities derived from the adiabatic potential energy surfaces. The evolution of the populations of the ground state and the two lowest excited adiabatic states is calculated up to 200 fs. The results agree well with quantum simulations available for the first 100 fs based on the same potential energy surfaces. Four different time scales are detected for the nuclear dynamics: Ultrafast Jahn--Teller dynamics between the excited states on a 5 fs time scale; fast transitions between the excited state and the ground state within a time scale of 20 fs; relatively slow partial conversion of a first-excited-state population to the ground state within a time scale of 100 fs; and nearly constant populations after roughly 120 fs due to a dynamical equilibrium between all three states. The latter provides a possible explanation of the experimental evidence that ammonia cation is nonfluorescent.

preprint2014arXiv

Quasi-classical description of molecular dynamics based on Egorov's theorem

Egorov's theorem on the classical propagation of quantum observables is related to prominent quasi-classical descriptions of quantum molecuar dynamics as the linearized semiclassical initial value representation (LSC-IVR), the Wigner phase space method or the statistical quasiclassical method. The error estimates show that different accuracies are achievable for the computation of expectation values and position densities. Numerical experiments for a Morse model of diatomic iodine and confined Henon-Heiles systems in various dimensions illustrate the theoretical results.

preprint2013arXiv

Propagation of Quantum Expectations with Husimi Functions

We analyse the dynamics of expectation values of quantum observables for the time-dependent semiclassical Schrödinger equation. To benefit from the positivity of Husimi functions, we switch between observables obtained from Weyl and Anti-Wick quantization. We develop and prove a second order Egorov type propagation theorem with Husimi functions by establishing transition and commutator rules for Weyl and Anti-Wick operators. We provide a discretized version of our theorem and present numerical experiments for Schrödinger equations in dimensions two and six that validate our results.

preprint2012arXiv

Wigner Measure Propagation and Conical Singularity for General Initial Data

We study the evolution of Wigner measures of a family of solutions of a Schrödinger equation with a scalar potential displaying a conical singularity. Under a genericity assumption, classical trajectories exist and are unique, thus the question of the propagation of Wigner measures along these trajectories becomes relevant. We prove the propagation for general initial data.