Researcher profile

Karsten Balzer

Karsten Balzer contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
9works
0followers
2topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

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

9 published item(s)

preprint2021arXiv

Neutralization dynamics of slow highly charged ions passing through graphene nanoflakes--an embedding self-energy approach

We study the time-dependent neutralization of a slow highly charged ion that penetrates a hexagonal hollow-centred graphene nanoflake. To compute the ultrafast charge transfer dynamics, we apply an effective Hubbard nanocluster model and use the method of nonequilibrium Green functions (NEGF) in conjunction with an embedding self-energy scheme which allows one to follow the temporal changes of the number of electrons in the nanoflake. We perform extensive simulations of the charge transfer dynamics for a broad range of ion charge states and impact velocities. The results are used to put forward a simple semi-analytical model of the neutralization dynamics that is in very good agreement with transmission experiments, in which highly charged xenon ions pass through sheets of single-layer graphene.

preprint2016arXiv

Stopping dynamics of ions passing through correlated honeycomb clusters

A combined nonequilibrium Green functions-Ehrenfest dynamics approach is developed that allows for a time-dependent study of the energy loss of a charged particle penetrating a strongly correlated system at zero and finite temperature. Numerical results are presented for finite inhomogeneous two-dimensional Fermi-Hubbard models, where the many-electron dynamics in the target are treated fully quantum mechanically and the motion of the projectile is treated classically. The simulations are based on the solution of the two-time Dyson (Keldysh-Kadanoff-Baym) equations using the second-order Born, third-order and T-matrix approximations of the self-energy. As application, we consider protons and helium nuclei with a kinetic energy between 1 and 500 keV/u passing through planar fragments of the two-dimensional honeycomb lattice and, in particular, examine the influence of electron-electron correlations on the energy exchange between projectile and electron system. We investigate the time dependence of the projectile's kinetic energy (stopping power), the electron density, the double occupancy and the photoemission spectrum. Finally, we show that, for a suitable choice of the Hubbard model parameters, the results for the stopping power are in fair agreement with ab-initio simulations for particle irradiation of single-layer graphene.

preprint2016arXiv

The auxiliary Hamiltonian approach and its generalization to non-local self-energies

The recently introduced auxiliary Hamiltonian approach [Balzer K and Eckstein M 2014 Phys. Rev. B 89 035148] maps the problem of solving the two-time Kadanoff-Baym equations onto a noninteracting auxiliary system with additional bath degrees of freedom. While the original paper restricts the discussion to spatially local self-energies, we show that there exists a rather straightforward generalization to treat also non-local correlation effects. The only drawback is the loss of time causality due to a combined singular value and eigen decomposition of the two-time self-energy, the application of which inhibits one to establish the self-consistency directly on the time step. For derivation and illustration of the method, we consider the Hubbard model in one dimension and study the decay of the Neel state in the weak-coupling regime, using the local and non-local second-order Born approximation.

preprint2015arXiv

Non-thermal melting of Neel order in the Hubbard model

We study the unitary time evolution of antiferromagnetic order in the Hubbard model after a quench starting from the perfect Néel state. In this setup, which is well suited for experiments with cold atoms, one can distinguish fundamentally different pathways for melting of long-range order at weak and strong interaction. In the Mott insulating regime, melting of long-range order occurs due to the ultra-fast transfer of energy from charge excitations to the spin background, while local magnetic moments and their exchange coupling persist during the process. The latter can be demonstrated by a local spin-precession experiment. At weak interaction, local moments decay along with the long-range order. The dynamics is governed by residual quasiparticles, which are reflected in oscillations of the off-diagonal components of the momentum distribution. Such oscillations provide an alternative route to study the prethermalization phenomenon and its influence on the dynamics away from the integrable (noninteracting) limit. The Hubbard model is solved within nonequilibrium dynamical mean-field theory, using the density matrix-renormalization group as an impurity solver.

preprint2014arXiv

Multiconfiguration time-dependent Hartree impurity solver for nonequilibrium dynamical mean-field theory

Nonequilibrium dynamical mean-field theory (DMFT) solves correlated lattice models by obtaining their local correlation functions from an effective model consisting of a single impurity in a self-consistently determined bath. The recently developed mapping of this impurity problem from the Keldysh time contour onto a time-dependent single-impurity Anderson model (SIAM) [C. Gramsch et al., Phys. Rev. B 88, 235106 (2013)] allows one to use wave function-based methods in the context of nonequilibrium DMFT. Within this mapping, long times in the DMFT simulation become accessible by an increasing number of bath orbitals, which requires efficient representations of the time-dependent SIAM wave function. These can be achieved by the multiconfiguration time-dependent Hartree (MCTDH) method and its multi-layer extensions. We find that MCTDH outperforms exact diagonalization for large baths in which the latter approach is still within reach and allows for the calculation of SIAMs beyond the system size accessible by exact diagonalization. Moreover, we illustrate the computation of the self-consistent two-time impurity Green's function within the MCTDH second quantization representation.

preprint2013arXiv

Auxiliary Hamiltonian representation of the nonequilibrium Dyson equation

The nonequilibrium Dyson (or Kadanoff-Baym) equation, which is an equation of motion with long-range memory kernel for real-time Green functions, underlies many numerical approaches based on the Keldysh formalism. In this paper we map the problem of solving the Dyson equation in real-time onto a noninteracting auxiliary Hamiltonian with additional bath degrees of freedom. The solution of the auxiliary model does not require the evaluation of a memory kernel and can thus be implemented in a very memory efficient way. The mapping is derived for a self-energy which is local in space and is thus directly applicable within nonequilibrium dynamical mean-field theory (DMFT). We apply the method to study the interaction quench in the Hubbard model for an optical lattice with a narrow confinement, using inhomogeneous DMFT in combination with second-order weak-coupling perturbation theory. We find that, although the quench excites pronounced density oscillations, signatures of the two-stage relaxation similar to the homogeneous system can be observed by looking at the time-dependent occupations of natural orbitals.

preprint2013arXiv

Hamiltonian-based impurity solver for nonequilibrium dynamical mean-field theory

We derive an exact mapping from the action of nonequilibrium dynamical mean-field theory (DMFT) to a single-impurity Anderson model (SIAM) with time-dependent parameters, which can be solved numerically by exact diagonalization. The representability of the nonequilibrium DMFT action by a SIAM is established as a rather general property of nonequilibrium Green functions. We also obtain the nonequilibrium DMFT equations using the cavity method alone. We show how to numerically obtain the SIAM parameters using Cholesky or eigenvector matrix decompositions. As an application, we use a Krylov-based time propagation method to investigate the Hubbard model in which the hopping is switched on, starting from the atomic limit. Possible future developments are discussed.

preprint2012arXiv

Quantum Breathing Mode of Interacting Particles in a One-dimensional Harmonic Trap

Extending our previous work, we explore the breathing mode---the [uniform] radial expansion and contraction of a spatially confined system. We study the breathing mode across the transition from the ideal quantum to the classical regime and confirm that it is not independent of the pair interaction strength (coupling parameter). We present the results of time-dependent Hartree-Fock simulations for 2 to 20 fermions with Coulomb interaction and show how the quantum breathing mode depends on the particle number. We validate the accuracy of our results, comparing them to exact Configuration Interaction results for up to 8 particles.

preprint2009arXiv

Finite elements and the discrete variable representation in nonequilibrium Green's function calculations. Atomic and molecular models

In this contribution, we discuss the finite-element discrete variable representation (FE-DVR) of the nonequilibrium Green's function and its implications on the description of strongly inhomogeneous quantum systems. In detail, we show that the complementary features of FEs and the DVR allows for a notably more efficient solution of the two-time Schwinger/Keldysh/Kadanoff-Baym equations compared to a general basis approach. Particularly, the use of the FE-DVR leads to an essential speedup in computing the self-energies. As atomic and molecular examples we consider the He atom and the linear version of H$_3^+$ in one spatial dimension. For these closed-shell models we, in Hartree-Fock and second Born approximation, compute the ground-state properties and compare with the exact findings obtained from the solution of the few-particle time-dependent Schrödinger equation.