Researcher profile

Othmar Koch

Othmar Koch contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 17 - UnverifiedVerification L1Unclaimed author
4works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

4 published item(s)

preprint2021arXiv

Efficient Magnus-type integrators for solar energy conversion in Hubbard models

Strongly interacting electrons in solids are generically described by Hubbardtype models, and the impact of solar light can be modeled by an additional time-dependence. This yields a finite dimensional system of ordinary differential equations (ODE)s of Schrödinger type, which can be solved numerically by exponential time integrators of Magnus type. The efficiency may be enhanced by combining these with operator splittings. We will discuss several different approaches of employing exponential-based methods in conjunction with an adaptive Lanczos method for the evaluation of matrix exponentials and compare their accuracy and efficiency. For each integrator, we use defect-based local error estimators to enable adaptive time-stepping. This serves to reliably control the approximation error and reduce the computational effort

preprint2020arXiv

Chaos-induced depletion of a Bose-Einstein condensate

The mean-field limit of a bosonic quantum many-body system is described by (mostly) non-linear equations of motion which may exhibit chaos very much in the spirit of classical particle chaos, i.e. by an exponential separation of trajectories in Hilbert space with a rate given by a positive Lyapunov exponent $λ$. The question now is whether $λ$ imprints itself onto measurable observables of the underlying quantum many-body system even at finite particle numbers. Using a Bose-Einstein condensate expanding in a shallow potential landscape as a paradigmatic example for a bosonic quantum many-body system, we show, that the number of non-condensed particles is subject to an exponentially fast increase, i.e. depletion. Furthermore, we show that the rate of exponential depletion is given by the Lyapunov exponent associated with the chaotic mean-field dynamics. Finally, we demonstrate that this chaos-induced depletion is accessible experimentally through the visibility of interference fringes in the total density after time of flight, thus opening the possibility to measure $λ$, and with it, the interplay between chaos and non-equilibrium quantum matter, in a real experiment.

preprint2020arXiv

Time adaptive Zassenhaus splittings for the Schrödinger equation in the semiclassical regime

Time dependent Schrödinger equations with conservative force field U commonly constitute a major challenge in the numerical approximation, especially when they are analysed in the semiclassical regime. Extremely high oscillations originate from the semiclassical parameter, and call for appropriate methods. We propose to employ a combination of asymptotic Zassenhaus splitting with time adaptivity. While the former turns the disadvantage of the semiclassical parameter into an advantage, leading to highly efficient methods with low error constants, the latter enables to choose an optimal time step and to speed up the calculations when the oscillations subside. We support the results with numerical examples.

preprint2019arXiv

Computable upper error bounds for Krylov approximations to matrix exponentials and associated $φ$-functions

An a posteriori estimate for the error of a standard Krylov approximation to the matrix exponential is derived. The estimate is based on the defect (residual) of the Krylov approximation and is proven to constitute a rigorous upper bound on the error, in contrast to existing asymptotical approximations. It can be computed economically in the underlying Krylov space. In view of time-stepping applications, assuming that the given matrix is scaled by a time step, it is shown that the bound is asymptotically correct (with an order related to the dimension of the Krylov space) for the time step tending to zero. This means that the deviation of the error estimate from the true error tends to zero faster than the error itself. Furthermore, this result is extended to Krylov approximations of $φ$-functions and to improved versions of such approximations. The accuracy of the derived bounds is demonstrated by examples and compared with different variants known from the literature, which are also investigated more closely. Alternative error bounds are tested on examples, in particular a version based on the concept of effective order. For the case where the matrix exponential is used in time integration algorithms, a step size selection strategy is proposed and illustrated by experiments.