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

by   Jason Kaye, et al.

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.


A high-order integral equation-based solver for the time-dependent Schrodinger equation

We introduce a numerical method for the solution of the time-dependent S...

Absorbing boundary conditions for the time-dependent Schrödinger-type equations in R^3

Absorbing boundary conditions are presented for three-dimensional time-d...

Numerical analysis for electromagnetic scattering from nonlinear boundary conditions

This work studies time-dependent electromagnetic scattering from obstacl...

A spectral deferred correction method for incompressible flow with variable viscosity

This paper presents a semi-implicit spectral deferred correction (SDC) m...

Time-dependent electromagnetic scattering from thin layers

The scattering of electromagnetic waves from obstacles with wave-materia...

A fourth-order compact time-splitting method for the Dirac equation with time-dependent potentials

In this paper, we present an approach to deal with the dynamics of the D...

Fast Ewald summation for Stokes flow with arbitrary periodicity

A fast and spectrally accurate Ewald summation method for the evaluation...

Please sign up or login with your details

Forgot password? Click here to reset