Optimal-continuum and multicentered Gaussian basis sets for high-harmonic generation spectroscopy

High-harmonic generation (HHG) is used to produce coherent XUV and soft X-ray radiation with atto-second resolution and is a sensitive tool for probing atomic and molecular structures. In this work, we have used time-dependent configuration interaction with a Gaussian basis set to compute the HHG spectrum of the hydrogen atom. To get a correct description of the HHG optical spectrum, the Gaussian basis set has to provide an accurate representation of the bound and the continuum states. Two strategies have been proposed: (1) multicentered (defining ghost atoms around the hydrogen) and (2) optimal-continuum Gaussian basis sets. We have systematically investigated these two approaches for the hydrogen atom, which permits a non-biased analysis of the basis set. Several basis sets have been constructed and tested by combining multicentered and optimal-continuum functions together in order to obtain a reliable and accurate Gaussian basis set to be used for HHG. We have studied the effect of changing the number of ghost atoms and the distance between the ghosts and the hydrogen atom, with and without optimal-continuum Gaussian functions. We conclude that multicentered basis sets are less efficient than basis sets using only optimal-continuum Gaussian functions for a proper description of HHG.