Source author record

Heiko Appel

Heiko Appel 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

22works
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

22 published item(s)

preprint2022arXiv

Frequency-dependent Sternheimer linear-response formalism for strongly coupled light-matter systems

The rapid progress in quantum-optical experiments especially in the field of cavity quantum electrodynamics and nanoplasmonics, allows to substantially modify and control chemical and physical properties of atoms, molecules and solids by strongly coupling to the quantized field. Alongside such experimental advances has been the recent development of ab-initio approaches such as quantum electrodynamical density-functional theory (QEDFT) that is capable of describing these strongly coupled systems from first-principles. To investigate response properties of relatively large systems coupled to a wide range of photon modes, ab-initio methods that scale well with system size become relevant. In light of this, we extend the linear-response Sternheimer approach within the framework of QEDFT to efficiently compute excited-state properties of strongly coupled light-matter systems. Using this method, we capture features of strong light-matter coupling both in the dispersion and absorption properties of a molecular system strongly coupled to the modes of a cavity. We exemplify the efficiency of the Sternheimer approach by coupling the matter system to the continuum of an electromagnetic field. We observe changes in the spectral features of the coupled system as Lorentzian line shapes turn into Fano resonances when the molecule interacts strongly with the continuum of modes. This work provides an alternative approach for computing efficiently excited-state properties of large molecular systems interacting with the quantized electromagnetic field.

preprint2021arXiv

Photoionization and transient Wannier-Stark ladder in silicon: First principle simulations versus Keldysh theory

Nonlinear photoionization of dielectrics and semiconductors is widely treated in the frames of the Keldysh theory whose validity is limited to small photon energies compared to the band gap and relatively low laser intensities. The time-dependent density functional theory (TDDFT) simulations, which are free of these limitations, enable to gain insight into non-equilibrium dynamics of the electronic structure. Here we apply the TDDFT to investigate photoionization of silicon crystal by ultrashort laser pulses in a wide range of laser wavelengths and intensities and compare the results with predictions of the Keldysh theory. Photoionization rates derived from the simulations considerably exceed the data obtained with the Keldysh theory within the validity range of the latter. Possible reasons of the discrepancy are discussed and we provide fundamental data on the photoionization rates beyond the limits of the Keldysh theory. By investigating the features of the Stark shift as a function of photon energy and laser field strength, a manifestation of the transient Wannier-Stark ladder states have been revealed which become blurred with increasing laser field strength. Finally, it is shown that the TDDFT simulations can potentially provide reliable data on the electron damping time that is of high importance for large-scale modeling.

preprint2020arXiv

Chemistry in Quantum Cavities: Exact Results, the Impact of Thermal Velocities and Modified Dissociation

In recent years tremendous progress in the field of light-matter interactions has unveiled that strong coupling to the modes of an optical cavity can alter chemistry even at room temperature. Despite these impressive advances, many fundamental questions of chemistry in cavities remain unanswered. This is also due to a lack of exact results that can be used to validate and benchmark approximate approaches. In this work we provide such reference calculations from exact diagonalisation of the Pauli-Fierz Hamiltonian in the long-wavelength limit with an effective cavity mode. This allows us to investigate the reliability of the ubiquitous Jaynes-Cummings model not only for electronic but also for the case of ro-vibrational transitions. We demonstrate how the commonly ignored thermal velocity of charged molecular systems can influence chemical properties, while leaving the spectra invariant. Furthermore, we show the emergence of new bound polaritonic states beyond the dissociation energy limit.

preprint2019arXiv

Ab Initio Bethe-Salpeter Equation Approach to Neutral Excitations in Molecules with Numeric Atom-Centered Orbitals

The Bethe-Salpeter equation (BSE) based on GW quasiparticle levels is a successful approach for calculating the optical gaps and spectra of solids and also for predicting the neutral excitations of small molecules. We here present an all-electron implementation of the GW+BSE formalism for molecules, using numeric atom-centered orbital (NAO) basis sets. We present benchmarks for low-lying excitation energies for a set of small organic molecules, denoted in the literature as "Thiel's set". Literature reference data based on Gaussian-type orbitals are reproduced to about one meV precision for the molecular benchmark set, when using the same GW quasiparticle energies and basis sets as the input to the BSE calculations. For valence correlation consistent NAO basis sets, as well as for standard NAO basis sets for ground state density-functional theory with extended augmentation functions, we demonstrate excellent convergence of the predicted low-lying excitations to the complete basis set limit. A simple and affordable augmented NAO basis set denoted "tier2+aug2" is recommended as a particularly efficient formulation for production calculations. We finally demonstrate that the same convergence properties also apply to linear-response time-dependent density functional theory within the NAO formalism.

preprint2019arXiv

Benchmarking Semiclassical and Perturbative Methods for Real-time Simulations of Cavity-Bound Emission and Interference

We benchmark a selection of semiclassical and perturbative dynamics techniques by investigating the correlated evolution of a cavity-bound atomic system to assess their applicability to study problems involving strong light-matter interactions in quantum cavities. The model system of interest features spontaneous emission, interference, and strong coupling behaviour, and necessitates the consideration of vacuum fluctuations and correlated light-matter dynamics. We compare a selection of approximate dynamics approaches including fewest switches surface hopping, multi-trajectory Ehrenfest dynamics, linearized semiclasical dynamics, and partially linearized semiclassical dynamics. Furthermore, investigating self-consistent perturbative methods, we apply the Bogoliubov-Born-Green-Kirkwood-Yvon hierarchy in the second Born approximation. With the exception of fewest switches surface hopping, all methods provide a reasonable level of accuracy for the correlated light-matter dynamics, with most methods lacking the capacity to fully capture interference effects.

preprint2019arXiv

Octopus, a computational framework for exploring light-driven phenomena and quantum dynamics in extended and finite systems

Over the last years extraordinary advances in experimental and theoretical tools have allowed us to monitor and control matter at short time and atomic scales with a high-degree of precision. An appealing and challenging route towards engineering materials with tailored properties is to find ways to design or selectively manipulate materials, especially at the quantum level. To this end, having a state-of-the-art ab initio computer simulation tool that enables a reliable and accurate simulation of light-induced changes in the physical and chemical properties of complex systems is of utmost importance. The first principles real-space-based Octopus project was born with that idea in mind, providing an unique framework allowing to describe non-equilibrium phenomena in molecular complexes, low dimensional materials, and extended systems by accounting for electronic, ionic, and photon quantum mechanical effects within a generalized time-dependent density functional theory framework. The present article aims to present the new features that have been implemented over the last few years, including technical developments related to performance and massive parallelism. We also describe the major theoretical developments to address ultrafast light-driven processes, like the new theoretical framework of quantum electrodynamics density-functional formalism (QEDFT) for the description of novel light-matter hybrid states. Those advances, and other being released soon as part of the Octopus package, will enable the scientific community to simulate and characterize spatial and time-resolved spectroscopies, ultrafast phenomena in molecules and materials, and new emergent states of matter (QED-materials).

preprint2019arXiv

Polaritonic Coupled-Cluster Theory

We develop coupled-cluster theory for systems of electrons strongly coupled to photons, providing a promising theoretical tool in polaritonic chemistry with a perspective of application to all types of fermion-boson coupled systems. We show benchmark results for model molecular Hamiltonians coupled to cavity photons. By comparing to full configuration interaction results for various ground-state properties and optical spectra, we demonstrate that our method captures all key features present in the exact reference, including Rabi splittings and multi-photon processes. Further, a path on how to incorporate our bosonic extension of coupled-cluster theory into existing quantum chemistry programs is given.

preprint2018arXiv

Real-time solutions of coupled Ehrenfest-Maxwell-Pauli-Kohn-Sham equations: fundamentals, implementation, and nano-optical applications

We present the theoretical foundations and the implementation details of a density-functional approach for coupled photons, electrons, and effective nuclei in non-relativistic quantum electrodynamics. Starting point of the formalism is a generalization of the Pauli-Fierz field theory for which we establish a one-to-one correspondence between external fields and internal variables. Based on this correspondence, we introduce a Kohn-Sham construction which provides a computationally feasible approach for ab-initio light-matter interactions. In the mean-field limit for the effective nuclei the formalism reduces to coupled Ehrenfest-Maxwell-Pauli-Kohn-Sham equations. We present an implementation of the approach in the real-space real-time code Octopus. For the implementation we use the Riemann-Silberstein formulation of classical electrodynamics and rewrite Maxwell's equations in Schrödinger form. This allows us to use existing time-evolution algorithms developed for quantum-mechanical systems also for Maxwell's equations. Overall, our approach is ideally suited for applications in nano-optics, nano-plasmonics, (photo) electrocatalysis, light-matter coupling in 2D materials, cases where laser pulses carry orbital angular momentum, or light-tailored chemical reactions in optical cavities to name but a few.

preprint2016arXiv

Correlated Electron-Nuclear Dynamics with Conditional Wave Functions

The molecular Schrödinger equation is rewritten in terms of non-unitary equations of motion for the nuclei (or electrons) that depend parametrically on the configuration of an ensemble of generally defined electronic (or nuclear) trajectories. This scheme is exact and does not rely on the tracing-out of degrees of freedom. Hence, the use of trajectory-based statistical techniques can be exploited to circumvent the calculation of the computationally demanding Born-Oppenheimer potential-energy surfaces and non-adiabatic coupling elements. The concept of potential-energy surface is restored by establishing a formal connection with the exact factorization of the full wave function. This connection is used to gain insight from a simplified form of the exact propagation scheme.

preprint2016arXiv

Systematic construction of density functionals based on matrix product state computations

We propose a systematic procedure for the approximation of density functionals in density functional theory that consists of two parts. First, for the efficient approximation of a general density functional, we introduce an efficient ansatz whose non-locality can be increased systematically. Second, we present a fitting strategy that is based on systematically increasing a reasonably chosen set of training densities. We investigate our procedure in the context of strongly correlated fermions on a one-dimensional lattice in which we compute accurate training densities with the help of matrix product states. Focusing on the exchange-correlation energy, we demonstrate how an efficient approximation can be found that includes and systematically improves beyond the local density approximation. Importantly, this systematic improvement is shown for target densities that are quite different from the training densities.

preprint2015arXiv

Exact Maps in Density Functional Theory for Lattice Models

In the present work, we employ exact diagonalization for model systems on a real-space lattice to explicitly construct the exact density-to-potential and for the first time the exact density-to-wavefunction map that underly the Hohenberg-Kohn theorem in density functional theory. Having the explicit wavefunction-to- density map at hand, we are able to construct arbitrary observables as functionals of the ground-state density. We analyze the density-to-potential map as the distance between the fragments of a system increases and the correlation in the system grows. We observe a feature that gradually develops in the density-to-potential map as well as in the density-to-wavefunction map. This feature is inherited by arbitrary expectation values as functional of the ground-state density. We explicitly show the excited-state energies, the excited-state densities, and the correlation entropy as functionals of the ground-state density. All of them show this exact feature that sharpens as the coupling of the fragments decreases and the correlation grows. We denominate this feature as intra-system steepening. We show that for fully decoupled subsystems the intra-system steepening transforms into the well-known inter-system derivative discontinuity. An important conclusion is that for e.g. charge transfer processes between localized fragments within the same system it is not the usual inter-system derivative discontinuity that is missing in common ground-state functionals, but rather the differentiable intra-system steepening that we illustrate in the present work.

preprint2015arXiv

Kohn-Sham Approach to Quantum Electrodynamical Density Functional Theory: Exact Time-Dependent Effective Potentials in Real Space

The density-functional approach to quantum electrodynamics is extending traditional density-functional theory and opens the possibility to describe electron-photon interactions in terms of effective Kohn-Sham potentials. In this work, we numerically construct the exact electron-photon Kohn-Sham potentials for a prototype system which consists of a trapped electron coupled to a quantized electromagnetic mode in an optical high-Q cavity. While the effective current that acts on the photons is known explicitly, the exact effective potential that describes the forces exerted by the photons on the electrons is obtained from a fixed-point inversion scheme. This procedure allows us to uncover important beyond-mean-field features of the effective potential which mark the breakdown of classical light-matter interactions. We observe peak and step structures in the effective potentials, which can be attributed solely to the quantum nature of light, i.e., they are real-space signatures of the photons. Our findings show how the ubiquitous dipole interaction with a classical electromagnetic field has to be modified in real-space in order to take the quantum nature of the electromagnetic field fully into account.

preprint2014arXiv

A quantum reactive scattering perspective on electronic nonadiabaticity

Based on quantum reactive-scattering theory, we propose a method for studying the electronic nonadiabaticity in collision processes involving electron-ion rearrangements. We investigate the state-to-state transition probability for electron-ion rearrangements with two comparable approaches. In the first approach the information of the electron is only contained in the ground-state Born-Oppenheimer potential-energy surface, which is the starting point of common reactive-scattering calculations. In the second approach, the electron is explicitly taken into account and included in the calculations at the same level as the ions. Hence, the deviation in the results between the two approaches directly reflects the electronic nonadiabaticity during the collision process. To illustrate the method, we apply it to the well-known proton-transfer model of Shin and Metiu (one electron and three ions), generalized by us in order to allow for reactive scattering channels. It is shown that our explicit electron approach is able to capture electronic nonadiabaticity and the renormalization of the reaction barrier near the classical turning points of the potential in nuclear configuration space. In contrast, system properties near the equilibrium geometry of the asymptotic scattering channels are hardly affected by electronic nonadiabatic effects. We also present an analytical expression for the transition amplitude of the asymmetric proton-transfer model based on the direct evaluation of integrals over the involved Airy functions.

preprint2014arXiv

Förster-induced energy transfer in functionalized graphene

Carbon nanostructures are ideal substrates for functionalization with molecules, since they consist of a single atomic layer giving rise to an extraordinary sensitivity to changes in their surrounding. The functionalization opens a new research field of hybrid nanostructures with tailored properties. Here, we present a microscopic view on the substrate-molecule interaction in the exemplary hybrid material consisting of graphene functionalized with perylene molecules. First experiments on similar systems have been recently realized illustrating an extremely efficient transfer of excitation energy from adsorbed molecules to the carbon substrate - a process with a large application potential for high-efficiency photovoltaic devices and biomedical imaging and sensing. So far, there has been no microscopically founded explanation for the observed energy transfer. Based on first-principle calculations, we have explicitly investigated the different transfer mechanisms revealing the crucial importance of Förster coupling. Due to the efficient Coulomb interaction in graphene, we obtain strong Förster rates in the range of 1/fs. We investigate its dependence on the substrate-molecule distance R and describe the impact of the momentum transfer q for an efficient energy transfer. Furthermore, we find that the Dexter transfer mechanism is negligibly small due to the vanishing overlap between the involved strongly localized orbital functions. The gained insights are applicable to a variety of carbon-based hybrid nanostructures.

preprint2014arXiv

Non-adiabatic and time-resolved photoelectron spectroscopy for molecular systems

We quantify the non-adiabatic contributions to the vibronic sidebands of equilibrium and explicitly time-resolved non-equilibrium photoelectron spectra for a vibronic model system of Trans-Polyacetylene. Using exact diagonalization, we directly evaluate the sum-over-states expressions for the linear-response photocurrent. We show that spurious peaks appear in the Born-Oppenheimer approximation for the vibronic spectral function, which are not present in the exact spectral function of the system. The effect can be traced back to the factorized nature of the Born-Oppenheimer initial and final photoemission states and also persists when either only initial, or final states are replaced by correlated vibronic states. Only when correlated initial and final vibronic states are taken into account, the spurious spectral weights of the Born-Oppenheimer approximation are suppressed. In the non-equilibrium case, we illustrate for an initial Franck-Condon excitation and an explicit pump-pulse excitation how the vibronic wavepacket motion of the system can be traced in the time-resolved photoelectron spectra as function of the pump-probe delay.

preprint2014arXiv

Optimized Effective Potential for Quantum Electrodynamical Time-Dependent Density Functional Theory

We propose a practical approximation to the exchange-correlation functional of (time-dependent) density functional theory for many-electron systems coupled to photons. The (time non-local) optimized effective potential (OEP) equation for the electron- photon system is derived. We test the new approximation in the Rabi model from weak to strong coupling regimes. It is shown that the OEP (i) improves the classical description, (ii) reproduces the quantitative behavior of the exact ground-state properties and (iii) accurately captures the dynamics entering the ultra-strong coupling regime. The present formalism opens the path to a first-principles description of correlated electron-photon systems, bridging the gap between electronic structure methods and quantum optics for real material applications.

preprint2014arXiv

Quantum Electrodynamical Density-Functional Theory: Bridging Quantum Optics and Electronic-Structure Theory

In this work we give a comprehensive derivation of an exact and numerically feasible method to perform ab-initio calculations of quantum particles interacting with a quantized electromagnetic field. We present a hierachy of density-functional-type theories that describe the interaction of charged particles with photons and introduce the appropriate Kohn-Sham schemes. We show how the evolution of a system described by quantum electrodynamics in Coulomb gauge is uniquely determined by its initial state and two reduced quantities. These two fundamental observables, the polarization of the Dirac field and the vector potential of the photon field, can be calculated by solving two coupled, non-linear evolution equations without the need to explicitly determine the (numerically infeasible) many-body wave function of the coupled quantum system. To find reliable approximations to the implicit functionals we present the according Kohn-Sham construction. In the non-relativistic limit this density-functional-type theory of quantum electrodynamics reduces to the density-functional reformulation of the Pauli-Fierz Hamiltonian, which is based on the current density of the electrons and the vector potential of the photon field. By making further approximations, e.g. restricting the allowed modes of the photon field, we derive further density functional-type theories of coupled matter-photon systems for the corresponding approximate Hamiltonians. In the limit of only two sites and one mode we deduce the according effective theory for the two-site Hubbard model coupled to one photonic mode. This model system is used to illustrate the basic ideas of a density-functional reformulation in great detail and we present the exact Kohn-Sham potentials for our coupled matter-photon model system.

preprint2013arXiv

Measuring excitation-energy transfer with a real-time time-dependent density functional theory approach

We investigate the time an electronic excitation travels in a supermolecular setup using a measurement process in an open quantum-system framework. The approach is based on the stochastic Schrödinger equation and uses a Hamiltonian from time-dependent density functional theory (TDDFT). It treats electronic-structure properties and intermolecular coupling on the level of TDDFT, while it opens a route to the description of dissipation and relaxation via a bath operator that couples to the dipole moment of the density. Within our study, we find that in supermolecular setups small deviations of the electronic structure from the perfectly resonant case have only minor influence on the pathways of excitation-energy transfer, thus lead to similar transfer times. Yet, sizable defects cause notable slowdown of the energy spread.

preprint2013arXiv

Reduced purities as measures of decoherence in many-electron systems

A hierarchy of measures of decoherence for many-electron systems that is based on the purity and the hierarchy of reduced electronic density matrices is presented. These reduced purities can be used to characterize electronic decoherence in the common case when the many-body electronic density matrix is not known and only reduced information about the electronic subsystem is available. Being defined from reduced electronic quantities, the interpretation of the reduced purities is more intricate than the usual (many-body) purity. This is because the nonidempotency of the $r$-body reduced electronic density matrix that is the basis of the reduced purity measures can arise due to decoherence or due to electronic correlations. To guide the interpretation, explicit expressions are provided for the one-body and two-body reduced purities for a general electronic state. Using them, the information content and structure of the one-body and two-body reduced purities is established, and limits on the changes that decoherence can induce are elucidated. The practical use of the reduced purities to understand decoherence dynamics in many-electron systems is exemplified through an analysis of the electronic decoherence dynamics in a model molecular system.

preprint2013arXiv

The time-dependent exchange-correlation functional for a Hubbard dimer: quantifying non-adiabatic effect

We address and quantify the role of non-adiabaticity ("memory effects") in the exchange-correlation (xc) functional of time-dependent density functional theory (TDDFT) for describing non-linear dynamics of many-body systems. Time-dependent resonant processes are particularly challenging for available TDDFT approximations, due to their strong non-linear and non-adiabatic character. None of the known approximate density functionals are able to cope with this class of problems in a satisfactory manner. In this work we look at the prototypical example of the resonant processes by considering Rabi oscillations within the exactly soluble 2-site Hubbard model. We construct the exact adiabatic xc functional and show that (i) it does not reproduce correctly resonant Rabi dynamics, (ii) there is a sizable non-adiabatic contribution to the exact xc potential, which turns out to be small only at the beginning and at the end of the Rabi cycle when the ground state population is dominant. We then propose a "two-level" approximation for the time-dependent xc potential which can capture Rabi dynamics in the 2-site problem. It works well both for resonant and for detuned Rabi oscillations and becomes essentially exact in the linear response regime. This new, fully non-adiabatic and explicit density functional constitutes one of the main results of the present work.

preprint2011arXiv

Stochastic quantum molecular dynamics for finite and extended systems

We present a detailed account of the technical aspects of stochastic quantum molecular dynamics, an approach introduced recently by the authors [H. Appel and M. Di Ventra, Phys. Rev. B 80 212303 (2009)] to describe coupled electron-ion dynamics in open quantum systems. As example applications of the method we consider both finite systems with and without ionic motion, as well as describe its applicability to extended systems in the limit of classical ions. The latter formulation allows the study of important phenomena such as decoherence and energy relaxation in bulk systems and surfaces in the presence of time-dependent fields.

preprint2009arXiv

Stochastic Quantum Molecular Dynamics

An approach to correlated dynamics of quantum nuclei and electrons both in dynamical interaction with external environments is presented. This stochastic quantum molecular dynamics rests on a theorem that establishes a one-to-one correspondence between the total ensemble-averaged current density of interacting nuclei and electrons and a given external vector potential. The theory allows for a first-principles description of phenomena previously inaccessible via standard quantum molecular dynamics such as electronic and nuclear relaxation in photochemistry, dissipative correlated electron-ion dynamics in intense laser fields, nuclear dephasing, etc. As a demonstration of the approach, we discuss the rotational relaxation of 4-(N,N-dimethylamino)benzonitrile in a uniform bath in the limit of classical nuclei.