Lorenz gauge formulation for time-dependent density functional theory

We describe the inclusion of electrodynamic fields in time-dependent density functional theory (TDDFT) by incorporating both the induced scalar and vector potentials within the time-dependent Kohn-Sham equation. The Hamiltonian is described in both the Coulomb and Lorenz gauges, and the advantages of the latter are outlined. Integral expressions are defined for the retarded potentials of each gauge and a methodological approach to evaluate these nontrivial expressions with low computational cost is adopted. Various molecular structures of relatively small sizes are studied, including water, benzene, and conductive carbon chains. Dipole strengths resulting from both pulse and boost excitations suggest a preserved gauge invariance of the proposed formal approach to TDDFT in the weak magnetic field limit.