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 Schr\"odinger equation for a hydrogen atom interacting with a linearly polarized laser pulse, whose electric field is given by E (t) =E₀f (t) cos (₀t+) with Gaussian carrier envelope f (t) =exp (-t^2∕^2). The carrier frequency ₀ is fixed to correspond to a wavelength of 8000. 3em{0ex}nm. Spectra for a selection of pulses, for which the intensity I₀=c₀E₀^2, duration T, and carrier-envelope phase are systematically varied, show that, depending on, all three forms are in good agreement for ``weak'' pulses with I₀I₁). 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₅) ⟩ or dipole velocity ⟨ \. {}z (t₅) ⟩ at the end (t₅) 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₅) ⟩0 and therefore ⟨z (t) ⟩ may increase without bound (i. e. , the atom may ionize), depending on, 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 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.