Development of Real-Time TDDFT Program with k-Point Sampling and DFT + U in a Gaussian and Plane Waves Framework

We developed a -point sampling real-time TDDFT (RT-TDDFT) program within the Gaussian and plane waves (GPW) framework of the CP2K software suite. In addition to standard real-time propagation of time-dependent Kohn-Sham orbitals, we make use of symmetry-based -point reduction and -point parallelization schemes so that our RT-TDDFT program in the GPW framework is feasible for practical large-scale calculations. We also implemented DFT + as a relevant extension for real-time simulations of systems with strong electron correlations. In particular, we extended the "tensorial" subspace representation approach for DFT + , following the formulation in [Chai, Z., et al. , 8984], to -point sampling RT-TDDFT. Our extension, which is, to our knowledge, the first application of the "tensorial" subspace representation approach to -point sampling RT-TDDFT, is found to be robust and efficient with small additional costs owing to the locality of Gaussian basis functions, indicating that it is a promising approach to RT-TDDFT + for solids. We show details of our implementation in CP2K and the results of our benchmark calculations.