We present a theoretical approach for the simulation of time-resolved harmonic spectra, including the effect of nuclear dynamics, which is applicable to complex systems involving many nuclear degrees of freedom. The method is based on the combination of our semiclassical field-induced surface hopping approach for the treatment of laser-induced nuclear dynamics with the time-dependent density functional theory for electron dynamics. We apply our method to the simulation of ultrafast nonadiabatic dynamics and time-resolved harmonic spectra in small silver clusters (Ag 2 and Ag 8 ), which exhibit discrete molecularlike electronic transitions. We demonstrate that the harmonic signal is highly sensitive to the nuclear dynamics and thus can be used as a probe of coupled electron-nuclear dynamics, which is complementary to common pump-probe methods such as time-resolved photoelectron spectroscopy. Our simulations allowed us also to determine the mechanism and the time scale of nonradiative relaxation in the "magic" Ag 8 cluster and have provided a fundamental insight into ultrafast dynamics of metal nanoclusters in the size regime where "each atom counts." The excited-state dynamics of Ag 8 involves an isomerization process from the initial structure with T d symmetry to the quadratic antiprism structure with D 4d symmetry which takes place on a time scale of ∼600 fs and is clearly identified in a time-resolved harmonic signal. Our theoretical approach is generally applicable for the prediction of time-resolved harmonic spectra in complex systems with many nuclear degrees freedom and should serve to stimulate new ultrafast experiments utilizing harmonic signals as a probe for nonadiabatic processes in molecular systems.PHYSICAL REVIEW A 83, 033408 (2011)where χ (N) I (R,t) and χ (N−1) J (R,E,t) represent the nuclear wave packet in the bound and continuum states, respectively. 033408-4 SIMULATION OF LASER-INDUCED COUPLED ELECTRON-. . . PHYSICAL REVIEW A 83, 033408 (2011) (N) I (r; R) are the eigenfunctions of the N-electron Hamiltonian, while the antisymmetrized product A[ (N−1) J (r; R)φ J (E)] represents the continuum eigenfunctions of the combined ion-free electron Hamiltonian. In this product (N−1) J (r; R) is the J th cationic state and φ J (E) is a free electron scattering state. The summation extends over the whole range of singly ionized states. The wave-function ansatz in Eq. (17) can be inserted in the full electron-nuclear TDSE including the coupling to the electric field and a set of equations for the time evolution of the continuum portion of the nuclear wave packet of the ionized system χ (N−1) J (R,E,t) can be derived as ihχ (N−1) J (R,E,t) = T + E (N−1) J + E χ (N−1) J (R,E,t) − I ε(t) · µ I J (R,E)χ (N) I (R,t). (18) The semiclassical limit of this equation can be obtained by reducing the wave packet χ (N−1) J (R,E,t) to swarms of trajectories, as described previously [18]. This yields a set of equations for the time evolution of the amplitudes c (N−1) J (E,t) associated with each trajectory of a continu...