Eliminating artificial boundary conditions in time-dependent density functional theory using Fourier contour deformation

Abstract

We present an efficient method for propagating the time-dependent Kohn-Sham equations in free space, based on the recently introduced Fourier contour deformation (FCD) approach. For potentials which are constant outside a bounded domain, FCD yields a high-order accurate numerical solution of the time-dependent Schrodinger equation directly in free space, without the need for artificial boundary conditions. Of the many existing artificial boundary condition schemes, FCD is most similar to an exact nonlocal transparent boundary condition, but it works directly on Cartesian grids in any dimension, and runs on top of the fast Fourier transform rather than fast algorithms for the application of nonlocal history integral operators. We adapt FCD to time-dependent density functional theory (TDDFT), and describe a simple algorithm to smoothly and automatically truncate long-range Coulomb-like potentials to a time-dependent constant outside of a bounded domain of interest, so that FCD can be used. This approach eliminates errors originating from the use of artificial boundary conditions, leaving only the error of the potential truncation, which is controlled and can be systematically reduced. The method enables accurate simulations of ultrastrong nonlinear electronic processes in molecular complexes in which the inteference between bound and continuum states is of paramount importance. We demonstrate results for many-electron TDDFT calculations of absorption and strong field photoelectron spectra for one and two-dimensional models, and observe a significant reduction in the size of the computational domain required to achieve high quality results, as compared with the popular method of complex absorbing potentials.

References in corpus (26)

• Correlation energy and spin polarization in the 2D electron gas
• Octopus, a computational framework for exploring light-driven phenomena and quantum dynamics in extended and finite systems
• Time-dependent complete-active-space self-consistent field method for multielectron dynamics in intense laser fields
• Notes on Perfectly Matched Layers (PMLs)
• Fast and oblivious convolution quadrature
• Exterior complex scaling as a perfect absorber in time-dependent problems
• Modeling electron dynamics coupled to continuum states in finite volumes
• Ab-initio angle and energy resolved photoelectron spectroscopy with time-dependent density-functional theory
• Ground state properties of the one dimensional Coulomb gas
• Density functional theory beyond the linear regime: Validating adiabatic LDA
• iSURF: A family of infinite-time surface flux methods
• Efficient and accurate modeling of electron photoemission in nanostructures with TDDFT
• Implementation of infinite-range exterior complex scaling to the time-dependent complete-active-space self-consistent-field method
• Harmonic generation in ring-shaped molecules
• Setting the photoelectron clock through molecular alignment
• Density Functional Resonance Theory of Unbound Electronic Systems
• tRecX -- an environment for solving time-dependent Schrödinger-like problems
• Multiple orbital effects in laser-induced electron diffraction of aligned molecules
• Perfectly Matched Layers versus discrete transparent boundary conditions in quantum device simulations
• First-principles simulations for attosecond photoelectron spectroscopy based on time-dependent density functional theory
• Stark Ionization of Atoms and Molecules within Density Functional Resonance Theory
• Exact transparent boundary conditions for the parabolic wave equations with linear and quadratic potentials
• A high-order integral equation-based solver for the time-dependent Schrodinger equation
• A study of one-dimensional transport of Bose-Einstein condensates using exterior complex scaling
• A quantum simulation of dissociative ionization of $H_2^{+}$ in full dimensionality with time dependent surface flux method
• Transparent Boundary Conditions for the Time-Dependent Schrödinger Equation with a Vector Potential