Quantum simulation of high-order harmonic spectra of the hydrogen atom

Three alternative forms of harmonic spectra, based on the dipole moment, dipole velocity, and dipole acceleration, are compared by a numerical solution of the Schrodinger equation for a hydrogen atom interacting with a linearly polarized laser pulse, whose electric field is given by E(t)=E(0)f(t)cos(omega(0)t+eta) with Gaussian carrier envelope f(t)=exp(-t(2)/delta(2)). The carrier frequency omega(0) is fixed to correspond to a wavelength of 800 nm. Spectra for a selection of pulses, for which the intensity I-0=c epsilon E-0(0)2, duration T proportional to delta, and carrier-envelope phase eta are systematically varied, show that, depending on eta, all three forms are in good agreement for "weak" pulses with I-0 < I-b, the over-barrier ionization threshold, but that marked differences among the three appear as the pulse becomes shorter and stronger (I-0>I-b). Except for scalings by powers of the harmonic frequency, the three forms differ from one another only by "limit contributions" proportional to the expectation values of the dipole moment < z(t(f))> or dipole velocity < z(t(f))> at the end (t(f)) of the pulse. For long, weak pulses the limit contributions are negligible, whereas for short, strong ones they are not. In the short, strong limit, where < z(t(f))>not equal 0 and therefore < z(t)> may increase without bound (i.e., the atom may ionize), depending on eta, an "infinite-time" spectrum based on the acceleration form provides a convenient computational pathway to the corresponding infinite-time dipole-velocity spectrum, which is related directly to the experimentally measured "harmonic photon number spectrum" (HPNS). For short, intense pulses the HPNS is quite sensitive to eta and exhibits not only the usual odd harmonics but also even ones. The analysis also reveals that most of the harmonic photons are emitted during the passage of the pulse. Because of the divergence of < z(t)> the dipole-moment form does not provide a numerically reliable route to the harmonic spectrum for very short (few-cycle), very intense laser pulses.